{"id":8804,"date":"2022-12-22T16:40:05","date_gmt":"2022-12-22T21:40:05","guid":{"rendered":"https:\/\/blogs.mathworks.com\/community\/?p=8804"},"modified":"2024-01-07T17:27:35","modified_gmt":"2024-01-07T22:27:35","slug":"cookies-code","status":"publish","type":"post","link":"https:\/\/blogs.mathworks.com\/community\/2022\/12\/22\/cookies-code\/","title":{"rendered":"Cookies & Code"},"content":{"rendered":"
\n
Seasonal greetings of the Greater Solstice Metropolitan Area! As I write this, the solstitial sun sets on the shortest day of my northern hemispheric year. Winter is here. So naturally, my mind turns to cookies.<\/div>\n
\"all_cookies.png\"<\/div>\n
The inspiration for this winter daydream comes from some code written by Eric Ludlam. Eric, in addition to being in charge of the MATLAB graphics team, is an inveterate tinkerer. He loves making fun graphics code. You may remember some of his pumpkin-inspired creations that I described in this space last year<\/a>. Or you may recall his many entries from the latest Mini Hack<\/a> (and the one before that<\/a>).<\/div>\n
He's back, this time with code for seasonal cookies and snowflakes. If you want to look at his code, you can find it in his Digital Art with MATLAB<\/a> repo. In fact, let's go look at it now.<\/div>\n
We can start by cloning our own copy in MATLAB. You can do this directly from the Current Folder Browser. Select \"Manage Files...\" from the Source Control context menu.<\/div>\n
\"cloning.png\"<\/div>\n
Enter the Repository Path (https:\/\/github.com\/zappo2\/digital-art-with-matlab<\/a>), and the files come right down.<\/div>\n
\"cloning2.png\"<\/div>\n
Eric made some cookie code using blobby implicit surfaces as described by Jim Blinn. We can use Eric's template as cookie-cutter (heh) to stamp out other designs. Shown below is the code from gingerbreadperson.m, only I've substituted my own design at the top: a candy cane.<\/div>\n
\n
\n
% Create a gingerbread person cookie<\/span><\/span><\/div>\n<\/div>\n
\n
K = [<\/span>' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxxxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
'xxx xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
'xxx xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
'xxx xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span><\/span><\/div>\n<\/div>\n
\n
' xxx '<\/span>]; <\/span><\/div>\n<\/div>\n
\n
% Convert K into array of blob centers<\/span><\/span><\/div>\n<\/div>\n
\n
[y,x]=find(~isspace(K));<\/span><\/div>\n<\/div>\n
\n
% Scale and offset centers into 3D cookie space<\/span><\/span><\/div>\n<\/div>\n
\n
C=[x*5,y*5]+5; C(:,3)=10;<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Compute a volume using a short implementation of Blinn's blobs<\/span><\/span><\/div>\n<\/div>\n
\n
nx=size(K,2)*5+10; ny=size(K,1)*5+10; nz=20; % Size of volume<\/span><\/span><\/div>\n<\/div>\n
\n
[y,x,z]=ndgrid(1:ny,1:nx,1:nz);<\/span><\/div>\n<\/div>\n
\n
vol=zeros(ny,nx,nz);<\/span><\/div>\n<\/div>\n
\n
for <\/span>i=1:size(C,1)<\/span><\/div>\n<\/div>\n
\n
vol=vol+exp(-.05*((C(i,1)-x).^2 + (C(i,2)-y).^2 + (C(i,3)-z).^2));<\/span><\/div>\n<\/div>\n
\n
end<\/span><\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Compute isosurface from the blobs<\/span><\/span><\/div>\n<\/div>\n
\n
S=isosurface(vol,.3);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Draw the isosurface using patch<\/span><\/span><\/div>\n<\/div>\n
\n
newplot<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
patch(S,'FaceColor'<\/span>,'#a56c3c'<\/span>,'EdgeColor'<\/span>,'none'<\/span>); % cookie<\/span><\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
axis equal ij off<\/span><\/span><\/div>\n<\/div>\n
\n
lighting gouraud<\/span><\/span><\/div>\n<\/div>\n
\n
camlight<\/span><\/div>\n<\/div>\n
\n
material([.6 .9 .3 ])<\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
It's not the most shapely candy cane, but it will do. Let's add some icing.<\/div>\n
\n
\n
newplot<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
patch(S,'FaceColor'<\/span>,'interp'<\/span>,'EdgeColor'<\/span>,'none'<\/span>,'FaceVertexCData'<\/span>,S.vertices(:,3));<\/span><\/div>\n<\/div>\n
\n
colormap(validatecolor({ '#a56c3c' '#f09090' <\/span>}, 'multiple'<\/span>));<\/span><\/div>\n<\/div>\n
\n
clim([10 20]);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
axis equal ij off<\/span><\/span><\/div>\n<\/div>\n
\n
lighting gouraud<\/span><\/span><\/div>\n<\/div>\n
\n
camlight<\/span><\/div>\n<\/div>\n
\n
material([.6 .9 .3 ])<\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
Obviously, it's easy to draw with X's into the variable K to make any cookie shape that pleases you. What happens if we make it random?<\/div>\n
\n
\n
% Create a gingerbread blob<\/span><\/span><\/div>\n<\/div>\n
\n
K = double(rand(10,10)>0.5);<\/span><\/div>\n<\/div>\n
\n
K(K==0) = char('*'<\/span>);<\/span><\/div>\n<\/div>\n
\n
K = char(K);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Convert K into array of blob centers<\/span><\/span><\/div>\n<\/div>\n
\n
[y,x]=find(~isspace(K));<\/span><\/div>\n<\/div>\n
\n
% Scale and offset centers into 3D cookie space<\/span><\/span><\/div>\n<\/div>\n
\n
C=[x*5,y*5]+5; C(:,3)=10;<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Compute a volume using a short implementation of Blinn's blobs<\/span><\/span><\/div>\n<\/div>\n
\n
nx=size(K,2)*5+10; ny=size(K,1)*5+10; nz=20; % Size of volume<\/span><\/span><\/div>\n<\/div>\n
\n
[y,x,z]=ndgrid(1:ny,1:nx,1:nz);<\/span><\/div>\n<\/div>\n
\n
vol=zeros(ny,nx,nz);<\/span><\/div>\n<\/div>\n
\n
for <\/span>i=1:size(C,1)<\/span><\/div>\n<\/div>\n
\n
vol=vol+exp(-.05*((C(i,1)-x).^2 + (C(i,2)-y).^2 + (C(i,3)-z).^2));<\/span><\/div>\n<\/div>\n
\n
end<\/span><\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
% Compute isosurface from the blobs<\/span><\/span><\/div>\n<\/div>\n
\n
newplot<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
S=isosurface(vol,.3);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
patch(S,'FaceColor'<\/span>,'interp'<\/span>,'EdgeColor'<\/span>,'none'<\/span>,'FaceVertexCData'<\/span>,S.vertices(:,3));<\/span><\/div>\n<\/div>\n
\n
colormap(validatecolor({ '#a56c3c' '#8090f0' <\/span>}, 'multiple'<\/span>));<\/span><\/div>\n<\/div>\n
\n
clim([10 20]);<\/span><\/div>\n<\/div>\n
\n
<\/div>\n<\/div>\n
\n
axis equal ij off<\/span><\/span><\/div>\n<\/div>\n
\n
lighting gouraud<\/span><\/span><\/div>\n<\/div>\n
\n
camlight<\/span><\/div>\n<\/div>\n
\n
material([.6 .9 .3 ])<\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
I like my random cookies, and I'm now inclined to try some in the kitchen.<\/div>\n
I'll close with Eric's snowflake code. Here's the default.<\/div>\n
\n
\n
snowflake<\/span><\/div>\n<\/div>\n
\n
axis off<\/span><\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
But, just as with the cookies (and as with real-life snowflakes), we can opt for a random flake.<\/div>\n
\n
\n
steps = randi(10,1,randi(8)+2)<\/span><\/div>\n
\n
\n
\n
\n
steps = <\/span>1\u00d710<\/span><\/div>\n<\/div>\n
\n
7 6 9 1 1 1 4 2 10 5<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n
\n
snowflake(steps)<\/span><\/div>\n<\/div>\n
\n
axis off<\/span><\/span><\/div>\n
\n
\n
<\/div>\n
<\/div>\n
<\/div>\n
<\/div>\n<\/div>\n<\/div>\n<\/div>\n<\/div>\n
Happy holidays, and thanks for the cookies, Eric!<\/div>\n<\/div>\n","protected":false},"excerpt":{"rendered":"
<\/div>\n

\nSeasonal greetings of the Greater Solstice Metropolitan Area! As I write this, the solstitial sun sets on the shortest day of my northern hemispheric year. Winter is here. So naturally, my mind... read more >><\/a><\/p>\n","protected":false},"author":69,"featured_media":0,"comment_status":"open","ping_status":"closed","sticky":false,"template":"","format":"standard","meta":[],"categories":[1],"tags":[],"_links":{"self":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8804"}],"collection":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/users\/69"}],"replies":[{"embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/comments?post=8804"}],"version-history":[{"count":3,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8804\/revisions"}],"predecessor-version":[{"id":8955,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/posts\/8804\/revisions\/8955"}],"wp:attachment":[{"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/media?parent=8804"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/categories?post=8804"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/blogs.mathworks.com\/community\/wp-json\/wp\/v2\/tags?post=8804"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}