其他
mathart：一个富有逼格的R包
我在《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 <- mollusc()
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, df4)
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")
利萨如曲线
set.seed(1)
df <- 1:100 %>% map_df(~lissajous(a = runif(1, 0, 10), A = runif(1, 0, 1)), .id = "id")
ggplot() +
geom_path(aes(x, y), df, size = 0.25, lineend = "round") +
facet_wrap(~id, nrow = 10) +
coord_equal() +
theme_blankcanvas(margin_cm = 1)
k-nearest neighbour graph
set.seed(1)
df <- heart3() %>%
sample_n(10000) %>%
k_nearest_neighbour_graph(40)
ggplot() +
geom_segment(aes(x, y, xend = xend, yend = yend), df, size = 0.03) +
coord_equal() +
theme_blankcanvas(margin_cm = 0)