作品名称：基于InSAR技术的海岸带城市沉降监测
作者单位：厦门理工学院计算机与信息工程学院
小组成员：张兰鑫、王棋、郑璐瑶、卢梅英
指导老师：李晖、满旺
获奖情况：三等奖
作品概述
为了克服传统差分干涉雷达在微小形变监测中受时空去相干以及大气效应等负面因素的影响，本项目应用永久散射体( Permanent Scatters，PS)。以厦门海岸带地区为例，选取2015至2018年间Sentinel 1卫星所获得的21幅SAR影像，进行ps反演，得到时间跨度久，形变信息丰富的检测结果。
本项目的工作共包括两个部分，第一部分是利用ENVI/SARscape平台SAR数据集。进行PS-InSAR处理。分别成对组合和做差分干涉，并产生相位缠绕干涉图；第二部分是对形变结果进行分析和实地验证。
作品制作流程
1.本研究拟采用欧洲空间局（ESA）发射的S1A（Sentinel-1A）卫星数据分析厦门市区的地表形变；
2.对哨兵数据和dem数据进行感兴趣区域的裁剪；
3.采用PS分别对SAR数据集进行地表形变反演；
4.得到地表形变时间序列后，确定沉降漏斗；
5.实地调查，验证形变结果。
图 作品制作流程
2.1 数据预处理
（1）哨兵数据导入
在input file 面板输入解压后的哨兵文件的safe文件。
再输入影像对应日期的eof精轨文件。
（2）数据裁剪
对导入后的哨兵数据进行裁剪，选择要进行ps处理的区域
对VV极化方式进行裁剪，所以最后的结果是VV极化方式的数据，同样包含pwr强度文件和slc相位信息的文件。
2.2 PS
PS 技术用于分析点目标，其结果与线性形变相关。输入图像的数量影响像素相干性的估算，像素相干性决定 PS 点的质量。因此推荐使用 20 景以上的数据参与计算，并且接收时 间要连续，最佳情况是一个月一景。PS 技术适合用于城市区域，或者干涉条件和辐射比较 稳定的区域。PS 技术可以探测到毫米级的位移，并推断在时间段上的变形速度。
图 PS流程图
（1）连接图生成
选择一个超级主影像，根据设定的基线阈值对所有图像建立主-从关系，所有的数据都 将配准到一个主图像，这个主图像可自动分配也可以用户自己选择。
图ps connection graph
（2）干涉工作流
这一步是根据像对的连接关系，对每一对像对进行干涉工作流处理。包括数据数据对配准和干涉图生成和干涉图取平。
（3） PS 反演Step1
第一次生成的形变结果，包括速率和高层改正值，这个结果没有去除大气影响的相位。结果用来对合成的干涉图进行取平。
主要方法是基于识别一定数量“相干雷达信号反射”（永久散射体），然后重点分析这些 可靠的单一目标（一个图像像素一个目标）的历史相位。永久散射体主要受两个限制：它必 须是稳定的(变动小于一毫米) 和它面向一个能被 SAR 天线探测到的方向。
（4）地理编码
所有 PS 相关结果（形变速率、高度残差、形变序列、KML、矢量文件等）投影到地图系统中。
图 ps geocoding（1）
在 Input Files 面板中，选择auxiliary.sml 文件。单击Optional File 面板，这是可选项，选择 GCP 文件（.shp 或者.xml）或者创建 GCP 文件。本操作不选择
图ps geocoding（2）
在 DEM/Cartographic System 面板中，选择srtm_dem 文件，参数保持默认，单击 Exec 按钮执行处理。
（5）处理结果
图ps结果
Ps 点的个数很多，所以分成了6个分块矢量文件。Ps点的个数568 065，厦门岛的面积158万平方公里。
2.3 结果分析与验证
在ArcGIS种对PS点进行克里金差值，得到结果如下：
由于在思明区分布较为集中的山体，山体地区的ps稀少，对该地区进行了掩模处理。
结果验证:
图 双子塔
图 财富中心
图 海湾公园
图 高崎码头
作品亮点
PS（永久散射体）处理工作量大，作品实现了基于Sentinel-1的厦门市地表沉降监测，并进行了实地验证，作品完成度较好。
