利用Sentinel-1A数据反演2016年台湾地震断层的同震滑动分布（一）

DInSAR技术获取的同震形变场，可为地震的受灾范围、地震机理研究等方面提供重要的信息，特别是同震形变数据作为约束条件，基于Okada弹性半空间位错模型反演断层的位置、走向、倾向等几何参数和断层滑动分布，对于深入理解断层活动特征和防震减灾具有重要的意义。
Okada模型是构建地下断层参数与地面形变数据之间的函数关系，主要是模拟观测到的干涉形变场和估计断层参数。反演断层的同震滑动分布主要步骤包括InSAR观测形变场降采样、非线性反演断层几何参数（经度、纬度、走向、倾向、滑动角、深度、断层的长度、宽度以及滑动量）、线性反演断层面上的精细滑动分布。因此，本文将分三个篇幅，分别讲解在SARscape中反演台湾地震断层的滑动分布。所有功能都在以下路径下操作（/SARscape/Interferometry/Displacement Modeling），本篇幅将要讲解InSAR观测形变场降采样。

1 准备工作
模型反演是在DInSAR技术结果的基础上进行的，所以需要对研究区进行DInSAR处理，在此不做赘述。需要DInSAR的最终结果包括：
Taiwan_S1A__reflat_disp
——传感器观测方向的形变，即LOS方向上的形变。
Taiwan _S1A__reflat_disp_ALOS
——视线方位角，正值是正北的顺时针方向，负值是正北的逆时针方向。
Taiwan _S1A__reflat_disp_ILOS
——视线入射角，视线与水平面垂线的夹角。
Taiwan _S1A__reflat_disp_dem
——重采样到制图输出分辨率上的参考DEM数据，经过地理编码，范围和输出的SAR产品一致。

2 反演区选取
创建矢量图层圈定形变场和模型反演范围。形变场采样间隔为1000米，模型反演范围采样间隔为2000米。

形变图进行linear拉伸显示，创建矢量图层圈定形变场（小区域）和模型反演范围（大区域）

编辑矢量文件的属性，分别设置大区域格网2000米，小区域格网1000米



3 InSAR观测形变场降采样
打开降采样工具/SARscape/Interferometry/Displacement Modeling/Image Subsampling，输入Subsampling.shp，默认输出工程区的shp文件


得到降采样矢量文件，可以查看采样值属性。

编辑矢量文件的属性，对形变值进行彩色渲染

得到采样结果：


注意事项：
DInSAR处理中应使用UTM投影下的DEM，即DInSAR所得产品结果，都是在UTM投影下。

为了便于操作练习附上数据下载路径：
链接：http://pan.baidu.com/s/1nuXC5CH 密码：tyxj
文章来源：http://blog.sina.com.cn/s/blog_764b1e9d0102xl2s.html