{"id":2061,"date":"2024-12-16T08:20:47","date_gmt":"2024-12-16T13:20:47","guid":{"rendered":"https:\/\/blogs.mathworks.com\/graphics-and-apps\/?p=2061"},"modified":"2024-12-16T10:08:52","modified_gmt":"2024-12-16T15:08:52","slug":"pinecode-creating-pinecones-with-fibonacci-spirals","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/graphics-and-apps\/2024\/12\/16\/pinecode-creating-pinecones-with-fibonacci-spirals\/","title":{"rendered":"Pinecode: Creating pinecones with Fibonacci spirals"},"content":{"rendered":"
\r\n
Today we explore phyllotactic spirals, a naturally emerging pattern found in densely packed plant structures such as pineapples, sunflowers, and pinecones. This generative spiral is based on tau<\/em>, a constant angular offset between successive seeds or scales. When tau<\/em> equals 2\u03c0 times the golden ratio, $\\frac{1+\\sqrt{5}}{2}\\approx 1\\ldotp 618$ the pattern is known as a Fibonacci spiral. Experiment with different values of tau<\/em> to see how the spirals transform.<\/div>\r\n
<\/div>\r\n
\"\"<\/div>\r\n
\r\n
<\/div>\r\n
<\/div>\r\n
tau = pi*(1+sqrt(5));<\/div>\r\n
n = 1:300;<\/div>\r\n
r = sqrt(n);<\/div>\r\n
theta = n * tau;<\/div>\r\n
polarscatter(theta,r,'bo')<\/div>\r\n<\/div>\r\n
<\/div>\r\n<\/div>\r\n
<\/div>\r\n
<\/div>\r\n
\r\n
Eric applies this principle to the 3D Cartesian domain in his parameterized pinecode.<\/div>\r\n
\r\n\r\n
\r\n\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n\r\n
\"\"<\/td>\r\nGuest Writer: Eric Ludlam<\/a><\/span><\/strong>\r\n\r\nEric, a seasoned MATLAB graphics expert, shares his approach to replicating the natural patterns of a pinecone through code. Learn about pinecone geometry and the efficacy of parameterized graphics. See Eric's full bio on our contributors page<\/a>.<\/td>\r\n<\/tr>\r\n<\/tbody>\r\n<\/table>\r\n
\r\n
\r\n
<\/div>\r\n
Let's explore a graphics algorithm for creating a pinecone. Here are a couple pictures of pinecones I found. This first was on the coast in Maine, and the second I found along a hiking trail near the MathWorks campus in Natick, Massachusetts.<\/div>\r\n
\"pc_ontree.jpg\" \"\"<\/div>\r\n
Pinecones are made of \"scales\", which you can think of as petals on a flower. In this article, \"scale\" specifically refers to these structures on the pinecone and not to the common mathematical transformation.<\/div>\r\n
The pinecone pictured on the left has a shape like an upside-down teardrop. The scales are wide and get smaller and more densely packed near the base. The base of the pinecone on the right shows how the scales radiate outward in a Fibonacci spiral, highlighted with yellow dots.<\/div>\r\n
The pinecone algorithm below attempts to replicate these patterns in a single surface object. It starts with a Fibonacci spiral which helps us compute the locations of the tips for each of the scales. Next, we compute the shape of one scale and replicate that shape for each location computed in our spiral.<\/div>\r\n
At the end, I add a little branch and some pine needles using a line and patch.<\/div>\r\n
If you open this blog post in MATLAB Online using the link below, all the parameters at the beginning of the script can be manipulated with sliders and color pickers making it easy for you to design your own pinecone!<\/div>\r\n
Click here to try it!<\/div>\r\n<\/div>\r\n
\"\"<\/a><\/div>\r\n
\r\n
\r\n
\r\n
% Pinecone Attributes<\/span><\/div>\r\n<\/div>\r\n
\r\n
nscales = 140; % Number of scales<\/span><\/div>\r\n<\/div>\r\n
\r\n
height = 2.5; % Height of inner part of pinecone<\/span><\/div>\r\n<\/div>\r\n
\r\n
swt = 0.3; % Scale width in theta<\/span><\/div>\r\n<\/div>\r\n
\r\n
sth = 0.2; % Scale thickness (as a curve in height)<\/span><\/div>\r\n<\/div>\r\n
\r\n
stilt = 1.2; % From inside to outside, tile of each scale<\/span><\/div>\r\n<\/div>\r\n
\r\n
nn = 18; % Number of needles<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Exponents for curves<\/span><\/div>\r\n<\/div>\r\n
\r\n
hexp = 2.8; % scale density over height<\/span><\/div>\r\n<\/div>\r\n
\r\n
rexp = 2; % radius shape across height<\/span><\/div>\r\n<\/div>\r\n
\r\n
scexp = 3; % scale rounded shape (end)<\/span><\/div>\r\n<\/div>\r\n
\r\n
suexp = 3; % scale roundness in height<\/span><\/div>\r\n<\/div>\r\n
\r\n
stexp = 2.1; % scale tilt for height<\/span><\/div>\r\n<\/div>\r\n
\r\n
stuexp = 2; % scale thickness of U by height<\/span><\/div>\r\n<\/div>\r\n
\r\n
cexp = 2; % color gradient<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Colors<\/span><\/div>\r\n<\/div>\r\n
\r\n
brown1 = [0 0 0]; % Dark center color<\/span><\/div>\r\n<\/div>\r\n
\r\n
brown2 = [0.470 0.240 0.010]; % Brigher outer color<\/span><\/div>\r\n<\/div>\r\n
\r\n
edge = [1 1 1]; % Scale edging<\/span><\/div>\r\n<\/div>\r\n
\r\n
green = [0.060 0.590 0.400]; % Pine Needles<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Resolution<\/span><\/div>\r\n<\/div>\r\n
\r\n
vps = 17; % verts per scale : An odd number so 1 vert is in center<\/span><\/div>\r\n<\/div>\r\n
\r\n
vpr = 20; % verts per radius<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Scale locations<\/span><\/div>\r\n<\/div>\r\n
\r\n
FIB = (1+sqrt(5))*(1:nscales); % Fibonacci golden ratio sequence as factor of pi\u00a0\u00a0<\/span><\/div>\r\n<\/div>\r\n
\r\n
H = linspace(0,1,nscales).^hexp*height; % Height of the root of each scale<\/span><\/div>\r\n<\/div>\r\n
\r\n
R = 1-abs(linspace(-1,1,nscales).^rexp); % Radius of pinecone over height<\/span><\/div>\r\n<\/div>\r\n
\r\n
U = 1-abs(linspace(-1,1,nscales).^stuexp); % Thickness of scale U shape by height<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Geometry of the scales at the locations<\/span><\/div>\r\n<\/div>\r\n
\r\n
ST = reshape((linspace(-.5,.5,vps)*swt+FIB')', 1,[]); % Scale Theta<\/span><\/div>\r\n<\/div>\r\n
\r\n
SR = reshape(((1-abs(linspace(-1,1,vps).^scexp)).*R')',1,[]); % Scale Radii<\/span><\/div>\r\n<\/div>\r\n
\r\n
SH = reshape((ones(1,vps).*H'+abs(linspace(-1,1,vps).^suexp.*U')*sth)',1,[]); % Scale Height<\/span><\/div>\r\n<\/div>\r\n
\r\n
MR = linspace(0,1,vpr)'.*SR;<\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Compute final geometry of the pinecone.<\/span><\/div>\r\n<\/div>\r\n
\r\n
X = cospi(ST).*MR;<\/div>\r\n<\/div>\r\n
\r\n
Y = sinpi(ST).*MR;<\/div>\r\n<\/div>\r\n
\r\n
Z = SH.*ones(vpr,1)+MR.*linspace(0,1,vps*nscales).^stexp*stilt;<\/div>\r\n<\/div>\r\n
\r\n
C = linspace(0,1,vpr).^cexp'.*ones(1,vps*nscales);<\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Plot pinecone<\/span><\/div>\r\n<\/div>\r\n
\r\n
set(gcf,'Color'<\/span>,'#fffafa'<\/span>); % Snow<\/span><\/div>\r\n<\/div>\r\n
\r\n
plot3(X(end,:),Y(end,:),-Z(end,:),'-'<\/span>,'Color'<\/span>,edge);<\/div>\r\n<\/div>\r\n
\r\n
surface(X,Y,-Z,C,'EdgeColor'<\/span>,'none'<\/span>);<\/div>\r\n<\/div>\r\n
\r\n
shading interp<\/span><\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Dark center to Brown colormap<\/span><\/div>\r\n<\/div>\r\n
\r\n
colormap(interp1([1 256],[brown1; brown2],1:256));<\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Branch<\/span><\/div>\r\n<\/div>\r\n
\r\n
bx = [ 0 .5 1 1.5 2 ];<\/div>\r\n<\/div>\r\n
\r\n
by = [ 0 0 0 0 0 ];<\/div>\r\n<\/div>\r\n
\r\n
bz = [ 0 .05 .2 .25 .2 ];<\/div>\r\n<\/div>\r\n
\r\n
line('XData'<\/span>,bx,'YData'<\/span>,by,'ZData'<\/span>,bz,'Color'<\/span>,brown1, 'LineWidth'<\/span>,5);<\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Pine Needles<\/span><\/div>\r\n<\/div>\r\n
\r\n
NT = linspace(0,5,nn)';<\/div>\r\n<\/div>\r\n
\r\n
NR = linspace(.1,.8,nn)';<\/div>\r\n<\/div>\r\n
\r\n
NV = [ 0 0 0<\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0-ones(nn,1)*2, cospi(NT).*NR, sinpi(NT).*NR+.9 ];<\/div>\r\n<\/div>\r\n
\r\n
NF = [ ones(nn,1) (1:nn)'+1 ];<\/div>\r\n<\/div>\r\n
\r\n
patch('Vertices'<\/span>,NV,'Faces'<\/span>,NF,'EdgeColor'<\/span>,green,'FaceColor'<\/span>,'none'<\/span>,...<\/span><\/div>\r\n<\/div>\r\n
\r\n
\u00a0\u00a0\u00a0\u00a0\u00a0\u00a0'LineWidth'<\/span>,1.5);<\/div>\r\n<\/div>\r\n
\r\n
<\/div>\r\n<\/div>\r\n
\r\n
% Setup Axes<\/span><\/div>\r\n<\/div>\r\n
\r\n
set(gca,'Position'<\/span>,[0 0 1 1],'Clipping'<\/span>,'off'<\/span>);<\/div>\r\n<\/div>\r\n
\r\n
daspect([1 1 1]);<\/div>\r\n<\/div>\r\n
\r\n
axis([-1.3 1.3 -1.3 1.3 -height+.2 .5],'off'<\/span>)<\/div>\r\n<\/div>\r\n
\r\n
material([.6 1 .3])<\/div>\r\n<\/div>\r\n
\r\n
lighting gouraud<\/span><\/div>\r\n<\/div>\r\n
\r\n
light('Position'<\/span>,[13.2 7.2 16.5]);<\/div>\r\n<\/div>\r\n<\/div>\r\n\r\n\r\n