CN105301589B - 高分辨宽测绘带sar地面运动目标成像方法 - Google Patents
高分辨宽测绘带sar地面运动目标成像方法 Download PDFInfo
- Publication number
- CN105301589B CN105301589B CN201510726789.3A CN201510726789A CN105301589B CN 105301589 B CN105301589 B CN 105301589B CN 201510726789 A CN201510726789 A CN 201510726789A CN 105301589 B CN105301589 B CN 105301589B
- Authority
- CN
- China
- Prior art keywords
- moving target
- echo
- channel
- signal
- doppler
- 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
- 238000003384 imaging method Methods 0.000 title claims abstract description 28
- 238000000034 method Methods 0.000 claims abstract description 24
- 238000001228 spectrum Methods 0.000 claims abstract description 22
- 239000013598 vector Substances 0.000 claims description 11
- 239000011159 matrix material Substances 0.000 claims description 4
- 238000012935 Averaging Methods 0.000 claims description 2
- 238000002592 echocardiography Methods 0.000 claims 2
- 230000003595 spectral effect Effects 0.000 abstract description 8
- 230000009466 transformation Effects 0.000 abstract 1
- 238000004088 simulation Methods 0.000 description 5
- 238000005457 optimization Methods 0.000 description 3
- 230000006835 compression Effects 0.000 description 2
- 238000007906 compression Methods 0.000 description 2
- 238000007796 conventional method Methods 0.000 description 2
- 239000000654 additive Substances 0.000 description 1
- 230000000996 additive effect Effects 0.000 description 1
- 238000004458 analytical method Methods 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 238000013508 migration Methods 0.000 description 1
- 230000005012 migration Effects 0.000 description 1
- 238000012545 processing Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 230000003068 static effect Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000017105 transposition Effects 0.000 description 1
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/9029—SAR image post-processing techniques specially adapted for moving target detection within a single SAR image or within multiple SAR images taken at the same time
-
- 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
-
- 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/02—Systems using reflection of radio waves, e.g. primary radar systems; Analogous systems
- G01S13/50—Systems of measurement based on relative movement of target
- G01S13/58—Velocity or trajectory determination systems; Sense-of-movement determination systems
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于稀疏空间谱估计的高分辨宽测绘带SAR地面运动目标成像方法,主要解决了现有方法在估计目标速度时引入的巨大运算量或谱分量增益损失的问题。本发明的实现步骤是:(1)对回波信号进行方位向傅里叶变换;(2)根据稀疏空间谱估计方法,利用某一多普勒单元的回波信号估计得到运动目标的速度;(3)根据估计得到的目标速度,得到目标的无模糊信号;(4)根据估计得到的目标速度,对无模糊信号进行传统高分辨运动目标成像。本发明将运动目标速度估计问题转换为稀疏空间谱估计问题,仅需较小的运算量即可精确估计得到运动目标速度,从而保证了运动目标多普勒频谱的精确重构,实现了高分辨宽测绘带运动目标SAR成像。
Description
技术领域
本发明属于雷达技术领域,更进一步涉及雷达成像技术领域中的高分辨率宽测绘带合成孔径雷达(Synthetic Aperture Radar,SAR)地面运动目标成像方法。
背景技术
合成孔径雷达因具有全天时、全天候、远距离成像的能力,在战略防御、地形测绘等军用和民用领域得到了广泛应用。其中,利用SAR获取运动目标成像结果已成为现代雷达的研究热点之一。
同时实现高分辨宽测绘带成像是雷达成像技术的追求目标。然而,传统星载 SAR***受最小天线面积限制,无法同时获取高分辨宽测绘带图像。结合数字波束形成技术,方位多通道SAR***通过发射低脉冲重复频率(Pulse Repetition Frequency,PRF)信号可克服此限制,但是低PRF将导致回波信号产生严重的多普勒模糊。针对地面静止目标,现有多种方法可实现多普勒模糊抑制从而实现高分辨宽测绘带SAR成像。但是,对于地面运动目标来说,由于其运动参数未知,无法直接精确构造运动目标的模糊重构函数,从而导致处理性能下降。
针对以上运动目标成像存在的问题,各国学者提出了多种方法。Li等人提出通过搜索空时平面内的运动目标区,得到运动目标所有可能的谱分量,但并未给出相应的判断准则。Baumgartner等人根据最大信杂噪比准则搜索运动目标的速度,但由于该方法需要对所有可能的运动目标速度都进行一次全孔径成像操作,因此运算量会大大增加。Yang和Zhang等人假设运动目标谱分量位于杂波谱分量中间,尽量可能地正确约束运动目标谱分量,虽然该方法可以避免繁重的搜索操作,但当运动目标的运动方向偏离该导向方向时,该方法将无法正确对运动目标的谱分量进行约束,从而无法对模糊分量进行抑制并导致运动目标增益的损失。
发明内容
本发明针对上述高分辨宽测绘带SAR地面运动目标成像方法面临的问题,提出了一种基于稀疏空间谱估计的高分辨宽测绘带SAR地面运动目标成像方法。与传统方法不同,本发明没有采用搜索或近似假设,并能以较小的运算量实现对地面运动目标的高质量成像。
为实现上述目的,本发明的主要步骤如下:
(1)对各通道接收的回波信号进行方位向傅里叶变换;
(2)根据稀疏空间谱估计方法,利用某一多普勒单元回波估计得到运动目标的速度;
(3)根据估计得到的运动目标速度,构造运动目标的频谱重构函数,得到运动目标的无模糊全带宽信号;
(4)根据估计得到的运动目标速度,对无模糊的全带宽回波信号进行传统高分辨运动目标成像,得到地面运动目标高分辨率宽测绘带图像。
本发明与现有的技术相比具有以下优点:
本发明利用稀疏空间谱估计的方法直接对运动目标的速度进行精确估计,从而避免了传统搜索操作引入的巨大运算量,另外运动目标速度的精确估计在保证运动目标的多普勒频谱精确重构的同时,克服了传统方法带来的运动目标谱分量增益的损失。
附图说明
图1为本发明的设计方法流程图;
图2为在斜距平面下的SAR***观测几何示意图;
图3为利用回波信号的某一多普勒单元进行空间谱估计的结果;
图4为利用传统方法得到的运动目标仿真成像结果;
图5为利用本发明方法的仿真成像结果。
具体实施方式
参照附图1,本发明的具体实施步骤如下:
步骤1,对各通道接收的回波信号进行方位向傅里叶变换。
假设SAR***共有M个沿航向均匀分布的接收通道,如图2所示,x轴为平台运动速度方向,y轴为斜距方向,雷达到地面目标的斜距矢量与零多普勒面的夹角φ称为锥角。不失一般性,本发明假设SAR***为正侧视成像。对于方位多通道高分辨宽测绘带SAR***来说,各通道接收回波均可等效为参考接收通道(如第1个通道)接收回波的时延。为保证宽距离测绘带,SAR***采用低PRF,这样回波信号将发生多普勒模糊。考虑加性白噪声,第m个通道接收的回波信号进行方位向傅里叶变换后可写为
其中,m=1,2,…,M,M为方位接收通道个数,τ为距离时间,fp为***PRF, fd∈[-fp/2,fp/2]为多普勒频率,vs为平台运动速度,a(τ,fd)为第1个通道的接收回波,Δxm为第m个通道相对第1个通道的距离,L=(N-1)/2,N是回波信号的多普勒模糊数,zm(τ,fd)为第m个通道的噪声。地面运动目标回波的多普勒频率与锥角φ的关系如下
其中,λ为载波频率,vt为运动目标的垂直航向速度,
忽略τ,式( 1 ) 可用矢量形式表示为
其中,
s(τ,fd)=[s1(τ,fd),s2(τ,fd),…,sM(τ,fd)]T, (4)
a(τ,fd)=[a(τ,fd-L·fp),a(τ,fd-(L-1)·fp),…,a(τ,fd+L·fp)]T, (5)
P(fd)=[p-L(fd),p-L+1(fd),…,pL(fd)], (6)
z(τ,fd)=[z1(τ,fd),z2(τ,fd),…,zM(τ,fd)]T, (8)
上标T表示矩阵转置。由此可知,各通道接收的多普勒模糊回波信号可等效看作为来自不同方向的信号叠加。与地面固定杂波不同,由于运动目标的速度vt未知,所以无法精确获知阵列流形矩阵P(fd),从而无法直接利用已有方法对回波信号进行频谱重构。
步骤2,根据稀疏空间谱估计方法,利用某个多普勒单元回波信号估计得到运动目标的速度。
根据式(2)所示的多普勒频率与锥角的关系,本发明将运动目标速度估计问题转换为信号的波达方向(direction-of-arrival,DOA)估计问题,也即空间谱估计问题。根据SAR***参数可计算出多普勒模糊数N,也即待估计DOA信号的个数。更进一步,由于运动目标空间谱在空间域上为稀疏分布,因此上述DOA 估计问题可转换为以下稀疏问题
其中,
Y=[s(τ1,fd),…,s(τK,fd)], (10)
Y为各通道接收的快拍信号,K为快拍个数,在实际中,可以利用相邻的K个距离单元回波替代。为所有可能方向的回波信号构成的阵列流形,φq (q=1,…,Q)的取值范围为雷达波束的照射范围。为Q×1维稀疏向量,若第n(n=1,…,N)个信号的方向角为φq,则的第q个值为 a(τ,fd-ln·fp),其中ln=n-L-1,否则其值为零。一般来说,Q>>M>N。由此可将代价函数写为
其中,||·||f为Frobenius范数,下标 q表示向量的第q个值,β为稀疏规则化参数。求解上述优化问题即可得到然后根据波峰位置得到信号的波达方向,再根据下式计算出运动目标的速度
其中,fd为所选取的多普勒单元,φq为估计得到的波达方向。在实际操作中,为了减小噪声的影响,可以选择多个多普勒单元估计得到多个vt,然后将其平均求得最终的vt。
式(15)所示的优化问题为凸优化问题,现已有多种成熟的求解方法。需要说明的是,对于SAR***来说,M和N的值都较小,而且在实际操作中,仅需要对式(15)进行若干几次求解,因此相比传统的搜索方法,本发明的运算量将大大降低。
步骤3,根据上一步估计得到的运动目标速度,构造运动目标的多普勒频谱重构函数,得到运动目标的无模糊全带宽信号。
利用步骤2估计得到的运动目标速度,即可正确构造出运动目标回波信号的阵列流形P(fd),然后可采用传统的模糊重构方法得到无模糊的运动目标全带宽信号。
步骤4,根据估计得到的运动目标速度,对无模糊的全带宽回波信号进行传统高分辨运动目标成像,得到地面运动目标高分辨率宽测绘带图像。
利用步骤2估计得到的运动目标速度,对回波信号进行距离压缩和距离徙动校正,然后再进行方位压缩,得到地面运动目标的高分辨宽测绘带SAR图像。
下面结合仿真数据实验对本发明的效果做进一步的说明。
1、仿真条件:
星载方位多通道SAR***的仿真参数如下表所示。多普勒模糊次数约为3,由于运动目标存在垂直航向速度,其多普勒频谱将偏离杂波谱约682Hz。
2、仿真数据包实验分析:
首先对各通道回波信号进行方位向傅里叶变换,然后选取某一多普勒单元回波进行稀疏空间谱估计。图3给出了利用第1360个多普勒单元回波进行空间谱估计的结果,相邻的10个距离单元作为快拍数。可以看出,利用本发明的方法可以准确估计得到运动目标的空间谱,由此估计得到的目标垂直航向速度为 10.53m/s,对应的多普勒偏移约为677Hz,与真实值非常接近。图4给出了假设运动目标谱偏离杂波谱fp/2所构造的频谱重构函数得到的运动目标成像结果的方位剖面图,可以看出成像结果存在较为严重的方位模糊。图5(a)给出了利用本发明所述的方法得到的运动目标成像结果,图5(b)为成像结果的局部放大图,可以看出运动目标的多普勒模糊分量得到了有效的抑制,并实现了很好的聚焦。
Claims (2)
1.高分辨宽测绘带合成孔径雷达地面运动目标成像方法,包括如下步骤:
(1)对各通道接收的回波信号进行方位向傅里叶变换;
(2)求解以下稀疏空间谱估计问题
其中,||·||f为Frobenius范数,||·||1为l1范数,Y为各通道接收的快拍信号,为所有可能方向的回波信号构成的阵列流形,为目标回波信号构成的稀疏矩阵,为各通道接收回波的快拍信号向量的l2范数组成的回波稀疏向量,β为稀疏规则化参数,求解上述稀疏空间谱估计问题即可得到然后根据波峰位置得到信号的波达方向,再根据下式计算出运动目标的速度
其中,λ为载波频率,fd为多普勒频率,l∈[-L,L],L=(N-1)/2,N是回波信号的多普勒模糊数,fp为***PRF,vs为平台运动速度,φq为估计得到的波达方向,在实际操作中,为了减小噪声的影响,可以选择多个多普勒单元估计得到多个vt,然后将其平均求得最终的vt;
(3)根据估计得到的运动目标速度,构造运动目标的频谱重构函数,得到运动目标的无模糊全带宽信号;
(4)根据估计得到的运动目标速度,对无模糊的全带宽信号进行传统高分辨运动目标成像,得到高分辨宽测绘带运动目标SAR图像。
2.根据权利要求1所述的高分辨宽测绘带合成孔径雷达地面运动目标成像方法,其特征在于:各通道接收的多普勒模糊回波信号可等效看作为来自不同方向的信号叠加,对各通道回波进行方位向傅里叶变换后可用矢量形式表示为
s(τ,fd)=P(fd)a(τ,fd)+z(τ,fd)
其中,τ为距离时间,fd为多普勒频率,s(τ,fd)为各通道接收回波构成的向量,P(fd)为各通道回波构成的阵列流形矩阵,a(τ,fd)为第1个通道接收的多普勒模糊回波信号构成的向量,z(τ,fd)为各通道噪声构成的向量。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510726789.3A CN105301589B (zh) | 2015-10-30 | 2015-10-30 | 高分辨宽测绘带sar地面运动目标成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510726789.3A CN105301589B (zh) | 2015-10-30 | 2015-10-30 | 高分辨宽测绘带sar地面运动目标成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN105301589A CN105301589A (zh) | 2016-02-03 |
CN105301589B true CN105301589B (zh) | 2018-10-23 |
Family
ID=55199068
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510726789.3A Expired - Fee Related CN105301589B (zh) | 2015-10-30 | 2015-10-30 | 高分辨宽测绘带sar地面运动目标成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN105301589B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105785366B (zh) * | 2016-04-13 | 2019-02-26 | 电子科技大学 | 基于稀疏空间谱估计的星载sar数字波束形成方法 |
US10539672B2 (en) * | 2016-12-09 | 2020-01-21 | GM Global Technology Operations LLC | Doppler ambiguity resolution at high signal to noise ratio |
CN106872954B (zh) * | 2017-01-23 | 2019-07-09 | 西安电子科技大学 | 一种高超声速平台杂波抑制和运动目标成像方法 |
CN107607952B (zh) * | 2017-10-16 | 2020-08-11 | 电子科技大学 | 基于电磁涡旋波的三维合成孔径雷达成像方法 |
CN111948653B (zh) * | 2020-07-31 | 2022-03-18 | 上海卫星工程研究所 | 基于p波段合成孔径雷达的林下目标检测方法及*** |
CN112946614A (zh) * | 2021-03-22 | 2021-06-11 | 哈尔滨工程大学 | 一种基于qga的多普勒频移目标运动参数的解算方法 |
Family Cites Families (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN100523865C (zh) * | 2005-12-28 | 2009-08-05 | 中国科学院电子学研究所 | 确定距离向多孔径接收宽测绘带合成孔径雷达子天线间距的方法 |
US7397418B1 (en) * | 2006-06-05 | 2008-07-08 | Sandia Corporation | SAR image formation with azimuth interpolation after azimuth transform |
DE102012021010B4 (de) * | 2012-10-26 | 2022-02-03 | Airbus Defence and Space GmbH | Synthetisches Apertur Radar zur gleichzeitigen Bildaufnahme und Bewegtzielerkennung |
CN104865568B (zh) * | 2015-06-02 | 2017-05-24 | 西安电子科技大学 | 基于稀疏重构的宽带雷达高速群目标分辨方法 |
-
2015
- 2015-10-30 CN CN201510726789.3A patent/CN105301589B/zh not_active Expired - Fee Related
Also Published As
Publication number | Publication date |
---|---|
CN105301589A (zh) | 2016-02-03 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105301589B (zh) | 高分辨宽测绘带sar地面运动目标成像方法 | |
CN109975807B (zh) | 一种适用于毫米波车载雷达的降维子空间测角方法 | |
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN106093870B (zh) | 高超声速飞行器下降段的sar-gmti杂波抑制方法 | |
CN105445704B (zh) | 一种sar图像中的雷达动目标抑制方法 | |
CN103383449B (zh) | 基于esprit算法的机载雷达近程杂波抑制方法 | |
CN108693511B (zh) | 时分复用mimo雷达的运动目标角度计算方法 | |
CN111007503B (zh) | 基于频率谱精确定位的运动目标聚焦和定位方法及*** | |
CN104698431B (zh) | 基于模糊分量doa估计的多通道sar方位解模糊方法 | |
CN106772253B (zh) | 一种非均匀杂波环境下的雷达杂波抑制方法 | |
CN110109113B (zh) | 一种基于级联对消的双基前视sar非平稳杂波抑制方法 | |
CN105785366B (zh) | 基于稀疏空间谱估计的星载sar数字波束形成方法 | |
CN114545411B (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN115128608B (zh) | 车载多普勒分集mimo合成孔径雷达时域成像方法 | |
CN111856461A (zh) | 基于改进pfa的聚束sar成像方法及其dsp实现 | |
CN106291547B (zh) | 基于天线方向图辅助的多普勒模糊分量自适应抑制方法 | |
CN109143236B (zh) | 适用于复杂飞行轨迹的双基聚束sar大场景成像方法 | |
CN109143235B (zh) | 一种双基前视合成孔径雷达地面运动目标检测方法 | |
CN106526544B (zh) | 基于高超声速平台的mimosar杂波抑制方法 | |
CN108732555A (zh) | 一种获取运动误差向量的方法以及自动驾驶阵列微波成像运动补偿的方法 | |
Fu et al. | ISAR translational motion compensation with simultaneous range alignment and phase adjustment in low SNR environments | |
CN113960597A (zh) | 一种多目标isar成像时间点自动选取方法 | |
CN104076361B (zh) | 无人机载战场监视雷达超分辨广域成像方法 | |
CN111638516B (zh) | 基于双频共轭处理技术的太赫兹频段sar运动补偿算法 | |
Han et al. | MIMO radar fast imaging algorithm based on sub-image combination |
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 |
Granted publication date: 20181023 Termination date: 20201030 |
|
CF01 | Termination of patent right due to non-payment of annual fee |