微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

4. 相位解缠

通过snaphu进行相位解缠。在虚拟机上新建一个文件夹,存放SNAP软件导出的数据,用于相位解缠。如图,单击右键,选择Create Folder。命名文件名为SNAPHU,并将第上一步中SnaphuExport导出的复制到这个文件夹中

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

在虚拟机中的SNAPHU文件夹下,点击鼠标右键,选择Open Terminal Here选项,如图所示。这个在linux系统上经常用到,叫Terminal(终端)。打开snaphu.conf文本文件,即运行snaphu软件进行解缠的命令,拷贝、粘贴到终端命令行中。粘贴好之后,点回车键,即运行该解缠命令,调用snaphu软件进行解缠。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

解缠的结果为两个文件,如图。将这两个文件拷贝出虚拟机,粘贴回上述数据处理文件夹中。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

将解缠结果重新导入SNAP软件中。如图,选择Radar->Interferometric->Unwrapping->Snaphu Import。在弹出的窗口中,1-Read-Phase选择缠绕(未解缠)的数据,在2-Read-Unwrapped-Phase选择解缠的数据,如图所示,选择解缠后的数据头文件(.hdr)。在3-SnaphuImport中选择勾选Do NOT save Wrapped Interferogram in the target product在4-Write把输出文件名加入_Unw,以表示该结果为解缠结果

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

点击Run运行,得到解缠结果如图所示

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

5. 相位转形变

此时,将干涉相位转换为形变图。 点击Radar - >Interferometric - >Products- >Phase to Displacement 设置如下参数然后点击Run运行。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

在Product Explorer中点击生成的形变图,点击displacement在视图窗口打开形变产品。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

6. 地理编码      

把雷达坐标系下的形变量转化为地理坐标系选择 Radar->Geometric->Terrain Correction->Range Doppler Terrain Correction。 

在Range Doppler Terrain Correction列表中,输入参数选择经过相位转为形变的数据,DEM选取准备好的External DEM(外部DEM)。参照图中所示。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

7.输出KMZ格式(Google Earth)

把形变结果导出到GoogleEarth。如图,在打开经过地理编码后的形变结果窗口中,选择Export View as Google Earth KMZ。生成的KMZ文件即可在Google Earth中打开。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码

打开Google Earth,点击File->Open,打开保存的KMZ文件,即可查看。

微波遥感SNAP(四)——检测地表沉降(2)相位解缠与地理编码