绘制图像的三维灰度曲面

im4=imread(img_name);   
[m4,n4] = size(im4);
[XX, YY] = meshgrid(1:1:n4,1:1:m4);
zz = im4;
figure,
surf(double(XX),double(YY),double(zz));title('3D_gray')  
xlabel('x');ylabel('y');
shading flat    %  interp   选择interp  或者flat  或者faceted 绘制出的图像的平滑程度不一样，如果把surf改为mesh，则绘制的是网格图像
axis ij   %设置绘制时的原点，ij为左上角，xy为左下角