InSAR形变结果可视化2(三维作图点线面)

SARculb——小孔径，画山河。

一、InSAR_Plot(3D)工具包(MatLab)：

1、
链接中为工具包的压缩包：链接：链接:https://pan.baidu.com/s/1tPuQLEx22ZcT6zeop3L5UA 
提取码:st8v

下载完成后会有两个文件夹，分别是示例数据（两个）与软件安装包。

安装过程较为简单，见上篇博客。

2、
点击InSAR_Plot(3D)出现工具包主界面，如下图所示，根据读者相关建议对Plot按键进行了一定更新，如下：

二、工具包的输入参数说明

DEM： 研究区的DEM，要尽量比形变点的范围大，格式为tif。注意：采用ArcGIS等裁剪时要按最大矩形框裁切，否则一些NAN值会自动赋值为1或者32768等值，会后面报错。也就是说要保证裁剪后的DEM跟原来的像素值一样，比如Arcgis的裁剪要注意如下选项：

Def： 研究区的形变速率或者累积形变量，格式为tif。注意：不管是DEM还是Def文件，其坐标系要统一，或者说两个文件的像素坐标都是经纬度（或者都是投影坐标）。举一个反例：如果读者采用了Hyp3的干涉图结果+SRTM/ALOS30m那么可能会产生错误，因为Hyp3的干涉图已经进行了投影处理，与DEM的坐标系不一样，因此如果要用Hyp3的形变图那就必须把DEM进行相同的投影转换。一般情况下是不用考虑这个问题的，请读者大胆尝试。

Thre： 该值为经纬度阈值，为了美化成图效果，使得形变点的高程更加贴合DEM，形变点的高程值采用的是DEM的高程，而不是时序算法反演的形变点高程。该阈值为通过形变点经纬度寻找在DEM中相同经纬度下的对应高程值。如果报错的话，下拉选择第二个阈值试试，一般情况下是不用考虑这个问题的，请读者大胆尝试。

三、作图效果说明

1、
上述的数据输入操作为点击Open_xx即可打开数据选择窗口，当数据导入完成时MatLab的命令行窗口会有导入成功的提示，如下所示：

2、
导入数据完成后点击Plot程序便会运行，两个示例数据的作图效果如下所示：

3、

注意： 请读者忽视命令窗口的报错!

作图成功后需要对图片进行相关修改，比如colorbar的标签“mm/y or mm” ，形变点的大小，也可以更改透明度、色带、视角等等。如下图为第二张图修改、放大之后的效果：

四、测试中的进阶功能

因为年关将近，时间紧迫，为尽早发布InSAR_Plot3/4D工具包，一些进阶的作图功能还未封装，且报错较多，正在修改当中，部分作图效果如下所示，后续将对InSAR_Plot3/4D工具包进行更新。在此，也提前敬祝各位读者新春快乐，万事如意！