- /
-
Misty moonshine
on 10 Oct 2021
- 39
- 205
- 1
- 0
- 268
b=508
image(zeros(b))
hold
X=linspace(-124,124,b)
a=-90:90
x=70*cosd(a);
y=70*sind(a);
Y=[x x*.7;y -y]
fill(Y(1,:)+270,Y(2,:)+120,'w')
axis([0 b 0 b])
colormap(bone)
[t,d]=cart2pol(X,X')
image(200*ones(b),'AlphaData',rescale(abs(ifft2((((d.^2).^-.8)).*cos(6*rand(b)))),0,1))
axis off