InSAR形变结果可视化1(二维作图点与线)

最新推荐文章于 2023-11-02 18:44:38 发布

SARClub

最新推荐文章于 2023-11-02 18:44:38 发布

阅读量3.7k

点赞数 11

文章标签： matlab 开发语言图像处理算法

本文链接：https://blog.csdn.net/qq_41341717/article/details/128564768

版权

SARculb——小孔径，画山河。

一、写在前面的话：

有一分热，发一分光，不必等候炬火（鲁迅）。

对于InSARer来说形变结果的可视化效果在很大程度上影响着实验的表达，在实际的处理中往往会碰见如下的困境：
在这里插入图片描述

当然这不仅与所使用的成像技术相关，也与处理策略，可视化方式相关密切。随着InSAR技术越来越受到更加广泛的关注，我国在科技领域中扮演的角色越来越重要，行业需求对InSARer的技术要求也越来越高，这无疑需要每一位InSARer付出更多的努力，不仅仅是InSAR结果的可视化，包括国产SAR卫星的技术迭代、SAR数据的公益化以及SAR影像处理软件的真正国产化。

二、InSAR形变结果二维作图工具包（Matlab）

链接中为工具包的压缩包：https://pan.baidu.com/s/1jHHyq4pYrvVch0nu9l4WEQ
提取码：9mz6

压缩包中共包含两个文件分别是示例数据（Data文件夹）和工具包（APP_package文件夹）。工具包的安装直接打开Matlab点击APP/安装APP/选中APP_package文件夹下的安装文件即可完成安装。安装完成后如下所示：
在这里插入图片描述

点击InSAR_plot(2D)，出现工具包的主界面，如下图。点击Open_xx即可打开数据选择窗口。

在这里插入图片描述

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

1、
本工具包暂时不能支持txt，excel，shp，tif等数据格式，只能支持MatLab的.mat格式，读者可以自行百度进行格式转换（如函数load、importdata等）。

lon：形变点经度，格式须转为matlab能够识别的.mat，类型为double。
lat：形变点纬度，格式同上。
vel：形变点形变速率，格式同上。
mm：形变点的时序沉降矩阵，格式同上。mm的横向代表日期，纵向代表不同的形变点。如SARscape、Mintpy、StamPS等结果格式。
date：影像的日期，格式为matlab能够识别的.mat，类型为cell。（如果是文本文件、excel或者其他的话可以百度一下将日期转为cell类型）。
point：需要作图的形变点个数，下拉选择。
style：作图的种类，一共四类，将在下一节叙述，需要注意的是如果选择了3或者4类型，point必须选择2。

2、
上述的数据输入操作为点击Open_xx即可打开数据选择窗口，当数据导入完成时MatLab的命令行窗口会有导入成功的提示并显示导入数据的维度，如下所示：
在这里插入图片描述

3、
参数输入、设置全部完成后点击Plot按钮，便开始运行了，首先会弹出选点交互界面如下图。此图为示例数据某采矿地区的形变速率图，从形变图中选取形变点，point参数是几就点几下，运行结束后选点窗口会出现黑色的三角形即为所选的点，等选点结束后程序会直接出图。
在这里插入图片描述

四、style参数说明

此部分主要写一下style参数，style一共具有4个参数。相关的作图效果如下：

1、style=1
style=1时为单点时序沉降图，一个点一张图，point可以设置为1~5，点击plot后，会画出每个点的历史沉降量（所选点的周围点平均沉降量）。如下图为style=1，point=3时出的图：
在这里插入图片描述

2、style=2
style=2时为多点的历史沉降量，在一张图上，point可以设置为1~5，点击plot后，会画出每个点的历史沉降量在一张图上（所选点的周围点平均沉降量）。如下图为style=2，point=5时出的图：
在这里插入图片描述

3、style=3
style=3时为形变速率的剖面线图，point只能设置为2，点击plot后，会画出线附近每个点的沉降速率，会出两张图，两张图的含义一样，只不过其中一个是加了颜色。如下图为style=3，point=2时出的图：
在这里插入图片描述

4、style=4
style=4时为历史形变量的剖面线图，point只能设置为2，点击plot后，会画出线附近每个点的历史沉降曲线。该曲线为非线性拟合的曲线，需要注意的是该拟合函数存在一定的局限性，对读者的数据点可能会存在较大的误差。如下图为style=4，point=2时出的图：
在这里插入图片描述

五、一些注意事项

1、MatLab的命令行窗口会显示选中点的经纬度。
2、MatLab命令行窗口弹出的信息为每条曲线拟合的标准差等。
3、如果颜色不喜欢的话可以在figure窗口下进行自定义编辑。也可以把拟合的直线删除、更改标记、编辑图例等等。
4、MatLab版本必须大于2016，尽量大于2018。
5、每次点击plot之前要关闭之前出的图包括选点窗口，不用关闭工具包窗口。
5、作图可能会出现错误，如果出现作图错误的情况重新选点出图，不要选取空白无点的地方，多试几次即可；如果还是不行的话，关闭工具包窗口重新运行，如果还是还是不行的话，那应该是不太可能的。
6、InSAR形变结果的3/4D可视化设计思路较为复杂，但已基本完成，希望能在年前发布InSAR_Plot(3/4D)工具包。