CN105785366B - 基于稀疏空间谱估计的星载sar数字波束形成方法 - Google Patents
基于稀疏空间谱估计的星载sar数字波束形成方法 Download PDFInfo
- Publication number
- CN105785366B CN105785366B CN201610229995.8A CN201610229995A CN105785366B CN 105785366 B CN105785366 B CN 105785366B CN 201610229995 A CN201610229995 A CN 201610229995A CN 105785366 B CN105785366 B CN 105785366B
- Authority
- CN
- China
- Prior art keywords
- echo
- signal
- aperture
- arrival
- evacuated space
- 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.)
- Expired - Fee Related
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
- G01S13/904—SAR modes
-
- 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)
Abstract
本发明公开了一种基于稀疏空间谱估计的星载SAR数字波束形成方法,主要解决了现有方法受小样本约束影响导致的接收增益损失的问题。本发明的实现步骤是:(1)对俯仰向各子孔径接收的回波信号进行距离压缩;(2)根据稀疏空间谱估计方法,利用各距离单元回波估计得到每个距离时刻对应的目标场景的回波方向;(3)根据估计得到的回波方向,计算每个距离时刻对应的子孔径的加权矢量,使得形成的接收波束中心准确地指向目标场景所在的位置。本发明将目标场景高程估计问题转换为稀疏空间谱估计问题,降低了现有方法极易受小样本约束影响的限制,实现了复杂地形下的正确波束指向,保证了回波信号的接收增益。
Description
技术领域
本发明属于雷达技术领域,更进一步涉及雷达成像技术领域中的星载合成孔径雷达(Synthetic Aperture Radar,SAR)数字波束形成(Digital Beamforming,DBF)方法。
背景技术
合成孔径雷达因具有全天时、全天候、远距离成像的能力,在战略防御、地形测绘等军用和民用领域得到了广泛应用。其中,结合数字波束形成(Digital Beamforming,DBF),利用俯仰向波束扫描(SCan-On-REceive,SCORE)技术来提高系统信噪比(Signal toNoise Ratio,SNR)已成为现代雷达的研究热点之一。
SCORE技术通过在俯仰向接收端设置多个子孔径,并对其进行时变加权处理,以形成一个等效的高增益窄波束来追踪接收回波。在理想地球模型的假设前提下,传统SCORE技术根据成像几何关系确定出回波时间与波达方向(direction-of-arrival,DOA)的关系,然后再对接收孔径进行时变加权,以使波束中心对准回波位置。但是,在实际场景中,特别是在山区等存在较大地形起伏的区域,若仍使用理想地球模型来计算波达方向将带来较大的偏差,从而导致目标回波接收增益的损失和SNR的恶化。
针对以上问题,各国学者提出了多种方法。F.Bordoni等人提出在星上对各通道接收回波分别进行距离压缩和配准后,再利用MUSIC等子空间类空间谱估计方法来获取回波信号等效的波达方向,从而缓解地形高程对波束指向的影响。冯帆等人和刘尧等人也提出了类似的方法,对接收回波信号进行距离压缩后,利用子空间谱估计方法来改善地形高程的影响。F.He等人提出将各通道接收回波下传至地面系统后,利用MUSIC等空间谱估计方法来估计回波信号的波达方向,从而缓解地形变化对波束形成的影响。这些方法均采用了子空间类测向方法来改善地形起伏的影响,然而子空间类测向方法受样本数的影响。当样本数较少,波达方向的估计性能将大大下降,从而造成SAR图像SNR的损失。
发明内容
本发明针对上述星载SAR DBF成像的问题,提出了一种新的基于稀疏空间谱估计的DBF方法。与传统方法不同,本发明在小样本的情况下,依然能够实现对复杂地形的正确波束指向,从而得到高质量SAR图像。
为实现上述目的,本发明的主要步骤如下:
(1)对俯仰向各子孔径接收的回波信号进行距离压缩;
(2)根据稀疏空间谱估计方法,利用各距离单元回波估计得到每个距离时刻对应的目标场景的波达方向;
(3)根据估计得到的波达方向,计算每个距离时刻对应的子孔径的加权矢量,使得形成的接收波束中心准确地指向目标场景所在的位置。
本发明与现有的技术相比具有以下优点:
本发明利用稀疏空间谱估计的方法对目标场景回波的波达方向进行估计,降低了传统方法极易受小样本约束影响的限制,实现了复杂地形下的正确波束指向,保证了回波信号的接收增益。
附图说明
图1为本发明的设计方法流程图;
图2为SAR系统在俯仰向的观测几何示意图;
图3为利用SCORE法和本发明方法估计得到的波束指向偏差和相应的增益损失随目标场景海拔高度的变化;
图4为利用传统自适应方法估计得到的波束指向偏差和相应的增益损失随目标场景海拔高度的变化;
具体实施方式
参照附图1,本发明的具体实施步骤如下:
步骤1,对俯仰向各子孔径接收的回波信号进行距离压缩。
假设SAR系统沿俯仰向划分为M个均匀分布的子孔径,如附图2所示,其中Hs为卫星平台高度,dm表示第m个子孔径到第1个子孔径的距离,α为天线水平倾角,r1和rm分别表示地面目标P至第1个和第m个子孔径的斜距,h为目标P的海拔高度,θ为目标至第1个子孔径的斜距连线与天线面板法线的夹角,其值离开天线法线向上为正,向下为负,可由下式计算得到
其中,Re为地球半径。根据以上几何关系可得
由于天线尺寸相比目标斜距来说非常小,因此上式可近似为
rm≈r1-dm sinθ (3)
不失一般性,本发明假设第1个子孔径发射信号,所有子孔径同时接收回波信号。
考虑加性白噪声且仅考虑距离向回波,第m个子孔径接收的回波信号进行距离压缩后可写为
其中,τ表示距离快时间,pr(τ)为压缩脉冲包络,对于矩形窗来说,pr(τ)为sinc函数,n(τ)为高斯白噪声。当dmsinθ的大小相对于分辨率来说不可忽略时,可对各子孔径接收回波相对第一个子孔径进行配准。这样上式可改写为
其中,
各子孔径接收回波可用矢量形式表示为
s(τ)=s1(τ)p(θ)+n(τ) (6)
其中,
s(τ)=[s1(τ),s2(τ),…,sM(τ)]T (7)
n(τ)=[n1(τ),n2(τ),…,nM(τ)]T (9)
上标T表示矩阵转置。
步骤2,根据稀疏空间谱估计方法,利用各距离单元回波估计得到每个距离时刻对应的目标场景的波达方向。
考虑回波信号的空域稀疏性,本发明将距离向回波信号构建为
其中,为所有可能波达方向的信号构成的阵列响应矩阵,即空域超完备字典集,θq(q=1,…,Q)的取值范围为雷达波束的照射范围,是稀疏向量,是回波信号s1(τ)的补零扩展形式,若回波信号的方向角为θq,则的第q个值为s1(τ),否则其值为零。
一般来说,Q>M,因此可将目标场景的波达方向估计问题转换为以下稀疏问题
其中,||·||2和||·||1分别为和范数,λ为稀疏规则化参数。求解上述优化问题即可得到然后根据波峰位置得到信号的波达方向。
上式所示的稀疏问题为凸优化问题,现已有多种高效的求解方法。与传统自适应DBF方法相比,本发明不需要利用多个训练样本来估计回波信号的协方差矩阵。当然,训练样本数的提高也有利于本发明方法精度的提高。
步骤3,根据估计得到的波达方向,计算每个距离时刻对应的子孔径的加权矢量,使得形成的接收波束中心准确地指向目标所在的位置。
利用估计得到的波达方向可得到时刻τ各子孔径的加权矢量ω(τ)为
其中,上标*表示取共轭。
在实际中,目标场景地形通常为连续缓慢变化。因此,为了降低运算量,可对距离向回波进行分块,选取每块场景中的某一距离单元估计得到波达方向,然后根据波达方向和斜距估计得到每块场景的地形高度,接着利用估计得到的各块距离向的地形高度拟合出整个场景的高度,最后再根据几何关系求解得到所有距离单元所对应的波达方向。
下面结合仿真数据实验对本发明的效果做进一步的说明。
1、仿真条件:
星载SAR系统的仿真参数如下表所示,并假设在场景中心有一点目标。
2、仿真数据包实验分析:
首先对各子孔径回波信号进行距离压缩,然后选取某一方位向回波进行波达方向估计。图3(a)和(b)分别给出了在不同目标海拔高度的情况下,利用SCORE法和本发明所提方法所引入的波束指向偏差和相应的接收增益损失,实线为利用SCORE法的结果,虚线为利用本发明所提方法的结果。从图3中可以看出,随着目标海拔高度的增加,SCORE法的波束指向偏差与增益损失也逐渐变大。而本发明所提方法的处理性能几乎不随目标海拔高度的变化而变化,也就是说本发明所提方法能适应复杂的地形。此外,需要说明的是本发明所提方法只利用了一个样本。
图4给出了利用传统自适应方法的估计结果,作为对比,这里也只采用了一个样本,可以看出自适应方法存在严重的估计误差。这是因为传统自适应方法要求样本数至少为2N+1,其中N为待估计的信号源个数,当只采用一个样本时,将导致严重的估计错误。
Claims (1)
1.基于稀疏空间谱估计的星载合成孔径雷达数字波束形成方法,其特征在于:包括如下步骤:
(1)对俯仰向各子孔径接收的回波信号进行距离压缩;
(2)根据稀疏空间谱估计方法,利用各距离单元回波估计得到每个距离时刻对应的目标场景的波达方向;
(3)根据估计得到的波达方向,计算每个距离时刻对应的子孔径的加权矢量,使得形成的接收波束中心准确地指向目标场景所在的位置;
其中,目标场景的波达方向与目标场景高程有关,即
其中θ为波达方向,h为目标场景高程;对各子孔径回波进行距离压缩后可用矢量形式表示为
s(τ)=s1(τ)p(θ)+n(τ)
其中,
s(τ)=[s1(τ),s2(τ),…,sM(τ)]T
n(τ)=[n1(τ),n2(τ),…,nM(τ)]T
上标T表示矩阵转置;
Re为地球半径;
Hs为卫星平台高度;
r1表示地面目标P至第1个子孔径的斜距;
α为天线水平倾角;
n(τ)为高斯白噪声;
所述稀疏空间谱估计方法按照以下方式进行:
将目标高程估计问题转换为回波信号的波达方向估计问题,更进一步地转换为求解以下稀疏空间谱估计问题
其中,为所有可能波达方向的信号构成的阵列响应矩阵,即空域超完备字典集,θq的取值范围为雷达波束的照射范围,q=1,…,Q;是稀疏向量,是回波信号s1(τ)的补零扩展形式,若回波信号的方向角为φq,则的第q个值为s1(τ),否则其值为零;||·||2和||·||1分别为l2和l1范数,λ为稀疏规则化参数;求解上述稀疏空间谱估计问题即可得到然后根据波峰位置得到信号的波达方向θ,最后再根据θ和h的关系得到目标高程;s表示距离向回波信号。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610229995.8A CN105785366B (zh) | 2016-04-13 | 2016-04-13 | 基于稀疏空间谱估计的星载sar数字波束形成方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610229995.8A CN105785366B (zh) | 2016-04-13 | 2016-04-13 | 基于稀疏空间谱估计的星载sar数字波束形成方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105785366A CN105785366A (zh) | 2016-07-20 |
CN105785366B true CN105785366B (zh) | 2019-02-26 |
Family
ID=56396532
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610229995.8A Expired - Fee Related CN105785366B (zh) | 2016-04-13 | 2016-04-13 | 基于稀疏空间谱估计的星载sar数字波束形成方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105785366B (zh) |
Families Citing this family (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN106291489B (zh) * | 2016-08-23 | 2019-03-01 | 中国人民解放军国防科学技术大学 | 适用于多种发射信号波形的合成孔径雷达回波仿真方法 |
CN106814343B (zh) * | 2017-01-22 | 2019-04-30 | 河海大学 | 一种分步估计的空域信号空间谱估计方法 |
CN107741586B (zh) * | 2017-09-29 | 2021-04-23 | 王辉 | 基于DBF-TOPS加权的星载Ka InSAR信号处理方法 |
CN108318879B (zh) * | 2017-12-29 | 2021-08-31 | 哈尔滨工业大学 | 基于iaa谱估计技术的isar图像横向定标方法 |
CN110488281B (zh) * | 2019-07-25 | 2021-07-09 | 中国科学院电子学研究所 | 一种大带宽dbf-sar色散校正方法 |
CN113567983B (zh) * | 2021-07-23 | 2024-01-16 | 中国人民解放军海军航空大学 | 一种利用滑轨的雷达合成孔径稀疏成像方法及系统 |
CN115856788B (zh) * | 2023-03-01 | 2023-05-16 | 中国科学院空天信息创新研究院 | 一种基于通道相消的俯仰向多通道sar干扰抑制方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996046A (zh) * | 2005-12-28 | 2007-07-11 | 中国科学院电子学研究所 | 距离向多孔径接收宽测绘带合成孔径雷达设计方法 |
US7397418B1 (en) * | 2006-06-05 | 2008-07-08 | Sandia Corporation | SAR image formation with azimuth interpolation after azimuth transform |
EP2725382A1 (de) * | 2012-10-26 | 2014-04-30 | Astrium GmbH | Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung |
CN104865568A (zh) * | 2015-06-02 | 2015-08-26 | 西安电子科技大学 | 基于稀疏重构的宽带雷达高速群目标分辨方法 |
CN105223573A (zh) * | 2015-10-10 | 2016-01-06 | 电子科技大学 | 宽域高分辨率多目标逆合成孔径雷达成像技术 |
CN105301589A (zh) * | 2015-10-30 | 2016-02-03 | 电子科技大学 | 高分辨宽测绘带sar地面运动目标成像方法 |
-
2016
- 2016-04-13 CN CN201610229995.8A patent/CN105785366B/zh not_active Expired - Fee Related
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN1996046A (zh) * | 2005-12-28 | 2007-07-11 | 中国科学院电子学研究所 | 距离向多孔径接收宽测绘带合成孔径雷达设计方法 |
US7397418B1 (en) * | 2006-06-05 | 2008-07-08 | Sandia Corporation | SAR image formation with azimuth interpolation after azimuth transform |
EP2725382A1 (de) * | 2012-10-26 | 2014-04-30 | Astrium GmbH | Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung |
CN104865568A (zh) * | 2015-06-02 | 2015-08-26 | 西安电子科技大学 | 基于稀疏重构的宽带雷达高速群目标分辨方法 |
CN105223573A (zh) * | 2015-10-10 | 2016-01-06 | 电子科技大学 | 宽域高分辨率多目标逆合成孔径雷达成像技术 |
CN105301589A (zh) * | 2015-10-30 | 2016-02-03 | 电子科技大学 | 高分辨宽测绘带sar地面运动目标成像方法 |
Non-Patent Citations (2)
Title |
---|
Performance Analysis for Multichannel HRWS SAR Systems Based on STAP Approach;Taoli Yang et.al;《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》;20131010;第10卷(第6期);第1409-1413页 * |
基于Capon谱估计的星载SAR自适应DBF研究;冯帆等;《雷达学报》;20140227;第3卷(第1期);第53-60页 * |
Also Published As
Publication number | Publication date |
---|---|
CN105785366A (zh) | 2016-07-20 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105785366B (zh) | 基于稀疏空间谱估计的星载sar数字波束形成方法 | |
CN107831480B (zh) | 弹载雷达和差通道稳健自适应杂波抑制方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN104316924B (zh) | 机载超高分辨率sar反投影图像的自聚焦运动补偿方法 | |
CN104615854B (zh) | 一种基于稀疏约束的波束展宽和旁瓣抑制方法 | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
CN104950305B (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN105629206B (zh) | 导向矢量失配下的机载雷达稳健空时波束形成方法及系统 | |
CN105445709B (zh) | 一种稀布阵列近场无源定位幅相误差校正方法 | |
CN103837861B (zh) | 基于特征子空间的子阵级线性约束自适应波束形成方法 | |
CN105301589B (zh) | 高分辨宽测绘带sar地面运动目标成像方法 | |
CN107607952B (zh) | 基于电磁涡旋波的三维合成孔径雷达成像方法 | |
CN111880180A (zh) | 一种高分辨运动舰船sar成像的自聚焦方法 | |
CN106597442A (zh) | 一种方位多通道脉内聚束sar成像方法 | |
CN107728112A (zh) | 在目标导向矢量严重失配情况下的稳健波束形成方法 | |
CN114594478B (zh) | 基于星载Ka波段SAR系统的船只目标干涉检测方法 | |
CN113671492A (zh) | 一种面向机动平台前视成像的samp重构方法 | |
CN111722225A (zh) | 基于先验相位结构信息的双基sar两维自聚焦方法 | |
CN108562898B (zh) | 一种前侧视sar的距离和方位两维空变自聚焦方法 | |
CN104698448B (zh) | 运动平台下基于流形分离的共形阵列稳健估角方法 | |
CN105824016A (zh) | 运动平台雷达检测超低空目标的稳健空时自适应处理方法 | |
CN107918112B (zh) | 基于有源-无源雷达协同的联合收发波束形成方法 | |
CN104020465B (zh) | 基于八单元小孔径圆阵天线的外辐射源雷达测角方法 | |
CN103248412B (zh) | 一种卫星导航多波束形成的方法 | |
CN116679306A (zh) | 一种基于系数修正的扩展孔径声纳成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
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: 20190226 Termination date: 20200413 |