A Time-Lapse of the Night Sky

on 15 Oct 2021

% This image appears much better on my

% machine than here, online 😕

% Dark sky

q=50:-.5:1;

hold% hold before imagesc to keep ydir normal

Current plot held

imagesc(q'.*(q./q),[20,80])

% Add star rotation

t=0:2.5E-5:.1;

n=4001; % n=numel(t);

[x,y]=pol2cart(6*rand(1,n)+t',rand(1,n).*(t*0+142));

m=99;

plot(x+m,y+m,'w')

% Add grassy hill

x=[1,1:.1:m,m];

y=[1,cos(0:.1/98:1)*30+randg(1,1,981),1];

fill(x,y,'k')

% Add house silhouette

fill([4 4 3 6 9 8 8]*5,[5 7 7 8 7 7 5]*5,'k')

axis([1 m 1 m])

camva(4)

colormap gray

Load full remix tree