• /
• Rainy evening in the MATropolis

on 8 Nov 2023
• 42
• 187
• 2
• 6
• 1460
drawframe(1);
function drawframe(f)
persistent handles
if f == 1 || isempty(handles)
rng('default')
axes(colormap=gray,View=[172,22],Color='k')
hold on
camva(3)
h=randg(1,25); % building heights
b=bar3(h);
axis equal
g=.17; % height and width of windows
[m,k]=ndgrid(1:25); % (x,y) building indices
handles.w = gobjects(0);
for i=1:625 % loop through each building
% Interpolate building color according to building height
% Ideally this should be outside the loop but that would require more characters
q=b(m(i));
set(q,'CData',q.ZData,'FaceColor','interp')
for j=1:2 % loop through each column of windows
L=k(i)-.27+(j-1)*.37; % in k-a+(j-1)*b, a is the signed distance from center of building in row k(i) to left edge of window in column j and b is the interval between left edges of windows.
x=L+[0,g,g,0]'; % [L;R;R;L] edges of windows in building i, window column j
z=ones(4,1)*(.2:.4:h(i)-.2)+[0;0;g;g]; % 4xn height of windows for n levels of windows; in a:b:c a is height of first row of windows, b is the interval between window-bottoms, c is height of tallest window.
y=z*0+m(i)+.45; % 4xn location of front face of building
% Plot windows for building i, column j
for w = 1:width(z)
handles.w(end+1)=patch(z(:,w)*0+x,y(:,w),z(:,w),'y');
if rand(1) < 0.25
handles.w(end).FaceColor = 'none';
end
end
end
end
set(gcf,'color','k')
axis off
camva(2.5)
campos([58.256 185.12 25.621])
view([165.27, 7.0091])
clim([0 8])
% lightning light
handles.light = light('style','infinite',...
'Position',[-0.21759 1.3841 0.1921], ...
'Visible', 'off');
% lightning background
ax = gca;
handles.Z = patch(ax.XLim([2,2,1,1]),...
ax.YLim([1 1 1 1]),...
ax.ZLim([1 2 2 1]),...
[6 7.5 -3 -6],...
'Visible','off');
handles.rain = gobjects(1);
xlim(xlim) % lazy way of freezing axis limits
ylim(ylim)
zlim(zlim)
end
% randomly toggle on/off some window lights
p = 0.001; % portion (0:1) of lights to toggle
n = numel(handles.w);
idx = randperm(n,ceil(n*p));
colors = {handles.w(idx).FaceColor};
isnone = cellfun(@ischar, colors);
set(handles.w(idx(isnone)), 'FaceColor', 'y')
set(handles.w(idx(~isnone)), 'FaceColor', 'none')
% Show lightning on select frames
lightningFrames = [24 26 28 29 31:34];
onoff = ismember(f,lightningFrames);
handles.light.Visible = onoff;
handles.Z.Visible = onoff;
% Rain
if mod(f,2)>0 % update every odd frame
handles.rain.delete
xyRain0 = rand(3,300).*[26;26;10]-[0;0;2];
xyRain1 = xyRain0 + [-1;0;-2]; % orientation and length of rain lines
handles.rain = plot3([xyRain0(1,:);xyRain1(1,:)], ...
[xyRain0(2,:);xyRain1(2,:)], ...
[xyRain0(3,:);xyRain1(3,:)], ...
'color',[.7 .7 .7 .4]);
end
end