将InSAR处理结果绘制为带有等高线的网格图

为了将InSAR处理的结果更加形象的显示，可以利用Matlab中自带的一些函数和方法对InSAR处理的结果图像进行绘制，绘制成为带有等高线的网格图

a=imread('E:\vel_los_utm.tiff');%%读取tiff影像数据，这里只能读取到值，其他的坐标等信息需要后面人为添加
b=a(6001:6500,6001:6500);%%选择绘制的区域，这里需要注意使用Matlab进行绘图时，数据量不能太大，我的电脑只能绘制1000*1000的格网，所以选择小一些的范围
b=smoothdata(b,'gaussian',80);%%这里由于数据有些许离散，继续高斯平滑，窗口范围为80，效果较好，也可以采用较小的数值
c=b(1:300,1:300);%%这里选择了一个更小的范围，以便于画图
meshc(c)%%绘制带有等高线的网格图，这里只是用一个变量，所以绘制结果的x，y范围为0-300，若需要显示经纬度则需定义与数据量相同的x，y向量，使用meshc(x,y,c)绘图
caxis([-0.025,0.015]);%%调整一下colorbar的范围

绘图完成

旋转查看

for i=-1:-1:-89    %%90°旋转查看图像
    %view(a,b):a是角度，b是仰视角
    view(i,18);
    pause(0.09);
end