CN112859077A - 一种多级合成孔径雷达干涉相位解缠方法 - Google Patents
一种多级合成孔径雷达干涉相位解缠方法 Download PDFInfo
- Publication number
- CN112859077A CN112859077A CN202110111087.XA CN202110111087A CN112859077A CN 112859077 A CN112859077 A CN 112859077A CN 202110111087 A CN202110111087 A CN 202110111087A CN 112859077 A CN112859077 A CN 112859077A
- Authority
- CN
- China
- Prior art keywords
- synthetic aperture
- aperture radar
- unwrapping
- phase
- differential
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
Images
Classifications
-
- G—PHYSICS
- G01—MEASURING; TESTING
- G01S—RADIO DIRECTION-FINDING; RADIO NAVIGATION; DETERMINING DISTANCE OR VELOCITY BY USE OF RADIO WAVES; LOCATING OR PRESENCE-DETECTING BY USE OF THE REFLECTION OR RERADIATION OF RADIO WAVES; ANALOGOUS ARRANGEMENTS USING OTHER WAVES
- G01S13/00—Systems using the reflection or reradiation of radio waves, e.g. radar systems; Analogous systems using reflection or reradiation of waves whose nature or wavelength is irrelevant or unspecified
- G01S13/88—Radar or analogous systems specially adapted for specific applications
- G01S13/89—Radar or analogous systems specially adapted for specific applications for mapping or imaging
- G01S13/90—Radar or analogous systems specially adapted for specific applications for mapping or imaging using synthetic aperture techniques, e.g. synthetic aperture radar [SAR] techniques
- G01S13/9021—SAR image post-processing techniques
- G01S13/9023—SAR image post-processing techniques combined with interferometric techniques
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种多级合成孔径雷达干涉相位解缠方法,包括步骤:(1)生成合成孔径雷达干涉影像;(2)应用外部数字高程模型生成合成孔径雷达差分干涉图;(3)利用最小费用流法对合成孔径雷达差分干涉图进行相位解缠,得到初始的合成孔径雷达解缠差分干涉图;(4)根据不同的相干性阈值,应用枝切法对初始的合成孔径雷达解缠差分干涉图进行迭代解缠,得到最终的合成孔径雷达解缠差分干涉图;(5)补偿地形相位,得到最终的合成孔径雷达解缠干涉图。所述方法能够提高现有的单一的合成孔径雷达干涉相位解缠的成功率,可在星载和机载干涉合成孔径雷达***对地表三维重建和地物高度提取中起到重要的作用,具有较佳的技术效果。
Description
技术领域
本发明涉及一种遥感影像的数字摄影测量处理方法,特别是涉及合成孔径雷达干涉相位解缠方法。
背景技术
干涉合成孔径雷达干涉测量(Interferometric Synthetic Aperture Radar,InSAR)作为一种主动的微波遥感技术在地形测绘已经显现出越来越大的优势。合成孔径雷达(Synthetic Aperture Radar,SAR)影像的干涉相位信息是提取地形和高度关键信息。在合成孔径雷达干涉测量中,合成孔径雷达干涉相位解缠(简称解缠)是把缠绕的干涉相位恢复到真实的干涉相位,是整个处理链中的关键和重要步骤。在地形和地物条件复杂的区域,由于受地形变化和植被的影响,极易造成干涉相位的去相干,在干涉影像中的某些区域表现出很低的相干性。这些低相干(相干系数小于0.2)区域在单一的解缠过程中解缠错误,造成了整幅合成孔径雷达干涉影像的解缠失败。为了保障在中高相干(相干系数大于0.3)区域解缠成功,需要发展一种多级的合成孔径雷达干涉相位解缠方法。
一般来说,合成孔径雷达干涉相位解缠方法主要有:(1)一类方法是经典的路径跟踪算法——枝切法,该方法用像元集形成的最小闭合路径积分检测出二维相位数据中的所有残数点,用枝切使得所有的积分路径不包含未平衡的残数点。这种方法的计算速度很快,解缠相位与缠绕相位只会相差整数个相位周期,不会破坏真实相位的主值;它最大的缺点是在相干性低,残数点密集的区域,很难正确地设置枝切,经常得不到准确的解,甚至枝切会形成一些孤立区域,在这些区域内无法得到解。(2)另一类方法称为网络流法,比较典型的是最小费用流法,它的主要思想是将解缠相位的导数和与缠绕相位的导数之间的差异最小化,它还可以将整个误差限制在一个小范围内,防止误差的再传递,解缠结果较精确,计算效率也相对较高。最小费用流法虽然解缠的成功率很高,但对于在地形和地物复杂的区域,也容易出现解缠错误,从而导致全局性的解缠失败。
发明内容
本发明提供了一种多级合成孔径雷达干涉相位解缠方法,通过控制相干系数阈值,根据二种合成孔径雷达干涉相位解缠方法(最小费用流法和枝切法)各自的特点,联合这两种解缠方法实施多次迭代解缠处理,在地形和地物复杂区域能够实现高精度的合成孔径雷达干涉相位解缠。
一种多级合成孔径雷达干涉相位解缠方法,包括以下步骤:
S1:生成合成孔径雷达干涉影像,包括干涉图和相干系数图;
S2:利用外部DEM生成合成孔径雷达差分干涉图;
S3:应用最小费用流法对合成孔径雷达差分干涉图解缠,得到初始的合成孔径雷达解缠差分干涉图;
S4:根据不同的相干系数阈值,利用枝切法对初始的合成孔径雷达解缠差分干涉图进行迭代解缠,得到最终的合成孔径雷达解缠差分干涉图;
S5:针对合成孔径雷达解缠差分干涉图补偿地形相位,得到合成孔径雷达解缠干涉图。
进一步的,步骤S1中,根据获取的合成孔径雷达影像对进行共轭相乘和相干系数计算,从而生成合成孔径雷达干涉影像,其中,所述合成孔径雷达影像对为:
其中,主影像和从影像共轭相乘生成干涉图为:
其中,通过相干性计算得到相干系数图为:
γ表示相干系数,E表示均值计算。
进一步的,步骤S2中,利用外部DEM生成高程相位,与生成的干涉图进行差分运算,得到合成孔径雷达差分干涉相位,进而生成合成孔径雷达差分干涉图,其中,外部DEM的高程相位:
其中,差分运算生成干涉图外部DEM的高程相位为:
进一步的,步骤S3中,选取设定的相干系数阈值,应用最小费用流法的相位解缠算法对合成孔径雷达差分干涉图进行解缠,得到初始的合成孔径雷达解缠差分干涉相位,从而形成初始的合成孔径雷达解缠差分干涉图;其中,初始的合成孔径雷达解缠差分干涉相位为:
ψMCF表示应用最小费用流法解缠后的初始的合成孔径雷达解缠差分干涉相位,W{.}代表解缠运算符,n表示解缠后的整周期数,γ表示相干系数,γh表示设定的相干系数阈值,所述合成孔径雷达解缠差分干涉相位ψMCF是针对相干系数大于γh的合成孔径雷达差分干涉相位所作的解缠运算结果。
进一步的,步骤S4中,以γ作为相干系数的步进数,逐级递减相干系数阈值,应用枝切法对初始的合成孔径雷达解缠差分干涉相位进行迭代解缠,直到迭代到一个设定的相干系数阈值γl,结束迭代,得到迭代解缠后的合成孔径雷达解缠差分干涉相位,从而形成最终的合成孔径雷达解缠差分干涉图;
其中,合成孔径雷达解缠差分干涉相位为:
表示应用枝切法每次迭代解缠后的合成孔径雷达解缠差分干涉相位,m表示迭代次数,γT表示每次迭代解缠所运用的不同的相干系数阈值,γ表示相干系数的步进数,γl表示设定的相干系数阈值,所述合成孔径雷达解缠差分干涉相位是针对相干系数大于γT的合成孔径雷达差分干涉相位所作的解缠运算结果,ψBR为每次迭代解缠后的合成孔径雷达解缠差分干涉相位的集合。
进一步的,步骤S5中,将不同相干系数阈值条件下的合成孔径雷达解缠差分相位ψMCF和ψBR进行地形相位补偿,
得到合成孔径雷达解缠干涉相位ψunw。
本发明提供的多级合成孔径雷达干涉相位解缠方法以后,实现了在地形和地物复杂区域的高精度的合成孔径雷达干涉相位解缠,保障了中高相干区域的解缠成功。多级合成孔径雷达干涉相位解缠方法在基于合成孔径雷达干涉测量技术的地表三维重建和地物高度提取起到重要作用,具有较佳的技术效果。
附图说明
下面结合附图对本发明进行进一步详细的描述,以使得本发明的上述优点更加明确。
图1是本发明所述的多级合成孔径雷达干涉相位解缠方法的流程图;
图2是应用单一的最小费用流法解缠的合成孔径雷达解缠差分干涉图;
图3是应用最小费用流法解缠的初始的合成孔径雷达解缠差分干涉图;
图4是应用枝切法迭代解缠后的合成孔径雷达解缠差分干涉图。
具体实施方式
如图1所示,本发明提供的多级合成孔径雷达干涉相位解缠方法,通过控制相干系数阈值,联合二种合成孔径雷达干涉相位解缠方法(最小费用流法和枝切法)多次迭代解缠,适用于星载和机载合成孔径雷达干涉影像,在地形和地物复杂区域能够实现高精度的合成孔径雷达干涉相位解缠。本发明包括以下步骤:
S1:首先根据获取的合成孔径雷达影像对,进行共轭相乘和相干系数计算,从而生成合成孔径雷达干涉影像,所述合成孔径雷达干涉影像包括干涉图和相干系数图。其中,合成孔径雷达影像对可表示为:
通过公式(2),主影像和从影像共轭相乘生成干涉图:
通过公式(3),通过相干性计算得到相干系数图。
γ表示相干系数,E表示均值计算。
S2:根据外部的DEM(DEM:全称为Digital Elevation Model,中文名称为:数字高程模型)生成高程相位,与生成的干涉图进行差分运算,得到合成孔径雷达差分干涉相位,进而生成合成孔径雷达差分干涉图。
通过公式(4),得到外部DEM的高程相位:
公式(4)中,表示外部DEM的高程相位,B⊥为合成孔径雷达***的垂直基线,λ表示合成孔径雷达***的波长,R表示合成孔径雷达***的斜距,θ表示合成孔径雷达***的入射角,Δz表示外部DEM中相邻像元的高度差。
通过公式(5),生成干涉图外部DEM的高程相位:
通过合成孔径雷达差分干涉相位,进而生成合成孔径雷达差分干涉图。
S3:选取较高的相干系数阈值,所述相干系数阈值为0.8以上,应用最小费用流法(MCF,The Minimum Cost Flow)的相位解缠算法对合成孔径雷达差分干涉图进行解缠,得到初始的合成孔径雷达解缠差分干涉相位,如公式(6)所示,从而形成初始的合成孔径雷达解缠差分干涉图。
公式(6)中,ψMCF表示应用最小费用流法解缠后的合成孔径雷达解缠差分干涉相位,W{.}代表解缠运算符,n表示解缠后的整周期数,γ表示相干系数,γh表示较高的相干系数阈值,公式(6)中计算的合成孔径雷达解缠差分干涉相位ψMCF是针对相干系数大于γh的合成孔径雷达差分干涉相位所作的解缠运算结果。
S4:以γ作为相干系数的步进数,逐级递减相干系数阈值,应用枝切法(BR,BranchCut)对初始的合成孔径雷达解缠差分干涉相位进行迭代解缠,如公式(7)所示,直到迭代到一个较低的相干系数阈值γl(0.2<γl<0.4),结束迭代,得到迭代解缠后的合成孔径雷达解缠差分干涉相位ψBR。ψMCF和ψBR共同组成了不同相干系数阈值条件下合成孔径雷达解缠差分干涉相位,从而形成最终的合成孔径雷达解缠差分干涉图。
公式(7)中,表示应用枝切法每次迭代解缠后的合成孔径雷达解缠差分干涉相位,m表示迭代次数,γT表示每次迭代解缠所运用的不同的相干系数阈值,γ表示相干系数的步进数,γl表示较低的相干系数阈值,(7)中计算的合成孔径雷达解缠差分干涉相位是针对相干系数大于γT的合成孔径雷达差分干涉相位所作的解缠运算结果,ψBR为每次迭代解缠后的合成孔径雷达解缠差分干涉相位的集合。
S5:针对合成孔径雷达解缠差分干涉图补偿地形相位,得到合成孔径雷达解缠干涉相位,形成合成孔径雷达解缠干涉图。
通过公式(8),将公式(6)和(7)计算得到的不同相干系数阈值条件下的合成孔径雷达解缠差分相位ψMCF和ψBR进行地形相位补偿,从而形成合成孔径雷达解缠干涉相位。
公式(8)中,ψunw表示合成孔径雷达解缠干涉相位。
利用星载合成孔径雷达干涉影像进行相位解缠实验,图2是实施例中的应用单一的最小费用流法解缠的合成孔径雷达解缠差分干涉图(相干系数阈值为0.3),从图中可以看出,该结果存在全局性的相位跳变,导致整幅影像解缠失败。
图3为应用最小费用流法的初始的合成孔径雷达解缠差分干涉图(相干系数阈值为0.85),图4为应用枝切法对初始的合成孔径雷达解缠差分干涉相位进行迭代解缠后的合成孔径雷达解缠差分干涉图(相干系数的步进数为0.05,相干系数阈值为0.3)。
可以看出,通过实施多级合成孔径雷达干涉相位解缠方法后,不存在全局性的解缠错误,在中高相干性区域,保留了正确的合成孔径雷达解缠干涉相位值,能够提升合成孔径雷达干涉测量技术在地形和高度测量中的精度。
上述具体实施例仅仅是示例性的,在本发明的上述教导下,本领域技术人员可以在上述实施例的基础上进行各种改进和变形,而这些改进或者变形落在本发明的保护范围内。本领域技术人员应该明白,上面的具体描述只是为了解释本发明的目的,并非用于限制本发明。本发明的保护范围由权利要求及其等同物限定。
Claims (6)
1.一种多级合成孔径雷达干涉相位解缠方法,包括以下步骤:
S1:生成合成孔径雷达干涉影像,包括干涉图和相干系数图;
S2:利用外部DEM生成合成孔径雷达差分干涉图;
S3:应用最小费用流法对合成孔径雷达差分干涉图解缠,得到初始的合成孔径雷达解缠差分干涉图;
S4:根据不同的相干系数阈值,利用枝切法对初始的合成孔径雷达解缠差分干涉图进行迭代解缠,得到最终的合成孔径雷达解缠差分干涉图;
S5:针对合成孔径雷达解缠差分干涉图补偿地形相位,得到合成孔径雷达解缠干涉图。
5.根据权利要求4所述的多级合成孔径雷达干涉相位解缠方法,其特征在于:步骤S4中,以γ作为相干系数的步进数,逐级递减相干系数阈值,应用枝切法对初始的合成孔径雷达解缠差分干涉相位进行迭代解缠,直到迭代到一个设定的相干系数阈值γl,结束迭代,得到迭代解缠后的合成孔径雷达解缠差分干涉相位,从而形成最终的合成孔径雷达解缠差分干涉图;
其中,合成孔径雷达解缠差分干涉相位为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110111087.XA CN112859077B (zh) | 2021-01-27 | 2021-01-27 | 一种多级合成孔径雷达干涉相位解缠方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110111087.XA CN112859077B (zh) | 2021-01-27 | 2021-01-27 | 一种多级合成孔径雷达干涉相位解缠方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112859077A true CN112859077A (zh) | 2021-05-28 |
CN112859077B CN112859077B (zh) | 2023-03-07 |
Family
ID=76009550
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110111087.XA Active CN112859077B (zh) | 2021-01-27 | 2021-01-27 | 一种多级合成孔径雷达干涉相位解缠方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112859077B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113624122A (zh) * | 2021-08-10 | 2021-11-09 | 中咨数据有限公司 | 融合GNSS数据与InSAR技术的桥梁变形监测方法 |
CN117633494A (zh) * | 2023-11-20 | 2024-03-01 | 中国矿业大学 | 一种基于awc-lstm模型的煤矿地表形变预测方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
CN105158761A (zh) * | 2015-08-31 | 2015-12-16 | 西安电子科技大学 | 基于枝切法和曲面拟合的雷达合成相位解缠方法 |
CN108663678A (zh) * | 2018-01-29 | 2018-10-16 | 西北农林科技大学 | 基于混合整数优化模型的多基线InSAR相位解缠算法 |
CN110109106A (zh) * | 2019-04-23 | 2019-08-09 | 中国电力科学研究院有限公司 | 一种InSAR地形复杂区域干涉相位解缠方法 |
CN110161501A (zh) * | 2019-05-24 | 2019-08-23 | 电子科技大学 | 一种多时序sar图像的目标区域地表起伏信息提取方法 |
CN112051571A (zh) * | 2020-08-14 | 2020-12-08 | 中冶(北京)交通科技发展有限公司 | 一种新型差分InSAR的LOS向形变量估计方法 |
-
2021
- 2021-01-27 CN CN202110111087.XA patent/CN112859077B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
JP2004191053A (ja) * | 2002-12-06 | 2004-07-08 | Mitsubishi Electric Corp | 合成開口レーダ装置及び数値標高モデル作成方法 |
CN105158761A (zh) * | 2015-08-31 | 2015-12-16 | 西安电子科技大学 | 基于枝切法和曲面拟合的雷达合成相位解缠方法 |
CN108663678A (zh) * | 2018-01-29 | 2018-10-16 | 西北农林科技大学 | 基于混合整数优化模型的多基线InSAR相位解缠算法 |
CN110109106A (zh) * | 2019-04-23 | 2019-08-09 | 中国电力科学研究院有限公司 | 一种InSAR地形复杂区域干涉相位解缠方法 |
CN110161501A (zh) * | 2019-05-24 | 2019-08-23 | 电子科技大学 | 一种多时序sar图像的目标区域地表起伏信息提取方法 |
CN112051571A (zh) * | 2020-08-14 | 2020-12-08 | 中冶(北京)交通科技发展有限公司 | 一种新型差分InSAR的LOS向形变量估计方法 |
Non-Patent Citations (6)
Title |
---|
李泓宇等: "引入质量图指导的改进的枝切法", 《国外电子测量技术》 * |
杨成生: "基于D-InSAR技术的煤矿沉陷监测", 《中国优秀硕士学位论文全文数据库工程科技Ⅰ辑》 * |
段建华等: "InSAR两种相位解缠算法的融合", 《微计算机信息》 * |
毕海霞等: "基于区域识别和区域扩展的相位解缠算法", 《电波科学学报》 * |
王紫燕等: "一种组合的InSAR数据的相位解缠算法", 《自然科学进展》 * |
王霖郁等: "一种枝切法和质量图相结合的InSAR相位解缠算法", 《应用科技》 * |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN113624122A (zh) * | 2021-08-10 | 2021-11-09 | 中咨数据有限公司 | 融合GNSS数据与InSAR技术的桥梁变形监测方法 |
CN117633494A (zh) * | 2023-11-20 | 2024-03-01 | 中国矿业大学 | 一种基于awc-lstm模型的煤矿地表形变预测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN112859077B (zh) | 2023-03-07 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN109633648B (zh) | 一种基于似然估计的多基线相位估计装置及方法 | |
CN107193005B (zh) | 一种无损卡尔曼滤波与粒子滤波相结合的相位展开算法 | |
CN111273293B (zh) | 一种顾及地形起伏的InSAR残余运动误差估计方法及装置 | |
CN112859077B (zh) | 一种多级合成孔径雷达干涉相位解缠方法 | |
CN109696675B (zh) | 基于迪杰斯特拉算法的InSAR时序影像集合配准方法 | |
CN106093939B (zh) | 一种基于相位差统计模型的InSAR图像相位解缠方法 | |
KR101712084B1 (ko) | 멀티간섭기법 기반의 이온왜곡 보정 방법 및 그 장치 | |
CN114114181B (zh) | 基于轨道误差相位基的星载sar干涉基线矫正方法 | |
CN112711021B (zh) | 一种多分辨率InSAR交互干涉时序分析方法 | |
CN110018476B (zh) | 时间差分基线集时序干涉sar处理方法 | |
CN109509219B (zh) | 基于最小生成树的InSAR时序影像集合的配准方法 | |
CN113311433B (zh) | 一种质量图和最小费用流结合的InSAR干涉相位两步解缠方法 | |
CN113589286B (zh) | 基于D-LinkNet的无迹卡尔曼滤波相位解缠方法 | |
CN111856459A (zh) | 一种改进的DEM最大似然约束多基线InSAR相位解缠方法 | |
CN115272683A (zh) | 一种基于深度学习的中心差分信息滤波相位解缠方法 | |
CN104849714A (zh) | 一种高效稳健的多基线ukf相位展开方法 | |
CN111965646B (zh) | 星载雷达数据处理方法、装置、设备及存储介质 | |
CN106707283B (zh) | 基于无味信息滤波的相位展开算法 | |
CN116224327B (zh) | 一种基于学习网络的矿区大梯度形变区相位解缠方法 | |
CN110146885B (zh) | 基于多卫星平台MT-InSAR融合的新成陆区超长形变时序提取方法 | |
CN115616575B (zh) | 一种星载sar立体测量辅助的干涉相位图解缠绕方法 | |
CN113341410A (zh) | 一种大范围林下地形估计方法、装置、设备及介质 | |
Lombardi et al. | Accuracy of high resolution CSK interferometric Digital Elevation Models | |
Hai et al. | Experimental DEM extraction from ASTER stereo pairs and 3D registration based on ICESat laser altimetry data in upstream area of Lambert Glacier, Antarctica | |
CN118011344B (zh) | 一种外部dem辅助的sar干涉定标方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |