%x1 = disp(:);
s = size(recon);
x1 = 1:s(1);
y1 = 1:s(2);

size(x1)
size(y1)
size(disp)

%smth = medfilt2(smth,[5,5]);

surf(x1,y1,recon,'EdgeColor','none')
%colormap(gray);
view(0,-90)
%rotate(img,[0,90],90)
%viewmtx(0,90,60)
axis off
%axis equal
axis square
%axis equal
axis vis3d
%axis image
axis fill
shading interp
%hold off
%pause;