system: [SYSTEM]

user: 2. 研究区概况及数据源 2.1. 研究区概况 研究区位于高海拔地区，属于高原大陆性气候，冬寒夏凉，雨热同季。既包含永久冻土也包含季节性冻土，区域内广泛发育着大面积的中、高纬度连续多年冻土，具有低温、高厚度的特点，对气温的变化非常敏感。冻土内部构造从上至下依次为季节冻土活动层与多年不融冻土层，地表形变受冻土变化的影响，热融湖、冻胀丘、石海、石河、斑土、冻胀草环等特征地貌形成。 2.2. 数据源 2.2.1. Sentinel-1数据 由于该地区植被较少，积雪相对较薄，相干性较好，很适合通过干涉合成孔径雷达(InSAR)去分析冻土的形变过程。Sentinel-1能够在较短的时间内获取长时间跨度的数据，这是长期挖掘缓慢变化的冻土的重要保证，极大地解决了冻土变形监测中雷达图像积累的瓶颈和困难，进一步降低了时间去相关的风险，使得产生具有良好空间覆盖的地表位移模型成为可能，使得其在长时间序列、短时间间隔的冻土变形监测中具有广阔的应用前景 [ 12 ]。Sentinel-1是一个全天时、全天候雷达成像系统，重访周期12天。Sentinel-1 基于C波段的成像系统采用4种成像模式来观测，工作模式如表1所示。可精确确定卫星位置和姿态角，可选用的轨道参数文件如表2所示，Sentinel-1具体的产品数据类型见表3。 本文研究的盐湖区域选用干涉宽幅模式的斜距单视复数产品，即通过IW模式获取Level-1的SLC产品来进行形变监测分析。本文选用的轨道文件为AUX_POEORB精确的轨道星历参数，分析得到的结果更加准确，该文件是在数据获取后21天内生产得到的，精度5 cm以内。通过获取121景Sentinel-1升轨影像数据对盐湖一带区域进行冻土形变研究，影像参数见表4所示。 Table 1 Sentinel-1 working mod 工作模式 扫描幅宽 空间分辨率 极化方式 SM (条带测绘模式) 80 km 5 m × 5 m HH-HV或VV-VH IW (干涉测量宽视场模式) 250 km 5 m × 20 m HH-HV或VV-VH EW (超宽视场模式) 400 km 25 m × 100 m HH-HV或VV-VH WV (波模式) 采样区域20 km × 20 km，采样区域间隔100 km 5 m × 20 m HH或VV 表1. 哨兵-1工作模式 Table 2 Track fil 轨道文件 获取时间 精度 AUX_POEORB (精确轨道星历参数) 影像数据可获取21天后 5 cm AUX_RESORB (修正轨道参数) 影像数据可获取3小时后 10 cm 表2. 轨道文件 Table 3 Product data typ 数据类型 成像模式 产品类型 SM SLC 条带模式 斜距单视复数产品 IW SLC 干涉宽幅模式 斜距单视复数产品 EW SLC 超宽幅模式 斜距单视复数产品 SM GRD 条带模式 地距多视产品 IW GRD 干涉宽幅模式 地距多视产品 EW GRD 超宽幅模式 地距多视产品 表3. 产品数据类型 Table 4 Image data parameter 日期 (Date) 光束模式 (Beam Mode) 极化方式 (Polarization) 轨道 (Path) 景 (Frame) 2014年10月22日 IW VV 70 110 2014年11月15日 IW VV 70 110 2014年12月9日 IW VV 70 110 2015年1月2日 IW VV 70 110 2015年1月26日 IW VV 70 110 2015年2月19日 IW VV 70 110 2015年3月3日 IW VV 70 110 2015年3月27日 IW VV 70 110 2015年5月14日 IW VV 70 110 2015年6月7日 IW VV 70 110 2015年7月1日 IW VV 70 112 2015年7月25日 IW VV 70 112 2015年8月18日 IW VV 70 112 2015年9月11日 IW VV 70 112 2015年10月5日 IW VV 70 110 2015年10月5日 IW VV 70 115 2015年10月29日 IW VV + VH 70 110 2015年10月29日 IW VV + VH 70 115 2015年11月22日 IW VV 70 109 2015年11月22日 IW VV 70 114 2015年12月16日 IW VV 70 109 2015年12月16日 IW VV 70 114 2016年1月9日 IW VV 70 109 2016年1月9日 IW VV 70 114 2016年2月2日 IW VV 70 109 2016年2月2日 IW VV 70 114 2016年2月26日 IW VV 70 109 2016年2月26日 IW VV 70 114 2016年3月21日 IW VV 70 109 2016年3月21日 IW VV 70 114 2016年4月14日 IW VV 70 109 2016年4月14日 IW VV 70 114 2016年5月8日 IW VV 70 109 2016年5月8日 IW VV 70 114 2016年6月1日 IW VV 70 109 2016年6月1日 IW VV 70 114 2016年7月19日 IW VV 70 109 2016年7月19日 IW VV 70 114 2016年8月12日 IW VV 70 109 2016年8月12日 IW VV 70 114 2016年9月5日 IW VV 70 109 2016年9月5日 IW VV 70 114 2016年9月29日 IW VV 70 109 2016年9月29日 IW VV 70 114 2016年10月23日 IW VV 70 109 2016年10月23日 IW VV 70 114 2016年11月16日 IW VV 70 109 2016年11月16日 IW VV 70 114 2016年12月10日 IW VV 70 109 2016年12月10日 IW VV 70 114 2017年1月3日 IW VV 70 109 2017年1月3日 IW VV 70 114 2017年1月27日 IW VV 70 109 2017年1月27日 IW VV 70 114 2017年2月8日 IW VV + VH 70 109 2017年2月8日 IW VV + VH 70 114 2017年2月20日 IW VV + VH 70 109 2017年2月20日 IW VV + VH 70 114 2017年3月4日 IW VV + VH 70 109 2017年3月4日 IW VV + VH 70 114 2017年3月16日 IW VV + VH 70 113 2017年3月28日 IW VV + VH 70 113 2017年4月9日 IW VV + VH 70 113 2017年4月21日 IW VV + VH 70 113 2017年5月3日 IW VV + VH 70 113 2017年5月15日 IW VV + VH 70 113 2017年5月27日 IW VV + VH 70 113 2017年7月2日 IW VV + VH 70 113 2017年7月14日 IW VV + VH 70 113 2017年7月26日 IW VV + VH 70 113 2017年8月7日 IW VV + VH 70 113 2017年8月19日 IW VV + VH 70 113 2017年8月31日 IW VV + VH 70 113 2017年9月12日 IW VV + VH 70 113 2017年9月24日 IW VV + VH 70 113 2017年10月6日 IW VV + VH 70 113 2017年10月18日 IW VV + VH 70 1296 2017年10月30日 IW VV + VH 70 1296 2017年11月11日 IW VV + VH 70 1296 2017年11月23日 IW VV + VH 70 1296 2017年12月5日 IW VV + VH 70 1296 2017年12月17日 IW VV + VH 70 1296 2017年12月29日 IW VV + VH 70 1296 2018年1月10日 IW VV + VH 70 1296 2018年1月22日 IW VV + VH 70 1296 2018年2月3日 IW VV + VH 70 1296 2018年2月15日 IW VV + VH 70 1296 2018年2月27日 IW VV + VH 70 1296 2018年3月11日 IW VV + VH 70 1296 2018年3月23日 IW VV + VH 70 1296 2018年4月4日 IW VV + VH 70 1296 2018年4月16日 IW VV + VH 70 1296 2018年4月28日 IW VV + VH 70 1296 2018年5月10日 IW VV + VH 70 1296 2018年5月22日 IW VV + VH 70 1296 2018年6月3日 IW VV + VH 70 1296 2018年6月15日 IW VV + VH 70 1296 2018年6月27日 IW VV + VH 70 1296 2018年7月9日 IW VV + VH 70 1296 2018年7月21日 IW VV + VH 70 1296 2018年8月2日 IW VV + VH 70 1296 2018年8月14日 IW VV + VH 70 1296 2018年9月7日 IW VV + VH 70 1296 2018年9月19日 IW VV + VH 70 1296 2018年10月1日 IW VV + VH 70 1296 2018年10月13日 IW VV + VH 70 1296 2018年10月25日 IW VV + VH 70 1296 2018年11月6日 IW VV + VH 70 1296 2018年11月18日 IW VV + VH 70 1296 2018年11月30日 IW VV + VH 70 1296 2018年12月24日 IW VV + VH 70 1296 2019年1月5日 IW VV + VH 70 1296 2019年1月17日 IW VV + VH 70 1296 2019年1月29日 IW VV + VH 70 1296 2019年2月22日 IW VV + VH 70 1296 2019年3月30日 IW VV + VH 70 1296 2019年4月11日 IW VV + VH 70 1296 2019年4月23日 IW VV + VH 70 1296 2019年6月22日 IW VV + VH 70 1296 2019年7月28日 IW VV + VH 70 1296 2019年8月21日 IW VV + VH 70 1296 表4. 影像数据参数 2.2.2. SRTM DEM数据 SRTM数据能获取60˚N~60˚S覆盖全球面积80%以上的三维地表地形信息。SRTM数据包括SRTM1和SRTM3两种，两种数据的空间分辨率分别为30米和90米，本文采用的是30 m空间分辨率的SRTM1数据。SRTM的主要技术参数如表5所示，本文实验区域为青藏高原玉树市治多县附近，地理覆盖范围为(34~38N，92~96E)，下载对应的SRTM1数据，共获得25个geotiff格式的DEM文件，影像数据三维显示如图2所示。 Table 5 Main technical parameters of SRT 轨道高度 233 km 扫描带宽 C波段225 km X波段50 km 雷达波长 C波段5.6 cm X波段3 cm 极化方式 HH， VV， HV和VH 覆盖范围 56˚S~60˚N之间 空间分辨率 SRTM1 30 m × 30 m SRTM3 90 m × 90 m 精度 C波段：水平和高程精度分别为20 m和10 m； X波段：水平和高程精度分别为20 m和4 m 表5. SRTM主要技术参数 图2. 研究区域SRTM数据三维显示 2.2.3. 轨道数据 利用SBAS-InSAR技术获取形变信息的过程中需要利用轨道姿态参数，由于雷达影像元数据中所包含的轨道姿态参数量较少，精度并不高，需要获取精度更高的轨道姿态量，通过结合外部卫星轨道数据数据来解决，对于Sentinel卫星轨道数据可在EARTHDATA平台上获取，本次研究共获取98个AUX_POEORB (精确轨道星历参数)轨道参数文件。