- /
-
Aolian ripples
on 26 Oct 2021
- 7
- 50
- 1
- 0
- 279
N=400;
H=randn(N);
y=1:N;
r=@round;
k=zeros(N);
k(r(rand(4)*N^2))=1;
k=conv2(k,ones(29),'same')<1;
for n=1:3e3
d=gradient(H);
g=(d+1)/6;
m=mod(r(y+randn(N).*k*4),N)+1;
o=mod(r(y'+(9-9*d).*k),N)+1;
H=conv2(H-g+accumarray([o(:),m(:)],g(:),[N,N]),ones(3)/9,'same');
end
imagesc(diff(H))
axis off
colormap(gray)