- /
-
Late summer evening
on 12 Oct 2021
- 29
- 33
- 0
- 0
- 273
a=508;
X=linspace(-125,124,a);
A=[18 16 13];
C=[10 8 8
1 .5 3]/10;
z=([0:.1:21].'./A).^(A/4);
y=rescale(z.*exp(-z),C(2,:),C(1,:));
[t,r]=cart2pol(X,X'-90);
colormap([1 1 .95
flip(y)]);
s=r;
s(s<12)=1;
imagesc(s);
hold
image(ones(a),'AlphaData',1E4*abs(ifft2(r.^-1.5.*cos(7*rand(a)))));
axis equal off