CN102621549A - 多基线/多频段干涉相位解缠频域快速算法 - Google Patents
多基线/多频段干涉相位解缠频域快速算法 Download PDFInfo
- Publication number
- CN102621549A CN102621549A CN2011103124394A CN201110312439A CN102621549A CN 102621549 A CN102621549 A CN 102621549A CN 2011103124394 A CN2011103124394 A CN 2011103124394A CN 201110312439 A CN201110312439 A CN 201110312439A CN 102621549 A CN102621549 A CN 102621549A
- Authority
- CN
- China
- Prior art keywords
- phase
- value
- phase gradient
- interferometric
- cloth
- 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
Landscapes
- Image Analysis (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明提出一种多基线/多频段相位解缠的频域算法。技术方案的思路是:首先,利用干涉相位图的干涉相位观测值计算在水平方向及垂直方向上的相位梯度值,并根据时域无旋条件计算相位梯度图边界上的相位梯度值。然后,对相位梯度值进行傅立叶变换得到相位梯度值的频域函数,在相位梯度估计值的频域函数与各基线/各频段相位梯度值的频域函数加权和之差的平方和最小约束条件下计算出满足无旋条件的傅立叶系数的近似值。其后,对傅立叶系数近似值进行傅立叶反变换,得到满足相位无旋条件的相位梯度估计值。最后,对得到的满足无旋条件的相位梯度估计值进行积分,从而得到相位解缠值。本发明在不影响精度的情况下,减少了相位解缠算法的计算量。
Description
技术领域
本发明属于遥感和信号处理的交叉技术领域,特别涉及一种利用干涉合成孔径雷达进行多基线/多频段条件下的干涉相位解缠的频域快速方法。
背景技术
一切将相位由主值或相位差值恢复为真实值的过程统称为相位解缠,除了干涉合成孔径雷达应用外,相位解缠在合成孔径声纳、自适应光学、核磁共振、地震处理等方面都有重要应用。利用多部干涉合成孔径雷达形成多基线或多频段干涉可提高相位解缠的性能。
本发明以干涉合成孔径雷达应用为例。传统的单基线或单频段干涉合成孔径雷达系统受干涉相位模糊和高程叠掩影响,在复杂地形区域相位解缠难度较大,极大地限制了单基线或单频段干涉合成孔径雷达系统的高精度全球测绘能力。多基线或多频段干涉合成孔径雷达系统的提出与实现则有效地提高了干涉合成孔径雷达对复杂地形的测量精度和测量覆盖能力。多基线或多频段干涉合成孔径雷达系统的最大优点就是可以充分利用其长短基线或高低频段获取疏密不同的干涉相位条纹来提高相位解缠的性能,短基线或低频段可以保证相位解缠的可靠性,长基线或高频段可以提高测量精度。因此多基线/多频段干涉合成孔径雷达系统更具吸引力,是未来发展的趋势。
目前多基线/多频段相位解缠方法主要有:中国余数定律法、投影法以及线性组合法、迭代法、时域最小二乘法、Kalman滤波法、最大似然法、最大后验法、空-像域联合子空间正交投影法和网络流法等,其中多基线/多频段时域最小二乘法的基本思想是使相位梯度估计值与多个相位梯度值加权和之差的平方和最小,等效于求解具有牛曼边界的泊松方程,这种方法本质是对误差进行平均, 特点是十分稳健,但效率不高。
发明内容
本发明提出一种计算量小的多基线/多频段相位解缠的频域算法,该算法的基本思想是在频域范围内使相位梯度估计值与多个相位梯度值加权和之差的平方和最小。
本发明技术方案的思路是:首先,针对各基线/各频段的干涉相位图,利用干涉相位图的干涉相位观测值计算在水平方向及垂直方向上的相位梯度值,并根据时域无旋条件计算相位梯度图边界上的相位梯度值。然后,对两个方向上的相位梯度值进行傅立叶变换得到相位梯度值的频域函数,在相位梯度估计值的频域函数与各基线/各频段相位梯度值的频域函数加权和之差的平方和最小约束条件下计算出满足无旋条件的傅立叶系数的近似值。其后,对计算得到的满足无旋条件的傅立叶系数近似值,进行傅立叶反变换,得到满足相位无旋条件的相位梯度估计值。最后,对得到的满足无旋条件的相位梯度估计值进行沿着任意路径的积分,从而得到相位解缠值。
本发明技术方案是:
假定通过多基线/多频段观测得到K幅干涉相位图,其中第k幅干涉相位图对应第k条基线/第k个频段,第k条基线/第k个频段对应的缠绕相位函数为: k=1,2,糑,;m=0,1,2,糓,-1;n=0,1,2,糔,-1,M和N分别表示干涉相位图的方位向和距离向点数。记第k幅干涉相位图的垂直有效基线为bk(或记第k幅干涉相位图的波长为lk),记第k幅干涉相位图与第一幅干涉相位图的垂直有效基线之比为ak=b1/bk(或记第k幅干涉相位图与第一幅干涉相位图的波长之比为ak=lk/l1)。利用上述的观测信息,完成以下步骤:
第一步:干涉相位图的相位梯度值计算。
本步骤对多基线/多频段干涉相位进行相位梯度值计算,同时根据相位的时域无旋条件计算相位梯度图边界上的相位梯度值。对每一幅干涉相位图进行下述计算:
第(1)步,计算非边界上的相位梯度值。
上式中W{×}表示取相位主值运算。
第(2)步:计算边界上的相位梯度值。
第二步,傅立叶系数近似值计算。
其中,p=0,1,K,M-1;q=0,1,K,N-1。
第三步:相位梯度估计值计算。
第四步:相位解缠值计算。
用下式计算第一幅干涉相位图的相位解缠值:
利用上面计算的相位解缠值,可进行干涉合成孔径雷达系统的后续高程反演和地形形变估计等应用。
采用本发明可取得以下技术效果:
本发明的第二步计算两个方向相位梯度的傅立叶变换系数近似值,计算公式是在频域范围内使相位梯度估计值与多个相位梯度值加权和之差的平方和最小这一约束条件推导得到的,无需镜像对称操作,相对时域最小二乘法而言大大减少了计算量,同时保持了相位解缠值的计算精度。
附图说明
图1为本发明提供的多基线/多频段相位解缠频域快速算法流程示意图;
图2为基线长度为100米情况下的去平地后的干涉相位图;
图3为基线长度为200米情况下的去平地后的干涉相位图;
图4为基线长度为300米情况下的去平地后的干涉相位图;
图5为利用多基线时域最小二乘法的解缠结果;
图6为利用本发明得到的解缠结果;
图7为图5和图6利用的两种解缠算法的差值图;
图8为计算图5和图6的处理速度对比图。
具体实施方式
图1为本发明提供的多基线相位解缠频域快速算法流程示意图。整个流程分为四步。第一步,干涉相位图的相位梯度值计算;这一步包括计算非边界上的相位梯度值和计算边界上的相位梯度值这两个步骤,可以得到边界满足无旋条件的相位梯度场。第二步,傅立叶系数近似值计算;基于第一步中建立的相位梯度场,利用最小二乘方法进行傅立叶系数计算,得到满足无旋条件的傅立叶系数近似值。第三步,相位梯度估计值计算;这一步利用第二步中的傅立叶系数估计值,进行傅立叶反变换,获得无旋条件下的相位梯度估计值。第四步,相位解缠值计算;对于第三步所得的满足无旋条件的相位梯度估计值进行沿着任意路径的积分,从而得到相位解缠值。
图2~图4是利用实验室的天基雷达仿真系统仿真得到的干涉基线分别为100米、200米和300米条件下的干涉相位图,仿真过程中用到的地形是Etna火山口,图像大小为512×768像素。
图5~图8是利用图2~图4提供的3幅干涉相位图进行仿真实验的结果,即K=3。
图5为利用多基线时域最小二乘法的解缠结果;图6为利用本发明得到的解缠结果;图7为图5和图6利用的两种解缠算法的差值图。由图7可看出,两种方法的差异在10的负3次方量级,说明两种算法的精度很接近。
图8为计算图5和图6的处理速度对比图,即利用本发明的解缠算法和传统的多基线时域最小二乘法的解缠方法处理速度对比表。
从理论上对比:时域最小二乘方法等效于求解具有牛曼边界的泊松方程,为快速求解也采用傅立叶变换到频域,但为了满足牛曼边界条件,需要作镜像对称操作,其傅立叶变换的总运算量为:
本发明的方法直接在频域上满足最小二乘约束,避免了牛曼边界条件约束,无需镜像对称操作,其傅立叶变换总运算量为:
利用上述两个公式,比较两种方法的傅立叶变换运算量,可知由于本发明不用做镜像对称操作,FFT数据处理点数小,因此可以更为有效快速地得到解缠相位。
对于两种算法,实验环境均为Inter Core2 Quad CPU 2.33GHz,内存2GB。对于上述多基线干涉相位,图8中结果表明,本发明提出的算法处理时间小于多基线时域最小二乘解缠方法的1/3,速度更快。
Claims (2)
1.一种多基线或多频段相位解缠的频域算法,其特征在于,包括下述步骤:
假定通过多基线或多频段观测得到K幅干涉相位图,K32,其中第k幅干涉相位图对应第k条基线或第k个频段,第k条基线或第k个频段对应的缠绕相位函数为: k=1,2,糑,;m=0,1,2,糓,-1;n=0,1,2,糔,-1,M和N分别表示干涉相位图的方位向和距离向点数;
记第k幅干涉相位图的垂直有效基线为bk或记第k幅干涉相位图的波长为lk;
记第k幅干涉相位图与第一幅干涉相位图的垂直有效基线之比为ak=b1/bk或记第k幅干涉相位图与第一幅干涉相位图的波长之比为ak=lk/l1;
利用上述的观测信息和ak,完成以下步骤:
第一步:干涉相位图的相位梯度值计算;
第(1)步,计算非边界上的相位梯度值;
上式中W{×}表示取相位主值运算;
第(2)步:计算边界上的相位梯度值;
第二步,傅立叶系数近似值计算;
第三步:相位梯度估计值计算;
第四步:相位解缠值计算;
用下式计算第一幅干涉相位图的相位解缠值:
上式中 m ii=0,1,2,K M-2;n ii=0,1,2,K N-2, 表示第一幅干涉相位图在m=0,n=0时缠绕相位函数的值。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110312439.4A CN102621549B (zh) | 2011-10-14 | 2011-10-14 | 多基线/多频段干涉相位解缠频域快速算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201110312439.4A CN102621549B (zh) | 2011-10-14 | 2011-10-14 | 多基线/多频段干涉相位解缠频域快速算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN102621549A true CN102621549A (zh) | 2012-08-01 |
CN102621549B CN102621549B (zh) | 2014-04-16 |
Family
ID=46561576
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201110312439.4A Expired - Fee Related CN102621549B (zh) | 2011-10-14 | 2011-10-14 | 多基线/多频段干涉相位解缠频域快速算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN102621549B (zh) |
Cited By (15)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103325092A (zh) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | 一种生成二维相位解缠质量图的方法及装置 |
CN104062658A (zh) * | 2014-06-09 | 2014-09-24 | 桂林电子科技大学 | 一种基于频率补偿的多基线sar干涉相位估计方法 |
CN104730519A (zh) * | 2015-01-15 | 2015-06-24 | 电子科技大学 | 一种采用误差迭代补偿的高精度相位解缠方法 |
CN105093226A (zh) * | 2015-08-31 | 2015-11-25 | 西安电子科技大学 | 基于全局最小均方算法的雷达相位解缠方法 |
CN103869290B (zh) * | 2014-03-20 | 2016-03-02 | 西安电子科技大学 | 基于大规模数据加权枝切线的相位解缠绕方法 |
CN107092022A (zh) * | 2017-04-21 | 2017-08-25 | 哈尔滨工业大学 | 基于InSAL的区域滤波质量引导相位解缠方法 |
CN107390218A (zh) * | 2017-08-28 | 2017-11-24 | 西安电子科技大学 | 基于最小无穷范数的二维相位解缠绕方法 |
CN107544069A (zh) * | 2017-08-28 | 2018-01-05 | 西安电子科技大学 | 基于平面近似模型的多基线相位解缠绕方法 |
CN108363056A (zh) * | 2018-01-25 | 2018-08-03 | 北京无线电测量研究所 | 一种多频点干涉成像方法及系统 |
CN108519588A (zh) * | 2018-04-12 | 2018-09-11 | 视缘(上海)智能科技有限公司 | 一种多频相位解缠方法及装置 |
CN109472834A (zh) * | 2018-10-23 | 2019-03-15 | 桂林电子科技大学 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
CN109633648A (zh) * | 2019-01-22 | 2019-04-16 | 北京航空航天大学 | 一种基于似然估计的多基线相位估计装置及方法 |
CN110082709A (zh) * | 2019-05-14 | 2019-08-02 | 中国电子科技集团公司第三十六研究所 | 一种均匀五元圆阵测向方法 |
CN111681299A (zh) * | 2020-06-16 | 2020-09-18 | 中国科学院空天信息创新研究院 | 基于InSAR解缠相位生成数字表面模型的方法及装置 |
CN113238227A (zh) * | 2021-05-10 | 2021-08-10 | 电子科技大学 | 一种结合深度学习的改进最小二乘相位解缠方法及系统 |
Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887121A (zh) * | 2010-06-22 | 2010-11-17 | 北京航空航天大学 | 基于半牛顿迭代法的星载干涉合成孔径雷达基线估计方法 |
CN102073035A (zh) * | 2010-12-13 | 2011-05-25 | 中国人民解放军国防科学技术大学 | 星载干涉合成孔径雷达系统仿真性能评估方法 |
-
2011
- 2011-10-14 CN CN201110312439.4A patent/CN102621549B/zh not_active Expired - Fee Related
Patent Citations (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101887121A (zh) * | 2010-06-22 | 2010-11-17 | 北京航空航天大学 | 基于半牛顿迭代法的星载干涉合成孔径雷达基线估计方法 |
CN102073035A (zh) * | 2010-12-13 | 2011-05-25 | 中国人民解放军国防科学技术大学 | 星载干涉合成孔径雷达系统仿真性能评估方法 |
Non-Patent Citations (1)
Title |
---|
向茂生等: "一种基于梯度估值的LMS相位解缠方法", 《遥感学报》, 31 May 1998 (1998-05-31), pages 94 - 96 * |
Cited By (27)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN103325092A (zh) * | 2013-03-13 | 2013-09-25 | 中国科学院电子学研究所 | 一种生成二维相位解缠质量图的方法及装置 |
CN103325092B (zh) * | 2013-03-13 | 2016-06-08 | 中国科学院电子学研究所 | 一种生成二维相位解缠质量图的方法及装置 |
CN103869290B (zh) * | 2014-03-20 | 2016-03-02 | 西安电子科技大学 | 基于大规模数据加权枝切线的相位解缠绕方法 |
CN104062658A (zh) * | 2014-06-09 | 2014-09-24 | 桂林电子科技大学 | 一种基于频率补偿的多基线sar干涉相位估计方法 |
CN104062658B (zh) * | 2014-06-09 | 2016-05-11 | 桂林电子科技大学 | 一种基于频率补偿的多基线sar干涉相位估计方法 |
CN104730519B (zh) * | 2015-01-15 | 2017-04-05 | 电子科技大学 | 一种采用误差迭代补偿的高精度相位解缠方法 |
CN104730519A (zh) * | 2015-01-15 | 2015-06-24 | 电子科技大学 | 一种采用误差迭代补偿的高精度相位解缠方法 |
CN105093226A (zh) * | 2015-08-31 | 2015-11-25 | 西安电子科技大学 | 基于全局最小均方算法的雷达相位解缠方法 |
CN107092022A (zh) * | 2017-04-21 | 2017-08-25 | 哈尔滨工业大学 | 基于InSAL的区域滤波质量引导相位解缠方法 |
CN107092022B (zh) * | 2017-04-21 | 2019-08-02 | 哈尔滨工业大学 | 基于InSAL的区域滤波质量引导相位解缠方法 |
CN107390218A (zh) * | 2017-08-28 | 2017-11-24 | 西安电子科技大学 | 基于最小无穷范数的二维相位解缠绕方法 |
CN107544069A (zh) * | 2017-08-28 | 2018-01-05 | 西安电子科技大学 | 基于平面近似模型的多基线相位解缠绕方法 |
CN107390218B (zh) * | 2017-08-28 | 2020-06-05 | 西安电子科技大学 | 基于最小无穷范数的二维相位解缠绕方法 |
CN108363056B (zh) * | 2018-01-25 | 2020-08-18 | 北京无线电测量研究所 | 一种多频点干涉成像方法及系统 |
CN108363056A (zh) * | 2018-01-25 | 2018-08-03 | 北京无线电测量研究所 | 一种多频点干涉成像方法及系统 |
CN108519588B (zh) * | 2018-04-12 | 2023-04-21 | 视缘(上海)交通科技有限公司 | 一种多频相位解缠方法及装置 |
CN108519588A (zh) * | 2018-04-12 | 2018-09-11 | 视缘(上海)智能科技有限公司 | 一种多频相位解缠方法及装置 |
CN109472834A (zh) * | 2018-10-23 | 2019-03-15 | 桂林电子科技大学 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
CN109472834B (zh) * | 2018-10-23 | 2023-04-14 | 桂林电子科技大学 | 一种基于小波变换的卡尔曼滤波相位展开方法 |
CN109633648B (zh) * | 2019-01-22 | 2022-04-29 | 北京航空航天大学 | 一种基于似然估计的多基线相位估计装置及方法 |
CN109633648A (zh) * | 2019-01-22 | 2019-04-16 | 北京航空航天大学 | 一种基于似然估计的多基线相位估计装置及方法 |
CN110082709B (zh) * | 2019-05-14 | 2021-04-13 | 中国电子科技集团公司第三十六研究所 | 一种均匀五元圆阵测向方法 |
CN110082709A (zh) * | 2019-05-14 | 2019-08-02 | 中国电子科技集团公司第三十六研究所 | 一种均匀五元圆阵测向方法 |
CN111681299A (zh) * | 2020-06-16 | 2020-09-18 | 中国科学院空天信息创新研究院 | 基于InSAR解缠相位生成数字表面模型的方法及装置 |
CN111681299B (zh) * | 2020-06-16 | 2024-03-15 | 中国科学院空天信息创新研究院 | 基于InSAR解缠相位生成数字表面模型的方法及装置 |
CN113238227A (zh) * | 2021-05-10 | 2021-08-10 | 电子科技大学 | 一种结合深度学习的改进最小二乘相位解缠方法及系统 |
CN113238227B (zh) * | 2021-05-10 | 2022-09-30 | 电子科技大学 | 一种结合深度学习的改进最小二乘相位解缠方法及系统 |
Also Published As
Publication number | Publication date |
---|---|
CN102621549B (zh) | 2014-04-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN102621549A (zh) | 多基线/多频段干涉相位解缠频域快速算法 | |
Sun et al. | Lightning VHF radiation location system based on short-baseline TDOA technique—Validation in rocket-triggered lightning | |
CN102175988B (zh) | 一种基于维度拆分的相关干涉仪测向方法 | |
CN106093939B (zh) | 一种基于相位差统计模型的InSAR图像相位解缠方法 | |
CN102866001B (zh) | 基于横向剪切干涉结构的成像系统像差的测量方法 | |
CN102955157B (zh) | 一种用于干涉合成孔径雷达图像精配准的快速相干系数法 | |
CN107632964A (zh) | 一种平面地磁异常场向下延拓递归余弦变换法 | |
CN103630900B (zh) | 3-d sar波数域快速成像的方法 | |
CN107992676A (zh) | 一种运动目标时域散射回波的快速仿真建模方法 | |
CN104730519B (zh) | 一种采用误差迭代补偿的高精度相位解缠方法 | |
CN109633648A (zh) | 一种基于似然估计的多基线相位估计装置及方法 | |
CN103456015A (zh) | 基于最优分数域Gabor谱特征的SAR目标检测方法 | |
CN104251990A (zh) | 合成孔径雷达自聚焦方法 | |
CN102902841A (zh) | 一种计算机载复杂结构雷达横截面积的方法 | |
CN108387880B (zh) | 一种基于多尺度有向Lyapunov指数的海杂波背景下弱小目标检测方法 | |
CN109725290A (zh) | 一种误差提取方法、装置、电子设备及可读存储介质 | |
Sinha et al. | Numerical modeling of tide-surge interaction along Orissa coast of India | |
CN102183761B (zh) | 星载干涉合成孔径雷达数字高程模型重建方法 | |
CN105572629B (zh) | 一种适用于任意阵列结构的低运算复杂度的二维测向方法 | |
CN103630901B (zh) | 机载下视阵列3-d sar成像的方法 | |
CN104123462A (zh) | 实多项式求根实现均匀线阵的谱music方法 | |
CN103728617B (zh) | 双基地合成孔径雷达时域快速成像方法 | |
CN105469368A (zh) | 一种干涉相位滤波方法 | |
CN103869315B (zh) | 临近空间圆周合成孔径雷达快速后向投影成像方法 | |
Zhou et al. | SAR imaging realization with FPGA based on VIVADO HLS |
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 |
Granted publication date: 20140416 Termination date: 20181014 |
|
CF01 | Termination of patent right due to non-payment of annual fee |