论文部分内容阅读
近年来,城市的地表沉降问题逐渐引起人们的关注。在对地表形变的测量技术上,除传统的GPS测量计算书和水准测量技术外,还可应用遥感领域的合成孔径雷达(Synthetic Aperture Radar,SAR)进行测量。SAR工作在微波波段,能够在短时间内获取到大面积的地物散射特性,如相关的频率、相位、幅度以及极化等信息。差分雷达干涉(Differential Interferometric SAR,D-InSAR)测量技术是在合成孔径雷达干涉(Synthetic Aperture Radar Interferometry,InSAR)测量技术上发展起来的地表形变测量技术。利用InSAR技术可以得到数字地面高程模型(Digital Elevation Model,DEM)。D-InSAR技术利用DEM数据来消除干涉图中由地形因素带来的影响,进而实现对地表微小形变的检测。Sentinel-1A卫星由欧空局(European Space Agency,ESA)于2014年4月发射,其中心波长为0.56 cm。研究区域为成都市中心城区,中心位置103°55’48"E,30°40’48"N。本文对Sentinel-1A卫星干涉宽幅(Interferometric Wide swath,IW)模式数据,使用双轨法D-InSAR技术对2017年2月12日至2017年12月21日期间的33景影像数据进行处理。研究内容有:(1)对升轨与降轨影像数据的处理效果进行比较;(2)对于Sentinel-1A数据VV极化与VH极化的处理效果进行比较;(3)对于研究区域比较不同的滤波方法、解缠方法的处理结果;(4)分析降水量和地震对研究区域地表形变的影响。针对研究区域的研究结果有:(1)升轨影像对于研究区域的覆盖情况较降轨影像好,处理效果上也稍好于降轨影像;(2)通过比较VV极化与VH极化的干涉处理效果,分析两种极化方式干涉图的相干系数,VV极化干涉处理结果的效果更好;(3)在滤波方法上,比较了Adaptive、Boxcar、固定幂指数Goldstein、利用相干系数调整幂指数的Goldstein四种滤波方法。通过比较噪声去除和边界条纹效果,最后一种滤波方法的处理效果最好。在解缠方法上,比较区域增长法和最小费用流法的解缠孤岛情况,最小费用流法的处理效果更好;(4)研究区域自2017年2月12日至2017年12月21日期间的地表总体形变量在–10~10 mm,其周边的总体形变量在–30~30 mm。形变结果与GPS测量数据及历年数据基本一致;(5)研究区域及其周边的降水量增加会引起研究区域的地表提升。因研究区域位于西高东低的成都平原,地下水含量有自西向东递减的特点。受到西侧的降水量影响较大,且影响时间较长。在2017年2月12日至12月21日期间,研究区域地表形变的范围在0~40 mm,可能为降水量增加影响导致;(6)受到周边地区地震次数的增多及震级提高的影响,研究区域的地表会产生提升。这种提升效果明显但具有短时性,在长时间段的体现微弱。在降水量少的时期抬该升量相对明显,在降水量多且集中的长时间时段内,其影响效果比降水量造成的影响弱。在2017年2月12日至12月21日期间,研究区域地表形变量的范围在0~30 mm,可能为地震增多及震级提高影响导致。