CN112630778A - 一种超高音速目标成像方法 - Google Patents
一种超高音速目标成像方法 Download PDFInfo
- Publication number
- CN112630778A CN112630778A CN202011274953.9A CN202011274953A CN112630778A CN 112630778 A CN112630778 A CN 112630778A CN 202011274953 A CN202011274953 A CN 202011274953A CN 112630778 A CN112630778 A CN 112630778A
- Authority
- CN
- China
- Prior art keywords
- target
- hypersonic
- speed
- frames
- dimensional
- 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
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)
- Measurement Of Velocity Or Position Using Acoustic Or Ultrasonic Waves (AREA)
Abstract
本发明公开一种超高音速目标成像方法,真实可靠、简单易行。本发明方法,包括:(10)亚音速目标数据处理:将实测亚音速目标总回波数据处理并保存到各帧的二维采样矩阵中;(20)亚音速目标运动一维距离成像:提取两帧回波数据分别到二维矩阵中,通过IFFT变换和图像拼接,得到两帧运动距离像;(30)亚音速目标静止成像:对目标回波进行速度补偿,再次通过IFFT变换和像拼接,得到亚音速目标静止一维距离距离像;(40)等效超高音速运动一维距离像:根据高速目标信息等效出两帧的超高音速目标运动一维距离像;(50)超高音速目标静止一维距离成像:对超高音速目标回波作速度补偿,IFFT变换和像拼接,得到静止一维距离像。
Description
技术领域
本发明属于制导引信技术领域,特别是一种超高音速目标成像方法。
背景技术
高分辨距离像是用宽带雷达信号获取的目标散射点回波在雷达射线上投影的向量和,它提供了目标散射点沿距离方向的分布信息,其特点是通过发出某一波长的高频信号,通过反射成像时间和位置,从而得出高分辨距离像,具有目标重要的结构特征,对目标识别与分类十分有价值,因而成为雷达自动目标识别的新技术。
然而,现有技术存在的问题是:对于大多数高分辨雷达进行实际测量目标时,绝大多数目标运动速度都处于亚音速状态,很难接触到超高音速目标的测试,使得人们在超高音速的高分辨距离像测试部分进展较慢了解较少,无法探究目标以超高音速运动状态下的成像情况以及和亚音速的目标像的区别和联系。
发明内容
本发明的目的在于提供一种超高音速目标成像方法,真实可靠、简单易行。
实现本发明目的的技术解决方案为:
一种超高音速目标成像方法,包括如下步骤:
(10)亚音速目标数据处理:将实测亚音速目标总回波数据处理为各帧回波数据点集并保存到各帧的二维采样矩阵中;
(20)亚音速目标运动一维距离成像:提取两帧回波数据分别到二维矩阵中,通过IFFT变换和图像拼接,得到两帧运动距离像;
(30)亚音速目标静止成像:设置回波距离门,根据两帧运动距离像及熵值测速法测出两帧下的目标速度值,并对目标回波进行速度补偿,再次通过IFFT变换和像拼接,得到亚音速目标静止一维距离距离像;
(40)等效超高音速运动一维距离像:设置超高音速的速度值,分别计算两帧静止距离像对应的超高音速目标相应位置及幅值大小,并将静止目标回波乘以高速等效因子,获得超高音速相位信息,根据高速目标位置、幅度和相位信息等效出两帧的超高音速目标运动一维距离像;
(50)超高音速目标静止一维距离成像:用熵值测速法测出两帧运动距离像下的超高音速的速度大小,并对超高音速目标回波作速度补偿,再通过IFFT变换和像拼接算法得到两帧的超高音速目标静止一维距离像。本发明与现有技术相比,其显著优点有:
1、等效真实:根据当前测出的亚音速目标的成像形状还有周边环境的噪声和背景杂波有依据地等效出超高音速环境的背景噪声和杂波以及成像质量。
2、实现简单:根据真实环境及亚音速目标仿真出超高音速对应的成像,实施简单,大大地减轻了实际测试的困难程度。
下面结合附图和具体实施方式对本发明作进一步的详细描述。
附图说明
图1为本发明超高音速目标成像方法的主流程图。
图2为图1中亚音速目标运动一维距离成像步骤的流程图。
图3为图1中亚音速目标静止成像步骤的流程图。
图4为图1中等效超高音速运动一维距离像步骤的流程图。
图5为超高音速目标静止一维距离成像步骤的流程图。
图6为IFFT和像拼接后一维距离像。
图7为亚音速目标静止成像。
图8为等效超高音速运动一维距离像。
图9为超高音速目标静止一维距离成像。
具体实施方式
如图1所示,本发明超高音速目标成像方法,包括如下步骤:
(10)亚音速目标数据处理:将实测亚音速目标总回波数据处理为各帧回波数据点集并保存到各帧的二维采样矩阵中;
(20)亚音速目标运动一维距离成像:提取两帧回波数据分别到二维矩阵中,通过IFFT变换和图像拼接,得到两帧运动距离像;
如图2所示,所述(20)亚音速目标回波信号处理步骤包括:
(21)逆傅里叶变换:将两帧的回波矩阵S(m,n)每一行N个数据做IFFT计算,从而获得目标的距离像冗余矩阵;
(22)距离像拼接:按最大值拼接算法将两帧距离像二维冗余矩阵拼接为一维距离像,以此获得两帧目标的运动一维距离像,如图6所示。
(30)亚音速目标静止成像:设置回波距离门,根据两帧运动距离像及熵值测速法测出两帧下的目标速度值,并对目标回波进行速度补偿,再次通过IFFT变换和像拼接,得到亚音速目标静止一维距离距离像;
如图3所示,所述(30)亚音速目标静止成像步骤包括:
(32)求熵值:先设置不同的速度估计值,求出两帧目标的运动一维距离像的熵,根据如下公式求出两帧像的最小的熵,其速度估计值即为测出的目标运动速度值:
V=min(Hn(A,Vn))
式中,Hn(A,Vn)为以估计速度Vn得到的熵,A(i)为目标的一维距离像序列,N为脉冲个数,V为当熵最小时的估计速度值。
(40)等效超高音速运动一维距离像:设置超高音速的速度值,分别计算两帧静止距离像对应的超高音速目标相应位置及幅值大小,并将静止目标回波乘以高速等效因子,获得超高音速相位信息,根据高速目标位置、幅度和相位信息等效出两帧的超高音速目标运动一维距离像;
如图4所示,所述(40)等效超高速运动一维距离像步骤包括:
(44)计算出距离像:根据求出的位移信息,相位信息,以及衰减信息,推算出两帧的目标回波,并作IFFT变换和像拼接,得到对应的超高音速目标的运动一维距离像,如图8所示。
(50)超高音速目标静止一维距离成像:用熵值测速法测出两帧运动距离像下的超高音速的速度大小,并对超高音速目标回波作速度补偿,再通过IFFT变换和像拼接算法得到两帧的超高音速目标静止一维距离像。
如图5所示,所述(50)超高速目标静止一维距离成像步骤包括:
(51)求熵值:先设置不同的超高音速速度估计值,求出两帧目标的运动一维距离像的熵,根据如下公式求出两帧像的最小的熵,其速度估计值即为测出的目标超高音速运动速度值:
V=min(Hn(A,Vn))
式中,Hn(A,Vn)为以估计速度Vn得到的熵,A(i)为目标的一维距离像序列,N为脉冲个数,V为当熵最小时的估计超高音速速度值;
使用本发明方案,用于实现不同目标在不同环境下以超高音速运动的成像情况,对进行超高音速实验成像测试提供了思路,解决了实验实施困难的问题,有效提升了实验效率。
Claims (5)
1.一种超高音速目标成像方法,其特征在于,包括如下步骤:
(10)亚音速目标数据处理:将实测亚音速目标总回波数据处理为各帧回波数据点集并保存到各帧的二维采样矩阵中;
(20)亚音速目标运动一维距离成像:提取两帧回波数据分别到二维矩阵中,通过LFFT变换和图像拼接,得到两帧运动距离像;
(30)亚音速目标静止成像:设置回波距离门,根据两帧运动距离像及熵值测速法测出两帧下的目标速度值,并对目标回波进行速度补偿,再次通过IFFT变换和像拼接,得到亚音速目标静止一维距离距离像;
(40)等效超高音速运动一维距离像:设置超高音速的速度值,分别计算两帧静止距离像对应的超高音速目标相应位置及幅值大小,并将静止目标回波乘以高速等效因子,获得超高音速相位信息,根据高速目标位置、幅度和相位信息等效出两帧的超高音速目标运动一维距离像;
(50)超高音速目标静止一维距离成像:用熵值测速法测出两帧运动距离像下的超高音速的速度大小,并对超高音速目标回波作速度补偿,再通过IFFT变换和像拼接算法得到两帧的超高音速目标静止一维距离像。
2.根据权利要求1所述的超高音速目标成像方法,其特征在于,所述(20)亚音速目标回波信号处理步骤包括:
(21)逆傅里叶变换:将两帧的回波矩阵S(m,n)每一行N个数据做IFFT计算,从而获得目标的距离像冗余矩阵;
(22)距离像拼接:按如下的最大值拼接算法将两帧距离像二维冗余矩阵拼接为一维距离像,以此获得两帧目标的运动一维距离像
Qs=fix(mLTs/Tu)
ne=N×mod(mLTs/Tu,1)
sndr=QsN+mod(ne-W+N,N)
Wp=NTP/Tu
m=mL,mL+1,...,mH-1,mH
q=fix[(mTs-dt)/Tu]+1
ne=N×mod(mTs/Tu,1)
ns=mod(ns-W+N,N)
qs=q×N
式中,Qs为第mL个采样点mLTs所在不模糊时间长度,N为脉冲串长度,Ts为采样时间,dt为时间分辨率,ne为该采样点对应的结束位置,ns为开始位置,sndr为开始位置,W为有效距离像宽度,Tu为不模糊距离窗,Z为去冗余后提取出的距离像存放数组,S(m,n)为原距离像冗余矩阵。
4.根据权利要求3所述的超高音速目标成像方法,其特征在于,所述(40)等效超高速运动一维距离像步骤包括:
(44)计算出距离像:根据求出的位移信息,相位信息,以及衰减信息,推算出两帧的目标回波,并作IFFT变换和像拼接,得到对应的超高音速目标的运动一维距离像。
5.根据权利要求4所述的超高音速目标成像方法,其特征在于,所述(50)超高速目标静止一维距离成像步骤包括:
(51)求熵值:先设置不同的超高音速速度估计值,求出两帧目标的运动一维距离像的熵,根据如下公式求出两帧像的最小的熵,其速度估计值即为测出的目标超高音速运动速度值:
V=min(Hn(A,Vn))
式中,Hn(A,Vn)为以估计速度Vn得到的熵,A(i)为目标的一维距离像序列,N为脉冲个数,V为当熵最小时的估计超高音速速度值;
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011274953.9A CN112630778B (zh) | 2020-11-16 | 2020-11-16 | 一种超高音速目标成像方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011274953.9A CN112630778B (zh) | 2020-11-16 | 2020-11-16 | 一种超高音速目标成像方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112630778A true CN112630778A (zh) | 2021-04-09 |
CN112630778B CN112630778B (zh) | 2022-08-16 |
Family
ID=75304375
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011274953.9A Active CN112630778B (zh) | 2020-11-16 | 2020-11-16 | 一种超高音速目标成像方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112630778B (zh) |
Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590812A (zh) * | 2012-01-13 | 2012-07-18 | 西安电子科技大学 | 基于调频连续波的sar实时成像方法 |
CN105068072A (zh) * | 2015-07-22 | 2015-11-18 | 北京环境特性研究所 | 一种运动目标的一维距离像的速度补偿方法 |
CN106872974A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 基于高超声速平台双通道雷达的高精度运动目标成像方法 |
CN111007503A (zh) * | 2019-09-03 | 2020-04-14 | 西安电子科技大学 | 基于频率谱精确定位的运动目标聚焦和定位方法及*** |
-
2020
- 2020-11-16 CN CN202011274953.9A patent/CN112630778B/zh active Active
Patent Citations (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102590812A (zh) * | 2012-01-13 | 2012-07-18 | 西安电子科技大学 | 基于调频连续波的sar实时成像方法 |
CN105068072A (zh) * | 2015-07-22 | 2015-11-18 | 北京环境特性研究所 | 一种运动目标的一维距离像的速度补偿方法 |
CN106872974A (zh) * | 2017-01-23 | 2017-06-20 | 西安电子科技大学 | 基于高超声速平台双通道雷达的高精度运动目标成像方法 |
CN111007503A (zh) * | 2019-09-03 | 2020-04-14 | 西安电子科技大学 | 基于频率谱精确定位的运动目标聚焦和定位方法及*** |
Non-Patent Citations (1)
Title |
---|
吴晓芳 等: "《频率步进雷达目标运动参数估计的新方法》", 《信号处理》 * |
Also Published As
Publication number | Publication date |
---|---|
CN112630778B (zh) | 2022-08-16 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Rao et al. | Adaptive sparse recovery by parametric weighted L $ _ {1} $ minimization for ISAR imaging of uniformly rotating targets | |
CN107229048B (zh) | 一种高分宽幅sar动目标速度估计与成像方法 | |
CN111142105B (zh) | 复杂运动目标isar成像方法 | |
CN108008385B (zh) | 基于稀疏贝叶斯学习的干扰环境isar高分辨成像方法 | |
CN104950306B (zh) | 一种海杂波背景下前视海面目标角超分辨成像方法 | |
CN107271993B (zh) | 一种基于最大后验的扫描雷达角超分辨成像方法 | |
CN103091674B9 (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
CN110244303B (zh) | 基于sbl-admm的稀疏孔径isar成像方法 | |
CN111505639B (zh) | 基于变重频采样模式的合成孔径雷达宽幅稀疏成像方法 | |
CN109507666B (zh) | 基于离网变分贝叶斯算法的isar稀疏频带成像方法 | |
CN108226928B (zh) | 基于期望传播算法的逆合成孔径雷达成像方法 | |
Sun et al. | High-resolution ISAR imaging of maneuvering targets based on sparse reconstruction | |
CN111580104B (zh) | 基于参数化字典的机动目标高分辨isar成像方法 | |
CN104950305A (zh) | 一种基于稀疏约束的实波束扫描雷达角超分辨成像方法 | |
CN113030972B (zh) | 基于快速稀疏贝叶斯学习的机动目标isar成像方法 | |
CN108646247B (zh) | 基于伽马过程线性回归的逆合成孔径雷达成像方法 | |
CN105137425A (zh) | 基于卷积反演原理的扫描雷达前视角超分辨方法 | |
CN107843875A (zh) | 基于奇异值分解降噪的贝叶斯压缩感知雷达数据融合方法 | |
CN114114260A (zh) | 一种基于深度学习的机动目标isar自聚焦成像方法 | |
CN111681272B (zh) | 一种基于奇异性功率谱的sar图像处理方法 | |
CN114252878B (zh) | 一种基于逆合成孔径雷达对运动目标进行成像及横向定标的方法 | |
CN113608217B (zh) | 基于增强矩阵补全的isar稀疏成像方法 | |
CN106680809B (zh) | 一种穿墙雷达自聚焦压缩感知成像方法 | |
Wu et al. | SAR imaging from azimuth missing raw data via sparsity adaptive StOMP | |
CN112630778B (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 |