- /
-
Seabed forms
on 18 Oct 2021
- 8
- 59
- 1
- 0
- 280
N=400;
H=rand(N);
[y,x]=meshgrid(1:N);
for n=1:999
d=H-circshift(H,1);
g=(1+erf(d))*0.5;
L=100.*(1-erf(d));
m=mod(round(x+randn(1)*7)-1,N)+1;
o=mod(round(y+L)-1,N)+1;
ig=accumarray([m(:),o(:)],g(:),[N,N]);
H=H-g+ig;
H=conv2(H,ones(3)/9,'same');
end
imshow(diff(H(25:N, 25:N)'))
colormap(gray)
caxis([-3,4])