- /
-
A Time-Lapse of the Night Sky
on 15 Oct 2021
- 72
- 531
- 5
- 0
- 280
% 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
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