显然,人们已经知道如何在Mathematica和LaTeX中制作xkcd风格的图形。我们可以用R来做吗?Ggplot2-ers吗?一个geom_xkcd和/或theme_xkcd?
我猜在基本图形,par(xkcd=TRUE)?我该怎么做?
作为ggplot2中的第一次尝试(下面将更优雅地展示),将jitter参数添加到一行中可以获得很棒的手绘外观。所以- - -
ggplot(mapping=aes(x=seq(1,10,.1), y=seq(1,10,.1))) +
geom_line(position="jitter", color="red", size=2) + theme_bw()
这是一个很好的例子-但是坐标轴和字体看起来更棘手。但是字体似乎已经解决了(如下)。处理斧头的唯一方法是把它们空出来,然后用手画出来吗?有没有更优雅的解决方案?特别是,在ggplot2中,是否可以修改新主题系统中的element_line以接受类似jitter的参数?
下面是基于xkcd论坛和extrafont包的链接对字体的尝试:
如上所述,在xkcd网站上有一个关于字体的论坛讨论:我抓住了我能找到的第一个,可能还有其他(更好的?)选项(@jebyrnes在上面的评论中发布了另一个可能的字体来源——TTF文件在这里;有人报告了404错误的来源,你可以在这里或这里,替代这些url适当xkcdFontURL下面;你可能需要更加努力地检索github上发布的链接)
xkcdFontURL <- "http://simonsoftware.se/other/xkcd.ttf"
download.file(xkcdFontURL,dest="xkcd.ttf",mode="wb")
(这是为了快速,一次性的使用:为了经常使用,你应该把它放在一些标准的系统字体目录中。)
library(extrafont)
关于字体最有用的信息是在extrafont github网站上——这是从那里截取的
font_import(".") ## because we downloaded to working directory
loadfonts()
从github网站上逐字逐句的例子:
library(ggplot2)
p <- ggplot(mtcars, aes(x=wt, y=mpg)) + geom_point() +
ggtitle("Fuel Efficiency of 32 Cars") +
xlab("Weight (x1000 lb)") + ylab("Miles per Gallon") +
theme(text=element_text(size=16, family="xkcd"))
ggsave("xkcd_ggplot.pdf", plot=p, width=4, height=4)
## needed for Windows:
## Sys.setenv(R_GSCMD = "C:/Program Files/gs/gs9.05/bin/gswin32c.exe")
embed_fonts("xkcd_ggplot.pdf")
下面是我使用上面的一些代码编写的ggplot2代码行:
ggplot()+geom_line(aes(x=seq(0,1,length.out=1000),y=sin(x)),position=position_jitter(width=0.02),lwd=1.5,col="white")+
geom_line(aes(x=seq(0,1,length.out=1000),y=sin(x)),position=position_jitter(width=0.004),lwd=1.4,col="red")+
geom_line(aes(x=seq(0,1,length.out=1000),y=cos(x)),position=position_jitter(width=0.02),lwd=1.5,col="white")+
geom_line(aes(x=seq(0,1,length.out=1000),y=cos(x)),position=position_jitter(width=0.004),lwd=1.4,col="blue")+
theme_bw()+theme(panel.grid.major=element_blank(),panel.grid.minor=element_blank())
不确定如何替换轴,但可以使用相同的方法与抖动。然后从XKCD导入字体,并使用geom_text进行分层。
这是一个非常,非常粗糙的开始,只覆盖(部分)手绘的外观和感觉的线条。这将需要一些工作来自动化,但在响应函数中添加一些AR(1)噪声可以使它看起来有点手绘
set.seed(551)
x <- seq(0, 1, length.out = 1000)
y <- sin(x)
imperfect <- arima.sim(n = length(y), model = list(ar = c(.9999)))
imperfect <- scale(imperfect)
z <- y + imperfect*.005
plot(x, z, type = "l", col = "blue", lwd = 2)
基本画线功能:
xkcd_line <- function(x, y, color) {
len <- length(x);
rg <- par("usr");
yjitter <- (rg[4] - rg[3]) / 1000;
xjitter <- (rg[2] - rg[1]) / 1000;
x_mod <- x + rnorm(len) * xjitter;
y_mod <- y + rnorm(len) * yjitter;
lines(x_mod, y_mod, col='white', lwd=10);
lines(x_mod, y_mod, col=color, lwd=5);
}
基本轴:
xkcd_axis <- function() {
rg <- par("usr");
yaxis <- 1:100 / 100 * (rg[4] - rg[3]) + rg[3];
xaxis <- 1:100 / 100 * (rg[2] - rg[1]) + rg[1];
xkcd_line(1:100 * 0 + rg[1] + (rg[2]-rg[1])/100, yaxis,'black')
xkcd_line(xaxis, 1:100 * 0 + rg[3] + (rg[4]-rg[3])/100, 'black')
}
和示例代码:
data <- data.frame(x=1:100)
data$one <- exp(-((data$x - 50)/10)^2)
data$two <- sin(data$x/10)
plot.new()
plot.window(
c(min(data$x),max(data$x)),
c(min(c(data$one,data$two)),max(c(data$one,data$two))))
xkcd_axis()
xkcd_line(data$x, data$one, 'red')
xkcd_line(data$x, data$two, 'blue')
生产: