作者中心
专家审稿
责编办公
编辑办公
High-resolution High-dimensional Imaging of Urban Building Based on GaoFen-3 SAR Data（in English）
-
摘要: 传统合成孔径雷达(SAR)只能获取方位-距离二维图像，无法准确反映目标的三维散射结构信息。层析合成孔径雷达(TomoSAR)是一种多基线干涉测量模式，它将合成孔径原理扩展至高程向，除了可对目标进行二维成像之外，还可以准确恢复目标的高度向散射信息，真正实现三维成像。差分层析合成孔径雷达(D-TomoSAR)将合成孔径原理延伸至高程和时间方向，不仅可以获得目标的三维散射结构，还可以高精度获取观测目标的形变速率，实现对目标形变的有效监测。高分三号是我国首颗1 m分辨率C频段多极化SAR卫星。它具有高分辨率、大成像幅宽、多成像模式等特点，对我国高分对地观测技术的发展具有重要意义。目前高分三号数据主要应用于目标识别等图像处理领域，没有充分利用SAR图像的相位信息。而且，由于设计之初未考虑后续高维成像应用，现有高分三号获取的SAR图像存在有一定的空间、时间去相干问题，对应用于后续干涉系列处理产生了一定影响。为解决上述问题，该文基于7景高分三号SAR复图像，开展了对北京雁栖湖周围建筑的三维、四维层析成像研究，在获取了建筑物三维散射结构信息的同时，实现了对建筑物形变的毫米级高精度监测。该初步实验结果证明了高分三号SAR数据的应用潜力，为后续进一步扩展高分三号SAR卫星在城市感知与监测中的应用提供了技术支撑。
-
关键词:
- 三维成像 /
- 四维成像 /
- 层析合成孔径雷达 /
- 差分层析合成孔径雷达 /
- 高分三号
Abstract: Conventional Synthetic Aperture Radar (SAR) can only obtain two-dimensional (2-D) azimuth-range images without accurately reflecting the three-Dimensional (3-D) scattering structure information of the targets. However, SAR Tomography (TomoSAR) is a multi-baseline interferometric measurement mode that extends the synthetic aperture principle into the elevation direction, making it possible to recover the true height of the target, thereby achieving 3-D imaging. Moreover, Differential SAR Tomography (D-TomoSAR) extends the synthetic aperture principle into the elevation and time directions simultaneously. Thus, it can obtain the target 3-D scattering structure along with the deformation speed of the observed target. GaoFen-3 (GF-3) is the first C-band multi-polarization 1 m resolution SAR satellite of China. It has several advantages, such as high-resolution, large swath width, and multiple imaging modes, which are crucial to the development of a high-resolution earth observation technology for China. Presently, GF-3 data are mainly used in the image processing field, such as target identification. However, the phase information of the SAR images is not yet fully utilized. Moreover, because of the high-dimensional imaging ability that was overlooked at the beginning of designing the system, existing SAR images acquired by GF-3 have spatial and temporal de-coherence problems. Thus, it is difficult to use the images in further interference series processing. To solve the above problems, this study achieved 3-D and four-Dimensional (4-D) imaging of buildings around Yanqi Lake, in Beijing, based on the data of seven SAR complex images. We obtained the 3-D scattering structure information of buildings and achieved millimeter-level high-precision monitoring of building deformation. The preliminary experimental results demonstrate the application potential of GF-3 SAR data and provide a technical support for the subsequent further application of the GF-3 SAR satellite in urban sensing and monitoring. -
1. 引言
与传统光学观测手段不同，合成孔径雷达(Synthetic Aperture Radar, SAR)具有全天时、全天候的工作能力，因而在国土资源勘测、自然灾害监测等领域得到广泛应用[1]。然而传统SAR成像只能获取目标的方位-距离二维图像，无法准确反映目标的三维散射特征，一定程度上影响了SAR图像的进一步应用。层析合成孔径雷达(Synthetic Aperture Radar Tomography, TomoSAR)是SAR成像技术进一步扩展。它将合成孔径原理延伸至高程向，可基于多景二维复图像获得目标的方位-距离-高程信息，进而实现三维成像[2]。1999年，Reigber等人[3,4]首次展示了机载TomoSAR成像结果，实现了基于谱估计方法的森林区域三维成像。2006年，Fornaro和Serafino[5]基于长时间基线星载SAR数据，开展了TomoSAR成像星载实验，证实了谱估计技术应用于高程向重构的可行性与有效性。2011年，Reale等人[6]证实了高分辨率数据和先进干涉处理技术相结合可以更好地重建建筑物三维结构。2013年，Shahzad和Zhu[7]提出了一种全新的建筑物立面重建方法，获取了拉斯维加斯百乐宫酒店的三维点云，证实了TomoSAR点云在构建动态城市模型方面的优越性。2017年，Wang等人[8]提出了迭代重加权的交替方向乘子算法，用于实现快速TomoSAR成像。2018年，Wang和Zhu[9]提出了一种基于核主成分分析的TomoSAR成像方法，使用极小代价即可分离同一方位-距离分辨单元中沿高程向分布的多个散射体。2019年，秦斐等人[10]针对TomoSAR成像中高程向分辨率较低、建筑物叠掩、提取建筑物目标特征效率较低等问题，提出了一种基于机器学习的建筑物目标识别和提取算法，提高了观测目标的特征提取效率，并通过机载阵列三维SAR实验数据验证了方法的有效性。差分层析合成孔径雷达(Differential Synthetic Aperture Radar Tomography, D-TomoSAR)是TomoSAR的进一步扩展。它基于多基线观测数据可实现对目标的方位-距离-高程-时间四维成像，不仅解决了TomoSAR成像中的高度错位与模糊问题，还可高精度获取目标的形变信息。D-TomoSAR概念由Lombardini于2003年首次提出[11]。2007年，Fornaro等人[12]证明了D-TomoSAR可作为传统永久散射体监测技术的有效替代方案，实现对大场景形变的有效监测。2008年，Fornaro等人[13]提出了一种可分离干扰散射体相关时间序列的技术，并应用于罗马地区的D-TomoSAR成像中，证明了该技术可用于城市复杂场景的形变监测。2009年，Zhu等人[14]基于多景TerraSAR-X数据，实现了对拉斯维加斯的D-TomoSAR成像，获取了城市四维雷达地图。2010年，Fornaro等人[15]使用ERS数据对罗马市Grotta Perfetta地区进行了D-TomoSAR成像，结果显示差分层析成像方法可以有效区分同一分辨单元中的不同散射体，克服了差分干涉的局限性，进一步提高了对城市基础设施形变的监测能力。2011年，Zhu和Bamler[16]提出了一种“时间扭曲”方法，并基于TerraSAR-X数据实现了城市区域的D-TomoSAR成像，有效获取了城市建筑的线性运动和季节性运动速率。2015年，Siddique等人[17]将D-TomoSAR方法和永久散射体干涉技术相结合，基于50幅TerraSAR-X图像获得了巴塞罗那中高层建筑物立面散射体的时空反演结果。2020年，Wang和Liu[18]提出了一种广义D-TomoSAR成像系统模型和一种基于拟极大似然的成像算法，同时反演出了建筑物线性运动、季节性运动等多个形变运动速率。目前，TomoSAR和D-TomoSAR成像技术已在城市建筑物和基础设施的三维重建和长期形变监测方面展现了极大的应用潜力。
压缩感知(Compressive Sensing, CS)是一种重要的稀疏信号处理技术，它可以使用较少样本实现对稀疏信号的高质量恢复[19-21]。CS-TomoSAR成像的前提是观测场景的高程向分布稀疏，而城市区域主要为人造建筑，其高程向分布都满足稀疏性条件。因此，CS在城市区域三维、四维成像中具有广阔的应用前景。起初，由于现代米级分辨率星载SAR系统轨道限制，TomoSAR成像高程向分辨率远低于方位向和距离向，因而迫切需要超分辨算法来解决这一问题，2010年，Zhu和Bamler[22]介绍了一种基于CS理论的TomoSAR成像方法，相比于传统谱估计方法，该算法实现了对目标高程向分布的超分辨重构。2012年，Zhu和Bamler[23]将所提出的SL1MMER算法应用到TerraSAR-X星载数据处理中，获得了拉斯维加斯百乐宫酒店的高分辨率TomoSAR成像结果，并证明了该算法具有超分辨能力。2015年，Weiss等人[24]提出了一种适用于TomoSAR的自适应CS算法，准确识别了同一分辨单元中的两个散射体位置。2017年，Li等人[25]研究了基于SPICE的TomoSAR成像方法，并利用8幅TerraSAR-X条带影像，实现了对内蒙古根河市某建筑的高精度三维重建。2010年，Zhu和Bamler[26]将CS技术应用在D-TomoSAR成像中，证明了在高程向多散射体分离上CS技术相比于传统谱估计方法的优越性，同时指出CS可自动识别散射体数量，非常适用于星载SAR系统的三维、四维成像。2010年，Zhu和Bamler[27]基于SL1MMER算法重构获得了拉斯维加斯会议中心建筑群的D-TomoSAR成像结果，展示了CS技术在四维成像方面的优势和能力。2014年，Leng等人[28]将最小绝对收缩和选择算子CS算法应用于建筑区域，并展示了巴塞罗那的高程重建和形变监测结果。
高分三号卫星是我国首颗分辨率达到1 m的C频段多极化SAR卫星，于2016年8月10日在中国太原卫星发射中心由长征四号丙运载火箭发射升空[29]。高分三号卫星是“国家高分辨率对地观测系统重大专项”中唯一的民用微波遥感成像卫星，具有高分辨率、大成像幅宽、多成像模式、长寿命运行等特点，可实现全天时、全天候的全球海洋与陆地监测[30]。目前高分三号已成功应用于高精度测绘、自然灾害监测等多个领域[31,32]。然而，由于设计之初未考虑后续高维成像应用，现有高分三号获取的SAR图像存在有一定的空间、时间去相干问题，对其进一步应用于干涉SAR、差分干涉SAR、TomoSAR、D-TomoSAR等存在一定挑战。2019年，余博等人[33]基于高分三号数据，对河南省登封市周围地区进行了干涉测量实验，通过与哨兵一号卫星数据结果进行对比分析，验证了高分三号的干涉能力以及可提取地表形变信息的能力。2021年7月，黄震等人[34]基于高分三号SAR数据进行了干涉测量实验，成功提取了观测区域的数字高程模型。
本文基于7景高分三号SAR复图像数据，利用CS技术，开展了TomoSAR和D-TomoSAR成像实验研究，获取了北京市雁栖湖周围建筑的高分辨率三维、四维SAR图像，实现了建筑物的高质量三维重建以及高精度形变监测，为后续基于高分三号SAR数据的干涉系列应用及多维高分辨率成像提供了技术支撑。
本文后续结构如下：第2节主要介绍TomoSAR, D-TomoSAR成像模型，并给出了上述两个模型的CS求解方案；第3节介绍了本文所使用的高分三号SAR复图像数据集；第4节基于仿真数据，开展了TomoSAR, D-TomoSAR成像实验，证明了CS技术在高分辨三维成像及高精度形变监测方面的有效性；第5节基于7景高分三号SAR复图像数据，对北京雁栖湖周围建筑进行了TomoSAR, D-TomoSAR成像研究，获取了代表性建筑和大观测区域的三维、四维雷达图像；第6节对文章进行了总结和展望。
2. 成像模型
2.1 TomoSAR成像模型
TomoSAR利用对同一场景观测获取的多幅配准的二维SAR复图像(多基线观测数据)在高程向上进行孔径合成，以获得高程向上的分辨能力，进而重构目标的三维散射信息[3,4,35,36]。TomoSAR成像几何如图1所示。设共有N条基线用于数据获取，令
bn(n=1,2,⋯,N) 表示高程向孔径分布，对于一个选定的方位-距离分辨单元，第n 幅SAR图像对应的聚焦测量值可以表示为gn=∫Δsγ(s)exp(−j2πξns)ds (1) 式中，
ξn=−2bn/(λr) 表示高程向频率，其中λ 为波长，r 为斜距；γ(s) 表示沿高程向s 的复反射函数；Δs 为高程向跨度。沿高程向s 对高程向复反射函数γ(s) 进行离散化，则式(1)中的成像模型可近似表示为gN×1=RN×L⋅γL×1 (2) 式中，L为高程向离散化点数；
g=[g1,g2,⋯,gN]T 表示测量值向量；R=exp(−j2πξnsl) 为根据TomoSAR成像几何所构建的观测矩阵；γ=[γ(s1),γ(s2),⋯,γ(sl)]T 表示高程向离散复反射函数，其中sl(l=1,2,⋯,L) 为离散高程向分布。从式(2)可以看出，TomoSAR成像模型可视为对γ(s) 不规则采样的离散傅里叶变换。因此，一个SAR测量值可以看作目标复反射函数沿高程向的一个谱参数。对于非参数化谱分析问题，高程向理论分辨率ρs 依赖于高程孔径大小Δb ，在高程向采样密集的情况下，ρs 可由式(3)进行计算。ρs=λr2Δb (3) 2.2 D-TomoSAR成像模型
相比于TomoSAR, D-TomoSAR在三维的基础上多了一个时间方向的维度。它沿高程向和形变速度向合成两个孔径，进而获取被观测目标的高程和形变速率的联合分辨率，实现四维成像[11]。对于N个复图像而言，当高程孔径位置为
bn 、时间基线为tn 时，第n 幅图像的聚焦测量值可以表示为gn=∫Δsγ(s)exp(−j2πξns)exp(−j2πηnV(s))ds (4) 式中，
ηn=−2πtn/λ 表示形变速度频率；V(s) 表示形变速率。式(4)中的模型也可写为gn=∫Δs∫Δvγ(s)δ(v−V(s))⋅exp(−j2π(ξns+ηnv))dvds (5) 式中，
Δv 表示被观测目标的形变速率跨度；δ(⋅) 是与形变项相关的谱分布。令aγ(s,v)=γ(s)δ(v−V(s)) ，则式(5)可以写成gn=∫Δs∫Δva(s,v)exp(−j2π(ξns+ηnv))dvds (6) 式(6)中的模型可以视为
a(s,v) 在高程-形变平面的二维傅里叶变换。因此，其在立面轴上的投影为反射率剖面γ(s) [37]。将式(6)中的s 和v 离散化后，D-TomoSAR成像模型可以表示为gN×1=RN×LQ⋅γLQ×1 (7) 其中，
g=[g1,g2,⋯,gN]T 表示测量向量；R=exp(−j2π⋅(ξnsl+ηnvq)) 为D-TomoSAR成像观测矩阵；sl(n=1,2,⋯,L) 为离散高程向分布；vq(q=1,2,⋯,Q) 为离散形变向分布；γ 由离散化的a(s,v) 组成。若时间孔径大小Δt ，则形变分辨率ρv 可由式(8)进行计算ρv=λ2Δt (8) 2.3 CS重建
城市中被观测目标主要为人造建筑，其高程向分布通常都是稀疏的，即每个方位-距离分辨单元中的散射体个数有限。因此，当测量矩阵
R 满足有限等距性质条件时，面向式(2)和式(7)中的模型，本文通过解决如下的最优化问题分别实现基于CS的TomoSAR和D-TomoSAR成像。ˆγ=argmin (9) 式中，
\beta 为正则化参数，与噪声水平和样本数目有关。CS算法可以在短时间内从获取的样本数据中实现高质量信号恢复[38,39]。基于该算法的优势，本文采用CS算法进行TomoSAR和D-TomoSAR成像。3. 高分三号数据集简介
本文所使用的高分三号数据集总共包含7景复图像，具体参数如表1所示。该数据集的7景图像的空间基线孔径大小约为1417 m，时间基线跨度是从2018年6月到2019年9月，共464 d，表2给出其时空基线具体参数，其时空基线分布情况如图2所示。本文以2019年3月1日获取的SAR图像为主影像，其余6景为辅图像。辅图像时空基线位置是相对于主影像计算得到的。本文第4节将基于表1的参数进行点目标仿真实验，第5节将基于该高分三号数据集进行实验。
表 1 高分三号数据集参数Table 1. Parameters of GF-3 dataset参数名称 数值 参数名称 数值 空间基线跨度 1417.4 m 数据景数 7景 时间基线跨度 464 d 方位向分辨率 0.3626 m 斜距 1052747 m 距离向分辨率 0.765692 m 波长 0.056 m 高程向理论分辨率 20.6174 m 入射角 47.2330015° 形变理论分辨率 21.8毫米/年 表 2 高分三号数据集时空基线参数Table 2. Spatial-temporal baseline parameters of GF-3 dataset编号 获取时间 空间基线(m) 时间基线(d) 1 2018.06.13 –459.108 –261 2 2019.01.31 –628.551 –29 3 2019.03.01 0 0 4 2019.03.30 –724.517 29 5 2019.07.24 692.863 145 6 2019.08.22 –38.211 174 7 2019.09.20 –510.491 203 4. 仿真实验
基于表1中的实验参数，本文设定高程向分布有两个散射体，模拟产生了7景仿真数据，添加信噪比为20 dB的噪声，并分别使用经典谱估计方法和CS算法对高程向进行TomoSAR和D-TomoSAR成像[40]。本节将展示3种经典谱估计算法即Beamforming (BF)[41], Adaptive beamforming (Capon)[42]和Multiple signal classification (MUSIC)[43,44]的仿真结果，用于与CS算法重构结果进行比较，以说明CS算法在TomoSAR, D-TomoSAR成像中的优势。图3为高程向两个散射点的TomoSAR成像结果，横坐标为高程向位置分布，纵坐标为散射点的幅度值。本文设置两个散射体的距离分别为11 m和50 m。由图3可以看出，当距离为11 m时，由于小于高程向理论分辨率，3种谱估计算法的成像结果均存在分辨率较低和模糊严重等问题，无法准确分离两个散射点，造成重建失败；而CS算法仍可有效识别两个散射体，实现高程向的超分辨成像。当距离为50 m时，可以看出，BF算法重构结果具有严重的、不规则的旁瓣；Capon和MUSIC算法相较于BF算法，重构结果具有更低的旁瓣，提升了TomoSAR高程向重构质量；而相较于3种谱估计算法，CS则可以更加有效地抑制旁瓣和噪声，进一步提升了高程向散射体的可分辨能力。图4为高程向3个散射点的TomoSAR成像结果，3个散射点的间隔是20 m，与两个散射体的成像结果相似，可以看出CS算法对于多散射体分离也具有很好效果。
图 3 高程向两个散射点的TomoSAR成像结果(左图：两个散射点之间的距离为11 m；右图：两个散射点之间的距离为50 m)Figure 3. TomoSAR reconstructed reflectivity profiles of two scattering points along the elevation direction (left image: the distance between two scattering points is 11 m; right image: the distance between two scattering points is 50 m)
图 4 高程向3个散射点的TomoSAR成像结果(3个散射点之间的间隔为20 m)Figure 4. TomoSAR reconstructed reflectivity profiles of three scattering points along the elevation direction (the distance between three scattering points is 20 m)随后，本文进一步开展了D-TomoSAR仿真实验研究。同样设所关注的方位-距离分辨单元中有两个散射体，位置分别处于–10 m和10 m处，二者之间的距离小于高程向理论分辨率，二者形变大小分别设为4毫米/年和–7毫米/年。基于表1中参数本文获取了相应的仿真数据，并添加了信噪比为20 dB的加性白噪声。图5为不同算法的D-TomoSAR仿真实验结果，横坐标为目标高程向位置，纵坐标为散射体的线性形变量。由图5可以看出，受困于算法局限性，3种谱估计方法均很难分辨出高程向两个散射体，尤其是BF算法，难以突破理论分辨率；而CS算法则可精确区分相距20 m的两个散射体，这与TomoSAR实验结果相吻合，再一次验证了CS技术在三维成像中的超分辨能力。由表1可知，仿真实验的形变理论分辨率约为21.8毫米/年，而仅有CS算法可以准确实现对4毫米/年和–7毫米/年形变量的精确估计，证实了其相较于谱估计技术在散射体形变监测中的优越性。
图 5 D-TomoSAR仿真结果(两个散射体高程位置为–10 m, 10 m；散射体形变速率分别为4毫米/年、–7毫米/年)Figure 5. D-TomoSAR simulation results (elevation position of two scatters are –10 m and 10 m; deformation velocity of two scatters are 4 mm/year and –7 mm/year, respectively)5. 高分三号实验结果
5.1 代表性建筑物反演
本节将展示基于此高分三号SAR数据集的TomoSAR和D-TomoSAR重建结果。首先选取观测场景内两处代表性建筑进行了TomoSAR和D-TomoSAR成像，分别为某生态农业公司和雁栖湖会展中心。图6(a)所示为北京某生态农业公司的光学图像，红色虚线框出的区域是所关注的目标。图6(b)为该生态农业公司的二维SAR图像。图7(a)为TomoSAR重建结果，获取了区域内5栋建筑的高度信息。由图7(a)可以看出，5个独立建筑的高度均为35 m，这与实际该5栋完全相同的建筑高度相符，验证了CS-TomoSAR成像技术的有效性，说明其可以用于高分三号SAR数据的三维高精度成像中。图7(b)为基于CS算法的D-TomoSAR重建结果，可以看出，5栋建筑均有不同程度的形变，这可能与建筑材料随着季节变化的热胀冷缩紧密相关，在日常监测中应该有所关注，防止出现相关危险。为更加直观地展示所关注场景，图8给出了生态农业公司的三维点云，更加准确地反映出了建筑物的三维散射结构。
图 7 生态农业公司高程图及形变速率图Figure 7. Elevation and deformation velocity maps of ecological agricultural company图9(a)为北京雁栖湖国际会展中心的光学影像，该建筑真实高度约30 m。图9(b)为雁栖湖国际会展中心的SAR图像。图10(a)给出了CS-TomoSAR成像结果，并在图11给出了其三维点云图。由图10(a)和图11可以看出，重建结果基本反映了建筑物的三维结构，重构的建筑物高度与实际相吻合，这也说明了本文基于高分三号SAR数据，已可实现对复杂建筑物较高精度的三维重建。图10(b)展示了基于CS技术的雁栖湖国际会展中心D-TomoSAR重建结果。由图10(b)可以看出，该建筑的形变大约在–10毫米/年～10毫米/年之间，且建筑的左半部分和右半部的线性形变速率正好相反，这反映了该建筑下方地面可能处于一边抬高一边塌陷的变化之中，应当着重关注。
图 9 北京雁栖湖国际会展中心Figure 9. Beijing Yanqi lake international convention and exhibition center
图 10 北京雁栖湖国际会展中心高程图及形变速率图Figure 10. Elevation and deformation velocity maps of Beijing Yanqi lake international convention and exhibition center
图 11 北京雁栖湖国际会展中心三维点云图Figure 11. 3-D point cloud of Beijing Yanqi lake international convention and exhibition center5.2 大场景三维、四维成像
本节展示了基于高分三号数据集的大场景TomoSAR和D-TomoSAR成像结果。图12分别给出了所选取的顶秀美泉小镇区域的光学影像和SAR图像。图13(a)给出了基于CS算法的TomoSAR成像结果，重构图像显示该区域内建筑高度均在15 m至20 m之间，准确反映了小镇内均为4～6层居民楼建筑的实际情况。图13(b)给出了基于CS算法的D-TomoSAR成像结果。其显示小镇左下角区域和右上角区域的建筑形变约为–10毫米/年，可能由于这两块区域周围的某些施工导致地面沉降，小镇中心的形变速率约为5毫米/年，总体来说比较稳定。该结果表明，基于高分三号SAR数据，已经可以实现大场景的高质量三维重建以及高精度形变监测，进一步验证了高分三号SAR卫星在城市感知与监测中的应用潜力。
图 13 顶秀美泉小镇区域的高程图及形变速率图Figure 13. Elevation and deformation velocity maps of Dingxiumeiquan town6. 结论
本文开展了基于高分三号SAR复图像数据的TomoSAR和D-TomoSAR成像研究，获取了北京雁栖湖地区两处代表性建筑的三维、四维成像结果，并给出了大面积观测场景的三维、四维雷达图像。研究结果显示了我国高分三号SAR卫星应用于TomoSAR和D-TomoSAR成像方面的潜力，为后续拓展高分三号的干涉系列应用领域提供了技术支撑。
后续本团队将继续开展相关研究工作，收集高分三号多个区域、多景数据信息，同时探索全新的高分辨TomoSAR, D-TomoSAR成像算法，实现面向复杂城市场景的大范围三维、四维成像，进一步挖掘高分三号在干涉系列应用中的巨大潜力。
-
图 3 高程向两个散射点的TomoSAR成像结果(左图：两个散射点之间的距离为11 m；右图：两个散射点之间的距离为50 m)
Figure 3. TomoSAR reconstructed reflectivity profiles of two scattering points along the elevation direction (left image: the distance between two scattering points is 11 m; right image: the distance between two scattering points is 50 m)
图 4 高程向3个散射点的TomoSAR成像结果(3个散射点之间的间隔为20 m)
Figure 4. TomoSAR reconstructed reflectivity profiles of three scattering points along the elevation direction (the distance between three scattering points is 20 m)
图 5 D-TomoSAR仿真结果(两个散射体高程位置为–10 m, 10 m；散射体形变速率分别为4毫米/年、–7毫米/年)
Figure 5. D-TomoSAR simulation results (elevation position of two scatters are –10 m and 10 m; deformation velocity of two scatters are 4 mm/year and –7 mm/year, respectively)
图 7 生态农业公司高程图及形变速率图
Figure 7. Elevation and deformation velocity maps of ecological agricultural company
图 9 北京雁栖湖国际会展中心
Figure 9. Beijing Yanqi lake international convention and exhibition center
图 10 北京雁栖湖国际会展中心高程图及形变速率图
Figure 10. Elevation and deformation velocity maps of Beijing Yanqi lake international convention and exhibition center
图 11 北京雁栖湖国际会展中心三维点云图
Figure 11. 3-D point cloud of Beijing Yanqi lake international convention and exhibition center
图 13 顶秀美泉小镇区域的高程图及形变速率图
Figure 13. Elevation and deformation velocity maps of Dingxiumeiquan town
图 3 TomoSAR reconstructed reflectivity profiles of two scattering points along the elevation direction (left image: the distance between two scatterers is 11 m; right image: the distance between two scatters is 50 m)
图 4 TomoSAR reconstructed reflectivity profiles of three scattering points along the elevation direction (the distance between three scatterers is 20 m)
图 5 D-TomoSAR simulation results (elevation position of two scatterers are –10 m and 10 m; deformation velocity of two scatterers are 4 mm/year and –7 mm/year, respectively)
图 10 Elevation and deformation velocity maps of Beijing Yanqi Lake International Convention and Exhibition Center
图 11 3-D point cloud of Beijing Yanqi Lake International Convention and Exhibition Center
表 1 高分三号数据集参数
Table 1. Parameters of GF-3 dataset
参数名称 数值 参数名称 数值 空间基线跨度 1417.4 m 数据景数 7景 时间基线跨度 464 d 方位向分辨率 0.3626 m 斜距 1052747 m 距离向分辨率 0.765692 m 波长 0.056 m 高程向理论分辨率 20.6174 m 入射角 47.2330015° 形变理论分辨率 21.8毫米/年 
下载: 导出CSV
表 2 高分三号数据集时空基线参数
Table 2. Spatial-temporal baseline parameters of GF-3 dataset
编号 获取时间 空间基线(m) 时间基线(d) 1 2018.06.13 –459.108 –261 2 2019.01.31 –628.551 –29 3 2019.03.01 0 0 4 2019.03.30 –724.517 29 5 2019.07.24 692.863 145 6 2019.08.22 –38.211 174 7 2019.09.20 –510.491 203
下载: 导出CSV
表 1 Parameters of the GF-3 dataset
Parameter Value Parameter Value Spatial baseline span 1417.4 m Number of scenes 7 Temporal baseline span 464 d Azimuth resolution 0.3626 m Slant range 1052747 m Range resolution 0.765692 m Wavelength 0.056 m Elevation resolution 20.6174 m Incident angle 47.2330015° Information resolution 21.8 mm/year
下载: 导出CSV
表 2 Spatio-temporal baseline parameters of the GF-3dataset
Number Time of acquisition Spatial baseline (m) Temporal baseline (d) 1 2018.06.13 –459.108 –261 2 2019.01.31 –628.551 –29 3 2019.03.01 0 0 4 2019.03.30 –724.517 29 5 2019.07.24 692.863 145 6 2019.08.22 –38.211 174 7 2019.09.20 –510.491 203
下载: 导出CSV
-
[1] CURLANDER J C and MCDONOUGH R N. Synthetic Aperture Radar: Systems and Signal Processing[M]. New York: Wiley, 1991. [2] BUDILLON A, JOHNSY A C, and SCHIRINZI G. Contextual information based SAR tomography of urban areas[C]. 2019 Joint Urban Remote Sensing Event (JURSE), Vannes, Italy, 2019: 1–4. doi: 10.1109/JURSE.2019.8809076. [3] REIGBER A, MOREIRA A, and PAPATHANASSIOU K P. First demonstration of airborne SAR tomography using multibaseline L-band data[C]. IEEE 1999 International Geoscience and Remote Sensing Symposium, Hamburg, Germany, 1999: 44–46. doi: 10.1109/IGARSS.1999.773395. [4] REIGBER A and MOREIRA A. First demonstration of airborne SAR tomography using multibaseline L-band data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2000, 38(5): 2142–2152. doi: 10.1109/36.868873 [5] FORNARO G and SERAFINO F. Imaging of single and double scatterers in urban areas via SAR tomography[J]. IEEE Transactions on Geoscience and Remote Sensing, 2006, 44(12): 3497–3505. doi: 10.1109/TGRS.2006.881748 [6] REALE D, FORNARO G, PAUCIULLO A, et al. Tomographic imaging and monitoring of buildings with very high resolution SAR data[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(4): 661–665. doi: 10.1109/LGRS.2010.2098845 [7] SHAHZAD M and ZHU Xiaoxiang. Reconstruction of building façades using spaceborne multiview TomoSAR point clouds[C]. 2013 IEEE International Geoscience and Remote Sensing Symposium - IGARSS, Melbourne, Australia, 2013: 624–627. doi: 10.1109/IGARSS.2013.6721234. [8] WANG Xiao, XU Feng, and JIN Yaqiu. The iterative reweighted alternating direction method of multipliers for separating structural layovers in SAR tomography[J]. IEEE Geoscience and Remote Sensing Letters, 2017, 14(11): 1883–1887. doi: 10.1109/LGRS.2017.2694825 [9] WANG Yuanyuan and ZHU Xiaoxiang. Robust nonlinear blind SAR tomography in urban areas[C]. 12th European Conference on Synthetic Aperture Radar, Aachen, Germany, 2018: 1–6. [10] 秦斐, 梁兴东, 张福博, 等. 基于机器学习的阵列层析SAR建筑物目标提取方法[J]. 信号处理, 2019, 35(2): 176–186. doi: 10.16798/j.issn.1003-0530.2019.02.003QIN Fei, LIANG Xingdong, ZHANG Fubo, et al. Building target extraction methods in array SAR tomography based on machine learning[J]. Journal of Signal Processing, 2019, 35(2): 176–186. doi: 10.16798/j.issn.1003-0530.2019.02.003 [11] LOMBARDINI F. Differential tomography: A new framework for SAR interferometry[C]. IGARSS 2003. 2003 IEEE International Geoscience and Remote Sensing Symposium, Toulouse, French, 2003: 1206–1208. doi: 10.1109/IGARSS.2003.1294059. [12] FORNARO G, PAUCIULLO A, and SERAFINO F. Multipass SAR processing for urbanized areas imaging and deformation monitoring at small and large scales[C]. 2007 Urban Remote Sensing Joint Event, Paris, 2007: 1–7. doi: 10.1109/URS.2007.371879. [13] FORNARO G, REALE D, and SERAFINO F. 4D SAR focusing: A tool for improved imaging and monitoring of urban areas[C]. IGARSS 2008 - 2008 IEEE International Geoscience and Remote Sensing Symposium, Boston, USA, 2008: V-475–V-478. doi: 10.1109/IGARSS.2008.4780132. [14] ZHU Xiaoxiang, ADAM N, and BAMLER R. Space-borne high resolution tomographic interferometry[C]. 2009 IEEE International Geoscience and Remote Sensing Symposium, Cape Town, South Africa, 2009: IV-869–IV-872. doi: 10.1109/IGARSS.2009.5417515. [15] FORNARO G, SERAFINO F, and REALE D. 4-D SAR imaging: The case study of Rome[J]. IEEE Geoscience and Remote Sensing Letters, 2010, 7(2): 236–240. doi: 10.1109/LGRS.2009.2032133 [16] ZHU Xiaoxiang and BAMLER R. Let’s do the time warp: Multicomponent nonlinear motion estimation in differential SAR tomography[J]. IEEE Geoscience and Remote Sensing Letters, 2011, 8(4): 735–739. doi: 10.1109/LGRS.2010.2103298 [17] SIDDIQUE M A, HAJNSEK I, AERSOSPACE G, et al. Investigating the combined use of differential SAR tomography and PSI for spatio-temporal inversion[C]. 2015 Joint Urban Remote Sensing Event (JURSE), Lausanne, Swiss, 2015: 1–4. doi: 10.1109/JURSE.2015.7120504. [18] WANG Zhigui and LIU Mei. Seasonal deformation and accelerated motion of infrastructure monitoring using a generalized differential SAR tomography[J]. IEEE Geoscience and Remote Sensing Letters, 2020, 17(4): 626–630. doi: 10.1109/LGRS.2019.2917324 [19] CANDÈS E J. Compressive sampling[C]. The International Congress of Mathematicians, Madrid, Spain, 2006: 1433–1452. [20] DONOHO D L. Compressed sensing[J]. IEEE Transactions on Information Theory, 2006, 52(4): 1289–1306. doi: 10.1109/TIT.2006.871582 [21] CANDÈS E J, ROMBERG J K, and TAO T. Stable signal recovery from incomplete and inaccurate measurements[J]. Communications on Pure and Applied Mathematics, 2006, 59(8): 1207–1223. doi: 10.1002/cpa.20124 [22] ZHU Xiaoxiang and BAMLER R. Tomographic SAR inversion by L1-Norm regularization—The compressive sensing approach[J]. IEEE Transactions on Geoscience and Remote Sensing, 2010, 48(10): 3839–3846. doi: 10.1109/TGRS.2010.2048117 [23] ZHU Xiaoxiang and BAMLER R. Super-resolution of sparse reconstruction for tomographic SAR imaging-demonstration with real data[C]. EUSAR 2012; 9th European Conference on Synthetic Aperture Radar, Nuremberg, USA, 2012: 215–218. [24] WEISS M, FORNARO G, and REALE D. Multi scatterer detection within tomographic SAR using a compressive sensing approach[C]. 2015 3rd International Workshop on Compressed Sensing Theory and its Applications to Radar, Sonar and Remote Sensing (CoSeRa), Pisa, 2015: 11–15. doi: 10.1109/CoSeRa.2015.7330254. [25] LI Xinwu, PENG Xing, and LIANG Lei. Sparse reconstruction-based SAR tomography and it’s application[C]. 2017 IEEE International Conference on Imaging Systems and Techniques (IST), Beijing, China, 2017: 1–5. doi: 10.1109/IST.2017.8261514. [26] ZHU Xiaoxiang and BAMLER R. Super-resolution for 4-D SAR tomography via compressive sensing[C]. 8th European Conference on Synthetic Aperture Radar, Aachen, Germany, 2010: 1–4. [27] ZHU Xiaoxiang and BAMLER R. Compressive sensing for high resolution differential SAR tomography - the SL1MMER algorithm[C]. 2010 IEEE International Geoscience and Remote Sensing Symposium, Honolulu, USA, 2010: 17–20. doi: 10.1109/IGARSS.2010.5654021. [28] LENG Ying, LI Fei, SHAO Yunfeng, et al. Elevation and deformation estimation based on TomoSAR[C]. EUSAR 2014; 10th European Conference on Synthetic Aperture Radar, Berlin, Germany, 2014: 1–4. [29] 姚天宇. 长征四号丙运载火箭成功发射高分三号卫星[J]. 中国航天, 2016(8): 8.YAO Tianyu. The Long March 4C carrier rocket successfully launched the GaoFen-3 satellite[J]. Aerospace China, 2016(8): 8. [30] 张庆君. 高分三号卫星总体设计与关键技术[J]. 测绘学报, 2017, 46(3): 269–277. doi: 10.11947/j.AGCS.2017.20170049ZHANG Qingjun. System design and key technologies of the GF-3 satellite[J]. Acta Geodaetica et Cartographica Sinica, 2017, 46(3): 269–277. doi: 10.11947/j.AGCS.2017.20170049 [31] 李胜阳, 许志辉, 陈子琪, 等. 高分3号卫星影像在黄河洪水监测中的应用[J]. 水利信息化, 2017(5): 22–26, 72.LI Shengyang, XU Zhihui, CHEN Ziqi, et al. Application of GF-3 satellite remote sensing image on Yellow River flood monitoring[J]. Water Resources Informatization, 2017(5): 22–26, 72. [32] 马建威, 孙亚勇, 陈德清, 等. 高分三号卫星在洪涝和滑坡灾害应急监测中的应用[J]. 航天器工程, 2017, 26(6): 161–166. doi: 10.3969/j.issn.1673-8748.2017.06.026MA Jianwei, SUN Yayong, CHEN Deqing, et al. Applications of GF-3 satellite in flood and landslide disasters emergency monitoring[J]. Spacecraft Engineering, 2017, 26(6): 161–166. doi: 10.3969/j.issn.1673-8748.2017.06.026 [33] 余博, 李如仁, 陈振炜, 等. 基于高分三号数据的干涉测量研究[J]. 测绘工程, 2019, 28(4): 19–23, 28.YU Bo, LI Ruren, CHEN Zhenwei, et al. Experimental research on interferometry based on GF-3 SAR data[J]. Engineering of Surveying and Mapping, 2019, 28(4): 19–23, 28. [34] 黄震, 赵双明, 常青. 国产高分三号卫星干涉测量试验[J]. 测绘地理信息, 2021, 46(4): 7–11. doi: 10.14188/j.2095-6045.2019083HUANG Zhen, ZHAO Shuangming, and CHANG Qing. Interferometric test using domestic GF-3 satellite[J]. Journal of Geomatics, 2021, 46(4): 7–11. doi: 10.14188/j.2095-6045.2019083 [35] FORNARO G, SERAFINO F, and SOLDOVIERI F. Three-dimensional focusing with multipass SAR data[J]. IEEE Transactions on Geoscience and Remote Sensing, 2003, 41(3): 507–517. doi: 10.1109/TGRS.2003.809934 [36] ZHU X. High-resolution spaceborne radar tomography[D]. Technische Universität München, 2008. [37] FORNARO G, REALE D, and SERAFINO F. Four-dimensional SAR imaging for height estimation and monitoring of single and double scatterers[J]. IEEE Transactions on Geoscience and Remote Sensing, 2009, 47(1): 224–237. doi: 10.1109/TGRS.2008.2000837 [38] PATI Y C, REZAIIFAR R, and KRISHNAPRASAD P S. Orthogonal matching pursuit: Recursive function approximation with applications to wavelet decomposition[C]. 27th Asilomar Conference on Signals, Systems and Computers, Pacific Grove, USA, 1993: 40–44. doi: 10.1109/ACSSC.1993.342465. [39] BI Hui and BI Guoan. Performance analysis of iterative soft thresholding algorithm for L1 regularization based sparse SAR imaging[C]. 2019 IEEE Radar Conference (RadarConf), Boston, USA, 2019: 1–6. doi: 10.1109/RADAR.2019.8835551. [40] 毕辉. 基于稀疏信号处理的SAR/TomoSAR成像方法研究[D]. [博士论文], 中国科学院大学, 2017: 33–41.BI Hui. Research on the imaging method of SAR/TomoSAR based on sparse signal processing[D]. [Ph. D. dissertation], University of Chinese Academy of Sciences, 2007: 33–41. [41] STOICA P and MOSES R. Spectral Analysis of Signals[M]. Upper Sadle River, New Jersey: Pearson Prentice Hall, 2005. [42] CAPON J. High-resolution frequency-wavenumber spectrum analysis[J]. Proceedings of the IEEE, 1969, 57(8): 1408–1418. doi: 10.1109/PROC.1969.7278 [43] SCHMIDT R. Multiple emitter location and signal parameter estimation[J]. IEEE Transactions on Antennas and Propagation, 1986, 34(3): 276–280. doi: 10.1109/TAP.1986.1143830 [44] STOICA P and NEHORAI A. MUSIC, maximum likelihood, and Cramer-Rao bound: Further results and comparisons[J]. IEEE Transactions on Acoustics,Speech,and Signal Processing, 1990, 38(12): 2140–2150. doi: 10.1109/29.61541 期刊类型引用(4)
1. Xiaolan Qiu，Zekun Jiao，Zhe Zhang，Qiancheng Yan，Chibiao Ding. Advances and prospects in SAR microwave vision three-dimensional imaging. National Science Open. 2024(05): 53-83 . 
必应学术2. 任子帅，张照，高雨欣，郭睿. 基于自适应高程约束的TomoSAR三维成像. 雷达学报. 2023(05): 1056-1068 . 
本站查看3. 郭华东，吴文瑾，张珂，李新武. 新型SAR对地环境观测. 测绘学报. 2022(06): 862-872 . 百度学术4. 林赟，张琳，韦立登，张汉卿，冯珊珊，王彦平，洪文. 无先验模型复杂结构设施SAR全方位三维成像方法研究. 雷达学报. 2022(05): 909-919 . 本站查看其他类型引用(13)
-
下载:
必应学术
计量
- 文章访问数: 3759
- HTML全文浏览量: 1253
- PDF下载量: 494
- 被引次数: 17
