CN113391309A - 一种火星探测器雷达径向下视成像方法 - Google Patents
一种火星探测器雷达径向下视成像方法 Download PDFInfo
- Publication number
- CN113391309A CN113391309A CN202110658684.4A CN202110658684A CN113391309A CN 113391309 A CN113391309 A CN 113391309A CN 202110658684 A CN202110658684 A CN 202110658684A CN 113391309 A CN113391309 A CN 113391309A
- Authority
- CN
- China
- Prior art keywords
- mars
- radar
- looking
- radial downward
- probe
- 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
Links
Images
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
-
- 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
-
- 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
-
- Y—GENERAL TAGGING OF NEW TECHNOLOGICAL DEVELOPMENTS; GENERAL TAGGING OF CROSS-SECTIONAL TECHNOLOGIES SPANNING OVER SEVERAL SECTIONS OF THE IPC; TECHNICAL SUBJECTS COVERED BY FORMER USPC CROSS-REFERENCE ART COLLECTIONS [XRACs] AND DIGESTS
- Y02—TECHNOLOGIES OR APPLICATIONS FOR MITIGATION OR ADAPTATION AGAINST CLIMATE CHANGE
- Y02A—TECHNOLOGIES FOR ADAPTATION TO CLIMATE CHANGE
- Y02A90/00—Technologies having an indirect contribution to adaptation to climate change
- Y02A90/10—Information and communication technologies [ICT] supporting adaptation to climate change, e.g. for weather forecasting or climate simulation
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
本发明公开一种火星探测器雷达径向下视成像方法,应用于雷达成像领域,为解决现有技术无法应用于火星探测器径向下视成像模式的问题,本发明首先对火星表面进行雷达径向下视回波精准表征,突破径向下视成像机理性瓶颈;然后构建描述天线方向图与火星表面回波信号关系的可解析二维加权最小二乘代价函数,将径向下视雷达分辨率提升问题转化为最优化求解问题;最后采用快速二维自相关迭代重构方法求取全局最优解,实现雷达径向下视超分辨成像。与传统深空探测器着陆指示方法相比,本发明方法能够全天时、全天候表征火星探测器径向下视方向的地形变化信息,为火星探测器的软着陆选址提供新的技术手段。
Description
技术领域
本发明属于雷达成像领域,特别涉及一种适用于火星探测器动力下降阶段雷达径向下视区域成像技术。
背景技术
在火星探测任务中,探测器的安全软着陆是关键一环。由于火星表面自然光照度变化大,且经常发生尘暴现象,基于光学传感器的着陆手段受天候影响大,难以为探测器提供全天时、全天候着陆引导。而雷达具有强穿透性、不受光照限制等优势,是实现探测器安全软着陆的有效手段。要想实现火星探测器的安全软着陆,就必须实现探测器雷达径向下视成像,为其提供实时的地形探测信息。但是,现有雷达成像技术,比如合成孔径雷达成像技术和多普勒波束锐化技术,在下视区域存在成像盲区,无法应用于火星探测器径向下视成像模式。
针对雷达下视成像难题,文献“Gierull C H.On a concept for an airbornedownward-looking imaging radar.AEU International Journal of Electronics andCommunications,1999,53(6):295-304.”提出一种雷达下视SAR成像方法。该方法在切航向和沿航向上结合了线性阵列的真实孔径与平台运动的合成孔径,形成二维孔径,获得了二维分辨力。但是该方法需要雷达平台进行平动,不适用于火星探测器软着陆阶段的垂直降落运动模式。文献“Foessel-Bunting A,Whittaker W.Mmw-scanning radar for descentguidance and landingsafeguard[C],Proceedings of the 6th InternationalSymposium on Artificial Intelligence,Robotics,and Automation in Space,2001”提出了一种二维实波束扫描雷达成像方法。该方法不依赖于特定的成像模式,但是存在分辨率低的缺点,无法为探测器的软着陆提供精确引导。
发明内容
针对火星探测器雷达径向下视成像难题,本发明提出一种火星探测器雷达径向下视成像方法,分别在成像机理和成像方法提出创新,能够解决传统光学着陆指示技术受火星表面复杂天时和天候影响而失效的问题。
本发明采用的技术方案为:一种火星探测器雷达径向下视成像方法,包括:
S1、对火星表面雷达径向下视回波进行精准表征;
S2、基于步骤S1的精准表征进行火星表面超分辨率重建。
步骤S1根据精准表征得到火星探测器雷达接收回波的表达式为:
其中,
步骤S2具体为:根据火星探测器雷达接收回波的表达式,将火星表面重建问题转化为加权最小二乘最优化问题,定义代价函数为:
通过最小化代价函数,第k个火星表面信源的重建结果为:
其中,t为迭代次数。
本发明的有益效果:本发明首先对火星表面进行雷达径向下视回波精准表征,突破径向下视成像机理性瓶颈;然后构建描述天线方向图与火星表面回波信号关系的可解析二维加权最小二乘代价函数,将径向下视雷达分辨率提升问题转化为最优化求解问题;最后采用快速二维自相关迭代重构方法求取全局最优解,实现雷达径向下视超分辨成像。与传统深空探测器着陆指示方法相比,本发明方法能够全天时、全天候表征火星探测器径向下视方向的地形变化信息,为火星探测器的软着陆选址提供新的技术手段。
附图说明
图1火星探测器雷达径向下视成像几何构型;
图2火星探测器雷达阵列分布图;
图3待成像火星表面原始仿真场景;
图4传统方法成像结果图;
图5本发明方法成像结果图。
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2019上验证正确。下面结合附图和具体实施例对本发明方法做进一步的阐述。
具体实施步骤如下:
步骤一:火星表面雷达径向下视回波精准表征
采用如图1所示的火星探测器雷达径向下视模型,火星探测器平台在着陆过程中,以一速度V前行,火星探测器雷达面阵阵列对火星探测器径向(速度方向)下视(下方)方向进行二维探测成像。
选取如表1所示的雷达仿真系统参数。仿真采用的探测器雷达阵列分布如图2所示,该雷达采用M×N=32×32阵元的均匀面阵,其阵元间距为d=0.005m,本实施例采用火星表面“阿拉伯高地”的光学图像灰度值作为仿真场景中各散射点的幅度,如图3所示。假设该场景有P×Q=100×100=K个非相干信源表示复数集合,第k个信源的俯仰角和方位角分别表示为θk和定义变量uk、vk在这里仅作为中间变量,用于简化后续公式书写,没有实际的含义。
表1仿真参数
阵面的第1个子阵列的接收信号表示为
y1(n)=Axs(n)+e1(n),n=1,…,Nr
其中,第1个子阵列的导向矢量Ax=[ax(v1),ax(v2),…,ax(vK)],而[·]T表示转置。为第1个子阵列的加性高斯白噪声,Nr=10表示快拍数。信源矢量n表示当前快拍,vec(S)表示信源矩阵S的列堆叠。考虑N=64个子阵列,则阵列的接收信号可以表示为
Ax可以理解为面阵阵列某行(某子阵列)的导向矩阵,其维数为M×K,Ay为面阵阵列某列的导向矩阵,维数为N×K。
根据上式,火星探测器雷达各阵元的接收信号可进一步简化为
步骤二:火星表面超分辨率快速重建方法
由步骤一中火星探测器雷达接收回波表达式,火星表面重建问题可转化为一加权最小二乘最优化问题,定义代价函数
其中表示A的第k列元素构成的列向量,Qk为与ak和sk(n)有关的加权矩阵,sk(n)表示s(n)在第k个方向上的信源幅度,根据Kronecker积特性和Hadamard积特性,最小化代价函数f,可得火星表面信源的重建结果为
Λl+1=Zl+ηlΛl
其中,1M为长度为M的全1矢量,1N为长度为N的全1矢量, η为引入的中间变量,l=10为迭代次数。最终可得W(t)=Wl+1。上述迭代式的初始化为W0=0,η0=0,Z0=Y,Λ0=0,其中满足y=vec(Y)。
对整个成像场景S进行10次本方法的迭代处理,二维火星表面地形重建的结果如图5所示。对比图4和图5,采用传统雷达成像方法所得的仿真结果分辨率低,难以为火星探测器选择合适的着陆位置。在参数不变的情况下,采用本发明方法得到的成像结果角分辨率有大幅度提升,可获得地形的准确分布,为火星探测器的着陆选址提供了更丰富的信息。
综上所述,本发明方法相较于传统成像方法,有效地提升了火星探测器雷达径向下视表面地形重建的分辨率。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。对于本领域的技术人员来说,本发明可以有各种更改和变化。凡在本发明的精神和原则之内,所作的任何修改、等同替换、改进等,均应包含在本发明的权利要求范围之内。
Claims (6)
1.一种火星探测器雷达径向下视成像方法,其特征在于,包括:
S1、对火星表面雷达径向下视回波进行精准表征;
S2、基于步骤S1的精准表征进行火星表面超分辨率重建。
2.根据权利要求1所述的一种火星探测器雷达径向下视成像方法,其特征在于,所述雷达采用M×N阵元的均匀面阵。
4.根据权利要求3所述的一种火星探测器雷达径向下视成像方法,其特征在于,步骤S2具体为:根据火星探测器雷达接收回波的表达式,将火星表面重建问题转化为加权最小二乘最优化问题,定义代价函数为:
通过最小化代价函数,第k个火星表面信源的重建结果为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110658684.4A CN113391309B (zh) | 2021-06-15 | 2021-06-15 | 一种火星探测器雷达径向下视成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202110658684.4A CN113391309B (zh) | 2021-06-15 | 2021-06-15 | 一种火星探测器雷达径向下视成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN113391309A true CN113391309A (zh) | 2021-09-14 |
CN113391309B CN113391309B (zh) | 2022-09-09 |
Family
ID=77620807
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202110658684.4A Active CN113391309B (zh) | 2021-06-15 | 2021-06-15 | 一种火星探测器雷达径向下视成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN113391309B (zh) |
Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1537157A (en) * | 1975-06-24 | 1978-12-29 | Thomson Csf | High repetition frequency side-looking pulse radar system |
US20030132873A1 (en) * | 1999-12-14 | 2003-07-17 | Jean-Jacques Berthelier | Method for obtaining underground imagery using a ground-penetrating radar |
US20050134440A1 (en) * | 1997-10-22 | 2005-06-23 | Intelligent Technolgies Int'l, Inc. | Method and system for detecting objects external to a vehicle |
CN102914773A (zh) * | 2012-09-11 | 2013-02-06 | 电子科技大学 | 一种多航过圆周sar三维成像方法 |
CN102967859A (zh) * | 2012-11-14 | 2013-03-13 | 电子科技大学 | 一种前视扫描雷达成像方法 |
JP2014013180A (ja) * | 2012-07-04 | 2014-01-23 | Mitsubishi Electric Corp | レーダ処理装置 |
EP2896971A1 (en) * | 2014-01-16 | 2015-07-22 | Institute of Electronics, Chinese Academy of Sciences | Spaceborne Multi-Channel Synthetic Aperture Radar Imaging Device |
CN105116382A (zh) * | 2015-09-17 | 2015-12-02 | 西安电子科技大学 | 多峰值低多普勒旁瓣相位编码信号的失配滤波器设计方法 |
CN106680817A (zh) * | 2016-12-26 | 2017-05-17 | 电子科技大学 | 一种实现前视雷达高分辨成像的方法 |
CN107255805A (zh) * | 2017-08-01 | 2017-10-17 | 西安电子科技大学 | 基于加权最小二乘的雷达目标rcs的预测方法 |
CN107607945A (zh) * | 2017-08-31 | 2018-01-19 | 电子科技大学 | 一种基于空间嵌入映射的扫描雷达前视成像方法 |
CN107678000A (zh) * | 2016-08-02 | 2018-02-09 | 北京遥感设备研究所 | 一种在轨星体地下探测雷达系统 |
CN111754392A (zh) * | 2020-05-30 | 2020-10-09 | 同济大学 | 一种火星轨道器高分辨率光学相机动态成像仿真方法 |
-
2021
- 2021-06-15 CN CN202110658684.4A patent/CN113391309B/zh active Active
Patent Citations (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
GB1537157A (en) * | 1975-06-24 | 1978-12-29 | Thomson Csf | High repetition frequency side-looking pulse radar system |
US20050134440A1 (en) * | 1997-10-22 | 2005-06-23 | Intelligent Technolgies Int'l, Inc. | Method and system for detecting objects external to a vehicle |
US20030132873A1 (en) * | 1999-12-14 | 2003-07-17 | Jean-Jacques Berthelier | Method for obtaining underground imagery using a ground-penetrating radar |
JP2014013180A (ja) * | 2012-07-04 | 2014-01-23 | Mitsubishi Electric Corp | レーダ処理装置 |
CN102914773A (zh) * | 2012-09-11 | 2013-02-06 | 电子科技大学 | 一种多航过圆周sar三维成像方法 |
CN102967859A (zh) * | 2012-11-14 | 2013-03-13 | 电子科技大学 | 一种前视扫描雷达成像方法 |
EP2896971A1 (en) * | 2014-01-16 | 2015-07-22 | Institute of Electronics, Chinese Academy of Sciences | Spaceborne Multi-Channel Synthetic Aperture Radar Imaging Device |
CN105116382A (zh) * | 2015-09-17 | 2015-12-02 | 西安电子科技大学 | 多峰值低多普勒旁瓣相位编码信号的失配滤波器设计方法 |
CN107678000A (zh) * | 2016-08-02 | 2018-02-09 | 北京遥感设备研究所 | 一种在轨星体地下探测雷达系统 |
CN106680817A (zh) * | 2016-12-26 | 2017-05-17 | 电子科技大学 | 一种实现前视雷达高分辨成像的方法 |
CN107255805A (zh) * | 2017-08-01 | 2017-10-17 | 西安电子科技大学 | 基于加权最小二乘的雷达目标rcs的预测方法 |
CN107607945A (zh) * | 2017-08-31 | 2018-01-19 | 电子科技大学 | 一种基于空间嵌入映射的扫描雷达前视成像方法 |
CN111754392A (zh) * | 2020-05-30 | 2020-10-09 | 同济大学 | 一种火星轨道器高分辨率光学相机动态成像仿真方法 |
Non-Patent Citations (4)
Title |
---|
HAN KUOYE 等: "《Statistical Analysis of the Effects of Virtual Element Position Errors on Airborne Down-Looking LASAR 3-D Imaging》", 《IEEE GEOSCIENCE AND REMOTE SENSING LETTERS》 * |
张永超: "《机载雷达前视虚拟阵列超分辨成像方法研究》", 《中国优秀博硕士学位论文全文数据库(博士)信息科技辑》 * |
王明远 等: "《"火星快车"下视雷达对火星电离层探测简介》", 《物理》 * |
黄旭华 等: "《一种交替发射模式的机载MIMO-SAR下视三维成像方法》", 《西安电子科技大学学报》 * |
Also Published As
Publication number | Publication date |
---|---|
CN113391309B (zh) | 2022-09-09 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108107431B (zh) | 一种圆柱扫描sar三维成像快速实现方法 | |
CN106249237B (zh) | 一种曲线轨迹下大斜视sar频域成像方法 | |
CN104950305B (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN106908787A (zh) | 一种实波束扫描雷达前视角超分辨率成像方法 | |
CN108594228B (zh) | 基于isar图像重聚焦的空间目标姿态估计方法 | |
CN109116311A (zh) | 基于知识辅助稀疏迭代协方差估计的杂波抑制方法 | |
CN103487803B (zh) | 迭代压缩模式下机载扫描雷达成像方法 | |
CN101430380A (zh) | 基于非均匀采样的大斜视角机载sar聚束模式成像方法 | |
CN106291543A (zh) | 一种运动平台扫描雷达超分辨成像方法 | |
CN107966686B (zh) | 一种基于线性调频的时间依赖性频控阵目标探测方法 | |
CN112612024A (zh) | 微波阵列快速成像方法 | |
CN110082764B (zh) | 基于稳健正则化层析方法的sar图像成像方法 | |
CN109298420B (zh) | 一种合成孔径雷达的运动目标迭代最小熵成像方法及装置 | |
CN104898119A (zh) | 一种基于相关函数的动目标参数估计方法 | |
Fortune et al. | Statistical autofocus of synthetic aperture sonar images using image contrast optimisation | |
CN106646466B (zh) | 一种基于主成分分析的加权后向投影算法的成像方法 | |
CN108107427A (zh) | 基于超分辨技术的机载/弹载阵列雷达前视成像方法 | |
CN105824016B (zh) | 运动平台雷达检测超低空目标的稳健空时自适应处理方法 | |
Yuan et al. | Range alignment in ISAR imaging based on deep recurrent neural network | |
CN113534065B (zh) | 一种雷达目标微动特征提取与智能分类方法及系统 | |
Zhang et al. | High‐resolution imaging for the multireceiver SAS | |
CN113391309B (zh) | 一种火星探测器雷达径向下视成像方法 | |
CN112684448B (zh) | 面向多参数组合的星载sar波位计算及生成方法及系统 | |
Wei et al. | Learning-based split unfolding framework for 3-D mmW radar sparse imaging | |
CN111190168B (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 | ||
GR01 | Patent grant | ||
GR01 | Patent grant |