Starting today, recreational weed is legal in Canada. This news has some how lead me to find Cannabis Curve, a mathematical equation to draw Cannabis….!!!

So to celebrate? being 2nd country in the world (1st was Uruguay) to legalize the green stuff for fun, I decided I’ll try drawing cannabis curve with ggplot. Here’s the final results.

Cannabis_Final

Here’s the step I took, because I couldn’t really understand the mathematical equation, so I’ve break it down step by step to sort of understand what each part of equation is doing.

library(tidyverse)

cannabis <- tibble(
  t = seq(-pi,pi, length.out=1000),
  r1 = (1+.9*cos(8*t)), ## this will draw 8 petals  ## this number determines number of leafs!
  r2 = r1 * (1+.1*cos(24*t)), ## this make it pointy
  r3 = r2 * (.9+0.5*cos(200*t)), ## this makes it jaggy
  r4 = r3 * (1+sin(t)), ## Hmm.. I think I want to rorate it 90 degree... 
  r4_alt = r3 * (1+sin(t-pi/2)), ## one way to do it...
  r = (1+.9*cos(8*t)) * (1+.1*cos(24*t)) * (.9+0.5*cos(200*t)) * (1+sin(t))  ## Put all in line line!
) 

cannabis %>% 
  ggplot(aes(x=t, y=r1)) + 
  geom_path(color="#7ABA71", size=2) +
  coord_polar() +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "(1+.9*cos(8*t) draws 8 petals")

cannabis %>% 
  ggplot(aes(x=t, y=r2)) + 
  geom_path(color="#7ABA71", size=2) +
  coord_polar() +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "(1+.9*cos(8*t) * * (1+.1*cos(24*t)) makes the tip pointy")

cannabis %>% 
  ggplot(aes(x=t, y=r3)) + 
  geom_path(color="#7ABA71", size=0.5) +
  coord_polar() +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "(1+.9*cos(8*t) * * (1+.1*cos(24*t)) * (.9+0.5*cos(200*t)) makes zaggy")

cannabis %>% 
  ggplot(aes(x=t, y=r4)) + 
  geom_path(color="#7ABA71", size=0.5) +
  coord_polar(start=pi/2) +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "(1+.9*cos(8*t) * * (1+.1*cos(24*t)) * (.9+0.5*cos(200*t)) * (1+sin(t)) - OK Cool, Now 2 leaves are small!", subcaption="Notice I used start=pi/2 to rotate!")

cannabis %>% 
  ggplot(aes(x=t, y=r)) + 
  geom_polygon(fill="#499b4a", color="#74Ba71", size=0.1) +
  coord_polar(theta="x", start=pi/2) +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "Instead of using geom_path, I used geom_polygon")

I couldn’t figure out how to “crop” the polar coordinate image, so there’s lots of white space on final image, but I like my little cannabis!

I couldn’t find way to crop the image, but I realized I just need to plot the above in Cartesian coordinate…!

# To translate polar coordinate to cartesian, I can use below 
## x = r * cos(theta)
## y = r * sin(theta)

cannabis <- cannabis %>% mutate(
  x = r * cos(t),
  y = r * sin(t)
)

cannabis %>% 
  ggplot(aes(x=x, y=y)) + 
  geom_polygon(fill="#499b4a", color="#74Ba71", size=0.1) +
  theme_void(base_family="Roboto Condensed") +
  labs(title = "Cannabis on Cartesian Coordinate")