CN106872977A - 一种基于分段弱正交匹配追踪的层析sar三维成像方法 - Google Patents
一种基于分段弱正交匹配追踪的层析sar三维成像方法 Download PDFInfo
- Publication number
- CN106872977A CN106872977A CN201611239581.XA CN201611239581A CN106872977A CN 106872977 A CN106872977 A CN 106872977A CN 201611239581 A CN201611239581 A CN 201611239581A CN 106872977 A CN106872977 A CN 106872977A
- Authority
- CN
- China
- Prior art keywords
- sar
- oblique distance
- vertical
- target
- imaging methods
- 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.)
- Pending
Links
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
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)
- Image Analysis (AREA)
Abstract
本发明提供一种基于分段弱正交匹配追踪的层析SAR三维成像方法,包括以下步骤:S1、多次航过采集目标的原始SAR数据,并对原始数据进行二维高分辨成像、配准以及相位补偿的预处理,从而得到可以进行层析SAR成像的观测数据集;S2、利用预处理的多景二维SAR数据在斜距垂向进行孔径合成,并重构出斜距垂向的信息;S3、获取目标高程信息,显示目标建筑物的三维点云结果。该方法能够在未知斜距垂向信号稀疏度的情况下很好的重构斜距垂向信息,实现目标建筑物的三维成像。
Description
技术领域
本发明涉及合成孔径雷达(Synthetic Aperture Radar,SAR)成像技术领域,更具体地,涉及一种基于分段弱正交匹配追踪的层析SAR三维成像方法。
背景技术
SAR是一种主动式的对地观测技术,与传统的光学传感器相比,SAR可以实现全天时、全天候对地实时观测,并具有一定的地表穿透能力。SAR技术的发展带来了数据量和计算量的急剧膨胀,而且随着实时成像需求的增强,对大规模SAR回波数据的实时成像处理越发重要。
压缩感知理论在层析SAR成像领域的应用正在快速的发展。各种基于压缩感知理论算法的提出拓展了层析SAR多维成像的方式。2012年12月26日出版的《电子测量与仪器学报》公开了“基于正则化正交匹配追踪的SAR层析成像”技术,其提出基于正则化正交匹配追踪(Regularized Orthogonal Matching Pursuit,ROMP)的层析SAR三维成像。文献中利用仿真数据验证了正交匹配追踪在层析SAR三维成像的准确性,基于ROMP的层析SAR三维成像算法能精确重构斜距垂向(文献中称为高度向)信号,实现高分辨率成像。虽然该算法保证了重构的稳健性以及运算的效率,但ROMP是非自适应的重构算法,在实际应用中利用ROMP进行层析SAR三维成像时需要已知信号的稀疏度K才能够准确的重构出目标的斜距垂向信息,真实地还原目标地物的真实三维场景。而在实际应用中很难准确确定斜距垂向信号的稀疏度,这大大增加了生产应用的难度。
发明内容
本发明提供一种克服上述问题或者至少部分地解决上述问题的基于分段弱正交匹配追踪的层析SAR三维成像方法,在未知稀疏度的情况下很好的进行自适应,在层析SAR三维成像中很好的重构斜距垂向,实现目标的三维成像。
根据本发明的一个方面,提供一种层析SAR三维成像方法,包括以下步骤:
S1、获取目标建筑物的二维SAR影像并进行预处理,构成观测数据集;
S2、基于分段弱正交匹配追踪的压缩感知算法重构斜距垂向信息;
S3、获取目标高程信息,显示目标建筑物的三维点云结果。
作为优选的,所述步骤S1具体包括:采集多次航过目标的原始SAR数据,并对原始数据进行二维高分辨成像、配准以及相位补偿的预处理,从而得到可以进行层析SAR成像的观测数据集
作为优选的,所述二维SAR影像数据包括轨道高度、重访周期、入射角范围、中心斜距、方位分辨率、距离分辨率,单通道场景大小、最大获取长度。
作为优选的,所述步骤S2具体包括:
S21、将多个场景的二维SAR影像数据沿斜距垂向进行孔径合成,将其中一幅做为主影像,其他影像以此影像进行配准,配准后主影像上每个像素沿斜距垂向的同一散射单元构成一个观测信号;
S22、对每个观测信号利用分段弱正交匹配追踪的方法进行重构,得到目标建筑物沿斜距垂向的信息。
作为优选的,所述步骤S21进一步包括:分别确定被截取目标影像的方位向和距离向的起点和终点,,将一个散射单元的信号构成一个时间序列的信息。
作为优选的,所述步骤S22具体包括:
S221、获取第m航迹影像表达式Sm;
S222、将Sm表达式进行离散转换,得到用于压缩感知算法的观测向量、测量矩阵和散射稀疏矩阵;
S223、进行斜距垂向信息重构。
作为优选的,所述步骤S221具体包括,主影像中的任意一个像素的复数值Sm,作为真实三维反射情况在x-r平面的投影,即反射信号沿斜距垂向的积分:
式中,γ(s)为散射体反射率沿斜距垂向的分布函数,λ为波长,b⊥m为航迹m、m-1的基线距离,ro为主影像中心斜距。
作为优选的,所述步骤S222具体包括:
对第m航迹影像表达式Sm进行离散表示:
式中,Sn为第n个目标,N为沿斜距垂向的散射体数目,m=1...M,M<<N;将上式通过矩阵表示:
y=Φγ+ξ
式中,y=Sm,为M维观测量,表示M条航迹的观测量;
ξ=[ξ1,ξ2,…,ξM]T,为独立分布的高斯白噪声;
Φ=[Φ1,Φ2,…Φn,…ΦN],为M×N的测量矩阵;
γ=[γ1,γ2…,γn,…,γN]T为稀疏矩阵;
对观测向量y、测量矩阵Φ、稀疏矩阵γ进行分段弱正交匹配追踪算法处理,得到重构后的斜距信息
与现有技术相比,本发明的有益效果在于:本发明通过自适应的分段弱正交匹配追踪(Stagewise Weak Orthogonal Matching Pursuit,SWOMP)层析SAR三维成像算法,能够在未知斜距垂向信号稀疏度的情况下很好的重构斜距垂向信息,实现目标建筑物的三维成像。这极大地扩大层析SAR三维成像的成像范围。
附图说明
图1为本发明实施例1的基于分段弱正交匹配追踪的层析SAR三维成像方法流程框图;
图2为本发明实施例图1的具体流程示意图;
图3是本发明实施例1的分段弱正交匹配追踪算法流程图;
图4是本发明实施例2的SAR成像结果示意图。
具体实施方式
下面结合附图和实施例,对本发明的具体实施方式作进一步详细描述。以下实施例用于说明本发明,但不用来限制本发明的范围。
最后,本申请的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
实施例1
图1示出了一种基于分段弱正交匹配追踪的层析SAR三维成像方法的流程框图,包括以下步骤:
S1、获取目标建筑物的二维SAR影像并进行预处理,构成观测数据集;
S2、基于分段弱正交匹配追踪的压缩感知算法重构斜距垂向信息;
S3、获取目标高程信息,显示目标建筑物的三维点云结果。
在本实施例步骤S1中,如图2所示,[S1,S2,…Sm-1,Sm,…SM]为传感器在不同时间、不同航迹位置对同一目标区域获取的二维SAR影像。Sm-1Sm之间的基线距离为b⊥m,Sm到s轴的垂直距离为δ(δ是由于航线不稳定造成的,计算时可通过校正消除)。斜距垂向s(方位-距离向所在平面的法线)平行于s’。L为斜距垂向合成孔径长度。P1,Pn为目标地物在斜距垂向上的两个待测点,r0为主图像中心斜距。作为优选的,所述二维SAR影像数据包括轨道高度、重访周期、入射角范围、方位分辨率、距离分辨率,单通道场景大小、最大获取长度,还包括光速c、波长λ、入射角向量InA、近距(Near range)向量、中心距(Center Range)向量、远距(Far Range)向量、方位分辨率、距离分辨率、斜距R0、水平距离H0=R0*cos(InA)、基线向量。
作为优选的,所述步骤S1具体包括:通过传感器在不同时间、不同航迹位置对同一目标获取二维SAR影像数据。
如图2所示,所述步骤S2具体包括:
S21、将多个场景的二维SAR影像数据沿斜距垂向进行孔径合成,将其中一幅做为主影像,其他影像以此影像进行配准,配准后主影像上每个像素沿斜距垂向的同一散射单元构成一个观测信号;
S22、对每个观测信号利用分段弱正交匹配追踪的方法进行重构,得到目标建筑物沿斜距垂向的信息。
作为优选的,所述步骤S21进一步包括:分别确定被截取目标影像的方位向和距离向的起点和终点,将一个散射单元的信号构成一个时间序列的信息。
作为优选的,所述步骤S22具体包括:
S221、获取第m航迹影像表达式Sm;
S222、将Sm表达式进行离散转换,得到用于压缩感知算法的观测向量、测量矩阵和散射稀疏矩阵;
S223、进行斜距垂向信息重构。
作为优选的,所述步骤S221具体包括:
主影像中的任意一个像素的复数值Sm,可以看成真实三维反射情况在x-r平面的投影,即反射信号沿斜距垂向的积分:
式中,γ(s)为散射体反射率沿斜距垂向的分布函数,λ为波长,b⊥m为航迹m、m-1的基线距离,ro为主影像中心斜距。
实际情况中,获取的数据往往基线分布不均匀或基线数目少,利用傅立叶变换方法进行三维成像无法得到理想的结果。根据压缩感知理论,如果信号在某个变换域稀疏或可压缩,可以利用低于奈奎斯特采样率的观测数据对该信号重构,而雷达目标可以由少数几个散射中心描述,即满足目标的稀疏性。因此可以将压缩感知应用到层析SAR三维成像中,对斜距垂向信号进行重构,获取目标建筑在斜距垂向的信息,从而提取出建筑物的高程信息。所述步骤S222具体包括,对Sm的航迹影像进行离散表示:
式中,Sn为第n个目标,N为沿斜距垂向的散射体数目,m=1...M,M<N;将上式通过矩阵表示:
y=Φγ+ξ
式中,y=Sm,为M维观测量,表示M条航迹的观测量;ξ=[ξ1,ξ2,…,ξM]T,为独立分布的高斯白噪声,Φ=[Ф1,Ф2…Фn…ФN],为M×N的测量矩阵, γ=[γ1,γ2…,γn,…,γN]T,为N维散射稀疏矩阵;如图3和图4所示,对观测向量y、测量矩阵Φ、N维散射稀疏矩阵γ进行分段弱正交匹配追踪算法处理,得到重构后的斜距信息利用稀疏矩阵可得到重构信号为γ,从而提取出地物的高程信息。分段弱正交匹配追踪(SWOMP)是OMP改进算法中的一种,每次迭代可以选择多个原子。同时它在选择原子时降低了对测量矩阵的要求,最重要的是此算法不需要输出信号稀疏度K,相比其他需要已知信号稀疏度K的算法具有独到的优势。
由压缩感知理论可知,如果信号在某个变换域稀疏或可压缩,可以利用低于奈奎斯特采样率的观测数据对该信号重构。由于SAR的成像原理,SAR影像获取的斜距垂向信号可以作为稀疏信号,因此可以将压缩感知应用到层析SAR三维成像中,获取目标地物在斜距垂向的信息,从而提取出地物的高程信息。
实施例2
本实施例中方案提供某区域的14景二维SAR影像数据。利用ENVI读取14景目标区域的数据,并将数据保存成MATLAB可以读取的格式来作为准备数据。也可以只选取区域中的单一建筑物的14景影像数据作为准备数据,选取的区域需要记录下它在SAR影像中的坐标位置,方案选取14景北京某酒店的二维SAR影像作为准备数据。
对目标区域的三维点云生成(基于压缩感知SWOMP的TomoSAR三维点云成像算法)
(1)读取和加载数据
将准备数据读取并加载到MATLB中。
(2)确定目标区域位置
分别确定被截取目标影像的方位向和距离向的起点和终点位置。读取以二维数组的形式表示的14景影像数据,再将这14个二维数组构成一个三维的数组。
(3)初始化参数
影像数据参数(光速c、载波频率f、波长lambda、入射角向量InA、Near range向量、Center Range向量、Far Range向量、方位分辨率(azimuth resolution)、距离分辨率(range resolution)、斜距R0、水平距离H0=R0*COS(InA)、基线向量baseline)
高程参数:设置高程预估值
(4)获取斜距垂向信息
对于层析SAR三维成像斜距垂向的信息获取主要有傅立叶变换方法、谱估计算法、压缩感知重构算法。压缩感知重构算法较其他二者有特别的优势,方案选取基于压缩感知的SWOMP对目标地物的斜距垂向信息进行重构。
(5)计算高程信息
利用步骤(4)获取目标地物的斜距信息,再经过进一步计算提取目标的高程信息。根据SWOPM算法流程图3,实验中需要输入的参数有传感矩阵A、观测向量y、迭代次数S,阈值α。其中,S默认值为10,α取值范围为0<α≤1。
(6)绘制三维点云成果
将重构形成的目标地物表示出来,并输出效果图,如图4所示。
最后,本申请的方法仅为较佳的实施方案,并非用于限定本发明的保护范围。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的保护范围之内。
Claims (8)
1.一种层析SAR三维成像方法,其特征在于,包括以下步骤:
S1、获取目标建筑物的二维SAR影像并进行预处理,构成观测数据集;
S2、基于分段弱正交匹配追踪的压缩感知算法重构斜距垂向信息;
S3、获取目标高程信息,显示目标建筑物的三维点云结果。
2.根据权利要求1所述的层析SAR三维成像方法,其特征在于,所述步骤S1具体包括:采集多次航过目标的原始SAR数据,并对原始数据进行二维高分辨成像配准相位补偿的预处理,从而得到可以进行层析SAR成像的观测数据集。
3.根据权利要求2所述的层析SAR三维成像方法,其特征在于,所述二维SAR影像数据包括轨道高度、重访周期、入射角范围、中心斜距、方位分辨率、距离分辨率,单通道场景大小、最大获取长度。
4.根据权利要求2所述的层析SAR三维成像方法,其特征在于,所述步骤S2具体包括:
S21、将多个场景的二维SAR影像数据沿斜距垂向进行孔径合成,将其中一幅做为主影像,其他影像以此影像进行配准,配准后主影像上每个像素沿斜距垂向的同一散射单元构成一个观测信号;
S22、对每个观测信号利用分段弱正交匹配追踪的方法进行重构,得到目标建筑物沿斜距垂向的信息。
5.根据权利要求4所述的层析SAR三维成像方法,其特征在于,所述步骤S21进一步包括:分别确定被截取目标影像的方位向和距离向的起点和终点,将一个散射单元的信号构成一个时间序列的信息。
6.根据权利要求4所述的层析SAR三维成像方法,其特征在于,所述步骤S22具体包括:
S221、获取第m航迹影像表达式Sm;
S222、将Sm表达式进行离散转换,得到用于压缩感知算法的观测向量、测量矩阵和散射稀疏矩阵;
S223、进行斜距垂向信息重构。
7.根据权利要求6所述的层析SAR三维成像方法,其特征在于,所述步骤S221具体包括,主影像中的任意一个像素的复数值Sm,作为真实三维反射情况在x-r平面的投影,即反射信号沿斜距垂向的积分:
式中,γ(s)为散射体反射率沿斜距垂向的分布函数,λ为波长,b⊥m为航迹m、m-1的基线距离,ro为主影像中心斜距。
8.根据权利要求4所述的层析SAR三维成像方法,其特征在于,所述步骤S222具体包括:
对第m航迹影像表达式Sm进行离散表示:
式中,Sn为第n个目标,N为沿斜距垂向的散射体数目,m=1...M,M<<N;将上式通过矩阵表示:
y=Φγ+ξ
式中,y=Sm,为M维观测量,表示M条航迹的观测量;
ξ=[ξ1,ξ2,…,ξM]T,为独立分布的高斯白噪声;
Φ=[Φ1,Φ2,…Φn,…ΦN],为M×N的测量矩阵;
γ=[γ1,γ2…,γn,…,γN]T为稀疏矩阵;
对观测向量y、测量矩阵Φ、稀疏矩阵γ进行分段弱正交匹配追踪算法处理,得到重构后的斜距信息
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611239581.XA CN106872977A (zh) | 2016-12-28 | 2016-12-28 | 一种基于分段弱正交匹配追踪的层析sar三维成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201611239581.XA CN106872977A (zh) | 2016-12-28 | 2016-12-28 | 一种基于分段弱正交匹配追踪的层析sar三维成像方法 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN106872977A true CN106872977A (zh) | 2017-06-20 |
Family
ID=59164188
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201611239581.XA Pending CN106872977A (zh) | 2016-12-28 | 2016-12-28 | 一种基于分段弱正交匹配追踪的层析sar三维成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106872977A (zh) |
Cited By (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108646244A (zh) * | 2018-03-28 | 2018-10-12 | 中科卫星应用德清研究院 | 测量建筑五维形变的分析方法及*** |
CN109212529A (zh) * | 2018-11-02 | 2019-01-15 | 国网四川省电力公司电力应急中心 | 一种输电铁塔监测方法与装置 |
CN109444886A (zh) * | 2018-12-30 | 2019-03-08 | 成都汇蓉国科微***技术有限公司 | 一种双基sar实时成像的并行计算处理方法、***及介质 |
CN111679277A (zh) * | 2020-05-28 | 2020-09-18 | 电子科技大学 | 一种基于sbrim算法的多基线层析sar三维成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US20110175770A1 (en) * | 2009-06-30 | 2011-07-21 | Petros Boufounos | High Resolution SAR Imaging Using Non-Uniform Pulse Timing |
CN102662171A (zh) * | 2012-04-23 | 2012-09-12 | 电子科技大学 | 一种sar层析三维成像方法 |
US8665132B2 (en) * | 2008-12-10 | 2014-03-04 | The United States Of America As Represented By The Secretary Of The Army | System and method for iterative fourier side lobe reduction |
CN105388476A (zh) * | 2015-12-28 | 2016-03-09 | 河南工业大学 | 一种基于联合稀疏模型的层析sar成像方法 |
-
2016
- 2016-12-28 CN CN201611239581.XA patent/CN106872977A/zh active Pending
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8665132B2 (en) * | 2008-12-10 | 2014-03-04 | The United States Of America As Represented By The Secretary Of The Army | System and method for iterative fourier side lobe reduction |
US20110175770A1 (en) * | 2009-06-30 | 2011-07-21 | Petros Boufounos | High Resolution SAR Imaging Using Non-Uniform Pulse Timing |
CN102662171A (zh) * | 2012-04-23 | 2012-09-12 | 电子科技大学 | 一种sar层析三维成像方法 |
CN105388476A (zh) * | 2015-12-28 | 2016-03-09 | 河南工业大学 | 一种基于联合稀疏模型的层析sar成像方法 |
Non-Patent Citations (2)
Title |
---|
THOMAS BLUMENSATH等: "Stagewise Weak Gradient Pursuits", 《IEEE TRANSACTIONS ON SIGNAL PROCESSING》 * |
林珲等: "SAR层析成像的基本原理、关键技术和应用领域", 《测绘地理信息》 * |
Cited By (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108646244A (zh) * | 2018-03-28 | 2018-10-12 | 中科卫星应用德清研究院 | 测量建筑五维形变的分析方法及*** |
CN109212529A (zh) * | 2018-11-02 | 2019-01-15 | 国网四川省电力公司电力应急中心 | 一种输电铁塔监测方法与装置 |
CN109212529B (zh) * | 2018-11-02 | 2021-06-22 | 国网四川省电力公司电力应急中心 | 一种输电铁塔监测方法与装置 |
CN109444886A (zh) * | 2018-12-30 | 2019-03-08 | 成都汇蓉国科微***技术有限公司 | 一种双基sar实时成像的并行计算处理方法、***及介质 |
CN111679277A (zh) * | 2020-05-28 | 2020-09-18 | 电子科技大学 | 一种基于sbrim算法的多基线层析sar三维成像方法 |
CN111679277B (zh) * | 2020-05-28 | 2022-05-03 | 电子科技大学 | 一种基于sbrim算法的多基线层析sar三维成像方法 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Alonso et al. | A novel strategy for radar imaging based on compressive sensing | |
CN103698763B (zh) | 基于硬阈值正交匹配追踪的线阵sar稀疏成像方法 | |
Wang et al. | Three-dimensional reconstruction from a multiview sequence of sparse ISAR imaging of a space target | |
CN103713288B (zh) | 基于迭代最小化稀疏贝叶斯重构线阵sar成像方法 | |
Costantini et al. | A general formulation for redundant integration of finite differences and phase unwrapping on a sparse multidimensional domain | |
CN105677942B (zh) | 一种重复轨道星载自然场景sar复图像数据快速仿真方法 | |
CN102662171B (zh) | 一种sar层析三维成像方法 | |
Fornaro et al. | A null-space method for the phase unwrapping of multitemporal SAR interferometric stacks | |
Gilman et al. | Transionospheric synthetic aperture imaging | |
US20100141508A1 (en) | Method and system for forming an image with enhanced contrast and/or reduced noise | |
CN106872977A (zh) | 一种基于分段弱正交匹配追踪的层析sar三维成像方法 | |
CN109061642A (zh) | 一种贝叶斯迭代重加权稀疏自聚焦阵列sar成像方法 | |
CN108226927A (zh) | 基于加权迭代最小稀疏贝叶斯重构算法的sar成像方法 | |
CN104391295A (zh) | 一种图像熵最优的压缩传感sar稀疏自聚焦成像方法 | |
CN102788979B (zh) | 一种基于后向投影InSAR成像配准的GPU实现方法 | |
CN103018740B (zh) | 一种基于曲面投影的InSAR成像方法 | |
Girard et al. | Sparse representations and convex optimization as tools for LOFAR radio interferometric imaging | |
CN103018741A (zh) | 一种基于后向投影的InSAR成像去平地一体化方法 | |
Yang et al. | Compressed sensing radar imaging with magnitude sparse representation | |
Ren et al. | 3D Imaging Algorithm for Down‐Looking MIMO Array SAR Based on Bayesian Compressive Sensing | |
Douxchamps et al. | Stereoscopic and velocimetric reconstructions of the free surface topography of antidune flows | |
Riel et al. | Uncertainties in forest canopy height estimation from polarimetric interferometric SAR data | |
Ndoye et al. | An MM-Based Algorithm for ${\ell} _ {1} $-Regularized Least-Squares Estimation With an Application to Ground Penetrating Radar Image Reconstruction | |
CN106646470A (zh) | 基于广义正交匹配追踪的层析sar三维点云生成方法 | |
CN113608218A (zh) | 一种基于后向投影原理的频域干涉相位稀疏重构方法 |
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 | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20170620 |
|
RJ01 | Rejection of invention patent application after publication |