CN106842197A - 一种机载isar舰船成像中心成像时刻选择方法 - Google Patents
一种机载isar舰船成像中心成像时刻选择方法 Download PDFInfo
- Publication number
- CN106842197A CN106842197A CN201510882973.7A CN201510882973A CN106842197A CN 106842197 A CN106842197 A CN 106842197A CN 201510882973 A CN201510882973 A CN 201510882973A CN 106842197 A CN106842197 A CN 106842197A
- Authority
- CN
- China
- Prior art keywords
- isar
- extreme point
- center
- imagings
- data
- 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
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
- G01S13/9064—Inverse SAR [ISAR]
-
- 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/9023—SAR image post-processing techniques combined with interferometric 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
本发明提供了一种机载ISAR舰船成像中心成像时刻选择方法,其特征在于,包括以下步骤:步骤1,对待处理的机载ISAR舰船成像的原始回波数据进行解包以及距离向脉冲压缩,再对整段ISAR数据进行距离对准和相位补偿处理,并在距离向对数据进行截取;步骤2,对距离对准和相位补偿处理过后的数据段进行短时快速傅里叶变换处理,数据段以可成图像的规则划分为若干子段,通过式(1)和式(2)对于每个子段的图像分别进行多普勒展宽频率fds和多普勒中心频率fdc估计,再分别进行插值和拟合,从而得到两条估计曲线Φds和Φdc。
Description
技术领域
本发明涉及一种结合多普勒展宽及多普勒中心估计的机载ISAR舰船成像最优时刻选择方法,属于ISAR舰船成像中最优成像时刻选择的技术领域。
背景技术
非合作舰船目标的成像一直是雷达成像领域关注的一个方面。近年来我国海军装备发展迫切,舰载和机载雷达系统急需具备舰船高分辨率雷达成像能力,舰船类目标的ISAR成像成为研究热点。ISAR成像中,舰船成像的难度在于其运动形式的多样化,对于机载ISAR舰船成像而言,不仅包括舰船和飞机的相对运动,还包括由于海面波浪起伏造成的横摇、纵摇以及偏航这三种摇摆运动。若对舰船的摇摆运动进行合理的利用,可以在较短的成像时间内获得高质量的舰船图像,但实际情况下舰船的各维摇摆近似于正弦函数特性,且摇摆的幅度和周期各不相同,使得舰船在三维空间中做非平稳运动,增加了舰船成像的复杂度。而载机与雷达之间的相对运动能够改变合成转动矢量的方向与大小,从而引起成像面的变化,更加增加了成像难度。
最优成像时间选择法首先选择最佳成像数据段,然后通过传统距离多普勒(RD)方法,即方位向通过快速傅里叶变换(FFT)即可成像。最优成像时间选择法比较典型的有:
1、彭石宝等人在“基于相位线性度的ISAR非平稳目标成像时间选择新算法[J].电子与信息学报,2010,32(12):2637-2640”提出基于相位线性度估计的最优成像时间选择方法。该方法基于目标的平均距离像提取两个特显点距离单元信号来估计转动相位,选择转动相位线性度最大的时间段作为成像时间,但运算较为复杂,存在一定的局限性。
2、程玉平等人在“一种改进的多普勒中心频率估计方法[J].西安电子科技大学学报,1999,26(1):44-48”提出基于多普勒中心估计的方法。但该方法在估计过程中可能会存在较大的估计误差,影响估计精度。
不同于上述几种方法,本发明通过一种数学概率密度模型来估计多普勒中心频率和多普勒展宽频率,不需要设置门限,并且结合两个估计曲线再次对有效转动矢量的变化趋势进行估计,两种方法的结合能够有效的地减少估计误差,具有良好的鲁棒性,且其过程简单、计算量小。根据本发明人的了解,关于结合多普勒展宽及多普勒中心估计的ISAR舰船成像最优时刻选择方法,目前尚没有公开发表的技术文献。
发明创造的目的
本发明所要解决的技术问题是针对上述背景技术的不足,提供了一种结合多普勒展宽及多普勒中心估计的机载ISAR舰船成像最优时刻选择方法。
技术方案
本发明为实现上述目的,本发明的技术方案包括如下步骤:
一种机载ISAR舰船成像中心成像时刻选择方法,其特征在于,包括以下步骤:
步骤1,对待处理的机载ISAR舰船成像的原始回波数据进行解包以及距离向脉冲压缩,再对整段ISAR数据进行距离对准和相位补偿处理,并在距离向对数据进行截取;
步骤2,对距离对准和相位补偿处理过后的数据段进行短时快速傅里叶变换处理,数据段以可成图像的规则划分为若干子段,通过式(1)和式(2)对于每个子段的图像分别进行多普勒展宽频率fds和多普勒中心频率fdc估计,再分别进行插值和拟合,从而得到两条估计曲线Φds和Φdc,
其中fds为多普勒展宽,fdc表示多普勒中心频率,·为点乘运算,PRF为脉冲重复频率,Na为方位向脉冲数,pn表示对Fn归一化结果,Fn表示一维方位向包络(下标n表示方位向坐标)。
步骤3,对多普勒中心估计曲线Φdc取绝对值得|Φdc|,对|Φdc|和Φds进行归一化处理,分别对|Φdc|、Φds极值点位置和零点位置排序,以|Φdc|极值点位置和零点位置为准,若Φds极值点位置与|Φdc|极值点位置相差±500以内,则保留该Φds极值点;若Φds零点位置与|Φdc|零点位置相差±500以内,则保留该Φds零点;若Φds极值点位置等于|Φdc|零点位置,则舍弃该Φds极值点,最终将所有保留的Φds极值点和零点进行插值拟合,获得估计曲线即为新的估计有效转动矢量变化趋势的曲线Φnew。
步骤4,对步骤3处理得到的曲线Φnew进行差分处理,从而选取曲线Φnew的极值点对应的方位向脉冲点数作为ISAR舰船成像的中心成像时刻。
步骤5,以步骤4中选择的中心成像时刻,选取合适的脉冲数进行积累成像,从而获得ISAR舰船图像。
本发明方法简单实用,主要适用于ISAR舰船成像的应用领域。与现有技术相比具有如下优点:
(1)与传统的多普勒中心及多普勒展宽估计方法不同,本发明采用基于数学概率密度模型的估计方法,不需要设置门限,实现简单。
(2)与传统的ISAR舰船最优成像时刻选择方法相比,本发明将多普勒中心估计及多普勒展宽估计方法进行有效结合,从而减小有效转动矢量的估计误差,提高估计精度。
附图说明
图1是本发明时刻选择方法的流程图。
具体实施方式
参阅附图1,本发明的一种机载ISAR舰船成像中心成像时刻选择方法,其特征在于,包括以下步骤:
步骤1,对待处理的机载ISAR舰船成像的原始回波数据进行解包以及距离向脉冲压缩,再对整段ISAR数据进行距离对准和相位补偿处理,并在距离向对数据进行截取;
步骤2,对距离对准和相位补偿处理过后的数据段进行短时快速傅里叶变换处理,数据段以可成图像的规则划分为若干子段,通过式(1)和式(2)对于每个子段的图像分别进行多普勒展宽频率fds和多普勒中心频率fdc估计,再分别进行插值和拟合,从而得到两条估计曲线Φds和Φdc,
其中fds为多普勒展宽,fdc表示多普勒中心频率,·为点乘运算,PRF为脉冲重复频率,Na为方位向脉冲数,pn表示对Fn归一化结果,Fn表示一维方位向包络(下标n表示方位向坐标)。
步骤3,对多普勒中心估计曲线Φdc取绝对值得|Φdc|,对|Φdc|和Φds进行归一化处理,分别对|Φdc|、Φds极值点位置和零点位置排序,以|Φdc|极值点位置和零点位置为准,若Φds极值点位置与|Φdc|极值点位置相差±500以内,则保留该Φds极值点;若Φds零点位置与|Φdc|零点位置相差±500以内,则保留该Φds零点;若Φds极值点位置等于|Φdc|零点位置,则舍弃该Φds极值点,最终将所有保留的Φds极值点和零点进行插值拟合,获得估计曲线即为新的估计有效转动矢量变化趋势的曲线Φnew。
步骤4,对步骤3处理得到的曲线Φnew进行差分处理,从而选取曲线Φnew的极值点对应的方位向脉冲点数作为ISAR舰船成像的中心成像时刻。
步骤5,以步骤4中选择的中心成像时刻,选取合适的脉冲数进行积累成像,从而获得ISAR舰船图像。
Claims (1)
1.一种机载ISAR舰船成像中心成像时刻选择方法,其特征在于,包括以下步骤:
步骤1,对待处理的机载ISAR舰船成像的原始回波数据进行解包以及距离向脉冲压缩,再对整段ISAR数据进行距离对准和相位补偿处理,并在距离向对数据进行截取;
步骤2,对距离对准和相位补偿处理过后的数据段进行短时快速傅里叶变换处理,数据段以可成图像的规则划分为若干子段,通过式(1)和式(2)对于每个子段的图像分别进行多普勒展宽频率fds和多普勒中心频率fdc估计,再分别进行插值和拟合,从而得到两条估计曲线Φds和Φdc,
其中fds为多普勒展宽,fdc表示多普勒中心频率,·为点乘运算,PRF为脉冲重复频率,Na为方位向脉冲数,pn表示对Fn归一化结果,Fn表示一维方位向包络(下标n表示方位向坐标)。
步骤3,对多普勒中心估计曲线Φdc取绝对值得|Φdc|,对|Φdc|和Φds进行归一化处理,分别对|Φdc|、Φds极值点位置和零点位置排序,以|Φdc|极值点位置和零点位置为准,若Φds极值点位置与|Φdc|极值点位置相差±500以内,则保留该Φds极值点;若Φds零点位置与|Φdc|零点位置相差±500以内,则保留该Φds零点;若Φds极值点位置等于|Φdc|零点位置,则舍弃该Φds极值点,最终将所有保留的Φds极值点和零点进行插值拟合,获得估计曲线即为新的估计有效转动矢量变化趋势的曲线Φnew。
步骤4,对步骤3处理得到的曲线Φnew进行差分处理,从而选取曲线Φnew的极值点对应的方位向脉冲点数作为ISAR舰船成像的中心成像时刻。
步骤5,以步骤4中选择的中心成像时刻,选取合适的脉冲数进行积累成像,从而获得ISAR舰船图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510882973.7A CN106842197B (zh) | 2015-12-03 | 2015-12-03 | 一种机载isar舰船成像中心成像时刻选择方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201510882973.7A CN106842197B (zh) | 2015-12-03 | 2015-12-03 | 一种机载isar舰船成像中心成像时刻选择方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106842197A true CN106842197A (zh) | 2017-06-13 |
CN106842197B CN106842197B (zh) | 2019-12-17 |
Family
ID=59149932
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201510882973.7A Active CN106842197B (zh) | 2015-12-03 | 2015-12-03 | 一种机载isar舰船成像中心成像时刻选择方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106842197B (zh) |
Cited By (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169746A (zh) * | 2017-12-21 | 2018-06-15 | 南京理工大学 | 线性调频脉冲半主动雷达导引头信号处理方法 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0544533B1 (en) * | 1991-11-26 | 1996-10-02 | Texas Instruments Incorporated | Improved ISAR imaging radar system |
CN101846741A (zh) * | 2010-05-07 | 2010-09-29 | 北京航空航天大学 | 一种逆合成孔径雷达成像数据段选择方法 |
CN103293527A (zh) * | 2013-05-15 | 2013-09-11 | 西安电子科技大学 | 基于置信框架的自适应isar成像方法 |
CN104122550A (zh) * | 2014-07-08 | 2014-10-29 | 上海无线电设备研究所 | 高分辨率isar实时成像系统 |
CN104931966A (zh) * | 2015-06-12 | 2015-09-23 | 北京航空航天大学 | 一种基于dcs算法的星载视频sar成像处理方法 |
-
2015
- 2015-12-03 CN CN201510882973.7A patent/CN106842197B/zh active Active
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0544533B1 (en) * | 1991-11-26 | 1996-10-02 | Texas Instruments Incorporated | Improved ISAR imaging radar system |
CN101846741A (zh) * | 2010-05-07 | 2010-09-29 | 北京航空航天大学 | 一种逆合成孔径雷达成像数据段选择方法 |
CN103293527A (zh) * | 2013-05-15 | 2013-09-11 | 西安电子科技大学 | 基于置信框架的自适应isar成像方法 |
CN104122550A (zh) * | 2014-07-08 | 2014-10-29 | 上海无线电设备研究所 | 高分辨率isar实时成像系统 |
CN104931966A (zh) * | 2015-06-12 | 2015-09-23 | 北京航空航天大学 | 一种基于dcs算法的星载视频sar成像处理方法 |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN108169746A (zh) * | 2017-12-21 | 2018-06-15 | 南京理工大学 | 线性调频脉冲半主动雷达导引头信号处理方法 |
CN108169746B (zh) * | 2017-12-21 | 2021-09-21 | 南京理工大学 | 线性调频脉冲半主动雷达导引头信号处理方法 |
Also Published As
Publication number | Publication date |
---|---|
CN106842197B (zh) | 2019-12-17 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN104574324B (zh) | 对地面建筑群遥感图像频谱特性约束的去噪方法和系统 | |
CN102540189A (zh) | 基于复数后向投影的自旋目标三维成像方法 | |
CN103091674B (zh) | 基于hrrp序列的空间目标高分辨成像方法 | |
CN102426360A (zh) | 空中微动旋转目标的二维isar成像方法 | |
CN103675759A (zh) | 一种改进的分数阶傅里叶变换机动弱目标检测方法 | |
CN107632304A (zh) | 一种大斜视条件下的合成孔径声呐成像处理方法 | |
CN107831473B (zh) | 基于高斯过程回归的距离-瞬时多普勒图像序列降噪方法 | |
CN106772375A (zh) | 基于参数估计的压缩感知成像方法 | |
CN102866260A (zh) | 非接触式河流表面流场成像量测方法 | |
CN107942314A (zh) | 基于lasso特征提取的多普勒穿墙雷达定位方法 | |
CN112764116B (zh) | 一种稀疏阵列稀疏频点平面扫描体制快速成像方法 | |
CN109308713A (zh) | 一种基于前视声纳的改进核相关滤波水下目标跟踪方法 | |
CN105068062A (zh) | 基于稀疏散射中心提取的距离像数据外推方法 | |
CN103630899B (zh) | 地面运动目标高分辨雷达压缩感知成像的方法 | |
CN111896957A (zh) | 基于小波变换及压缩感知的舰船目标前视三维成像方法 | |
CN114167423A (zh) | 基于深度回归网络的雷达海浪参数测量方法 | |
CN105005048B (zh) | 基于显著图的Laplacian协同压缩雷达成像方法 | |
CN102944876A (zh) | 一种高效的宽带单脉冲雷达三维成像方法 | |
CN112882016A (zh) | 基于改进的稳健自适应波束形成的多人生命体征检测方法 | |
CN106842197A (zh) | 一种机载isar舰船成像中心成像时刻选择方法 | |
CN103487807A (zh) | 一种isar最优成像时间选择方法 | |
CN105116410B (zh) | 基于线性模型匹配的干涉相位图自适应滤波算法 | |
CN106908782B (zh) | 基于水面状态连续成像系统的波浪传播方向的提取方法 | |
CN109444843A (zh) | 一种单频雷达成像方法和装置 | |
CN113096171B (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 |