CN1680826A - 生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法 - Google Patents
生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法 Download PDFInfo
- Publication number
- CN1680826A CN1680826A CNA2004100230989A CN200410023098A CN1680826A CN 1680826 A CN1680826 A CN 1680826A CN A2004100230989 A CNA2004100230989 A CN A2004100230989A CN 200410023098 A CN200410023098 A CN 200410023098A CN 1680826 A CN1680826 A CN 1680826A
- Authority
- CN
- China
- Prior art keywords
- window
- imaginary part
- real
- width
- aperture radar
- 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
本发明的方法是:对于两幅复数图像中对应的每一像素点,以当前点为中心,取一大小为m×n(m、n为大于0的数)的等值线窗口或其它形式窗口,在窗口内对两幅复图像数据的实部或虚部进行相关运算,求出其相关系数C1;再用上述相同方法对第一幅复图像数据的实部和第二幅复图像数据的虚部进行相关运算,求出其相关系数C2,或对第一幅复图像数据的虚部和第二幅复图像数据的实部进行相关运算,求出其相关系数-C2,进而得到C2;然后对C2与C1的比值求反正切即可获得干涉相位图像。本发明在干涉合成孔径雷达数据处理过程中,由两幅精确配准后的复图像生成完全免除相干斑噪声的干涉相位图。从而极大地提高相位解缠的可靠性及生成数字高程图的精度。
Description
技术领域
本发明涉及图像处理技术、遥感测量等领域,进一步是指用于生成干涉合成孔径雷达干涉相位图的条纹等值线窗口实虚部相关方法、矩形等窗口实虚部相关方法。
背景技术
干涉合成孔径雷达(Interferometric Synthetic Aperture Radar,缩写:InSAR或IFSAR)三维成像技术是新近发展起来的空间观测技术,它是传统的合成孔径雷达(Synthetic Aperture Radar,缩写:SAR)遥感技术与射电天文干涉技术结合的产物。干涉合成孔径雷达技术通过两幅天线同时观测或一副天线两次平行观测,接收地面目标的回波信号,经成像处理后得到同一观测区域两幅具有相干性的合成孔径雷达复数图像(包括强度信息和相位信息),再经过干涉处理后得到这两幅复图像的干涉相位图像,按照一定的几何关系进行变换,进而得到观测目标的三维信息。
目前干涉相位图的生成主要采用共轭相乘的方法。记V1(r,x)和V2(r,x)分别为精确配准后的两幅复图像的复值,则:
式中φ1和φ2分别为两根天线的回波相位,Δφ=φ2-φ1为回波相位差,即干涉相位。将两幅复图像对应像素进行复共轭相乘得:
式中*表示复共轭。由式(3)可得:
图1所示为用这种方法得到的干涉相位图。由于干涉合成孔径雷达系统中斑点噪声、空间失相关、时间失相关、数据处理噪声等因素的影响,使得干涉相位图的信噪比较低。这种高噪声的存在严重地影响了相位解缠的进行和高精度数字高程图的获取,已经成为干涉合成孔径雷达数据处理中一个最主要的瓶颈。
发明内容
本发明要解决的技术问题是,针对现有技术存在的缺陷,在干涉合成孔径雷达数据处理过程中,用条纹等值线窗口实虚部相关方法、矩形等窗口实虚部相关方法代替传统的共轭相乘方法,由两幅精确配准后的复图像生成完全免除相干斑噪声的干涉相位图。从而极大地提高相位解缠的可靠性及生成数字高程图的精度,解决了干涉合成孔径雷达数据处理中一个最主要的困难。
本发明的技术方案是,所述生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法是:
(1)对于两幅复数图像中对应的每一像素点,以当前点为中心,取一大小为m×n(m、n为大于0的数)的等值线窗口或其它形式窗口,在窗口内对两幅复图像数据的实部或虚部进行相关运算,求出其相关系数C1;
(2)用上述相同方法对第一幅复图像数据的实部和第二幅复图像数据的虚部进行相关运算,求出其相关系数C2,或对第一幅复图像数据的虚部和第二幅复图像数据的实部进行相关运算,求出其相关系数-C2,进而得到C2;
(3)对C2与C1的比值求反正切即可获得干涉相位图像。
所述窗口可以是条纹等值线窗口,也可以是矩形及其它形式窗口。
以下对本发明做出进一步说明。
关于本发明原理的公式推导。
使用的相关公式可以是协方差(均值归一化)相关、标准化协方差相关等多种数学表达形式。下面以标准化协方差相关公式为例来加以推导。
标准化协方差相关公式如下:
<·>m×n算符代表在m×n像素范围内对某一变量求均值。
对于两幅复图像V1(r,x)和V2(r,x)分别取其实部数据,记为:
f1=A1cosφ1 (6)
f2=A2cos(φ1+Δφ) (7)
由于图像精确配准,可以认为A1=A2=A,φ1和φ2为随机分布变量,由散斑统计理论,可以认为在满足一定尺度(m×n)的窗口上:
<cosφ1>m×n=<cos(φ1+Δφ)>m×n=0 (8)
并假设相位变化量Δφ在m×n个像素窗口区域内保持不变,由式(6)和式(7)可得:
<f1>m×n=0 (9)
<f2>m×n=0 (10)
(11)
将(12)(13)(14)代入(5)中,得
同理,分别取V1(r,x)和V2(r,x)的实部与虚部,记为:
f1=A1cosφ1 (16)
f2=A2sin(φ1+Δφ) (17)
用上面同样的方法进行推导可以得到:
由式(15)与式(18),可以得到:
Δφ=arctan(C2/C1) (19)
对比式(4)与式(19)可知二者得到干涉相位图的相位分布完全一致。
另外,将第一幅图的虚部与第二幅图的虚部进行相关运算也可以得到C1,将第一幅图的虚部与第二幅图的实部进行相关运算可以得到-C2,进而得到C2。
关于本发明的相关窗口选取。
由上述推导可知,上述结论是在条纹等值线窗口上才严格成立。因此,相关运算的最优窗口是条纹等值线窗口,用它可以得到质量和精度最好的干涉相位图。也可以选取矩形及其它形式窗口,仍可得到比较好的近似结果。条纹等值线窗口是一种曲线窗口,它是沿干涉相位图的条纹走向取一定长度和宽度所形成的窗口,可以从共轭相乘方法生成的干涉相位图得到,也可以从矩形等窗口相关方法生成的干涉相位图得到。可以用不同算法获取条纹等值线窗口,它不属本发明的内容,此处不赘述。
由以上可知,本发明为生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,它是一种用基于条纹等值线窗口、矩形等窗口的复图像实虚部相关方法生成无干涉斑干涉合成孔径雷达干涉相位图的方法;本发明方法在干涉合成孔径雷达数据处理过程中,用条纹等值线窗口实虚部相关方法、矩形等窗口实虚部相关方法代替传统的共轭相乘方法,由两幅精确配准后的复图像生成完全免除相干斑噪声的干涉相位图。从而极大地提高相位解缠的可靠性及生成数字高程图的精度,解决了干涉合成孔径雷达数据处理中一个最主要的困难。
附图说明
图1是精确配准后两幅复图像用共轭相乘方法生成的干涉相位图,图像大小为2048×2048像素。
图2是精确配准后两幅复图像用矩形窗口相关方法生成的干涉相位图,所用矩形窗口大小为19×19。
图3是精确配准后两幅复图像用条纹等值线窗口相关方法生成的干涉相位图,所用条纹等值线窗口大小为41×5。
具体实施方式
1.获取具有相干性的两幅SAR复图像数据,并将其精确配准;
2.对于两幅复数图像中对应的每一像素点,以当前点为中心,取一大小为m×n的窗口(可以是条纹等值线窗口,也可以是矩形及其它形式窗口),在窗口内对两幅复图像数据的实部(或虚部)进行相关运算,求出其相关系数C1;
3.用相同的方法对第一幅复图像数据的实部和第二幅复图像数据的虚部进行相关运算,求出其相关系数C2,或是对第一幅复图像数据的虚部和第二幅复图像数据的实部进行相关运算,求出其相关系数-C2,进而得到C2;
4.对C2与C1的比值求反正切即可获得干涉相位图像。
Claims (5)
1.生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,其特征是,该方法为:
(1)对于两幅复数图像中对应的每一像素点,以当前点为中心,取一大小为m×n(m、n为大于0的数)的窗口,在窗口内对两幅复图像数据的实部或虚部进行相关运算,求出其相关系数C1;
(2)用上述相同方法对第一幅复图像数据的实部和第二幅复图像数据的虚部进行相关运算,求出其相关系数C2,或对第一幅复图像数据的虚部和第二幅复图像数据的实部进行相关运算,求出其相关系数-C2,进而得到C2;
(3)对C2与C1的比值求反正切即可获得干涉相位图像。
2.根据权利要求1所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,其特征是,所述窗口是条纹等值线窗口或矩形及其它形式窗口。
3.根据权利要求2所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,其特征是,所述条纹等值线窗口为曲线窗口,它是沿干涉相位图的条纹走向取一定长度和宽度所形成的窗口。
4.根据权利要求3所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,其特征是,所述条纹等值线窗口由传统复图像共轭相乘方法生成的干涉相位图得到。
5.根据权利要求3所述的生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法,其特征是,所述条纹等值线窗口由矩形及其它形式窗口实虚部相关方法生成的干涉相位图得到。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100230989A CN1307428C (zh) | 2004-04-09 | 2004-04-09 | 生成无干涉斑合成孔径雷达干涉相位图的方法 |
PCT/CN2005/000346 WO2005098473A1 (fr) | 2004-04-09 | 2005-03-21 | Procede de correlation d'imagerie reelle permettant la generation d'un interferogramme sans tache de rso |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CNB2004100230989A CN1307428C (zh) | 2004-04-09 | 2004-04-09 | 生成无干涉斑合成孔径雷达干涉相位图的方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN1680826A true CN1680826A (zh) | 2005-10-12 |
CN1307428C CN1307428C (zh) | 2007-03-28 |
Family
ID=35067481
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CNB2004100230989A Expired - Fee Related CN1307428C (zh) | 2004-04-09 | 2004-04-09 | 生成无干涉斑合成孔径雷达干涉相位图的方法 |
Country Status (2)
Country | Link |
---|---|
CN (1) | CN1307428C (zh) |
WO (1) | WO2005098473A1 (zh) |
Cited By (9)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007071113A1 (fr) * | 2005-12-22 | 2007-06-28 | Qifeng Yu | Procede interferometrique de correlation d'isolignes en trois images pour la generation d'un interferogramme de radar |
CN102073992A (zh) * | 2010-12-09 | 2011-05-25 | 国网电力科学研究院 | 一种高分辨率sar卫星图像相干斑去噪方法 |
CN102323570A (zh) * | 2011-05-24 | 2012-01-18 | 中国人民解放军国防科学技术大学 | 一种雷达目标回波信号模拟器的幅相特性估计方法 |
CN103308914A (zh) * | 2012-12-06 | 2013-09-18 | 北京理工大学 | 一站固定双基地干涉sar处理方法 |
CN103871030A (zh) * | 2014-02-17 | 2014-06-18 | 中国科学院电子学研究所 | 一种干涉图像的滤波方法及设备 |
CN104463874A (zh) * | 2014-12-12 | 2015-03-25 | 中国人民解放军装备学院 | 一种自适应复数词典和稀疏编码的干涉相位图像降噪方法 |
CN108663678A (zh) * | 2018-01-29 | 2018-10-16 | 西北农林科技大学 | 基于混合整数优化模型的多基线InSAR相位解缠算法 |
CN109708584A (zh) * | 2018-12-30 | 2019-05-03 | 南京航空航天大学 | 一种油膜干涉法条纹间距图像识别方法 |
CN112098997A (zh) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | 三维全息成像安检雷达图像异物检测方法 |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110109106A (zh) * | 2019-04-23 | 2019-08-09 | 中国电力科学研究院有限公司 | 一种InSAR地形复杂区域干涉相位解缠方法 |
CN111461999B (zh) * | 2020-03-13 | 2023-02-14 | 西安工程大学 | 一种基于超像素相似性测量的sar图像相干斑抑制方法 |
US11933883B2 (en) | 2021-09-24 | 2024-03-19 | Aloft Sensing, Inc. | System and method for self-contained high-precision navigation |
CN114910909B (zh) * | 2022-06-17 | 2024-06-21 | 苏州大学 | 数据驱动下的合成孔径雷达干涉图相位恢复方法及系统 |
Family Cites Families (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US5677693A (en) * | 1994-05-26 | 1997-10-14 | Hughes Aircraft Company | Multi-pass and multi-channel interferometric synthetic aperture radars |
JP2877106B2 (ja) * | 1996-11-18 | 1999-03-31 | 日本電気株式会社 | アロングトラック・インターフェロメトリsar |
ITMI991154A1 (it) * | 1999-05-25 | 2000-11-25 | Milano Politecnico | Procedimento per misure radar di spostamento di aere urbane e zone franose |
JP4583619B2 (ja) * | 2000-09-13 | 2010-11-17 | 富士フイルム株式会社 | 縞画像解析誤差検出方法および縞画像解析誤差補正方法 |
CN1148067C (zh) * | 2001-06-01 | 2004-04-28 | 中国科学院计算技术研究所 | 合成孔径雷达复数图像数据压缩方法 |
CN1173291C (zh) * | 2001-06-26 | 2004-10-27 | 中国科学院遥感应用研究所 | 利用计算机重建含有斑点噪声的图像的方法 |
GB0219645D0 (en) * | 2002-08-23 | 2002-10-02 | Secr Defence | Radar imaging |
-
2004
- 2004-04-09 CN CNB2004100230989A patent/CN1307428C/zh not_active Expired - Fee Related
-
2005
- 2005-03-21 WO PCT/CN2005/000346 patent/WO2005098473A1/zh active Application Filing
Cited By (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2007071113A1 (fr) * | 2005-12-22 | 2007-06-28 | Qifeng Yu | Procede interferometrique de correlation d'isolignes en trois images pour la generation d'un interferogramme de radar |
CN102073992A (zh) * | 2010-12-09 | 2011-05-25 | 国网电力科学研究院 | 一种高分辨率sar卫星图像相干斑去噪方法 |
CN102323570A (zh) * | 2011-05-24 | 2012-01-18 | 中国人民解放军国防科学技术大学 | 一种雷达目标回波信号模拟器的幅相特性估计方法 |
CN102323570B (zh) * | 2011-05-24 | 2013-03-13 | 中国人民解放军国防科学技术大学 | 一种雷达目标回波信号模拟器的幅相特性估计方法 |
CN103308914A (zh) * | 2012-12-06 | 2013-09-18 | 北京理工大学 | 一站固定双基地干涉sar处理方法 |
CN103871030A (zh) * | 2014-02-17 | 2014-06-18 | 中国科学院电子学研究所 | 一种干涉图像的滤波方法及设备 |
CN104463874A (zh) * | 2014-12-12 | 2015-03-25 | 中国人民解放军装备学院 | 一种自适应复数词典和稀疏编码的干涉相位图像降噪方法 |
CN108663678A (zh) * | 2018-01-29 | 2018-10-16 | 西北农林科技大学 | 基于混合整数优化模型的多基线InSAR相位解缠算法 |
CN109708584A (zh) * | 2018-12-30 | 2019-05-03 | 南京航空航天大学 | 一种油膜干涉法条纹间距图像识别方法 |
CN112098997A (zh) * | 2020-09-18 | 2020-12-18 | 欧必翼太赫兹科技(北京)有限公司 | 三维全息成像安检雷达图像异物检测方法 |
Also Published As
Publication number | Publication date |
---|---|
CN1307428C (zh) | 2007-03-28 |
WO2005098473A1 (fr) | 2005-10-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN1680826A (zh) | 生成无干涉斑合成孔径雷达干涉相位图的实虚部相关方法 | |
CN100573189C (zh) | 基于多特显点的干涉式逆合成孔径雷达成像方法 | |
CN106249236B (zh) | 一种星载InSAR长短基线图像联合配准方法 | |
CN104101864B (zh) | 基于eof分解的导航x波段雷达海浪参数反演算法 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN110531355B (zh) | 一种星载合成孔径雷达实时统一成像方法 | |
CN103605119B (zh) | 一种星载合成孔径雷达条带模式方位模糊抑制方法 | |
CN107942327A (zh) | 基于脉冲相位编码的单通道hrws‑sar成像方法 | |
CN109270527B (zh) | 圆迹sar子孔径图像序列联合相关dem提取方法 | |
CN113204023B (zh) | 联合ps目标与ds目标的双极化相位优化地表形变监测方法 | |
CN112764116A (zh) | 一种稀疏阵列稀疏频点平面扫描体制快速成像方法 | |
Sletten et al. | Maritime signature correction with the NRL multichannel SAR | |
CN1808172A (zh) | 机载干涉合成孔径雷达原始回波生成方法 | |
CN103675815A (zh) | 一种在大斜视sar成像模式下对多普勒调频率进行精确估计的方法 | |
CN111273218A (zh) | 一种基于多层互质阵列的相干源波达方向估计方法 | |
Rao et al. | Satellite image fusion using fast discrete curvelet transforms | |
CN110554377B (zh) | 基于多普勒中心偏移的单通道sar二维流场反演方法及系统 | |
CN110618409B (zh) | 顾及叠掩及阴影的多通道InSAR干涉图仿真方法及系统 | |
Zhao et al. | Improved maximum likelihood estimation for optimal phase history retrieval of distributed scatterers in InSAR stacks | |
CN1415975A (zh) | 单次飞行全极化合成孔径雷达图像反演地面数字高程的方法 | |
Yang et al. | Phase imbalance estimation for azimuth multi-channel ScanSAR system | |
CN108919265B (zh) | 一种基于互信息量最大的星载TomoSAR空间基线序列设计方法 | |
CN113835089A (zh) | 一种基于差频InSAR的缓变高程反演方法 | |
CN111638515A (zh) | 基于双频联合处理技术的太赫兹频段sar运动补偿算法 | |
CN115015928B (zh) | 基于总功率同散射机制的高效极化时序InSAR方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
C14 | Grant of patent or utility model | ||
GR01 | Patent grant | ||
CF01 | Termination of patent right due to non-payment of annual fee | ||
CF01 | Termination of patent right due to non-payment of annual fee |
Granted publication date: 20070328 Termination date: 20200409 |