CN102662171A - 一种sar层析三维成像方法 - Google Patents

一种sar层析三维成像方法 Download PDF

Info

Publication number
CN102662171A
CN102662171A CN2012101213400A CN201210121340A CN102662171A CN 102662171 A CN102662171 A CN 102662171A CN 2012101213400 A CN2012101213400 A CN 2012101213400A CN 201210121340 A CN201210121340 A CN 201210121340A CN 102662171 A CN102662171 A CN 102662171A
Authority
CN
China
Prior art keywords
sar
height
centerdot
signal
phi
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
Application number
CN2012101213400A
Other languages
English (en)
Other versions
CN102662171B (zh
Inventor
闵锐
曹宗杰
黄伟
窦宏林
皮亦鸣
范录宏
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
University of Electronic Science and Technology of China
Original Assignee
University of Electronic Science and Technology of China
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by University of Electronic Science and Technology of China filed Critical University of Electronic Science and Technology of China
Priority to CN 201210121340 priority Critical patent/CN102662171B/zh
Publication of CN102662171A publication Critical patent/CN102662171A/zh
Application granted granted Critical
Publication of CN102662171B publication Critical patent/CN102662171B/zh
Expired - Fee Related legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Landscapes

  • Radar Systems Or Details Thereof (AREA)

Abstract

本发明公开了一种SAR层析三维成像方法。本发明的方法采用将正则化正交匹配追踪(ROMP)与SAR层析三维成像系统结合,通过准确重构高度向稀疏信号来实现高度维聚焦。当测量矩阵严格满足约束等距性条件(RIP)时,该方法可对任意稀疏信号进行精确重构。利用该成像方法,对稀疏度为k的稀疏信号经过k次迭代即可获得一个原子数|I|≥2k的支撑集I来精确重构高度向稀疏信号,运算量小,运算效率高,而且每次迭代时对观测向量进行排序,保证了迭代的最优性,具有较高的重构稳健性。

Description

一种SAR层析三维成像方法
技术领域
本发明属于电子信号处理技术领域,涉及空间遥感和空对地观测信息处理技术领域,特别涉及机载合成孔径雷达(Synthetic Aperture Radar,SAR)层析三维成像技术。
背景技术
SAR是二十世纪雷达技术发展的重要里程碑,它利用雷达回波信号的相关性,累积雷达运动过程中回波信号的多普勒频移,在雷达的运动方向上合成等效的雷达孔径,实现方位向的高分辨成像。因SAR采用主动工作模式,对自身发射电磁波的反射回波进行成像处理,不受光照、温度等外界环境的限制,可实现全天时、全天候的区域检测成像,且对植被、测绘和军事等领域得到了广泛的应用。
传统SAR是对三维场景的二维成像且由于距离向成像位置由目标与雷达的距离决定,SAR成像结果实际上是具有相同距离(不同地面位置、不同高程的综合结果)的所有散射体散射系数的叠加,因此在二维SAR成像中存在着圆柱对称模糊(即成像结果无法分辨具有相同距离的不同散射体)、叠掩现象(雷达接收到目标上部反射回波先于下部反射回波)等问题,难以满足越来越高的成像精度和复杂环境侦查的要求。而问题的根源在于二维SAR成像的本质是三维场景空间到二维成像平面沿距离为半径的圆柱面的投影。因此,解决问题的最直接方法就是获得目标的第三维信息,避免三维空间到二维平面投影的信息损失。
层析SAR(TomoSAR)将合成孔径原理扩展至第三维即高度向,从而得到三维图像。TomoSAR的基本思想是通过不同时刻和不同轨道位置对同一地区采集的多次SAR回波数据,在高度向上合成大孔径。国内外机载SAR层析三维成像的主要形式为:首先利用机载二维SAR系统在不同视角的轨迹上获得同一区域的多幅SAR二维成像结果,再将所有的SAR图像进行配准处理以使对图像的处理对应于同一散射体的处理;然后按照轨迹分布的顺序逐像素的将所有图像构建高度向序列信号,并对高度向信号进行调频校正后通过傅里叶变换得到散射体的高度向分布,加之二维成像的斜距-方位向共同实现三维无模糊成像。
然而,由于航迹数量和分布的稀疏性,高度向难以得到大孔径,因此高度向分辨率远低于方位向分辨率,因此在不影响方位-距离分辨率的条件下,高度向需使用超高分辨率算法。鉴于高度向信号的稀疏性,即在方位-距离分辨单元内的散射可以看作少量散射中心散射的合成,提出了基于压缩感知(Compress Sensing)理论的高度向超高分辨率算法。
对于稀疏基线的SAR成像过程,可以看作是稀疏信号的重构,目前,主要的稀疏信号表示方法有基追踪(Basis Pursuit,BP)算法、贪婪算法。BP算法有很好的稳建性,但是运算效率很差;另一方面,贪婪算法虽然运算效率高,但是稳健性得不到保证。
发明内容
本发明的目的是为了解决现有的SAR层析三维成像方法存在的上述问题,提出了一种SAR层析三维成像方法。
本发明的技术方案为:一种SAR层析三维成像方法,具体包括如下步骤:
步骤1,对SAR层析系统各基线接收的回波信号进行二维成像处理,将所有成像结果按基线顺序进行配准处理,逐像素构建高度向信号,进行相位校正后得到观测向量G:
G=X+W=Φ.Γ+W
其中,G=[x11,…,xmm,…,xMM]T,X=[x1,…,xm,…,xM]T为实测数据向量,W=[ω1,·,ωm,·,ωM]T为独立同分布的高斯白噪声;Φ=[Φ0,…,Φn,…,ΦN-1]为测量矩阵,其中, Φ n = [ φ n 1 , · · · , φ n m , · · · , φ n M ] T , φ n m = exp { - j 2 π λ [ ( nΔs ) 2 - 2 nΔs · s m + s m 2 r 0 ] } ,
Γ=[γ0,·,γn,·,γN-1]T为散射系数矩阵,其中,γn=γ(x0,r0,nΔs),M表示基线数目,N表示高度向的采样点数,λ表示雷达发射信号的波长,x0为方位向坐标,r0为距离向坐标,Δs表示高度向的采样间隔,sm代表第m条航迹的高度。
步骤2,构造索引值集合
Figure BDA0000156351050000023
Figure BDA0000156351050000024
初始化高度向信号稀疏度k;残差r=G,迭代次数p=1;
步骤3,计算相关系数向量y=Φ*r,找出y中最大的k个元素,将其下标构成集合J;若y中的非零元素个数小于k,那么就将所有这些非零元素的下标构成集合J;
步骤4,对步骤3得到的集合J中所有元素进行正则化处理,即从集合J中选择满足下面条件的子集
Figure BDA0000156351050000025
J0中可比较的坐标对所有i,j∈J0满足
2|y(i)|≥|y(j)|
步骤5,从步骤4得到的J0组合中,选择能量最大的,即
Figure BDA0000156351050000026
最大,将其添加进索引值集合I中,即I←I∪J0,得到匹配原子系数集合;
步骤6,利用步骤5得到的匹配原子系数集合更新残差r,具体过程为:计算重构向量其中ΦI表示通过步骤5中的集合I对矩阵各列的索引来限制测量矩阵Φ,即测量矩阵Φ由步骤1中的向量{Φn|n∈I}构成,然后更新残差r=G-Φx;
步骤7,如果p>k或者集合I中的元素个数|I|≥2k,则停止迭代,x即为最终所求的高度向信号的重构向量,否则,令p=p+1,转到步骤3继续迭代。
本发明的有益效果:本发明针对现有的BP算法稳健性好,但运算效率差,而贪婪算法运算效率高,但稳健性差的缺陷,设计了一种SAR层析三维成像方法,克服了现有的BP算法稳健性好,但运算效率差,而贪婪算法运算效率高,但稳健性差的缺陷。本发明采用将正则化正交匹配追踪(ROMP)与SAR层析三维成像系统结合,通过准确重构高度向稀疏信号来实现高度维聚焦。当测量矩阵严格满足约束等距性条件(RIP)时,该方法可对任意稀疏信号进行精确重构。利用该成像方法,对稀疏度为k的稀疏信号经过k次迭代即可获得一个原子数|I|≥2k的支撑集I来精确重构高度向稀疏信号,运算量小,运算效率高,而且每次迭代时对观测向量进行排序,保证了迭代的最优性,具有较高的重构稳健性。因此,这种基于正则化正交匹配追踪的SAR层析三维成像方法在运算效率和重构精度两方面都具有较好的性能,极大地拓展了SAR层析三维成像平台的应用前景。
附图说明
图1为本发明的SAR层析三维成像方法的流程图。
图2为本发明具体实施方案采用的层析SAR成像的几何模型示意图。
具体实施方式
下面结合附图和具体实施方式对本发明做进一步的说明。
本发明的SAR层析三维成像方法,如图1所示,具体包括如下步骤:
步骤1,SAR平台在不同时刻和不同航迹位置对同一被测区域采集得到m次回波数据。然后运用现有的保留相位信息的经典二维SAR成像方法得到每次飞行的二维复图像。二维复图像中方位-距离像素点(x0,r0)的第m次飞行得到的灰度值如下式:
g ( x 0 , r 0 , s m ) = ∫ γ ( x 0 , r 0 , s ) · exp [ - j 2 π λ · 2 R m ( r 0 , s ) ] ds
其中,γ(·)为高度向的三维反射系数,x0为方位向坐标,r0为距离向坐标,s为高度向坐标,λ为雷达发射信号的波长,Rm(·)为第m次航迹与散射点的最短距离,sm代表第m条航迹的高度。
第m次航迹与散射点的最短距离可以近似表达为:Rm(r0,s)≈r0+(sm-s)/2r0。令Δs表示高度向采样间隔,经过相位矫正后将上式近似表示为离散形式:
g m = Σ n = 0 N - 1 γ m ( x 0 , r 0 , nΔs ) · exp { - j 2 π λ [ ( s m - nΔs ) 2 r 0 ] }
其中,N为采样点数。将步骤1中得到的信号离散形式表示为矩阵形式:
G=X+W=Φ.Γ+W
其中,G=[x11,…,xmm,…,xMM]T为观测向量,X=[x1,…,xm,…,xM]T为实测数据向量,W=[ω1,…,ωm,…,ωM]T为独立同分布的高斯白噪声;Φ=[Φ0,…,Φn,…,Φn-1]为测量矩阵,其中 Φ n = [ φ n 1 , · · · , φ n m , · · · , φ n M ] T , φ n m = exp { - j 2 π λ [ ( nΔs ) 2 - 2 nΔs · s m + s m 2 r 0 ] } ; Γ=[γ0,…,γn,…,γN-1]T为散射系数矩阵,其中,γn=γ(x0,r0,nΔs)。M表示基线数目,N表示高度向的采样点数。
步骤1的实质:对SAR层析系统各基线接收的回波信号进行二维成像处理,将所有成像结果按基线顺序进行配准处理,逐像素构建高度向信号,进行相位校正后得到观测向量G;
步骤2,构造索引值集合
Figure BDA0000156351050000044
Figure BDA0000156351050000045
初始化信号稀疏度为k;残差r=G,迭代次数p=1。
这里,初始化信号稀疏度具体为:令k=T,T为场景目标个数,
Figure BDA0000156351050000046
表示空集合。
步骤3,计算相关系数向量y=Φ*r,找出y中最大的k个元素,将其下标构成集合J;若y中的非零元素个数小于k,那么就将所有这些非零元素的下标构成集合J;
步骤4,对步骤3得到的集合J中所有元素进行正则化处理,即从集合J中选择满足下面条件的子集
Figure BDA0000156351050000047
J0中可比较的坐标对所有i,j∈J0满足
2|y(i)|≥|y(j)|
步骤5,从步骤4得到的J0组合中,选择能量最大的,
Figure BDA0000156351050000048
最大,将其添加进索引值集合I中,即I←I∪J0,得到更新后的索引值集合I;
步骤6,利用步骤5得到的索引值集合I来更新残差r,具体过程为:计算重构向量
Figure BDA0000156351050000049
其中ΦI表示通过步骤5中的集合I对矩阵各列的索引来限制测量矩阵Φ,即测量矩阵Φ由步骤1中的向量{Φn|n∈I}构成,然后更新残差r=G-Φx;
步骤7,如果p>k或者集合I中的元素个数|I|≥2k,则停止迭代,x即为最终所求的高度向信号的重构向量,否则,令p=p+1,转到步骤3继续迭代。
需要说明的是:
1、步骤3中严格满足RIP条件保证了对稀疏度为k的信号,测量矩阵Φ的任意k列向量近似构成正交系统,每次迭代选择观测向量中的k个坐标值而不是选择一个坐标,降低了运算量,提高了运算效率。
2、步骤4中为避免因选择错误坐标对迭代结果产生不利影响,引入正则化处理以确保所选坐标均含有所需的信号信息。因此,即使选择了错误的坐标,只要其数量较少,对迭代结果也不会产生太大影响。
本实施方案采用的层析SAR成像几何模型如图2所示,x,r,z分别表示方位向、斜距向以及高度向,v为垂直于x-r平面的法线方向。这里对点目标进行正侧视情况的飞行模式成像,假设距离-方位向的聚焦已经完成,且是理想的,只对高度向点目标回波进行仿真实验,以证明本发明的方法应用于层析SAR高度维聚焦的有效性。
本发明的成像方法当测量矩阵满足等距约束性条件RIP(Restricted Isometry Property)条件时,可以精确重建任意稀疏信号,而且重建速度较快。重建稀疏度为k的信号,该算法首先根据相关原则进行原子的一次筛选,通过求残差r与测量矩阵Φ中各个原子之间内积的绝对值,来计算相关系数y,并将按照此方法筛选出的k个原子的索引值存到候选集J中以便进行原子的二次筛选,然后采用正则化过程进行原子的二次筛选,即
Figure BDA0000156351050000051
按照如下条件将J中索引值对应的原子的相关系数分成若干组:2|y(i)|≥|y(j)|,然后选择能量最大的一组相关系数对应的原子索引值存入J0中,通过正则化处理可以使得本发明的方法最多经过k次迭代便可以得到一个原子数|I|≥2k的支撑集I用于精确重建信号,对于没有选入支撑集的原子,正则化过程能保证它们的能量一定远小于被选入原子的能量,经过k次迭代得到用于信号重建的支撑集后,再采用最小二乘法进行信号逼近以及余量更新,最终得到用于信号重构的观测向量以重构高度向信号,实现高度向高分辨成像。因此,本发明的方法在运算效率和重构精度两方面的良好性能极大地拓展了SAR层析三维成像平台的应用前景。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。

Claims (2)

1.一种SAR层析三维成像方法,其特征在于,具体包括如下步骤:
步骤1,对SAR层析系统各基线接收的回波信号进行二维成像处理,将所有成像结果按基线顺序进行配准处理,逐像素构建高度向信号,进行相位校正后得到观测向量G:
G=X+W=Φ·Γ+W
其中,G=[x11,…,xmm,…,xMM]T,X=[x1,…,xm,…,xM]T为实测数据向量,W=[ω1,…,ωm,…,ωM]T为独立同分布的高斯白噪声;Φ=[Φ0,…,Φn,…,ΦN-1]为测量矩阵,其中 Φ n = [ φ n 1 , · · · , φ n m , · · · , φ n M ] T , φ n m = exp { - j 2 π λ [ ( nΔs ) 2 - 2 nΔs · s m + s m 2 r 0 ] } ,
Γ=[γ0,…,γn,…,γN-1]T为散射系数矩阵,其中,γn=γ(x0,r0,nΔs),M表示基线数目,N表示高度向的采样点数,λ表示雷达发射信号的波长,x0为方位向坐标,r0为距离向坐标,Δs表示高度向的采样间隔,sm代表第m条航迹的高度。
步骤2,构造索引值集合
Figure FDA0000156351040000013
Figure FDA0000156351040000014
初始化高度向信号稀疏度k;残差r=G,迭代次数p=1;
步骤3,计算相关系数向量y=Φ*r,找出y中最大的k个元素,将其下标构成集合J;若y中的非零元素个数小于k,那么就将所有这些非零元素的下标构成集合J;
步骤4,对步骤3得到的集合J中所有元素进行正则化处理,即从集合J中选择满足下面条件的子集J0中可比较的坐标对所有i,j∈J0满足
2|y(i)|≥|y(j)|
步骤5,从步骤4得到的J0组合中,选择能量最大的,即
Figure FDA0000156351040000016
最大,将其添加进索引值集合I中,即I←I∪J0,得到更新后的索引值集合I;
步骤6,利用步骤5得到的索引值集合I更新残差r,具体过程为:计算重构向量
Figure FDA0000156351040000017
其中ΦI表示通过步骤5中的集合I对矩阵各列的索引来限制测量矩阵Φ,即测量矩阵Φ由步骤1中的向量{Φn|n∈I}构成,然后更新残差r=G-Φx;
步骤7,如果p>k或者集合I中的元素个数|I|≥2k,则停止迭代,x即为最终所求的高度向信号的重构向量,否则,令p=p+1,转到步骤3继续迭代。
2.根据权利要求1所述的SAR层析三维成像方法,其特征在于,初始化信号稀疏度具体为:k=T,T为场景目标个数。
CN 201210121340 2012-04-23 2012-04-23 一种sar层析三维成像方法 Expired - Fee Related CN102662171B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN 201210121340 CN102662171B (zh) 2012-04-23 2012-04-23 一种sar层析三维成像方法

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN 201210121340 CN102662171B (zh) 2012-04-23 2012-04-23 一种sar层析三维成像方法

Publications (2)

Publication Number Publication Date
CN102662171A true CN102662171A (zh) 2012-09-12
CN102662171B CN102662171B (zh) 2013-11-06

Family

ID=46771699

Family Applications (1)

Application Number Title Priority Date Filing Date
CN 201210121340 Expired - Fee Related CN102662171B (zh) 2012-04-23 2012-04-23 一种sar层析三维成像方法

Country Status (1)

Country Link
CN (1) CN102662171B (zh)

Cited By (21)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076608A (zh) * 2013-01-27 2013-05-01 西安电子科技大学 轮廓增强的聚束式合成孔径雷达成像方法
CN103344949A (zh) * 2013-06-18 2013-10-09 中国人民解放军海军航空工程学院 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法
CN103698763A (zh) * 2013-12-12 2014-04-02 电子科技大学 基于硬阈值omp的线阵sar稀疏成像方法
CN103809205A (zh) * 2014-02-27 2014-05-21 中国矿业大学(北京) 一种地质雷达波速层析探测的走时数据快速采集方法
CN105388476A (zh) * 2015-12-28 2016-03-09 河南工业大学 一种基于联合稀疏模型的层析sar成像方法
CN105528792A (zh) * 2016-01-11 2016-04-27 南京邮电大学 一种医学图像配准混合算法
CN106355151A (zh) * 2016-08-30 2017-01-25 电子科技大学 一种基于深度置信网络的三维sar图像目标识别方法
CN106872977A (zh) * 2016-12-28 2017-06-20 北京建筑大学 一种基于分段弱正交匹配追踪的层析sar三维成像方法
CN107064933A (zh) * 2017-03-10 2017-08-18 中国科学院遥感与数字地球研究所 基于循环谱估计的sar层析建筑物高度的方法
CN107817492A (zh) * 2017-09-25 2018-03-20 中国科学院电子学研究所 宽角合成孔径雷达的成像方法及装置
CN107977680A (zh) * 2017-12-15 2018-05-01 武汉艾孚达信息技术有限公司 一种信号稀疏重构的方法
CN108627835A (zh) * 2018-06-29 2018-10-09 中国科学院电子学研究所 全极化差分sar层析的目标重构方法
CN108983229A (zh) * 2018-05-03 2018-12-11 电子科技大学 基于sar层析技术的高压输电铁塔高度及形变提取方法
CN109917389A (zh) * 2019-04-16 2019-06-21 中国人民解放军国防科技大学 一种机载全息sar成像中的相位校正方法
CN110082764A (zh) * 2019-04-26 2019-08-02 西安电子科技大学 基于稳健正则化层析方法的sar图像成像方法
CN110275164A (zh) * 2018-03-13 2019-09-24 中国人民解放军国防科技大学 一种用于多发多收合成孔径雷达的三维成像方法
CN110717949A (zh) * 2018-07-11 2020-01-21 天津工业大学 基于tromp的干涉高光谱图像稀疏重建
CN111025246A (zh) * 2019-11-28 2020-04-17 北京遥测技术研究所 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法
CN111679277A (zh) * 2020-05-28 2020-09-18 电子科技大学 一种基于sbrim算法的多基线层析sar三维成像方法
CN111948654A (zh) * 2020-08-12 2020-11-17 中国科学院空天信息创新研究院 机载层析sar三维点云生成方法
CN116400356A (zh) * 2023-06-07 2023-07-07 中山大学 一种基于同质区域联合的层析sar三维成像方法

Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221697A (zh) * 2011-03-25 2011-10-19 电子科技大学 一种机载多天线sar层析三维成像系统及其成像方法

Patent Citations (1)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN102221697A (zh) * 2011-03-25 2011-10-19 电子科技大学 一种机载多天线sar层析三维成像系统及其成像方法

Non-Patent Citations (6)

* Cited by examiner, † Cited by third party
Title
任笑真 等: "一种SAR层析成像的RELAX改进算法", 《数据采集与处理》, vol. 25, no. 3, 31 May 2010 (2010-05-31), pages 302 - 306 *
张红 等: "SAR层析技术的研究与应用", 《遥感技术与应用》, vol. 25, no. 2, 30 April 2010 (2010-04-30), pages 282 - 287 *
张郭 等: "基于层析成像模型的聚束式SAR图像复原", 《计算机科学》, vol. 36, no. 4, 30 April 2009 (2009-04-30), pages 19 - 22 *
王彦平 等: "长序列星载合成孔径雷达数据层析处理技术", 《测试技术学报》, vol. 22, no. 6, 31 December 2008 (2008-12-31), pages 471 - 477 *
王金峰 等: "一种机载SAR层析三维成像算法", 《电子与信息学报》, vol. 32, no. 5, 31 May 2010 (2010-05-31), pages 1029 - 1033 *
王金峰 等: "基于CLEAN算法的稀疏基线SAR三维层析成像", 《信号处理》, vol. 25, no. 8, 31 August 2009 (2009-08-31), pages 491 - 494 *

Cited By (34)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN103076608A (zh) * 2013-01-27 2013-05-01 西安电子科技大学 轮廓增强的聚束式合成孔径雷达成像方法
CN103076608B (zh) * 2013-01-27 2014-07-23 西安电子科技大学 轮廓增强的聚束式合成孔径雷达成像方法
CN103344949A (zh) * 2013-06-18 2013-10-09 中国人民解放军海军航空工程学院 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法
CN103344949B (zh) * 2013-06-18 2015-03-18 中国人民解放军海军航空工程学院 基于Radon-线性正则模糊函数的雷达微弱动目标检测方法
CN103698763A (zh) * 2013-12-12 2014-04-02 电子科技大学 基于硬阈值omp的线阵sar稀疏成像方法
CN103698763B (zh) * 2013-12-12 2016-01-13 电子科技大学 基于硬阈值正交匹配追踪的线阵sar稀疏成像方法
CN103809205A (zh) * 2014-02-27 2014-05-21 中国矿业大学(北京) 一种地质雷达波速层析探测的走时数据快速采集方法
CN105388476B (zh) * 2015-12-28 2017-12-29 河南工业大学 一种基于联合稀疏模型的层析sar成像方法
CN105388476A (zh) * 2015-12-28 2016-03-09 河南工业大学 一种基于联合稀疏模型的层析sar成像方法
CN105528792A (zh) * 2016-01-11 2016-04-27 南京邮电大学 一种医学图像配准混合算法
CN106355151A (zh) * 2016-08-30 2017-01-25 电子科技大学 一种基于深度置信网络的三维sar图像目标识别方法
CN106355151B (zh) * 2016-08-30 2019-10-01 电子科技大学 一种基于深度置信网络的三维sar图像目标识别方法
CN106872977A (zh) * 2016-12-28 2017-06-20 北京建筑大学 一种基于分段弱正交匹配追踪的层析sar三维成像方法
CN107064933A (zh) * 2017-03-10 2017-08-18 中国科学院遥感与数字地球研究所 基于循环谱估计的sar层析建筑物高度的方法
CN107064933B (zh) * 2017-03-10 2020-12-11 中国科学院遥感与数字地球研究所 基于循环谱估计的sar层析建筑物高度的方法
CN107817492A (zh) * 2017-09-25 2018-03-20 中国科学院电子学研究所 宽角合成孔径雷达的成像方法及装置
CN107977680A (zh) * 2017-12-15 2018-05-01 武汉艾孚达信息技术有限公司 一种信号稀疏重构的方法
CN110275164A (zh) * 2018-03-13 2019-09-24 中国人民解放军国防科技大学 一种用于多发多收合成孔径雷达的三维成像方法
CN110275164B (zh) * 2018-03-13 2022-03-18 中国人民解放军国防科技大学 一种用于多发多收合成孔径雷达的三维成像方法
CN108983229A (zh) * 2018-05-03 2018-12-11 电子科技大学 基于sar层析技术的高压输电铁塔高度及形变提取方法
CN108983229B (zh) * 2018-05-03 2022-04-19 电子科技大学 基于sar层析技术的高压输电铁塔高度及形变提取方法
CN108627835A (zh) * 2018-06-29 2018-10-09 中国科学院电子学研究所 全极化差分sar层析的目标重构方法
CN108627835B (zh) * 2018-06-29 2021-07-27 中国科学院电子学研究所 全极化差分sar层析的目标重构方法
CN110717949A (zh) * 2018-07-11 2020-01-21 天津工业大学 基于tromp的干涉高光谱图像稀疏重建
CN109917389A (zh) * 2019-04-16 2019-06-21 中国人民解放军国防科技大学 一种机载全息sar成像中的相位校正方法
CN110082764A (zh) * 2019-04-26 2019-08-02 西安电子科技大学 基于稳健正则化层析方法的sar图像成像方法
CN111025246A (zh) * 2019-11-28 2020-04-17 北京遥测技术研究所 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法
CN111025246B (zh) * 2019-11-28 2021-09-07 北京遥测技术研究所 一种静止轨道sar对海面及船只复合场景成像仿真系统及方法
CN111679277A (zh) * 2020-05-28 2020-09-18 电子科技大学 一种基于sbrim算法的多基线层析sar三维成像方法
CN111679277B (zh) * 2020-05-28 2022-05-03 电子科技大学 一种基于sbrim算法的多基线层析sar三维成像方法
CN111948654A (zh) * 2020-08-12 2020-11-17 中国科学院空天信息创新研究院 机载层析sar三维点云生成方法
CN111948654B (zh) * 2020-08-12 2022-03-18 中国科学院空天信息创新研究院 机载层析sar三维点云生成方法
CN116400356A (zh) * 2023-06-07 2023-07-07 中山大学 一种基于同质区域联合的层析sar三维成像方法
CN116400356B (zh) * 2023-06-07 2023-08-18 中山大学 一种基于同质区域联合的层析sar三维成像方法

Also Published As

Publication number Publication date
CN102662171B (zh) 2013-11-06

Similar Documents

Publication Publication Date Title
CN102662171A (zh) 一种sar层析三维成像方法
CN102645651B (zh) 一种sar层析超分辨成像方法
US9417323B2 (en) SAR point cloud generation system
JP6106285B2 (ja) 衛星レーダー干渉度のイオン歪み補正方法及びその装置
CN107238824B (zh) 基于先验dem数据的星载sar图像几何精校正方法
JP6320577B2 (ja) 合成開口レーダ信号処理装置
CN104007440A (zh) 一种加速分解后向投影聚束合成孔径雷达成像方法
CN100545676C (zh) 基于相关加权的干涉合成孔径雷达干涉相位估计方法
CN103018740B (zh) 一种基于曲面投影的InSAR成像方法
KR102151362B1 (ko) 극좌표변환을 이용한 항공기기반 영상복원장치 및 이를 이용한 영상복원방법
CN112882030B (zh) InSAR成像干涉一体化处理方法
CN103630905B (zh) 阵列天线sar极坐标交叠子孔径成像方法
CN102788979A (zh) 一种基于后向投影InSAR成像配准的GPU实现方法
Magnard et al. Analysis of a maximum likelihood phase estimation method for airborne multibaseline SAR interferometry
CN102004250A (zh) 基于频域展开的星机联合双基地合成孔径雷达成像方法
CN103616682A (zh) 一种基于曲面投影的多基线InSAR处理方法
RU2661941C1 (ru) Способ построения радиолокационного изображения с помощью радиолокационной станции с синтезированной апертурой
KR100441590B1 (ko) 간섭측정용 합성 개구 레이다의 기하학적 특성을 이용하여지형고도를 측정하기 위한 디지털 고도모형 생성방법
CN115685200A (zh) 一种高精度大前斜视sar成像运动补偿与几何校正方法
CN103630898A (zh) 对多基线干涉sar相位偏置进行估计的方法
CN106872977A (zh) 一种基于分段弱正交匹配追踪的层析sar三维成像方法
CN103728617B (zh) 双基地合成孔径雷达时域快速成像方法
CN103792534A (zh) 一种基于先验相位结构知识的sar两维自聚焦方法
CN102890277B (zh) 移不变双基地合成孔径雷达距离徙动成像方法
Li et al. An advanced DSS-SAR InSAR terrain height estimation approach based on baseline decoupling

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: 20131106

Termination date: 20160423

CF01 Termination of patent right due to non-payment of annual fee