我在《Use ggplot2》一文中画过一个「囧」字，也画过一只蝴蝶，这个R包，可以让你画出更多的数学艺术。
安装
install.packages(c('devtools', 'mapproj', 'tidyverse', 'ggforce'))
devtools::install_github('marcusvolz/mathart')
devtools::install_github('marcusvolz/ggart')
加载相应的包
library(mathart)
library(ggart)
library(ggforce)
library(tidyverse)
燃烧吧火鸟
df <- rbind(parrot()="" %="">% mutate(id = 1),
stork() %>% mutate(id = 2),
magpie() %>% mutate(id = 3))
ggplot() +
geom_point(aes(x, y, alpha = r), df, size = 0.03) +
facet_wrap(~id, nrow = 2, scales = 'free') +
scale_alpha_continuous(range = c(0.03, 0.06)) +
theme_blankcanvas(margin_cm = 1)
画植物
df <- rbind(olive_branch()="" %="">% mutate(id = 1),
palm_branch() %>% mutate(id = 2),
branch() %>% mutate(id = 3))
ggplot() +
geom_circle(aes(x0 = x, y0 = y, r = r), df, size = 0.03, alpha = 0.1) +
coord_equal() +
facet_wrap(~id, nrow = 3) +
theme_blankcanvas(margin_cm = 1)
画贝壳
df <->
df1 <- df="" %="">% mutate(id = 1)
df2 <- df="" %="">% mutate(id = 2)
df3 <- df="" %="">% mutate(id = 3)
ggplot() +
geom_point(aes(x, y), df1, size = 0.03, alpha = 0.03) +
geom_path( aes(x, y), df1, size = 0.03, alpha = 0.03) +
geom_point(aes(x, z), df2, size = 0.03, alpha = 0.03) +
geom_path( aes(x, z), df2, size = 0.03, alpha = 0.03) +
geom_point(aes(y, z), df3, size = 0.03, alpha = 0.03) +
geom_path( aes(y, z), df3, size = 0.03, alpha = 0.03) +
facet_wrap(~id, nrow = 2, scales = 'free') +
theme_blankcanvas(margin_cm = 0.5)
画蝴蝶
与我在《Use ggplot2》一文中画的蝴蝶不一样哦。
df <- rbind(butterfly1()="" %="">% mutate(id = 1),
butterfly2() %>% mutate(id = 2),
butterfly3() %>% mutate(id = 3))
ggplot() +
geom_point(aes(x, y, alpha = r), df, size = 0.03) +
coord_equal() +
facet_wrap(~id, nrow = 2) +
scale_alpha_continuous(range = c(0.03, 0.16)) +
theme_blankcanvas(margin_cm = 0)
鲜肉们最爱的心
df <- rbind(heart1()="" %="">% mutate(id = 1),
heart2() %>% mutate(id = 2),
heart3() %>% mutate(id = 3),
heart4() %>% mutate(id = 4))
ggplot() +
geom_path(aes(x, y), df, size = 0.25, lineend = 'round') +
facet_wrap(~id, nrow = 2, scales = 'free') +
theme_blankcanvas(margin_cm = 1)
各种小动物
df1 <- ant()="" %="">% mutate(id = 1)
df2 <- spider()="" %="">% mutate(id = 2)
df3 <- millipede()="" %="">% mutate(id = 3)
df <- rbind(df1="" %="">% mutate(x = (x - min(df1$x))*0.9 + 0.2, y = (y - min(df1$y))*0.9 + 1, r = r * 0.75),
df2 %>% mutate(x = (x - min(df2$x))*1.5 + 2.4, y = (y - min(df2$y))*1.5 + 2.2),
df3 %>% mutate(x = (x - min(df3$x))*1.75 + 2.4, y = (y - min(df3$y))*1.75))
ggplot() +
geom_circle(aes(x0 = x, y0 = y, r = r), df, size = 0.03, alpha = 0.03, n = 100) +
geom_path(aes(x, y, group = id), df, size = 0.5, alpha = 0.5) +
coord_equal() +
theme_blankcanvas(margin_cm = 1)
harmonographs
df1 <- harmonograph(a1="">1, A2 = 1, A3 = 1, A4 = 1,
d1 = 0.004, d2 = 0.0065, d3 = 0.008, d4 = 0.019,
f1 = 3.001, f2 = 2, f3 = 3, f4 = 2,
p1 = 0, p2 = 0, p3 = pi/2, p4 = 3*pi/2) %>% mutate(id = 1)
df2 <- harmonograph(a1="">1, A2 = 1, A3 = 1, A4 = 1,
d1 = 0.0085, d2 = 0, d3 = 0.065, d4 = 0,
f1 = 2.01, f2 = 3, f3 = 3, f4 = 2,
p1 = 0, p2 = 7*pi/16, p3 = 0, p4 = 0) %>% mutate(id = 2)
df3 <- harmonograph(a1="">1, A2 = 1, A3 = 1, A4 = 1,
d1 = 0.039, d2 = 0.006, d3 = 0, d4 = 0.0045,
f1 = 10, f2 = 3, f3 = 1, f4 = 2,
p1 = 0, p2 = 0, p3 = pi/2, p4 = 0) %>% mutate(id = 3)
df4 <- harmonograph(a1="">1, A2 = 1, A3 = 1, A4 = 1,
d1 = 0.02, d2 = 0.0315, d3 = 0.02, d4 = 0.02,
f1 = 2, f2 = 6, f3 = 1.002, f4 = 3,
p1 = pi/16, p2 = 3*pi/2, p3 = 13*pi/16, p4 = pi) %>% mutate(id = 4)
df <- rbind(df1,="" df2,="" df3,="">
p <- ggplot()="">
geom_path(aes(x, y), df, alpha = 0.25, size = 0.5) +
coord_equal() +
facet_wrap(~id, nrow = 2) +
theme_blankcanvas(margin_cm = 0)
ggsave('harmonograph01.png', p, width = 20, height = 20, units = 'cm')
