function H = view3Dmask(I,mask);
  
[m,n]=size(I);
%cols = n:-1:1;
cols=1:n;
%rows=m:-1:1;
rows = 1:m;
view(-45+180,30);
m = (mask*128+128);
m(1,1) = 1;
surf(0.4*double(I(rows,cols)),m(rows,cols));
cm = [gray(128); copper(128)];
colormap(cm);
light;
%shading interp;
shading flat
%lighting gouraud;
lighting phong
material dull
axis equal;
