EN中文

自然资源遥感 >

2025 , Vol. 37 >Issue 6: 128 - 137

DOI: https://doi.org/10.6046/zrzyyg.2024346

技术方法

双极化优化的时序InSAR形变监测研究

  • 玄甲斌 , 1 ,
  • 李如仁 1 ,
  • 傅文学 , 2, 3
展开
  • 1.沈阳建筑大学交通与测绘工程学院, 沈阳 110168
  • 2.中国科学院空天信息创新研究院数字地球重点实验室, 北京 100094
  • 3.可持续发展大数据国际研究中心, 北京 100094
傅文学(1977-),男,博士,副研究员,主要从事微波遥感方向研究。Email: 。

玄甲斌(2000-),男,硕士研究生,主要从事合成孔径雷达形变监测研究。Email: 。

Copy editor: 李瑜

收稿日期: 2024-10-19

  修回日期: 2025-01-23

  网络出版日期: 2026-06-03

基金资助

广西科技重大专项“多模态信息驱动的稀土矿区安全要素遥感监测与预警”(桂科AA24206025)

收起

Deformation monitoring using time-series InSAR with dual-polarization optimization

  • XUAN Jiabin , 1 ,
  • LI Ruren 1 ,
  • FU Wenxue , 2, 3
Expand
  • 1. School of Transportation and Geomatics Engineering, Shenyang Jianzhu University, Shenyang 110168, China
  • 2. Key Laboratory of Digital Earth Science, Aerospace Information Research Institute, Chinese Academy of Sciences, Beijing 100094, China
  • 3. International Research Center of Big Data for Sustainable Development Goals, Beijing 100094, China

Received date: 2024-10-19

  Revised date: 2025-01-23

  Online published: 2026-06-03

Fold

摘要

高质量监测点的空间密度及其干涉相位质量是时序合成孔径雷达干涉测量(interferometric synthetic aperture radar,InSAR)开展形变监测的重要指标,为了进一步提高InSAR技术在非城市区域的形变监测能力,使用双极化Sentinel-1数据,提出了一种顾及分布式散射体(distributed scatterer, DS)的极化时序InSAR技术方法。该方法根据分布式散射体的特点,并将振幅离差(dispersion of amplitude, DA)作为相位质量评价指标,使用不同方法分别对时序SAR数据的强度信息和相位信息进行极化优化处理,对所优化前后数据进行地表形变监测。以浙江省宁波市为例,采用40景双极化(VV-VH)Sentinel-1数据进行实验。结果表明,所提方法能够显著提高监测点的密度和干涉相位质量。与单极化相比,永久散射体(persistent scatterer, PS)数量提高了约20%,DS点数量提高了约57.5%,干涉图相位质量有明显提升,平均相干性能够提升15%以上,可以更加详细地反映区域形变状况。

关键词： 极化时序InSAR技术; 双极化Sentinel-1数据; 复杂地形的形变监测

本文引用格式

玄甲斌 , 李如仁 , 傅文学 . 双极化优化的时序InSAR形变监测研究[J]. 自然资源遥感, 2025 , 37(6) : 128 -137 . DOI: 10.6046/zrzyyg.2024346

Abstract

The spatial density and interferometric phase quality of high-quality monitoring points serve as key indicators for deformation monitoring using the time-series interferometric synthetic aperture radar (InSAR) technique. To further enhance the deformation monitoring ability of the InSAR technique for non-urban areas, this study proposed a polarization time-series InSAR method that takes into account distributed scatterers (DSs) using dual-polarization images from Sentinel-1. Specifically, polarization processing of the intensity and phase information of time-series SAR data was conducted using various methods based on the characteristics of DSs and taking the dispersion of amplitude (DA) as an indicator for the phase quality assessment. Then, surface deformation monitoring was performed using the data before and after optimization. This study carried out experiments on Ningbo City in Zhejiang Province using 40 scenes of dual-polarization (VV-VH) images from Sentinel-1. The results indicate that the proposed method can significantly increase the density of monitoring points and the interferometric phase quality. Compared to single polarization, the proposed method increased the quantities of persistent scatterers (PSs) and DSs by about 20% and 57.5%, respectively. Furthermore, the interferometric phase quality was also significantly improved, with the average coherence increasing by more than 15%. The proposed method allows for a more detailed reflection of regional deformations.

Key words： polarization time-series InSAR technique; dual-polarization images from Sentinel-1; deformation monitoring of complex terrain

0 引言

地表形变监测在预防地表沉降、滑坡、地震等相关地质灾害方面发挥着至关重要的作用,可有效降低灾害风险。具体而言,通过对地表形变的持续精准监测,可以及时发现潜在的地质灾害隐患,为灾害预警和应急处置提供科学依据,从而降低灾害发生的可能性和危害程度[1-3]。合成孔径雷达干涉测量(interferometric synthetic aperture radar, InSAR)技术能够全天候、全天时的以低成本、高空间分辨率和高效率监测大范围的地表形变,具有传统形变监测技术无可比拟的优越性[4-7]。随着对InSAR技术研究的不断深入,国内外学者先后提出了一系列的时序InSAR方法[6,8 -10]。其基本原理是通过对长时间序列覆盖同一区域的SAR影像进行处理,提取后向散射特性较强的高相干点目标的形变信息,理论上的监测精度可达毫米级[11-12]。
但是,针对人工目标稀少的非城市区域,由于受到去相关的强烈影响,这些传统的时序InSAR技术在实际应用中常常面临难以提取到数量充足且分布密度较高的高相干点目标这一问题[13]。相位解缠是InSAR技术中的关键环节,相干目标不足会导致解缠过程出现错误,进而影响整个形变监测的有效性和可靠性。与永久散射体(persistent scatterer, PS)相比,分布式散射体(dispersion of amplitude, DS)具有更为广泛的空间分布特点,其通常出现在裸地、荒漠以及稀疏植被覆盖区域等[14-15]。因此,分布式散射体干涉测量技术(DS-InSAR)可以很好地弥补传统时序InSAR技术在非城市区域监测点数量不足的问题,可以为非城市区域的地表形变监测提供更有力的技术支持[16-18]。
另一方面,受极化SAR干涉测量技术(polarimetric InSAR, PolInSAR)[19]的启发,一些研究人员提出利用多极化SAR影像中包含的极化信息对干涉相位进行优化。因此,Pipia等[20]于2009年基于全极化地基SAR数据成功将极化优化算法应用于PSI技术。随着卫星技术的发展,近年来多极化SAR卫星(TerraSAR-X,ALOS-2,Sentinel-1A,GF-3等)陆续发射升空,获取了大量的多极化SAR数据,为极化SAR干涉测量技术提供了最为坚实的数据基础[21-22]。利用多极化SAR数据,国内外学者开展了大量的研究工作,提出了多种基于PSI技术的极化优化技术(polarimetric persistent scatterer interferometry, PolPSI)[23]。研究表明,极化通道的极化信息蕴含着丰富的目标特性。所以,通过利用不同极化通道的极化信息,搜索目标的最佳散射机制(scattering mechanism,SM),可以降低SAR像元因去相关引起的干涉相位噪声,优化散射体相位质量,提高监测点的空间密度和相位质量,极大助力在中低相干区的时序形变监测能力[24-27]。
目前最常用的极化优化方法有最优极化通道(BEST)和遍历搜索极化优化算法(exhaustive search polarimetric optimization, ESPO)等[28]。BEST方法仅在HH,VV和VH极化通道中选择最佳极化通道,该方法简单高效,但并没有充分发挥极化通道的潜力。ESPO则是遍历极化空间的最优解,但计算负担非常高,会限制其在大型场景中的实际应用[29]。通常将相位质量指标振幅离差法(dispersion of amplitude, DA)作为最佳极化通道的判断准则。由于极化优化会导致相邻像元之间的极化方式存在差异,而DS-InSAR中同质像元识别的原理是基于幅度或强度影像以统计推断为手段度量邻域像元与中心像元的相似度进行判断[30-32],因此现有的极化优化方法会降低DS点的选取数量。那么如何利用极化信息,提高DS点的数量,将极化优化算法应用于DS-InSAR技术值得进一步的研究。
目前,学术界提出了多种PolPSI技术,但是由于DS-InSAR中同质像元识别原理的限制,常用的PolPSI并不适用于DS-InSAR技术,并且计算非常耗时。因此,为了更好地提高InSAR技术在复杂地形的形变监测效果,本文在永久散射体干涉测量技术(PS-InSAR)和DS-InSAR技术的基础上,结合极化优化算法及一种基于群体智能的优化算法,提出了一种顾及DS的高效的极化时序InSAR技术方法,该方法解决同质像元识别原理对极化优化的限制,并大幅提高了极化优化效率。基于本文提出的方法,选取浙江省宁波市部分区域为研究区,使用双极化Sentinel-1数据进行实验,并对极化优化前后的相关结果进行讨论。

1 研究区概况及数据源

浙江省宁波市地处宁波平原,地势呈现西南高、东北低的特点。该市地貌主要涵盖山脉、丘陵、盆地和平原等类型,其具体的地形分布情况如图1所示。目前,宁波市以现代化滨海大都市为总体框架,积极建设“全球海洋中心城市”,构建“一核、三湾、六片”发展格局,涵盖对前湾新区、镇海北仑片区、象山港片区和三门湾片区的产业开发。这些地区有不少沿海围垦区,需采取填海造地等方式解决土地资源问题。在未固结的泥沙填土上开展工程建设会导致地面沉降广泛出现[33]。
图1 研究区域地形图

Fig.1 Topographic map of study area

本文使用2022-01-20—2023-12-17间的40景Sentinel-1 IW模式的双极化(VV+VH)数据,首先对数据进行预处理,包括配准、裁剪等操作,得到研究区域的VV,VH这2种极化方式的单视复数(single look complex, SLC)影像,之后基于本研究的方法对宁波市进行时序形变监测。具体研究区为图1所示的2 000×2 500像元的沿海区域,该区域包含了人造设施及大量裸地、植被等自然地表,地形地貌较为复杂,研究区地表覆盖分类及卫星影像如图2所示。
图2 研究区地表覆盖分类及卫星影像图

Fig.2 Land cover classification and satellite image map of the study area

2 极化时序InSAR技术

2.1 时序极化优化

在全极化数据下,对于每个SAR像元,散射矩阵S表示为一个2×2的复矩阵,通常情况下它包含了散射体的完整信息,其目标散射矩阵S[19]如下:
S= $\left[\begin{array}{ll}{S}_{HH}& {S}_{HV}\\ {S}_{VH}& {S}_{VV}\end{array}\right]$,
式中: SHH和SVV为同极化分量; SVH和SHV为交叉极化分量。假设满足互异性,即SHV=SVH,散射矩阵S经Pauli基矢量化后可得极化散射向量k[4],计算公式为:
k= $\frac{1}{\sqrt{2}}{\left[{S}_{HH}+{S}_{VV},{S}_{HH}-{S}_{VV},2{S}_{HV}\right]}^{\mathrm{T}}$,
式中: T表示转置。当全极化向双极化转变时,式(2)可转化为:
kVV/VH= $\frac{1}{\sqrt{2}}{\left[{S}_{VV},2{S}_{HV}\right]}^{\mathrm{T}}$,
kHH/VH= $\frac{1}{\sqrt{2}}{\left[{S}_{HH},2{S}_{HV}\right]}^{\mathrm{T}}$,
kHH/VV= $\frac{1}{\sqrt{2}}{\left[{S}_{HH}+{S}_{VV},{S}_{HH}-{S}_{VV}\right]}^{\mathrm{T}}$。
极化优化是一种通过选择受去相关影响最小的极化状态来提高高质量相干点的数量及其相位质量的方法[23],即通过表征散射机制的复数投影向量ω[25],在极化空间中遍历使像元质量最佳的极化方式μ,计算公式为:

μ=ω†k,

式中: †表示共轭转置。
对于双极化数据[4,29],其投影向量ω为:
$ \omega=\left[\begin{array}{c} \cos \alpha \\ \sin \alpha \cdot \mathrm{e}^{j \vartheta} \end{array}\right] \begin{array}{c} (0 \leqslant \alpha \leqslant \pi / 2) \\ (-\pi \leqslant \vartheta \leqslant \pi) \end{array},$
式中: α,ϑ分别为不同散射机制的Pauli参数。
振幅离差指数(DA)是一种常用的相位质量评价指标[11,15],表示为时间序列上同一像元的振幅A的标准差σA与均值$\stackrel{-}{A}$之比。在多极化条件下,DA表示为:
DA= $\frac{{\sigma }_{A}}{\stackrel{-}{A}}$= $\frac{1}{\overline{\left|{\omega }^{†}k\right|}\sqrt{N}}\sqrt{\stackrel{N}{\sum _{i=1}}{(\left|{\omega }^{†}{k}_{i}\right|-\overline{\left|{\omega }^{†}k\right|})}^{2}}$,
式中: N为影像数量; $\left|\stackrel{-}{·}\right|$为求均值。
因此时序极化优化的主要目的就是在二维空间中搜索使DA最小的投影向量ω,求解参数α和ϑ。需要指出的是,由于极化优化需要遍历每个SAR像元的所有极化空间,计算成本较高,会限制其在大型场景中的实际应用。

2.2 双极化时序InSAR形变监测

不同于以往的PolPSI技术,为了进一步提高监测点的密度,本文将DS-InSAR引入时序极化优化中。根据DS点选取受邻域像元影响的问题,本文将时序SAR数据的强度信息和相位信息分别使用不同的极化优化方法,以此得到极化优化后的SLC影像,并使用优化后的影像联合PS点和DS点进行地表形变监测,具体流程如图3所示。
图3 双极化时序InSAR形变监测处理流程

Fig.3 Processing flow of dual-polarization time-series InSAR deformation monitoring

首先,为不影响DS点的选取,根据同质像元识别的原理可知,规定所有SAR像元必须在同一极化状态下,因此针对时序SAR影像的强度信息,本文采用全局优化策略,即最佳全局散射机制(best global scattering mechanism, BGSM)法[34]。该方法的原理是针对整个SAR影像数据对极化空间进行搜索,并通过最小化整个影像的的平均值找到对于全局最优的投影向量ω。该方法的优势在于规定了整幅影像都在同一极化状态下,不会改变邻域像元的相似度,能够提高同质像元识别数量。
其次,为最大程度降低SAR数据的相位噪声,提高相位质量,针对时序SAR影像的相位数据,逐像元进行极化优化处理。ESPO方法[35]对每个像元在整个极化空间上执行二维搜索,通过使每个像元的DA最小找到最优极化状态,图4展示了某SAR像元的极化空间参数搜索结果。但是,该方法的搜索耗时是一个不可否认的缺点。因此,为了减轻负担,本文在ESPO算法的基础上,将粒子群优化(particle sswarm optimization, PSO)算法加入其中,使其能够快速找到近似最优的极化空间,而不必遍历所有可能的情况,以提高算法的运行效率。
图4 SAR像元的极化空间参数搜索结果

Fig.4 Search results of polarization space parameters of SAR pixels

最后基于双极化优化得到的SLC影像,联合PS和DS技术,进行时序地表形变监测,并对监测结果进行讨论与分析。

3 结果与分析

为分析极化优化对PS点和DS点选取数量的提升效果,实验分别对比分析了2种方法在识别PS点的过程中不同阈值下的点数量,以及同质像元识别结果和最终DS点的选点结果。

3.1 算法性能评估

为了验证加入PSO算法对运行效率的提升,分别针对单个像元及研究区整体运行时间及运算结果进行分析。电脑配置为12th Gen Intel(R) Core(TM) i7-12 700单线程运算,对于单一像元[566,1 151],在设置ESPO方法搜索步长为3°情况下,运行时间为0.08 s,最小DA值为0.4 198; 加入PSO算法后,在设置粒子数5,迭代次数20的情况下,同一像元的优化时间为0.003 s,最小DA值为0.4 196,其具体迭代情况如图5所示。对于研究区整体来说,其大小为2 000×2 500像元,在设置参数不变的前提下,ESPO方法预计需要120 h,而PSO算法仅需2.5 h。由此可以看出PSO算法能够极大提高极化优化效率,且求出的最小DA值也要优于ESPO方法。为了更合理地设置PSO算法的参数,以提高运行效率,经过多次实验得出,在设置粒子数5,迭代次数10的情况下,绝大多数像元均能达到收敛状态,因此在接下来的实验中,PSO算法采用该参数设置。
图5 PSO算法收敛性分析

Fig.5 Convergence analysis of the PSO algorithm

3.2 选点结果分析

为了深入分析VV极化与极化优化后这2种不同数据对PS点识别所产生的影响,本实验基于振幅离差(DA)对PS点的选取结果进行分析。通常情况下,相干目标的后向散射性越强,其在时序上的表现就会越稳定,相应的DA值也就越小。所以,那些DA低于特定阈值的点,便可以被选作PS点。通过这样的方式,可以更加准确地识别出具有稳定后向散射特性的PS点,从而为后续的形变监测和分析提供可靠的数据基础。同时,不同的极化方式和极化优化处理会对相干目标的后向散射性产生不同的影响,进而影响到PS点的选取结果。因此,对这2种不同数据下的PS点识别进行分析,有助于更好地理解极化和极化优化在InSAR形变监测中的作用,为进一步提高形变监测的精度和可靠性提供有益的参考。
研究表明,经极化优化后,强度影像的旁瓣效应与噪声均在一定程度上有所减轻。图6呈现了VV,VH,极化优化后(OPT)这3种方式获取的时序平均强度图,从红色方框内放大的子区域能够更清晰地看出,VV极化因人工建筑等强散射体对周围弱散射体回波信号产生严重干扰,致使建筑物边界模糊; VH极化因强弱散射体回波信号均较弱,使得地物难以辨别; 而极化优化处理有力地抑制了旁瓣效应,提升了不同地物间的对比度。
图6 不同极化方式下的平均强度图

Fig.6 Plot of average intensity for different polarisation models

为定性分析极化优化对提高PS点数量的效果,选取研究区中的2种典型区域进行研究,如图7所示,分别为裸地、农田区域(白框1)和工厂厂区(白框2),其中底图为研究区的平均强度图,红色点表示PS点,框1中PS点个数分别为1 033和1 167个,框2中分别为15 078和15 830个。根据图7分析: 针对裸地、农田区域,通过黄色圆圈对比分析可知,极化优化能够提高PS点的数量,这是由于该区域在SAR影像分辨单元内存在大量的散射体导致散射回波信号较弱,通过双极化优化寻找最优散射机制,有利于改善SAR影像分辨单元的散射回波信号,因此极化优化能够作用于该区域。但是对于工厂厂区,极化优化并没有体现出很好的效果,这是因为在人工建筑物较多的区域,散射回波信号强,SAR分辨单元本来就具有较高的稳定性。其次,为定量分析极化优化对提高PS点数量的效果,不同阈值下点数量的统计结果如表1所示。在不同的DA指标范围下极化优化后选取的点数量均多于VV极化,并且优化效果明显,在阈值设为0.4时,PS点数量能够提升约20%。
图7 2种典型区域的PS点数量对比

Fig.7 Comparison of the number of PS points in 2 typical areas

表1 不同阈值下选取的PS点数量

Tab.1 Number of PS points selected under different thresholds

方法 指标范围
0~0.2 0~0.3 0~0.4 0~0.5
VV 1 969 10 797 35 874 92 982
OPT 2 592 13 573 42 976 107 726
(OPT-VV)/VV (↑) 31.6% 25.7% 19.8% 15.9%
为分析VV极化与极化优化2种不同数据对DS点识别的影响,实验分别对同质像元识别结果和最终DS点的选点结果分析讨论。图8分别显示了2种方法得到的同质像元识别数量分布图。结合图2进行分析可以看出,地表覆盖类型会对同质像元的识别产生影响。以城市区域为例,城市的建筑物散射特性较强,导致选取的同质像元数量相对较少。相比之下,裸地或低矮植被覆盖区域等属于典型的分布式目标集中区域,所以在这些地方选取的同质像元数量会比较多。由此可见,不同的地表覆盖类型具有不同的散射特性,进而影响了同质像元的识别结果。为了更直观地对比优化的效果,实验选择红框区域进行对比分析,通过图8对比分析可得: 极化优化能够显著提高同质像元识别的数量,并且能够提高同质像元识别的准确性。如图8中白圈区域所示,(b)中该区域的同质像元数量明显高于(a),并且地物框架则十分明显,说明此时同质像元和异质像元能够被很好地区分。
图8 同质像元识别数量对比

Fig.8 Comparison of the number of identified homogeneous pixels

之后基于同质像元识别结果,进行DS点选择。DS识别分布对比如图9所示,其中底图为研究区的平均强度图,黄色点表示DS点。经统计得出,VV极化下DS点数量为358 869个,而极化优化处理后识别出的DS点数量为565 117,增长了约57.5%。
图9 DS点识别分布对比图

Fig.9 Comparison diagram of DS point identification distribution

3.3 差分干涉图结果分析

选取一幅长时间基线(20230103—20221128)的差分干涉图进行质量评价。分析区域位于研究区东部沿海码头区域。图10(a)(b)分别显示了VV极化差分干涉图和极化优化后的差分干涉图,图10(c)为对应区域的平均强度图,红色虚线显示该码头的范围。通过对比优化前后的差分干涉图能够明显看出,极化优化后的干涉图条纹要比原始干涉图清晰得多。在原始干涉图中,码头内部构筑物的边缘区域干涉相位十分模糊,经过优化之后,能够降低差分干涉图的噪声水平,从而使得构筑物的轮廓更加清晰可见。
图10 局部区域干涉图优化前后对比

Fig.10 Comparison of local area interferograms before and after optimization

为定量分析极化优化对差分干涉图的优化效果,选用定量指标对差分干涉图质量进行评价。基于优化前后的差分干涉图,分别计算其残差点数(residue point number, RPN)、相位差之和(sum of phase differences, SPD)及干涉相干性(coherence, COH)。RPN和SPD越小,表明干涉图的质量越好。采用[15×15]规则窗口分别对优化前后的差分干涉图进行相干性估计,COH越大,对干涉图的优化效果越好。表2显示了在不同时间基线下,极化优化对差分干涉图的优化效果。结果表明: 极化优化对于长时间基线,使干涉图的RPN和SPD分别降低了11.1%和4.9%,使干涉图的COH提高了11.4%; 对于短时间基线,使干涉图的RPN和SPD分别降低了13.5%和6.2%,使干涉图的COH提高了15.2%。
表2 干涉图质量评价结果

Tab.2 Interferogram quality evaluation results

干涉图 长时间基线(20230103—20220201) 短时间基线(20230103—20221128)
RPN SPD COH RPN SPD COH
VV 925 008 6.2e+07 0. 289 872 105 6.5e+07 0. 322
OPT 822 700
(11.1%↓)		 5.9e+07
(4.9%↓)		 0. 322
(11.4%↑)		 754 483
(13.5%↓)		 6.1e+07
(6.2%↓)		 0. 371
(15.2%↑)

3.4 形变监测结果分析

图11为利用VV极化数据和极化优化后数据进行时序InSAR形变监测获得的研究区地表形变结果,方框中为宁波市鄞州区部分区域的地表形变速率结果。对比2种方法获得的整体形变趋势结果与已有研究数据[33,36]的平均沉降速率基本一致,且极化优化后的形变监测结果的监测点的密度更高,说明在该区域极化优化技术能够更加详细地反映区域形变状况。
图11 研究区地表形变速率图

Fig.11 Surface deformation rate map of the study area

从图11的3处红圈区域更加清晰地看出,在极化优化之后,由于形变监测点的密度较大,所以所反演出来的形变细节也更为清晰。此外,在最左侧红圈区域中2种方法的监测形变结果呈现出显著的差异,产生这种差异的主要原因是在时间序列InSAR相位解缠的过程中,相干目标需事先构建空间网络,之后再进行模糊度解算,而空间网络的构建方式会直接影响到后续模糊度解算的准确性,进而对形变监测的结果产生很大的影响。目前Delaunay三角网是当前时序相位解缠的主流构网方法,评价构网质量的主要指标有三角网的密度以及三角网边的相位梯度大小,若网络形态包含高相位梯度边缘,导致违背相位连续性假设,会对形变监测结果产生较大影响[37]。通过前文的分析可知,经极化优化处理后得到的形变监测点的数量及其相位质量都有明显的提升,因此优化后使得相干目标的连接更加紧密,构建的三角网密度更大,效果更好,能够更加准确地捕捉到细微的形变变化,可以更详细地反演真实地表形变,从而提高了地表形变监测能力。

4 结论

为了进一步提高InSAR技术在复杂区域的形变监测能力,针对提高高相干点数量和干涉相位质量进行研究,提出了一种估计DS的极化时序InSAR技术方法,并以浙江省宁波市为例,对比分析单极化方法与本文方法的形变监测效果,结论如下:
1)通过加入PSO算法,大幅提高极化优化的运算效率,比ESPO方法时间减少约48倍。
2)本文所提的方法能很好地保持PS目标分辨率,提高PS点数量,并提升DS点识别的准确性和数量,相较单极化数据,DS点增幅达57.5%。
3)通过对相位数据的优化,能够明显提升差分干涉相位图的质量,干涉平均相干性能够提升15%以上。
4)通过提高监测点密度及其干涉相位质量,可以得到更详细的研究区地表形变监测结果,为提高InSAR技术在复杂区域监测效果的问题提供了一种新的思路和方法。

参考文献

原文顺序 | 文献年度倒序 | 文中引用次数倒序
[1]
Zhu C H, Wang C S, Zhang B C, et al. Differential interferometric synthetic aperture radar data for more accurate earthquake catalogs[J]. Remote Sensing of Environment, 2021,266:112690.

[2]
Cigna F, Tapete D. Present-day land subsidence rates,surface faulting hazard and risk in Mexico City with 2014-2020 Sentinel-1 IW InSAR[J]. Remote Sensing of Environment, 2021,253:112161.

[3]
Liu Y F, Yang H L, Wang S Z, et al. Monitoring and stability analy-sis of the deformation in the woda landslide area in Tibet,China by the DS-InSAR method[J]. Remote Sensing, 2022, 14(3):532.

DOI

[4]
Zhao F, Wang T, Zhang L X, et al. Polarimetric persistent scatterer interferometry for ground deformation monitoring with VV-VH sentinel-1 data[J]. Remote Sensing, 2022, 14(2):309.

DOI

[5]
Blanco-Sánchez P, Mallorquí J J, Duque S, et al. The coherent pixels technique (CPT):An advanced DInSAR technique for nonlinear deformation monitoring[J]. Pure and Applied Geophysics, 2008, 165(6):1167-1193.

DOI

[6]
Ferretti A, Prati C, Rocca F. Permanent scatterers in SAR interferometry[C]// Proceedings of SAR Image Analysis,Modeling,and Techniques II. SPIE-International Society for Optical Engineering, 1999:132-139.

[7]
安炳琪, 罗海滨, 丁海勇, 等. 基于SBAS-InSAR技术的西宁地表形变监测[J]. 遥感技术与应用, 2021, 36(4) :838-846.

DOI

An B Q, Luo H B, Ding H Y, et al. Monitoring of surface deformation in Xining based on SBAS-InSAR[J]. Remote Sensing Technology and Application, 2021, 36(4):838-846.

[8]
Ferretti A, Prati C, Rocca F. Nonlinear subsidence rate estimation using permanent scatterers in differential SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5):2202-2212.

DOI

[9]
Zhang L, Lu Z, Ding X L, et al. Mapping ground surface deformation using temporarily coherent point SAR interferometry:Application to Los Angeles Basin[J]. Remote Sensing of Environment, 2012,117:429-439.

[10]
Fan H D, Gao X X, Yang J K, et al. Monitoring mining subsidence using a combination of phase-stacking and offset-tracking methods[J]. Remote Sensing, 2015, 7(7):9166-9183.

DOI

[11]
何倩. 联合PS和DS的时序InSAR地表沉降监测方法研究[D]. 徐州: 中国矿业大学, 2022.

He Q. Study on time-series InSAR method for monitoring surface subsidence with persistent and distributed scatterers[D]. Xuzhou: China University of Mining and Technology, 2022.

[12]
Huang J L, Bai Y, Lei S G, et al. Time-series SBAS pixel offset tracking method for monitoring three-dimensional deformation in a mining area[J]. IEEE Access, 2020,8:118787-118798.

[13]
张路, 廖明生, 董杰, 等. 基于时间序列InSAR分析的西部山区滑坡灾害隐患早期识别——以四川丹巴为例[J]. 武汉大学学报(信息科学版), 2018, 43(12):2039-2049.

Zhang L, Liao M S, Dong J, et al. Early detection of landslide hazards in mountainous areas of West China using time series SAR interferometry—a case study of Danba,Sichuan[J]. Geomatics and Information Science of Wuhan University, 2018, 43(12):2039-2049.

[14]
Ferretti A, Fumagalli A, Novali F, et al. A new algorithm for processing interferometric data-stacks:SqueeSAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2011, 49(9):3460-3470.

DOI

[15]
赵立峰. 基于DS-InSAR的矿区地表形变监测方法研究——以沛北矿区为例[D]. 徐州: 中国矿业大学, 2022.

Zhao L F. Study on surface deformation monitoring method based on DS-InSAR in mining area[D]. Xuzhou: China University of Mining and Technology, 2022.

[16]
Goel K, Adam N. A distributed scatterer interferometry approach for precision monitoring of known surface deformation phenomena[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(9):5454-5468.

DOI

[17]
Li T T, Zhang H Z, Fan H D, et al. Position inversion of goafs in deep coal seams based on DS-InSAR data and the probability integral methods[J]. Remote Sensing, 2021, 13(15):2898.

DOI

[18]
Chen B Q, Yang J L, Li Z H, et al. A new sequential homogeneous pixel selection algorithm for distributed scatterer InSAR[J]. GIScience & Remote Sensing, 2023, 60(1):31.

[19]
Cloude S R, Papathanassiou K P. Polarimetric SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 1998, 36(5):1551-1565.

DOI

[20]
Pipia L, Fabregas X, Aguasca A, et al. Polarimetric differential SAR interferometry:First results with ground-based measurements[J]. IEEE Geoscience and Remote Sensing Letters, 2009, 6(1):167-171.

DOI

[21]
赵峰, 张雷昕, 王腾, 等. 城市地表形变的双极化Sentinel-1数据极化时序InSAR技术监测[J]. 武汉大学学报(信息科学版), 2022, 47(9):1507-1514.

Zhao F, Zhang L X, Wang T, et al. Polarimetric persistent scatterer interferometry for urban ground deformation monitoring with sentinel-1 dual polarimetric data[J]. Geomatics and Information Science of Wuhan University, 2022, 47(9):1507-1514.

[22]
Iglesias R, Fabregas X, Aguasca A. A new technique for Atmospheric Phase Screen(APS) compensation based on a 2D MultiplernRegression Model in Ground-Based SAR[C]// 9th European Conference on Synthetic Aperture Radar.Nuremberg,Germany:VDE Verlag GMBH, 2012:12-15.

[23]
冯瀚, 赵峰, 汪云甲, 等. 双极化Sentinel-1数据极化时序InSAR技术地表形变监测——以上海市浦东机场为例[J]. 遥感学报, 2022, 26(12):2531-2541.

Feng H, Zhao F, Wang Y J, et al. Dual-polarization Sentinel-1 data polarization time series InSAR technology surface deformation monitoring—Taking Shanghai Pudong airport as an example[J]. National Remote Sensing Bulletin, 2022, 26(12):2531-2541.

DOI

[24]
Navarro-Sanchez V D, Lopez-Sanchez J M, Vicente-Guijalba F. A contribution of polarimetry to satellite differential SAR interferometry:Increasing the number of pixel candidates[J]. IEEE Geoscience and Remote Sensing Letters, 2010, 7(2):276-280.

DOI

[25]
Iglesias R, Monells D, Fabregas X, et al. Phase quality optimization in polarimetric differential SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(5):2875-2888.

DOI

[26]
Iglesias R, Monells D, López-Martínez C, et al. Polarimetric optimization of temporal sublook coherence for DInSAR applications[J]. IEEE Geoscience and Remote Sensing Letters, 2015, 12(1):87-91.

DOI

[27]
赵峰. 面向PSI技术的高质量SAR像元识别与干涉相位极化优化方法研究[J]. 地理与地理信息科学, 2020, 36(3):141.

Zhao F. Advanced pixel selection and optimization algorithms for persistent scatterer interferometry (PSI)[J]. Geography and Geo-Information Science, 2020, 36(3):141.

[28]
Navarro-Sanchez V D, Lopez-Sanchez J M. Improvement of persistent-scatterer interferometry performance by means of a polarimetric optimization[J]. IEEE Geoscience and Remote Sensing Letters, 2012, 9(4):609-613.

DOI

[29]
桂阳归, 汪长城, 卢丽君, 等. 顾及极化最优的时序双极化InSAR形变监测方法[J]. 测绘科学, 2021, 46(12):75-82.

Gui Y G, Wang C C, Lu L J, et al. A time-series dual-polarization InSAR method based on polarimetric optimization for monitoring deformation[J]. Science of Surveying and Mapping, 2021, 46(12):75-82.

[30]
蒋弥, 丁晓利, 何秀凤, 等. 基于快速分布式目标探测的时序雷达干涉测量方法:以Lost Hills油藏区为例[J]. 地球物理学报, 2016, 59(10):3592-3603.

DOI

Jiang M, Ding X L, He X F, et al. FaSHPS-InSAR technique for distributed scatterers:A case study over the lost hills oil field,California[J]. Chinese Journal of Geophysics, 2016, 59(10):3592-3603.

[31]
蒋弥, 丁晓利, 李志伟. 时序InSAR同质样本选取算法研究[J]. 地球物理学报, 2018, 61(12):4767-4776.

DOI

Jiang M, Ding X L, Li Z W. Homogeneous pixel selection algorithm for multitemporal InSAR[J]. Chinese Journal of Geophysics, 2018, 61(12):4767-4776.

[32]
Jiang M, Yong B, Tian X, et al. The potential of more accurate InSAR covariance matrix estimation for land cover mapping[J]. ISPRS Journal of Photogrammetry and Remote Sensing, 2017,126:120-128.

[33]
方海伦, 胡荣荣, 邱雨欣. 基于PS-InSAR技术的宁波市地面沉降监测研究[J]. 地质灾害与环境保护, 2023, 34(4):50-55.

Fang H L, Hu R R, Qiu Y X. Study on land subsidence monitoring basing on ps-insar technique[J]. Journal of Geological Hazards and Environment Preservation, 2023, 34(4):50-55.

[34]
Azadnejad S, Maghsoudi Y, Perissin D. Evaluation of polarimetric capabilities of dual polarized Sentinel-1 and TerraSAR-X data to improve the PSInSAR algorithm using amplitude dispersion index optimization[J]. International Journal of Applied Earth Observation and Geoinformation, 2020,84:101950.

[35]
Iglesias R, Mallorqui J J, López-Dekker P. DInSAR pixel selection based on sublook spectral correlation along time[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(7):3788-3799.

DOI

[36]
温浩, 高峰, 胡在凰, 等. 地面沉降InSAR监测数据融合方法——以宁波市为例[J]. 测绘通报, 2024(s2):12-16.

Wen H, Gao F, Hu Z H, et al. Data fusion method of ground subsidence InSAR monitoring—a case study of Ningbo City[J]. Bulletin of Surveying and Mapping, 2024(s2):12-16.

[37]
马张烽, 蒋弥, 李桂华, 等. 空间网络对时序InSAR相位解缠的影响——以Delaunay与Dijkstra网络为例[J]. 测绘学报, 2022, 51(2):248-257.

DOI

Ma Z F, Jiang M, Li G H, et al. Effects of spatial network on time series InSAR phase unwrapping:Take the Delaunay and Dijkstra networks for example[J]. Acta Geodaetica et Cartographica Sinica, 2022, 51(2):248-257.

Options
文章导航

/

〈 〉