CN107607952B - 基于电磁涡旋波的三维合成孔径雷达成像方法 - Google Patents
基于电磁涡旋波的三维合成孔径雷达成像方法 Download PDFInfo
- Publication number
- CN107607952B CN107607952B CN201710975293.9A CN201710975293A CN107607952B CN 107607952 B CN107607952 B CN 107607952B CN 201710975293 A CN201710975293 A CN 201710975293A CN 107607952 B CN107607952 B CN 107607952B
- Authority
- CN
- China
- Prior art keywords
- dimensional
- image
- synthetic aperture
- course
- aperture radar
- 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.)
- Active
Links
- 238000003384 imaging method Methods 0.000 title claims abstract description 26
- 238000012937 correction Methods 0.000 claims abstract description 13
- 238000013508 migration Methods 0.000 claims abstract description 13
- 230000005012 migration Effects 0.000 claims abstract description 12
- 230000006835 compression Effects 0.000 claims abstract description 7
- 238000007906 compression Methods 0.000 claims abstract description 7
- 238000000034 method Methods 0.000 abstract description 9
- 238000006243 chemical reaction Methods 0.000 abstract description 3
- 230000033001 locomotion Effects 0.000 abstract description 3
- 230000012292 cell migration Effects 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 230000000694 effects Effects 0.000 description 2
- 238000011160 research Methods 0.000 description 2
- 241000791900 Selene vomer Species 0.000 description 1
- 208000004350 Strabismus Diseases 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000004891 communication Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000007123 defense Effects 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000002592 echocardiography Methods 0.000 description 1
- 238000005516 engineering process Methods 0.000 description 1
- 238000002474 experimental method Methods 0.000 description 1
- 239000012467 final product Substances 0.000 description 1
- 238000013507 mapping Methods 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
Images
Landscapes
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明公开了一种基于电磁涡旋波的三维合成孔径雷达成像方法,主要解决了现有二维成像方法仅考虑目标与雷达相对静止的问题。本发明的实现步骤是:(1)发射具有不同轨道角动量的电磁涡旋波;(2)对接收到的回波信号进行距离压缩和距离徙动校正;(3)沿航向进行去斜操作,得到沿航向聚焦图像;(4)以场景中心为参考,对沿航向聚焦后的图像进行多普勒单元徙动校正;(5)在轨道角动量域进行傅里叶变换,得到三维图像,再根据角度转换关系,得到最终的三维合成孔径雷达图像。本发明通过利用电磁涡旋波携带的轨道角动量,结合目标与雷达的相对运动,实现了三维合成孔径雷达成像。
Description
技术领域
本发明属于雷达技术领域,更进一步涉及雷达成像技术领域中的合成孔径雷达(Synthetic Aperture Radar,SAR)三维成像方法。
背景技术
合成孔径雷达因具有全天时、全天候、远距离成像的能力,在战略防御、地形测绘等军用和民用领域得到了广泛应用。传统SAR仅能获取斜平面上的两维图像,且受叠掩和透视收缩等的影响。为了克服两维成像的缺陷,国内外学者提出了多种三维SAR成像的方法。目前三维SAR成像已成为现代雷达的研究热点之一。
下视线性阵列三维SAR系统作为实现三维SAR成像的主要技术之一,目前已受到各国学者的关注。通过利用垂直航向的线性阵列,下视线性阵列三维SAR可实现除传统SAR成像可得到的沿航向和距离向分辨以外,还能获得垂直航向的分辨能力,从而实现对照射场景的三维成像。但是,下视线性阵列三维SAR存在两个主要的问题:垂直航向分辨率受阵列尺寸的影响,与沿航向和距离向分辨率相比通常较低;为避免模糊图像,垂直航向上的阵元间距必须足够小,从而导致阵元数的大大增加。
电磁涡旋波携带的轨道角动量(Orbital Angular Momentum,OAM)受到了无线通信,射电天文学等越来越多领域的重视。近年来,电磁涡旋波被成功运行于二维雷达目标成像。但是目前的研究局限于雷达与目标相对静止,且仅能获取目标的两维图像。
发明内容
本发明针对上述三维SAR成像的问题,提出了一种新的基于电磁涡旋波的三维SAR成像方法。与传统方法不同,本发明充分利用雷达与目标的相对运动,在传统两维SAR成像的基础上,结合电磁涡旋波携带的OAM获取目标的第三维图像,从而得到三维SAR图像。
为实现上述目的,本发明的主要步骤如下:
(1)发射具有不同OAM态的电磁涡旋波;
(2)对接收到的回波信号进行距离压缩和距离徙动校正;
(3)沿航向进行去斜操作,得到沿航向聚焦图像;
(4)以场景中心为参考,对沿航向聚焦后的图像进行多普勒单元徙动校正;
(5)在OAM域进行傅里叶变换,得到三维图像,再根据角度转换关系,得到最终的三维SAR图像。
本发明与现有的技术相比具有以下优点:
本发明利用电磁涡旋波携带的OAM和雷达与场景间的相对运动,克服了传统方法仅能实现二维成像的局限,最终可得到三维SAR图像。
附图说明
图1为本发明的设计方法流程图;
图2为SAR系统的成像几何示意图;
图3为利用本发明方法得到的三维SAR图像;
图4为利用本发明方法得到的第二个点目标的沿锥角方向和沿方位角方向的成像剖面图;
具体实施方式
参照附图1,本发明的具体实施步骤如下:
步骤1,发射具有不同OAM态的电磁涡旋波。
目前有若干种OAM产生和发射的方式,本发明以均匀圆阵为例,经过适当的辐射波形调整,第n个阵元发射的线性调频信号为
步骤2,将接收到的回波信号进行距离压缩和距离徙动校正。
系统成像几何示意图如附图2所示,雷达以速度vs沿x轴方向运动,r为雷达至场景中某一目标的斜距,β,θ和φ分别表示锥角,入射角和方位角。假设观测场景中有一点目标,其极坐标位置为(r0,θ0,φ0),相应的直角坐标位置为(x,y,z),零时刻雷达的坐标位置(0,0,0),时刻t为(vst,0,0)。为简单起见,忽略噪声,对所有阵元接收回波进行累加解调后为
其中,σ为目标的后向散射系数,Jl(·)为第一类l阶贝塞尔函数,k=2π/λ为波数,λ=c/fc为波长,c为电磁波传播速度,al为生成OAM态为l的圆阵半径,wa为沿航迹方向包络,t0=x/vs。r(t)为雷达至目标的瞬时斜距,其值为
其中,β(t)为瞬时锥角,且有β(0)=β0。对于低斜视成像说,瞬时方位角也可写为
由此可将接收回波重写为
对接收回波信号进行传统距离压缩和距离徙动校正后为
其中,pr(τ)为距离压缩后的脉冲包络,τc=2rc/c。
步骤3,沿航向进行去斜操作,得到方位聚焦图像。
对距离压缩后的信号沿时间域t进行去斜操作,可得
再对上式沿时间域t进行傅里叶变换,可得
其中,ft为多普勒频率,pa(ft)为沿航向压缩脉冲包络。至此,已获取距离向和沿航向的压缩图像。由上式可知,OAM态l与多普勒频率ft相互耦合,目标在多普勒域的聚焦位置随着l的不同而不同。
步骤4,以场景中心为参考,对沿航向聚焦后的图像进行多普勒单元徙动校正。
与距离单元徙动类似,这里将多普勒聚焦位置随OAM的变化称为多普勒单元徙动,其徙动量为可以看出,多普勒单元徙动量由OAM,平台速度,斜距和入射角决定。对于给定的l和rc,仅θ0未知。通常情况下,场景中心的入射角θc可以根据系统参数估计得到,而与多普勒分辨率相比,可以忽略,因此可以利用场景中心为参考进行多普勒单元徙动校正。多普勒单元徙动校正后可得
步骤5,在OAM域进行傅里叶变换,得到三维图像,再根据角度转换关系,得到最终的三维SAR图像。
对多普勒单元徙动校正后的回波信号沿OAM域进行傅里叶变换可得
由此便得到了目标场景的三维SAR图像。
下面结合仿真数据实验对本发明的效果做进一步的说明。
1、仿真条件:
三维SAR成像系统的仿真参数如下表所示,并假设在场景中有三个点目标。
三个点目标对应的锥角和方位角分别为(2.66°,86.8°),(0.87°,88.9°)和(-0.85°,90.97°)
2、仿真数据包实验分析:
利用本发明提出的方法,对接收回波进行三维SAR成像后的结果如附图3所示,可以看出三个点目标得到了很好的区分,且聚焦位置与真实位置一致。为进一步给出聚焦效果说明,附图4(a)和(b)分别给出了第二个点目标的沿锥角方向和沿方位角的剖面图,可以看出,目标实现了很好的聚焦。
Claims (1)
1.一种基于电磁涡旋波的三维合成孔径雷达成像方法,其特征在于包括如下步骤:
(1)将接收回波信号进行距离压缩和距离徏动校正;
(2)对距离压缩后的回波信号沿航向进行去斜操作,得到方位聚焦图像;
(3)以场景中心为参考,利用系统参数估计得到场景中心入射角,在此基础上构建多普勒单元徙动量,然后再对沿航向聚焦后的图像进行多普勒单元徙动校正;
(4)对多普勒单元徙动校正后的回波信号在OAM域进行傅里叶变换,得到三维SAR图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710975293.9A CN107607952B (zh) | 2017-10-16 | 2017-10-16 | 基于电磁涡旋波的三维合成孔径雷达成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710975293.9A CN107607952B (zh) | 2017-10-16 | 2017-10-16 | 基于电磁涡旋波的三维合成孔径雷达成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107607952A CN107607952A (zh) | 2018-01-19 |
CN107607952B true CN107607952B (zh) | 2020-08-11 |
Family
ID=61078574
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710975293.9A Active CN107607952B (zh) | 2017-10-16 | 2017-10-16 | 基于电磁涡旋波的三维合成孔径雷达成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107607952B (zh) |
Families Citing this family (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109541594A (zh) * | 2018-11-12 | 2019-03-29 | 中国人民解放军国防科技大学 | 基于涡旋电磁波的条带sar三维成像方法 |
CN109490881A (zh) * | 2018-11-21 | 2019-03-19 | 中国科学院电子学研究所 | 基于涡旋电磁波的干涉sar高程测量系统及测量方法 |
CN110412571B (zh) * | 2019-07-19 | 2023-03-14 | 西安电子科技大学 | 基于电磁涡旋波的合成孔径雷达三维成像方法 |
CN110426707B (zh) * | 2019-08-12 | 2021-05-07 | 中国人民解放军国防科技大学 | 基于轨道角动量的涡旋sar成像方法及成像系统 |
CN111796279B (zh) * | 2020-07-13 | 2022-03-04 | 内蒙古科技大学 | 无源电磁涡旋sar方位向超分辨率成像方法及装置 |
CN112327299B (zh) * | 2021-01-04 | 2021-03-16 | 中国人民解放军国防科技大学 | 电磁涡旋前视目标三维成像方法 |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223573A (zh) * | 2015-10-10 | 2016-01-06 | 电子科技大学 | 宽域高分辨率多目标逆合成孔径雷达成像技术 |
CN105301589A (zh) * | 2015-10-30 | 2016-02-03 | 电子科技大学 | 高分辨宽测绘带sar地面运动目标成像方法 |
CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
-
2017
- 2017-10-16 CN CN201710975293.9A patent/CN107607952B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN105223573A (zh) * | 2015-10-10 | 2016-01-06 | 电子科技大学 | 宽域高分辨率多目标逆合成孔径雷达成像技术 |
CN105301589A (zh) * | 2015-10-30 | 2016-02-03 | 电子科技大学 | 高分辨宽测绘带sar地面运动目标成像方法 |
CN106526589A (zh) * | 2016-12-26 | 2017-03-22 | 浙江大学 | 一种基于涡旋电磁波的雷达目标二维成像方法 |
CN106886020A (zh) * | 2017-03-06 | 2017-06-23 | 中国人民解放军国防科学技术大学 | 一种单天线接收条件下的电磁涡旋成像方法 |
Non-Patent Citations (3)
Title |
---|
"Electromagnetic Vortex Imaging Using Uniform Concentric Circular Arrays";Tiezhu Yuan et al.;《IEEE Antennas and Wireless Progration Letters》;20151012;全文 * |
"Orbital-Augular-Momentum-Based Electromagnetic Vortex Imaging";Kang Liu et al.;《IEEE Antennas and Wireless Propagation Letters》;20141204;全文 * |
"基于电磁涡旋的雷达目标成像";郭桂蓉 等;《国防科技大学学报》;20130630(第06期);全文 * |
Also Published As
Publication number | Publication date |
---|---|
CN107607952A (zh) | 2018-01-19 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107607952B (zh) | 基于电磁涡旋波的三维合成孔径雷达成像方法 | |
Wang et al. | High-resolution imaging using a wideband MIMO radar system with two distributed arrays | |
Long et al. | A DBS Doppler centroid estimation algorithm based on entropy minimization | |
CN109581352B (zh) | 一种基于毫米波雷达的超分辨测角系统 | |
CN109738894B (zh) | 一种用于大视场合成孔径雷达的高斜视多角度成像方法 | |
CN105044719A (zh) | 一种基于圆周sar的太赫兹高精度垂直曲面成像方法 | |
CN105785366B (zh) | 基于稀疏空间谱估计的星载sar数字波束形成方法 | |
CN103454638B (zh) | 一种圆迹合成孔径雷达三维层析成像方法 | |
CN111796279B (zh) | 无源电磁涡旋sar方位向超分辨率成像方法及装置 | |
CN110967692A (zh) | 一种成像方法 | |
CN109490881A (zh) | 基于涡旋电磁波的干涉sar高程测量系统及测量方法 | |
CN105372657A (zh) | 基于回波数据的视频合成孔径雷达运动补偿成像方法 | |
Zhang et al. | Influence of the stop-and-hop assumption on synthetic aperture sonar imagery | |
CN104122552B (zh) | 一种滑动式双站圆周合成孔径雷达成像方法 | |
CN109298417B (zh) | 一种基于雷达信号处理的建筑物内部结构探测方法及装置 | |
Yarlequé et al. | FMCW GPR radar mounted in a mini-UAV for archaeological applications: First analytical and measurement results | |
CN106019280A (zh) | 基于距离多普勒校正的fmcw sar成像方法及装置 | |
CN110879391B (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN111208514A (zh) | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 | |
CN110568410B (zh) | 一种空间频率色散的微波雷达超分辨方法 | |
JP2004198275A (ja) | 合成開口レーダ装置および画像再生方法 | |
Liang et al. | A two-step processing method for diving-mode squint SAR imaging with subaperture data | |
CN109870685B (zh) | 基于改进的rd算法的室内距离向移动sar成像方法 | |
Wang et al. | Study on motion compensation method for W-band UAV MISAR real-time imaging | |
Han et al. | MIMO radar fast imaging algorithm based on sub-image combination |
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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |