为了将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