如何用R创建Phyllotaxis螺旋

Question

如何在R中创建Phyllotaxis螺旋？此类图像的示例显示在本页底部：http://www.mathrecreation.com/2008/09/phyllotaxis-spirals.html图像附在下面：

这些也类似于费马的螺旋：http://en.wikipedia.org/wiki/Fermat%27s_spiral

我搜索并找到了这种方法：Spiral Graph in R 但它不会创建上述模式之一。

Answer 1

使用http://algorithmicbotany.org/papers/abop/abop-ch4.pdf

中的公式

golden.ratio = (sqrt(5) + 1)/2
fibonacci.angle=360/(golden.ratio^2)
c=1
num_points=630
x=rep(0,num_points)
y=rep(0,num_points)

for (n in 1:num_points) {
    r=c*sqrt(n)
    theta=fibonacci.angle*(n)
    x[n]=r*cos(theta)
    y[n]=r*sin(theta)
}
plot(x,y,axes=FALSE,ann=FALSE,pch=19,cex=1)