Back to Gallery
  • Remix
  • Share
  • New Entry

  • Josh

  • /

  • lots of random ripples

on 26 Nov 2023
  • 21
  • 21
  • 0
  • 0
  • 811
drawframe(10);
Write your drawframe function below
function drawframe(f)
% adding f random pebbles at each frame iteration we see patterns emerge at
% higher orders within the chaos
% set tablesize nxn (ts), output image (o), timer matrix (tm), and seed
% random number generator (rng) so randomized pebbles remain in consistent
% locations as each frame is generated
ts=100;
o=zeros(ts);
tm=o;
rng(9);
% use starting image (im) for pebble locations
imr=[27 33 66];
imc=[30 70 50];
for i1=1:numel(imr)
tm(imr(i1),imc(i1))=1;
end
% sf (scale factor) is a multiplier allowing for deeper iterations at the cost of animation
% resolution (skips frames), due to limitation of f=1:48 (48 frame per animation)
sf=1;
% loop to iterate through previous frames (1 -> f-1) to generate current frame (f)
for i1=1:sf*f
if f==1
break;
end
% add random pebble every nth iteration
n=1;
if mod(i1,n)==0
nP=randi(ts^2,i1);
if tm(nP)==0
tm(nP)=1;
end
end
% pad output and timer with extra space for ripple overflow to be trimmed later
o=paddata(o,[ts+2*sf*f ts+2*sf*f],Side="both");
tm=paddata(tm,[ts+2*sf*f ts+2*sf*f],Side="both");
[rP,cP]=find(tm>0);
% core loop for evaluating locations > 0 in timer matrix (value equates to how
% many turns have lapsed since pebble dropped at this location) and drawing new ripples at
% these locations while also erasing old ripples
for i2=1:numel(rP)
t=tm(rP(i2),cP(i2));
newrip=ones(2*t+1);
newrip(2:end-1,2:end-1)=0;
o(rP(i2)-t:rP(i2)+t,cP(i2)-t:cP(i2)+t)=o(rP(i2)-t:rP(i2)+t,cP(i2)-t:cP(i2)+t)+newrip;
oldrip=ones(2*t-1);
if t>1
oldrip(2:end-1,2:end-1)=0;
end
o(rP(i2)-(t-1):rP(i2)+(t-1),cP(i2)-(t-1):cP(i2)+(t-1))=o(rP(i2)-(t-1):rP(i2)+(t-1),cP(i2)-(t-1):cP(i2)+(t-1))-oldrip;
tm(rP(i2),cP(i2))=t+1;
end
% trim overflow padding so output and timer both back to ts
c=size(o/2)/2;
o=o(c-ts/2:c+ts/2,c-ts/2:c+ts/2);
tm=tm(c-ts/2:c+ts/2,c-ts/2:c+ts/2);
end
o=uint16(o);
cmp=colormap(turbo(1+max(max(o))));
o=ind2rgb(o,cmp);
imshow(o)
end
Animation
Remix Tree