地基大视场SAR形变监测的非均匀大气相位校正方法

白泽朝; 王彦平; 王振海; 胡俊; 李洋; 林赟

doi:10.12000/JR22120

地基大视场SAR形变监测的非均匀大气相位校正方法

DOI: 10.12000/JR22120

白泽朝^{1, 2},
王彦平^2, ,,
王振海³,
胡俊³,
李洋²,
林赟²

1.
北方工业大学电气与控制工程学院北京 100144
2.
北方工业大学信息学院雷达监测技术实验室北京 100144
3.
中南大学地球科学与信息物理学院长沙 410083

基金项目: 国家自然科学基金重点国际合作研究项目(61860206013)，国家重点研发计划资助(2018YFC1505103)

详细信息

作者简介:
白泽朝，博士生，主要研究方向为星载/地基InSAR技术理论和应用

王彦平，博士，教授，博士生导师，主要研究方向为雷达三维成像、地基SAR系统研制及其应用

王振海，博士生，主要研究方向为星地联合三维形变监测

胡　俊，博士，教授，博士生导师，主要研究方向为InSAR大地测量技术及其在地质灾害监测中的应用

李　洋，博士，副教授，硕士生导师，主要研究方向为极化SAR、混合极化SAR信息处理与应用

林　赟，博士，副教授，硕士生导师，主要研究方向为合成孔径雷达三维成像技术、多角度SAR图像基础理论与方法

通讯作者:
王彦平 wangyp@ncut.edu.cn

责任主编：胡程 Corresponding Editor: HU Cheng
中图分类号: TN95
计量
- 文章访问数: 1122
- HTML全文浏览量: 733
- PDF下载量: 183
- 被引次数: 10
出版历程
- 收稿日期: 2022-06-22
- 修回日期: 2022-07-11
- 网络出版日期: 2022-07-28
- 刊出日期: 2023-02-28

A Non-homogenous Atmospheric Compensation Method for Deformation Monitoring of Wide-field Ground-based SAR

BAI Zechao^{1, 2},
WANG Yanping^{2
, ,},
WANG Zhenhai³,
HU Jun³,
LI Yang²,
LIN Yun²

1.
School of Electrical and Control Engineering, North China University of Technology, Beijing 100144, China
2.
Radar Monitoring Technology Laboratory, School of Information Science and Technology, North China University of Technology, Beijing 100144, China
3.
School of Geoscience and Info-physics, Central South University, Changsha 410083, China

Funds: The National Natural Science Foundation of China (61860206013), The National Key R&D Program of China (2018YFC1505103)

More Information

Corresponding author: WANG Yanping, wangyp@ncut.edu.cn

摘要

摘要: 大气变化是地基干涉合成孔径雷达(GB-InSAR)形变监测的主要干扰因素。由于监测环境的地形复杂，水汽、湿度和温度的空间异质性，基于均匀大气介质假设的校正方法可能导致大气校正精度较低。该文提出了一种两阶段半经验模型，用于估计和校正复杂大气条件下特大滑坡GB-InSAR监测过程中出现的大气相位误差。该方法兼顾线性大气相位和非线性大气相位，首先根据测区的范围和高程对观测到的大气相位进行建模，校正与地形相关的线性大气相位。然后，考虑复杂大气条件和方位向大视场角度场景下出现的空间域非均匀大气情况，选取稳定永久散射体(PS)通过插值的方式获取所有PS点的大气相位，校正非线性大气相位。采用该方法对三峡库区新铺和藕塘特大滑坡的地基大视场雷达图像进行处理，相比于常规方法减小大气相位误差最大约2 mm。能有效校正特大滑坡监测场景下出现的非均匀大气相位，满足特大滑坡广域范围监测需求。
- 地基干涉合成孔径雷达 /
- 大气相位校正 /
- 线性变化 /
- 非线性变化 /
- 稳定永久散射体
Abstract: Atmospheric influence is the main interference factor in Ground-Based Interferometric Synthetic Aperture Radar (GB-InSAR) deformation monitoring. Due to the complex terrain and various environmental factors, the correction method based on the assumption of a uniform atmospheric influence may lead to low atmospheric correction accuracy. In this paper, a two-stage semi-empirical model is proposed to correct the atmospheric phase screen during the GB-InSAR monitoring of a super large slope under complex atmospheric conditions. First, the observed atmospheric phase is modeled according to the height and range of the terrain structure to correct the linear atmospheric phase. Then, considering the complex atmospheric conditions and the spatially nonuniform atmosphere with a large azimuth field of view, stable Persistent Scatterers (PS) are selected to obtain the atmospheric phase of all PS by interpolation to correct the nonlinear atmospheric phase. This method is used to process a large field of view radar image of the foundation of the Xinpu and Outang landslides in the Three Gorges Reservoir area. Compared with the conventional method, the atmospheric phase error is reduced by approximately 2 mm. This method effectively corrects the nonuniform atmospheric phase under the landslide monitoring scene and meets the wide-area monitoring needs of the landslide.
- Ground-Based Interferometric Synthetic Aperture Radar (GB-InSAR) /
- Atmospheric phase correction /
- Linear variation /
- Nonlinear variation /
- Stable persistent scatterers

HTML全文

1. 引言

滑坡是我国主要地质灾害之一，其中特大型滑坡，其滑坡体量巨大，往往会造成严重的财产损失和人员伤亡，滑坡地表变形是实现地质灾害监测和预警的重要辅助信息。地基干涉合成孔径雷达(Ground-Based Interferometric Synthetic Aperture Radar, GB-InSAR)是一种高精度的地表形变监测设备，雷达天线在线性轨道上移动，形成方位合成孔径，有效提升了影像方位向分辨率，通过对同一场景的重复扫描和成像，采用干涉测量方法提取目标区域的形变信息^[1,2]。与星载干涉合成孔径雷达(Interferometric Synthetic Aperture Radar, InSAR)^[3,4]相比，GB-InSAR具有观测角度灵活、测量精度高(可达亚毫米)的优点，是局部区域高精度形变监测的重要监测手段。对于特大型滑坡监测，不仅要求GB-InSAR平台实现远距离覆盖，而且在方位向上具有大视场角度成像能力，来满足特大滑坡广域范围监测需求。

此外，观测范围在方位向上的扩大，使得场景内环境更为复杂，特别是观测中大气变化的影响，大气会引起电磁波折射，从而造成电磁波传播路径和方向的变化。在两次成像时刻大气条件(温度、气压和湿度)的变化会导致不同的传播延迟，并且在相位中将存在大气相位分量(通常称为“大气相位屏”)^[5-7]。目前的GB-InSAR大气相位校正方法主要分为3类，一类基于气象数据的方法^[5]，依据电磁波传播的折射率与场景内温度、气压和湿度数据之间的关系来校正大气相位。一类基于函数模型方法^[8-12]，通常大气影响在空间上具有强相关性，当满足大气在空间均匀性假设时，稳定的永久散射体点相位差在距离方向上的投影呈现线性函数关系，通过求取函数关系即可校正大气相位。此外，考虑了陡峭地形条件下大气折射率会随高程发生改变，有学者进一步考虑了高程对大气延迟系数的影响，给出了距离-高程的多元回归模型^[13,14]或者兼顾水平位置和高程等信息的多参数函数模型^[15,16]。另一类空间插值的方法^[17-21]，当大气在空间非均匀变化时，依据稳定的特征点，通过插值的方式进行大气相位校正。

基于稳定永久散射体(Persistent Scatterers, PS)线性函数模型和距离-高程的多元回归模型的方法易于实现，并且已广泛应用在GB-InSAR形变监测领域。然而，在特大滑坡变形监测场景下，需要在方位向上具有较大的视场角度(>100°)，GB-InSAR系统方位向范围的扩大，使得观测场景内环境更为复杂，可能存在明显的湍流或者局部水汽变化^[13]，导致大气介质在方位方向呈现非均匀的特性^[21]。因此，传统的基于稳定PS点函数模型的方法无法在复杂大气条件下准确校正大气相位。

在方位向大视场角度山区复杂大气条件下，大气参数的随机变化，可能在空间上存在线性和非线性两个特定特征的异常相位波动。针对上述问题，本文提出一种两阶段半经验模型校正线性和非线性大气相位。选择具有方位向较大视场角度和跨越长江观测的三峡库区新铺与藕塘滑坡实验区，详细介绍了大气校正的处理流程，并采用本文提出的两阶段半经验模型和常规方法进行大气相位校正实验，验证了算法的有效性。

2. 两阶段半经验模型大气相位校正方法

2.1 PS点选择

PS点选择是进行大气相位校正的必要前提，需要根据实际观测场景，避免植被波动和水域等造成误选点情况，选择合适的选点策略。一般根据点的幅度离散和相位标准差的统计特性，在高信噪比的情况下，目标点的稳定性直接由幅度离差指数 ${D_{\rm{A}}}$ 表示为^[22]

${D_{\rm{A}}} = \frac{{{\delta _{\rm{A}}}}}{{{m_{\rm{A}}}}}$

(1)

其中， ${\delta _{\rm{A}}}$ 是幅度标准偏差， ${m_{\rm{A}}}$ 是幅度平均值。幅度离差指数越小，幅度信息越稳定，然而场景中存在水域时，往往会出现误选等情况，可以通过估算相干系数来估计每个像素点的干涉相位质量，滤除水域的误选点，公式为^[23]

$\gamma = |\gamma | \cdot \exp \left\{ {{\rm{j}}\phi } \right\} = \frac{{{\rm{E}}\left\{ {{S_1} \cdot S_2^*} \right\}}}{{\sqrt {{\rm{E}}\left\{ {{{\left| {{S_1}} \right|}^2}} \right\} \cdot \left\{ {{{\left| {{S_2}} \right|}^2}} \right\}} }}$

(2)

其中， ${S_1}$ 和 ${S_2}$ 是与干涉图相同像素所对应的复值，|·|为绝对值运算符， ${\rm{E}}\left\{ {} \right\}$ 为期望值。将符合幅度离差阈值和相干性阈值的点选作PS点。

2.2 距离-高程模型初补偿

GB-InSAR目标点的干涉相位可以表示为^[1]

$\phi = {\phi _{{\text{def}}}} + {\phi _{{\text{atm}}}} + {\phi _{{\text{noise}}}} + 2k\pi$

(3)

其中，k为整数表示模糊度， ${\phi _{{\text{def}}}}$ 表示形变相位分量， ${\phi _{{\text{atm}}}}$ 表示大气相位分量， ${\phi _{{\text{noise}}}}$ 表示误差相位分量。为了获取高精度的形变相位 ${\phi _{{\text{def}}}}$ ，需要有效的大气相位 ${\phi _{{\text{atm}}}}$ 校正方法。

通过沿路径L积分得到雷达传感器与目标点之间大气相位为^[10]

${\phi _{{\text{atm}}}} = \frac{{4\pi }}{\lambda }\int_L {n(\gamma ,t)} {\text{d}}L$

(4)

其中，n表示折射率的变化，它随时间t和距离 $\gamma$ 而变化，L表示信号的传输路径， $\lambda$ 是波长。一般情况下，假设大气在空间上是匀质的，n在距离 $\gamma$ 上不发生变化，在时间t上随机变化，因此可以将干涉图中大气相位分量表示为随距离变化的公式^[8]

$\phi {(t)_{{\text{atm}}}} = \frac{{4\pi }}{\lambda }({\beta _0} + {\beta _1}\gamma )$

(5)

其中， ${\beta _0}$ 为常数项， ${\beta _1}$ 表示与距离相关的线性系数， $\gamma$ 表示传感器与目标点之间的距离。另外，考虑到地形对大气相位的影响，兼顾距离和地形的大气相位可以表示为^[16]

$\phi {(t)_{{\text{atm}}}} = \frac{{4\pi }}{\lambda }({\beta _0} + {\beta _1}\gamma + {\beta _2}\gamma h)$

(6)

其中， ${\beta _0}$ 为常数项， ${\beta _1}$ 表示与距离相关的线性系数， ${\beta _2}$ 表示与距离、地形相关的系数， $\gamma$ 表示传感器与目标点之间的距离，h表示目标点的高程值。可以建立大气相位与距离和地形的方程组

$\begin{split} & {\boldsymbol{\phi}} = {\boldsymbol{X\beta}} + {\boldsymbol{\varepsilon}}, {\boldsymbol{\phi}} = \left[ {\begin{array}{*{20}{c}} {{\phi _1}} \\ {{\phi _2}} \\ \vdots \\ {{\phi _n}} \end{array}} \right],\\ & {\boldsymbol{X}}= \dfrac{{4\pi }}{\lambda }\left[ {\begin{array}{*{20}{c}} 1&{{\gamma _1}}&{{\gamma _1}{h_1}} \\ 1&{{\gamma _2}}&{{\gamma _2}{h_2}} \\ \vdots & \vdots & \vdots \\ 1&{{\gamma _n}}&{{\gamma _n}{h_n}} \end{array}} \right] , {\boldsymbol{\beta}} = \left[ \begin{gathered} {\beta _0} \\ {\beta _1} \\ {\beta _2} \\ \end{gathered} \right], \\ & {\boldsymbol{\varepsilon}} = \left[ {\begin{array}{*{20}{c}} {{\varepsilon _1}} \\ {{\varepsilon _2}} \\ \vdots \\ {{\varepsilon _n}} \end{array}} \right] \end{split}$

(7)

其中， ${\boldsymbol{\phi}}$ 为n个PS点干涉相位构成的 $n \times 1$ 维向量，X为常数1和n个PS点距离与高程值构成的 $n \times 3$ 维向量， ${\boldsymbol{\beta}}$ 为待估计的 $3 \times 1$ 维向量， ${\boldsymbol{\varepsilon}}$ 为相位误差构成的 $n \times 1$ 维向量。采用最小二乘对 ${\boldsymbol{\beta}}$ 进行求解

$\hat {\boldsymbol{\beta}} = {({{\boldsymbol{X}}^{\bf{T}}}{\boldsymbol{X}})^{ - 1}}{{\boldsymbol{X}}^{\bf{T}}}{\boldsymbol{\phi}}$

(8)

则大气相位的估计值为

${\boldsymbol{\phi}} _{{\text{atm}}} = {\boldsymbol{X}}\hat {\boldsymbol{\beta}}$

(9)

将干涉相位 ${\boldsymbol{\phi}}$ 减掉估计的大气相位 ${{\boldsymbol{\phi}} _{{\text{atm}}}}$ 即为大气校正后相位。在线性大气相位校正阶段，基于相干性和幅度离差进行PS点选择，设置较高的阈值选作高质量PS点，采用距离-高程函数模型估算大气相位。此外，高质量PS点中可能存在一些不可靠点，比如形变区的点和相位跳变的点，为了确保估算结果的准确性，将大于二倍相位标准差 $\sigma$ 的点剔除，迭代求解距离-高程函数模型参数，重新校正大气相位。

$\left| {{\boldsymbol{\phi}} - {{\boldsymbol{\phi}} _{{\text{atm}}}}} \right| < 2\sigma$

(10)

2.3 空间维插值补偿

在校正随高程变化的大气相位后，计算观测场景内高质量PS点形变，设定形变阈值掩模选择出稳定PS点。稳定PS点的相位主要包含大气相位和噪声相位，大气相位在空间上表现出非匀质特征，但在较小的距离范围内可以视为匀质特征，噪声相位在空间上表现为随机特征，不具备空间相关性。在一定的空间范围内进行均值滤波，则可以有效滤除噪声相位影响，同时基于稳定的PS点，采用反距离权重插值算法，插值出采用较低阈值选取的PS点大气相位，从而有效校正非线性大气相位，本文所研究的高质量和低阈值设置在4.2节详细介绍。其中反距离权重插值算法是利用已知点与待插值点的距离来定义权重，根据周边的已知点加权计算待插值点的相位值，距离越近权重越大，其公式为

${Z_{(x,y)}} = {{\sum\limits_{i = 1}^n {\frac{{{Z_i}}}{{{{\left| {{d_i}} \right|}^\mu }}}} } \mathord{\left/ {\vphantom {{\sum\limits_{i = 1}^n {\frac{{{Z_i}}}{{{{\left| {{d_i}} \right|}^\mu }}}} } {\sum\limits_{i = 1}^n {\frac{1}{{{{\left| {{d_i}} \right|}^\mu }}}} }}} \right. } {\sum\limits_{i = 1}^n {\frac{1}{{{{\left| {{d_i}} \right|}^\mu }}}} }}$

(11)

其中， ${Z_{(x,y)}}$ 为插值点结果， $(x,y)$ 为待插值点的空间坐标， ${Z_i}$ 为第i个已知点相位值， $\left| {{d_i}} \right|$ 为第i个已知点与待插值点的空间距离， $\mu$ 为权重的幂指数，一般取2^[20]。此外，考虑到大气相位的空间相关性和计算效率，只选取待插值点周边最近的3个已知点参加计算。

本文算法的大气相位校正流程如图1所示，首先，利用幅度离差阈值和相干性阈值选取PS点，设置较高的阈值选作高质量PS点，参与后续的线性大气相位估计。然后，根据经验判断干涉图中大气误差影响范围和特点，基于高质量PS点采用距离-高程的函数模型对干涉相位校正高程相关的大气影响，在计算中剔除大于二倍相位标准差PS点，迭代求解距离-高程函数模型参数。最后计算序列图像形变，设定阈值选择出稳定PS点，采用反距离权重插值算法估计低阈值选取的PS点大气相位，从而有效校正非线性大气相位。

图 1 本文算法大气相位校正流程

Figure 1. The workflow for atmospheric phase screen processing steps

下载: 全尺寸图片幻灯片

理论上，在地形复杂区域距离-高程的函数模型可以提高大气校正的准确性。然而，在实践中干涉图可能同时存在不同特征的大气相位，简单的函数模型可能校正效果无法满足监测需求。为了在复杂大气条件下保持形变测量的可靠性，拓展GB-InSAR在大视场观测的应用场景，本文方法同时兼顾了线性和非线性大气相位，实现不同特征大气相位的有效校正。

3. 研究区描述和系统参数

3.1 研究区描述

研究区为新铺和藕塘滑坡实验区，位于长江三峡库区重庆市奉节县安坪镇。图2为研究区范围，其中雷达布设在长江北岸，跨江观测长江南岸的橘色半椭圆区，观测区长度约为3500 m，宽度约为2500 m，雷达在方位向上的观测视角为120°。新铺滑坡位于长江南岸(右岸)斜坡地带，北抵长江河床，南至泰山庙山脊，总体斜长约2 km，相对高差约300 m，斜坡坡度15°～20°。藕塘滑坡位于新铺滑坡的西侧，属浅中切割单斜低山河谷地貌，岩层倾向与坡向近于一致，长江从藕塘北边由西流向北东，与岩层走向夹角为10°～15°。

图 2 研究区范围

Figure 2. The study area

下载: 全尺寸图片幻灯片

3.2 地基大视场SAR系统

实验采用自主研制的地基大视场SAR系统对新铺和藕塘滑坡进行了观测，连续获取了29幅雷达图像，时间从2021年7月27日17:44～22:56。该GB-InSAR系统采用角度关联合成大视场成像技术，实现方位向120°以上大视场角度范围的形变监测。如图3所示，其工作在Ku波段，波长为1.74 cm，其方位分辨率在1 km处为5.4 m，距离分辨率为0.37 m。

图 3 GB-InSAR系统

Figure 3. The GB-InSAR system

下载: 全尺寸图片幻灯片

4. 实验分析

4.1 大气相位分析

本文提出的大气相位校正模型应用于三峡库区新铺和藕塘滑坡形变观测中，并对大气校正效果进行了评估。首先，对相邻图像进行干涉处理，每个干涉图的时间基线约为10 min。根据观测区情况，选择同时满足幅度离差指数和相干系数的点作为PS点，其中设置幅度离差指数0.15，相干系数0.9，筛选出25837个高质量PS点。图4(a)和图4(b)分别为平均幅度图和平均相干系数图，图4(c)为高质量PS点的幅度离差图，可以看出坡体上大部分像素点的幅度稳定性和相干系数都很高，不存在水域和植被区的误选点。

图 4 GB-InSAR数据参数

Figure 4. GB-InSAR data parameters

下载: 全尺寸图片幻灯片

图5显示了受大气影响小的干涉图，图6和图7分别显示了受到随高程变化的大气和非匀质大气影响的干涉图，使用不同的色标表示– $\pi$ ～ $\pi$ 范围内观察到的相位变化。干涉对A的影像采集时间分别为19:10与19:20，干涉对B的影像采集时间分别为20:04与20:14，干涉对C的影像采集时间分别为19:53与20:04。在数据采集的间隔时间内，图5(a)中的干涉相位没有超过色标范围，也没有显示出重复的色标变化，图5(b)的高质量PS相位散点也都在0 rad附近波动，表明在2500 m的最大范围内没有观察到相位缠绕现象。另外，由于GB-InSAR的成像几何，可能在雷达近处存在有强目标的干扰，会在部分干涉图零方位角上出现一条相位变化十分均匀的条纹，后续通过PS点选择可以避免对大气相位校正的影响。

图 5 干涉对A相位

Figure 5. Interferogram A

下载: 全尺寸图片幻灯片

图 6 干涉对B相位

Figure 6. Interferogram B

下载: 全尺寸图片幻灯片

图 7 干涉对C相位

Figure 7. Interferogram C

下载: 全尺寸图片幻灯片

大气引起的相位变化如图6(a)和图7(a)所示，可以将大气根据空间分布分为两种不同情况。图6(a)为受随高程变化大气影响的干涉相位，从其相位散点图6(b)中可以看出大气在1000～1500 m间具有恒定的相位递减，由于大气条件变化存在明显大气相位分量，该距离范围内干涉图中高质量PS点相位值位于正相位区域的0～1 rad。相位变化沿距离和地形具有恒定的递减，可以通过应用距离-高程模型对大气进行校正。

如图7(a)所示，干涉对C存在更为复杂的非匀质大气影响。由于地基大视场SAR观测场景的扩大，在受到风或者湍流大气影响时，场景内匀质大气的假设不再成立，从图7(b)相位散点可以看出，大气在1000～1500 m处呈现0～1 rad的相位递减，而在1500～2500 m处存在明显的非线性相位波动，因此采用传统的距离-高程的大气相位校正方法建模将难以有效去除该区域大气影响。

传统的基于函数模型的大气校正方法仅适用于小范围的匀质大气，即当气压、温度和湿度沿距离或者地形均匀变化时。因此，传统的大气校正方法不能作为所有干涉图大气校正的通用方法，必须设计一套能够处理复杂大气条件下的算法，使得可以适应地基大视场SAR观测场景，尤其是跨江等水汽复杂的观测中。

4.2 大气相位校正效果分析

由于观测范围广且地形复杂，大气条件变化明显，相应的大气相位误差也较大(如图6和图7)。在第1阶段采用距离-高程模型校正大气相位，图8和图9为在1000～1500 m范围校正随高程变化大气后的干涉图和对应高质量PS点相位散点图，PS点相位在0 rad波动。与图6和图7对比，干涉对B和干涉对C都有效校正了距离-高程相关的线性大气相位，最大的校正相位分量约为1 rad。但是，由于环境复杂，大气在空间上呈现不均匀变化情况，即使校正线性大气相位后，两个干涉图仍然存在非线性大气相位变化。如图9(b)所示，随着距离的变化，高质量PS点的相位呈现非线性变化趋势。

图 8 距离-高程校正后干涉对B相位

Figure 8. Interferogram B after distance-elevation correction

下载: 全尺寸图片幻灯片

图 9 距离-高程校正后干涉对C相位

Figure 9. Interferogram C after distance-elevation correction

下载: 全尺寸图片幻灯片

在第2阶段采用非线性大气相位校正方法，对选取的高质量PS点计算累计形变量，设定形变阈值为5 mm，选择出稳定PS点，然后选取同时满足幅度离差指数0.25和相干系数0.8的75104个较低阈值的PS点，基于稳定的PS点通过插值的方式估算所有PS点的大气相位。如图10和图11为校正后的干涉对B和干涉对C，在校正前干涉对B在沿江边距离向约1000 m位置存在约1 rad相位波动，在距离向约1500 m位置存在约–0.7 rad相位波动(图10(a))。经过非线性大气校正后，补偿了明显的大气相位，干涉对B整体相位在0 rad附近波动。另外干涉对C存在更为复杂的大气相位，在校正前，不仅沿江边约1000 m, 1500 m位置存在明显大气相位，在距离向约2500 m最远距离位置也存在约–0.5 rad相位波动(图11(a))，非线性大气校正后，干涉对C得到了有效相位补偿。

图 10 非匀质大气校正后干涉对B相位

Figure 10. Interferogram B after correcting nonlinear atmospheric

下载: 全尺寸图片幻灯片

图 11 非匀质大气校正后干涉对C相位

Figure 11. Interferogram C after correcting nonlinear atmospheric

下载: 全尺寸图片幻灯片

与小范围观测不同，GB-InSAR在远距离、大视场的观测中显示出更为复杂的大气相位变化，但可以通过选取稳定的有效测量点，提高形变信息的准确性。图12为获取的观测区累计形变量，如图12(a)所示，未校正大气相位时在江沿岸约1000 m和远距离2000~2500 m位置处存在明显大气相位引起的约5 mm形变误差，如图12(b)所示，仅使用距离-高程函数方法在沿江和远距离处大气相位有所削弱，但沿江边区域仍然存在大气相位，造成了约–3 mm的形变误差，如图12(c)所示，使用本文方法有效补偿了沿江和远距离区域大气相位，表明本文算法的大气相位校正结果相位误差显著降低。理论上，在复杂的场景根据大气相位的特征本文方法可以提高大气估计的准确性和形变提取的鲁棒性。从累计形变量图分析，监测区较稳定，在距离向约2000 m位置处，如图12(c)中红色矩形框所示，存在约4.2 mm的形变区，结合已有研究^[24]，推测在监测时间段内，滑坡的后缘区存在明显的形变，即在藕塘滑坡存在向下滑动的情况。此外，在形变区外存在零星分布的少量形变点，推测可能是植被扰动造成的。

图 12 累计形变量

Figure 12. Cumulative deformation

下载: 全尺寸图片幻灯片

为验证大气相位校正的准确性，挑选了在不同距离向和不同滑坡位置具有代表性的4个像元点变形曲线进行分析，图12(c)中标记了像元点位置，图13为像元点变形曲线图。P1位于新铺滑坡的高程较大位置，经过距离-高程函数校正后反而误差更大，最大误差约为2.4 mm。P2位于藕塘滑坡的高程较大位置，在未进行大气相位校正前，形变的波动幅度较大，最大误差约为–1.6 mm，经过距离-高程函数校正后，形变波动有所变小，但仍然存在约为–1.1 mm误差，本文方法校正后，变形序列变化更为随机且在0.5 mm范围内波动。P3位于监测到的藕塘形变区，总体呈现向下滑动的趋势，由于距离-高程函数方法主要校正沿江和远距离区域随高程变化的大气相位，该区域未进行大气校正，而本文方法大气校正后没有将真实的形变补偿掉，仍然保留了原有的形变趋势。P4位于新铺滑坡靠近长江的沿岸，经过距离-高程函数校正后最大误差由1.5 mm减小为1.1 mm，本文方法校正后，变形序列变化更为随机且在0.5 mm范围内波动。因此，本文方法不仅保留了真实的形变区，而且在非形变区实现了误差的有效降低。

图 13 典型点变形曲线对比

Figure 13. Cumulative deformation map

下载: 全尺寸图片幻灯片

通过三峡库区新铺和藕塘特大滑坡的地基大视场雷达图像实验，验证和分析了本文方法的有效性与准确性。在校正江沿岸和远距离位置随高程变化的大气相位，本文方法可以达到与距离-高程函数方法相同甚至更优的精度。在校正非均匀大气相位，函数模型方法往往会失效，采用大量的稳定PS点作为参考点，空间维插值的方式，可以有效保证非均匀大气相位的校正精度。

5. 结语

本文采用地基大视场SAR系统识别新铺和藕塘滑坡的形变区，观测到了由大气条件变化引起的干涉相位变化。由于该地区地形和大气条件复杂，传统的大气相位校正方法通常会提供不准确的形变结果。复杂大气环境下的干涉相位不仅存在随距离-高程变化的线性大气相位成分，还存在非线性大气相位成分。因此，提出了一种两阶段半经验模型估计各层大气相位，以消除空间域中的大气相位变化。该方法基于空间模型和统计补偿方法，在不使用温度、气压或湿度等大气参数的情况下，对空间非均匀性大气相位进行补偿，对于场景环境信息较少的GB-InSAR观测，可以有效地消除大气效应和准确地获取形变信息。本文方法不需要在站点和实验阶段收集与处理气象数据，拓展了GB-InSAR的适用性且能更好地对大范围下的大气误差进行校正。

本文方法还存在着一些不足，由于地基SAR系统工作在Ku波段，容易受到植被扰动影响，在形变区外仍然存在零星分布的由于植被扰动造成的误差点，需要在后续工作中进一步研究如何避免植被扰动带来的误差。

图 1 本文算法大气相位校正流程

Figure 1. The workflow for atmospheric phase screen processing steps

下载: 全尺寸图片幻灯片

图 2 研究区范围

Figure 2. The study area

下载: 全尺寸图片幻灯片

图 3 GB-InSAR系统

Figure 3. The GB-InSAR system

下载: 全尺寸图片幻灯片

图 4 GB-InSAR数据参数

Figure 4. GB-InSAR data parameters

下载: 全尺寸图片幻灯片

图 5 干涉对A相位

Figure 5. Interferogram A

下载: 全尺寸图片幻灯片

图 6 干涉对B相位

Figure 6. Interferogram B

下载: 全尺寸图片幻灯片

图 7 干涉对C相位

Figure 7. Interferogram C

下载: 全尺寸图片幻灯片

图 8 距离-高程校正后干涉对B相位

Figure 8. Interferogram B after distance-elevation correction

下载: 全尺寸图片幻灯片

图 9 距离-高程校正后干涉对C相位

Figure 9. Interferogram C after distance-elevation correction

下载: 全尺寸图片幻灯片

图 10 非匀质大气校正后干涉对B相位

Figure 10. Interferogram B after correcting nonlinear atmospheric

下载: 全尺寸图片幻灯片

图 11 非匀质大气校正后干涉对C相位

Figure 11. Interferogram C after correcting nonlinear atmospheric

下载: 全尺寸图片幻灯片

图 12 累计形变量

Figure 12. Cumulative deformation

下载: 全尺寸图片幻灯片

图 13 典型点变形曲线对比

Figure 13. Cumulative deformation map

下载: 全尺寸图片幻灯片

参考文献(24)

[1]	WANG Yanping, HONG Wen, ZHANG Yuan, et al. Ground-based differential interferometry SAR: A review[J]. IEEE Geoscience and Remote Sensing Magazine, 2020, 8(1): 43–70. doi: 10.1109/MGRS.2019.2963169
[2]	刘斌, 葛大庆, 李曼, 等. 地基合成孔径雷达干涉测量技术及其应用[J]. 国土资源遥感, 2017, 29(1): 1–6. doi: 10.6046/gtzyyg.2017.01.01 LIU Bin, GE Daqing, LI Man, et al. Ground-based interferometric synthetic aperture radar and its applications[J]. Remote Sensing for Land &Resources, 2017, 29(1): 1–6. doi: 10.6046/gtzyyg.2017.01.01
[3]	BAI Zechao, WANG Yanping, and BALZ T. Beijing land subsidence revealed using PS-InSAR with long time series TerraSAR-X SAR data[J]. Remote Sensing, 2022, 14(11): 2529. doi: 10.3390/rs14112529
[4]	HU Jun, LIU Jihong, LI Zhiwei, et al. Estimating three-dimensional coseismic deformations with the SM-VCE method based on heterogeneous SAR observations: Selection of homogeneous points and analysis of observation combinations[J]. Remote Sensing of Environment, 2021, 255: 112298. doi: 10.1016/j.rse.2021.112298
[5]	RODELSPERGER S. Real-Time Processing of Ground Based Synthetic Aperture Radar (GB-SAR) Measurements[M]. Darmstadt: Technische Universität Darmstadt, 2011.
[6]	IGLESIAS R, AGUASCA A, FABREGAS X, et al. Ground-based polarimetric SAR interferometry for the monitoring of terrain displacement phenomena–Part I: Theoretical description[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2015, 8(3): 980–993. doi: 10.1109/JSTARS.2014.2360040
[7]	PIERACCINI M and MICCINESI L. Ground-based radar interferometry: A bibliographic review[J]. Remote Sensing, 2019, 11(9): 1029. doi: 10.3390/rs11091029
[8]	LUZI G, PIERACCINI M, MECATTI D, et al. Ground-based radar interferometry for landslides monitoring: Atmospheric and instrumental decorrelation sources on experimental data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2004, 42(11): 2454–2466. doi: 10.1109/TGRS.2004.836792
[9]	NOFERINI L, PIERACCINI M, MECATTI D, et al. Permanent scatterers analysis for atmospheric correction in ground-based SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2005, 43(7): 1459–1471. doi: 10.1109/TGRS.2005.848707
[10]	PIPIA L, FABREGAS X, AGUASCA A, et al. Atmospheric artifact compensation in ground-based DInSAR applications[J]. IEEE Geoscience and Remote Sensing Letters, 2008, 5(1): 88–92. doi: 10.1109/LGRS.2007.908364
[11]	徐亚明, 周校, 王鹏, 等. GB-SAR构建永久散射体网改正气象扰动方法[J]. 武汉大学学报:信息科学版, 2016, 41(8): 1007–1012, 1020. doi: 10.13203/j.whugis20140507 XU Yaming, ZHOU Xiao, WANG Peng, et al. A method of constructing permanent scatterers network to correct the meteorological disturbance by GB-SAR[J]. Geomatics and Information Science of Wuhan University, 2016, 41(8): 1007–1012, 1020. doi: 10.13203/j.whugis20140507
[12]	黄其欢, 岳建平. 基于稳定点加权的GBSAR大气扰动校正方法[J]. 西南交通大学学报, 2017, 52(1): 202–208. doi: 10.3969/j.issn.0258-2724.2017.01.028 HUANG Qihuan and YUE Jianping. GBSAR atmospheric turbulence calibration based on weighted stable points[J]. Journal of Southwest Jiaotong University, 2017, 52(1): 202–208. doi: 10.3969/j.issn.0258-2724.2017.01.028
[13]	IANNINI L and GUARNIERI A M. Atmospheric phase screen in ground-based radar: Statistics and compensation[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(3): 537–541. doi: 10.1109/LGRS.2010.2090647
[14]	IGLESIAS R, FABREGAS X, AGUASCA A, et al. Atmospheric phase screen compensation in ground-based SAR with a multiple-regression model over mountainous regions[J]. IEEE Transactions on Geoscience and Remote Sensing, 2014, 52(5): 2436–2449. doi: 10.1109/TGRS.2013.2261077
[15]	KARUNATHILAKE A and SATO M. Atmospheric phase compensation in extreme weather conditions for ground-based SAR[J]. IEEE Journal of Selected Topics in Applied Earth Observations and Remote Sensing, 2020, 13: 3806–3815. doi: 10.1109/JSTARS.2020.3004341
[16]	LIU Jie, YANG Honglei, XU Linlin, et al. Novel model-based approaches for non-homogenous atmospheric compensation of GB-InSAR in the azimuth and horizontal directions[J]. Remote Sensing, 2021, 13(11): 2153. doi: 10.3390/rs13112153
[17]	HU Cheng, DENG Yunkai, TIAN Weiming, et al. A compensation method for a time-space variant atmospheric phase applied to time-series GB-SAR images[J]. Remote Sensing, 2019, 11(20): 2350. doi: 10.3390/rs11202350
[18]	DENG Yunkai, HU Cheng, TIAN Weiming, et al. A grid partition method for atmospheric phase compensation in GB-SAR[J]. IEEE Transactions on Geoscience and Remote Sensing, 2021, 60: 5206713. doi: 10.1109/TGRS.2021.3074161
[19]	DENG Yunkai, HU Cheng, TIAN Weiming, et al. 3-D deformation measurement based on three GB-MIMO radar systems: Experimental verification and accuracy analysis[J]. IEEE Geoscience and Remote Sensing Letters, 2021, 18(12): 2092–2096. doi: 10.1109/LGRS.2020.3014342
[20]	胡程, 邓云开, 田卫明, 等. 地基干涉合成孔径雷达图像非线性大气相位补偿方法[J]. 雷达学报, 2019, 8(6): 831–840. doi: 10.12000/JR19073 HU Cheng, DENG Yunkai, TIAN Weiming, et al. A compensation method of nonlinear atmospheric phase applied for GB-InSAR images[J]. Journal of Radars, 2019, 8(6): 831–840. doi: 10.12000/JR19073
[21]	IZUMI Y, ZOU Lilong, KIKUTA K, et al. Iterative atmospheric phase screen compensation for near-real-time ground-based InSAR measurements over a mountainous slope[J]. IEEE Transactions on Geoscience and Remote Sensing, 2020, 58(8): 5955–5968. doi: 10.1109/TGRS.2020.2973533
[22]	FERRETTI A, PRATI C, and ROCCA F. Permanent scatterers in SAR interferometry[J]. IEEE Transactions on Geoscience and Remote Sensing, 2001, 39(1): 8–20. doi: 10.1109/36.898661
[23]	ROSEN P A, HENSLEY S, ZEBKER H A, et al. Surface deformation and coherence measurements of Kilauea Volcano, Hawaii, from SIR-C radar interferometry[J]. Journal of Geophysical Research:Planets, 1996, 101(E10): 23109–23125. doi: 10.1029/96JE01459
[24]	肖捷夫. 库水涨落和降雨条件下藕塘滑坡变形演化机制及其预测模型研究[D]. [博士论文], 中国地质大学, 2021. XIAO Jiefu. Deformation evolution mechanism and displacement prediction model of Outang landslide under water level fluctuation and rainfall[D]. [Ph. D. dissertation], China University of Geosciences, 2021.

施引文献

期刊类型引用(2)

1.	林娜，谭力兵，张迪，丁凯，李双桃，肖茂池，张精平，王小华. 时序InSAR和LSTM结合的滑坡形变时空分析与预测方法. 地球信息科学学报. 2024(12): 2772-2787 . 百度学术
2.	徐甫，王政，李振洪，李永生. 复杂环境下的地基雷达大气改正方法. 武汉大学学报(信息科学版). 2023(12): 2069-2081 . 百度学术

其他类型引用(8)

资源附件(0)

访问统计

图(13)

计量

文章访问数: 1122
HTML全文浏览量: 733
PDF下载量: 183
被引次数: 10

1. 引言
2. 两阶段半经验模型大气相位校正方法
2.1 PS点选择
2.2 距离-高程模型初补偿
2.3 空间维插值补偿
3. 研究区描述和系统参数
3.1 研究区描述
3.2 地基大视场SAR系统
4. 实验分析
4.1 大气相位分析
4.2 大气相位校正效果分析
5. 结语

地基大视场SAR形变监测的非均匀大气相位校正方法

DOI: 10.12000/JR22120

通讯作者:
王彦平 wangyp@ncut.edu.cn

计量

A Non-homogenous Atmospheric Compensation Method for Deformation Monitoring of Wide-field Ground-based SAR

Corresponding author: WANG Yanping, wangyp@ncut.edu.cn

1. 引言

2. 两阶段半经验模型大气相位校正方法

2.1 PS点选择

2.2 距离-高程模型初补偿

2.3 空间维插值补偿

3. 研究区描述和系统参数

3.1 研究区描述

3.2 地基大视场SAR系统

4. 实验分析

4.1 大气相位分析

4.2 大气相位校正效果分析

5. 结语

期刊类型引用(2)

其他类型引用(8)

计量

目录

1. 引言

2. 两阶段半经验模型大气相位校正方法

2.1 PS点选择

2.2 距离-高程模型初补偿

2.3 空间维插值补偿

3. 研究区描述和系统参数

3.1 研究区描述

3.2 地基大视场SAR系统

4. 实验分析

4.1 大气相位分析

4.2 大气相位校正效果分析

5. 结语

期刊介绍

联系我们

地基大视场SAR形变监测的非均匀大气相位校正方法

DOI: 10.12000/JR22120

通讯作者: 王彦平 wangyp@ncut.edu.cn

计量

出版历程

A Non-homogenous Atmospheric Compensation Method for Deformation Monitoring of Wide-field Ground-based SAR

Corresponding author: WANG Yanping, wangyp@ncut.edu.cn

1. 引言

2. 两阶段半经验模型大气相位校正方法

2.1 PS点选择

2.2 距离-高程模型初补偿

2.3 空间维插值补偿

3. 研究区描述和系统参数

3.1 研究区描述

3.2 地基大视场SAR系统

4. 实验分析

4.1 大气相位分析

4.2 大气相位校正效果分析

5. 结语

期刊类型引用(2)

其他类型引用(8)

计量

出版历程

目录

1. 引言

2. 两阶段半经验模型大气相位校正方法

2.1 PS点选择

2.2 距离-高程模型初补偿

2.3 空间维插值补偿

3. 研究区描述和系统参数

3.1 研究区描述

3.2 地基大视场SAR系统

4. 实验分析

4.1 大气相位分析

4.2 大气相位校正效果分析

5. 结语

期刊介绍

联系我们

通讯作者:
王彦平 wangyp@ncut.edu.cn