CN108427115B - 合成孔径雷达对运动目标参数的快速估计方法 - Google Patents
合成孔径雷达对运动目标参数的快速估计方法 Download PDFInfo
- Publication number
- CN108427115B CN108427115B CN201810083374.2A CN201810083374A CN108427115B CN 108427115 B CN108427115 B CN 108427115B CN 201810083374 A CN201810083374 A CN 201810083374A CN 108427115 B CN108427115 B CN 108427115B
- Authority
- CN
- China
- Prior art keywords
- signal
- moving target
- distance
- azimuth
- time
- 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
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
- 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
- G01S13/904—SAR modes
-
- 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/9047—Doppler beam sharpening mode
Abstract
本发明提供了一种合成孔径雷达对运动目标参数的快速估计方法,属于雷达技术领域。本发明通过将脉冲压缩后的信号进行二阶Keystone变换操作,去除了距离频率与方位时间的二阶耦合。其次构造相位补偿信号,并通过乘法器将相位补偿信号与原信号相乘,降低了距离频率与方位时间的耦合阶数,即去除了一阶耦合。最后对信号进行方位向傅里叶变换和距离向傅里叶逆变换得到最终运动目标参数域的积累结果,从而可以估计运动目标的参数。本发明与现有技术相比,可以同时估计多普勒质心和多普勒调频率,并且无需距离徙动校正后对一维信号的复杂操作,仅存在距离徙动校正后的傅里叶变换操作,形式简单,运算效率高。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种合成孔径雷达对运动目标参数的快速估计方法。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种全天时、全天候、高分辨率的微波遥感成像雷达,被广泛应用于地形测绘、灾害监测、植被分析、海洋观测、城市规划等领域。
SAR成像对静止目标成像具有操作简单、精度高、多普勒参数易计算等优点。但是,对动目标成像结果中会出现运动目标散焦、移位等现象,主要由于运动参数未知,导致运动目标无法进行高分辨聚焦。
传统的运动目标参数估计方法有Hough变换、Radon变换等,其中,Radon变换可对多普勒质心进行估计,但是该算法仅能估计多普勒质心,不能实现对多普勒调频率的估计。并且该算法为非相干积累方法,抗噪性能差,运算效率低。
现有的技术中,还有一种相干积累的估计参数算法,但是该算法是对一维信号进行处理,二维回波需要进行校正后再进行复杂参数估计过程,操作复杂、运算效率较低。
发明内容
为了解决上述不能同时估计多普勒质心和多普勒调频率,并且估计过程复杂的问题,本发明提供一种合成孔径雷达对运动目标参数的快速估计方法。
本发明的一种合成孔径雷达对运动目标参数的快速估计方法,所述合成孔径雷达装载于飞行平台上,所述方法包括以下步骤:
步骤1,建立飞行平台和运动目标的几何模型并初始化参数;
步骤2,获取所述运动目标的回波信号,并对所述回波信号进行距离向匹配滤波操作;
步骤3,对匹配滤波后的信号进行二阶Keystone变换;
步骤4,将经过二阶Keystone变换后的信号通过方位向定量延时系统,构造相位补偿信号;
步骤5,对相位补偿信号进行共轭操作,并利用乘法器将共轭后的相位补偿信号和经过二阶Keystone变换后的信号相乘;
步骤6,对相乘后的信号进行距离向傅里叶逆变换和方位向傅里叶变换,得到参数估计结果。
进一步地,所述步骤1的具体流程包括:
以O为坐标原点建立直角坐标系,所述飞行平台的初始位置坐标设为(0,0,H0),地面运动目标的初始位置设为(X0,0,0),其中,H0为所述飞行平台的Z轴坐标,X0为所述运动目标的X轴坐标,所述飞行平台的飞行速度为V,所述运动目标的距离向速度为Vx,所述运动目标的距离向加速度为ax,所述运动目标的方位向速度为Vy,所述运动目标的方位向加速度为ay,得到所述飞行平台到所述运动目标的最短斜距R0,
得到所述运动目标到所述飞行平台的距离历史RM(t)为
其中,t为方位向的时间变量,
将所述距离历史RM(t)在中心时刻t=0沿泰勒级数展开,得到
进一步地,所述步骤2的具体流程包括:
获取下变频后的所述运动目标回波,所述运动目标的回波信号S(τ,t)在距离时域、方位时域的表达式为
其中,rect(·)为距离向和方位向矩形窗,Ta为合成孔径时间,τ为距离向的时间变量,Δτ=2RM(t)/c为线性调频信号的时延,Tp为发射信号时宽,为虚数单位,kτ为距离向调频率,c为光速,fc为载波频率;
对所述回波信号进行距离向匹配滤波处理,在距离频域、方位时域的信号记为
其中,fτ为距离频率变量,Bτ为发射信号带宽。
进一步地,所述步骤3的具体流程包括:
对匹配滤波后的信号进行二阶Keystone变换,二阶Keystone变换表示为
其中,tm为新的方位向时间变量,
将二阶Keystone变换后的信号相位沿方位时间tm=0进行泰勒展开,得到的信号表达式为
进一步地,所述步骤4的具体流程包括:
将信号通过方位向定量延时系统,构造相位补偿信号
其中,t0为常量时延。
进一步地,所述步骤5的具体流程包括:
对相位补偿信号进行共轭操作,并利用乘法器将共轭后的信号和二阶Keystone处理后的信号相乘,得到所述运动目标的参数估计结果
进一步地,所述步骤6的具体流程包括:
对相乘后的信号进行距离向傅里叶逆变换和方位向傅里叶变换,得到所述运动目标的参数估计结果
本发明的有益效果:首先将脉冲压缩后的信号进行二阶Keystone变换操作,去除了距离频率与方位时间的二阶耦合。其次构造相位补偿信号,并通过乘法器将相位补偿信号与原信号相乘,降低了距离频率与方位时间的耦合阶数,即去除了一阶耦合。最后对信号进行方位向傅里叶变换和距离向傅里叶逆变换得到最终运动目标参数域的积累结果,从而可以估计运动目标的参数。本发明和传统的运动目标参数方法相比,可以同时估计多普勒质心和多普勒调频率,并且无需距离徙动校正后对一维信号的复杂操作,仅存在距离徙动校正后的傅里叶变换操作,形式简单,运算效率高。本发明可以应用于机载SAR运动目标参数估计与成像领域。
附图说明
图1为本发明的流程图。
图2为本发明步骤1中建立的几何模型。
图3为经步骤2中距离向匹配滤波后的回波信号图像。
图4为经步骤3中二阶Keystone变换处理后的信号图像。
图5为经步骤5后的信号在距离时域-方位时域的图像。
图6为实施例中运动目标在参数域的图像。
具体实施方式
下面结合附图对本发明的实施例做进一步的说明。
本发明实施例以仿真实验的方法展现,所有步骤和结果都可以在Mattlab2017上验证正确。
如图1所示,本发明的合成孔径雷达对运动目标参数的快速估计方法,包括以下步骤:
步骤1,建立飞行平台和运动目标的几何模型并初始化参数。
本发明中,合成孔径雷达(SAR)装载于飞行平台上,对地面上的运动目标进行参数估计。数据见表一。建立如图2所示的飞行平台和运动目标的直角坐标系,以O为坐标原点建立直角坐标系,飞行平台的初始位置坐标设为(0m,0m,2000m),地面运动目标的初始位置设为(500m,0m,0m),飞行平台的飞行速度V为150m/s,运动目标在距离向的速度Vx为32m/s,运动目标在距离向的加速度ax为1m/s2,在方位向的速度Vy为8m/s,方位向的加速度ay为2m/s2,得到飞行平台到运动目标的最短斜距R0为2061m。
根据上述参数数据,得到运动目标到飞行平台的距离历史RM(t)为
其中,t为方位向的时间变量,
将距离历史RM(t)在中心时刻t=0沿泰勒级数展开,得到
步骤2,获取运动目标的回波信号,并对回波信号进行距离向匹配滤波操作。
获取经过脉冲压缩下变频的运动目标的回波,回波信号S(τ,t)在距离时域、方位时域的表达式为
其中,rect(·)为距离向和方位向矩形窗,Ta为合成孔径时间,τ为距离向的时间变量,Δτ=2RM(t)/c为线性调频信号的时延,Tp为发射信号时宽,为虚数单位,kτ为距离向调频率,c为光速,fc为载波频率。
对回波信号经过距离向匹配滤波处理,在距离频域、方位时域的信号记为
其中,fτ为距离频率变量,Bτ为发射信号带宽。
获取的信号如图3所示,信号是一条弯曲倾斜的线,其中包含距离频率和方位时域的一阶和二阶耦合。
步骤3,对匹配滤波后的信号进行二阶Keystone变换。
对匹配滤波后的信号进行二阶Keystone变换,二阶Keystone变换表示为
其中,tm为新的方位向时间变量,
将二阶Keystone变换后的信号相位沿方位时间tm=0进行泰勒展开,得到的信号表达式为
将步骤2后的脉冲压缩后的信号进行二阶Keystone变换操作,得到如图4所示的信号图,信号轨迹变成一条斜线,即二阶耦合被去除。
步骤4,将经过二阶Keystone变换后的信号通过方位向定量延时系统,构造相位补偿信号。
将信号通过方位向定量延时系统,构造相位补偿信号
其中,t0为常量时延。
步骤5,对相位补偿信号进行共轭操作,并利用乘法器将共轭后的相位补偿信号和经过二阶Keystone变换后的信号相乘。
对相位补偿信号进行共轭操作,并利用乘法器将共轭后的信号和二阶Keystone处理后的信号相乘,得到所述运动目标的参数估计结果
步骤4和5,构造相位补偿信号,并通过乘法器将相位补偿信号与原信号相乘,得到如图5所示信号图,图5经补偿方程补偿后的信号,信号的一阶耦合量被去除,降低了距离频率与方位时间的耦合阶数。
步骤6,对相乘后的信号进行距离向傅里叶逆变换和方位向傅里叶变换,得到参数估计结果。
对相乘后的信号进行距离向傅里叶逆变换和方位向傅里叶变换,得到所述运动目标的参数估计结果
图6为最终的参数域聚焦结果,可以看出该方法能够实现信号能量在参数域上聚焦并最终实现对运动目标的参数估计。
数据见表一。
表(一)
本实施例仅以一组数据为例,对运动目标参数进行了快速估计,本发明适用于机载SAR动目标参数估计与成像领域。
本领域的普通技术人员将会意识到,这里所述的实施例是为了帮助读者理解本发明的原理,应被理解为本发明的保护范围并不局限于这样的特别陈述和实施例。本领域的普通技术人员可以根据本发明公开的这些技术启示做出各种不脱离本发明实质的其它各种具体变形和组合,这些变形和组合仍然在本发明的保护范围内。
Claims (3)
1.一种合成孔径雷达对运动目标参数的快速估计方法,其特征在于,所述合成孔径雷达装载于飞行平台上,所述方法包括以下步骤:
步骤1,建立飞行平台和运动目标的几何模型并初始化参数;具体流程包括:
以O为坐标原点建立直角坐标系,所述飞行平台的初始位置坐标设为(0,0,H0),地面运动目标的初始位置设为(X0,0,0),其中,H0为所述飞行平台的Z轴坐标,X0为所述运动目标的X轴坐标,所述飞行平台的飞行速度为V,所述运动目标的距离向速度为Vx,所述运动目标的距离向加速度为ax,所述运动目标的方位向速度为Vy,所述运动目标的方位向加速度为ay,得到所述飞行平台到所述运动目标的最短斜距R0,
得到所述运动目标到所述飞行平台的距离历史RM(t)为
其中,t为方位向的时间变量,
将所述距离历史RM(t)在中心时刻t=0沿泰勒级数展开,得到
步骤2,获取所述运动目标的回波信号,并对所述回波信号进行距离向匹配滤波操作;具体流程包括:
获取下变频后的所述运动目标回波,所述运动目标的回波信号S(τ,t)在距离时域、方位时域的表达式为
其中,rect(·)为距离向和方位向矩形窗,Ta为合成孔径时间,τ为距离向的时间变量,Δτ=2RM(t)/c为线性调频信号的时延,Tp为发射信号时宽,为虚数单位,kτ为距离向调频率,c为光速,fc为载波频率;
对所述回波信号进行距离向匹配滤波处理,在距离频域、方位时域的信号记为
其中,fτ为距离频率变量,Bτ为发射信号带宽;
步骤3,对匹配滤波后的信号进行二阶Keystone变换;具体流程包括:
对匹配滤波后的信号进行二阶Keystone变换,二阶Keystone变换表示为
其中,tm为新的方位向时间变量,
将二阶Keystone变换后的信号相位沿方位时间tm=0进行泰勒展开,得到的信号表达式为
步骤4,将经过二阶Keystone变换后的信号通过方位向定量延时系统,构造相位补偿信号;具体流程包括:
将信号通过方位向定量延时系统,构造相位补偿信号
其中,t0为常量时延;
步骤5,对相位补偿信号进行共轭操作,并利用乘法器将共轭后的相位补偿信号和经过二阶Keystone变换后的信号相乘;
步骤6,对相乘后的信号进行距离向傅里叶逆变换和方位向傅里叶变换,得到参数估计结果。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810083374.2A CN108427115B (zh) | 2018-01-29 | 2018-01-29 | 合成孔径雷达对运动目标参数的快速估计方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201810083374.2A CN108427115B (zh) | 2018-01-29 | 2018-01-29 | 合成孔径雷达对运动目标参数的快速估计方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN108427115A CN108427115A (zh) | 2018-08-21 |
CN108427115B true CN108427115B (zh) | 2020-06-02 |
Family
ID=63156075
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201810083374.2A Active CN108427115B (zh) | 2018-01-29 | 2018-01-29 | 合成孔径雷达对运动目标参数的快速估计方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN108427115B (zh) |
Families Citing this family (8)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN109669183B (zh) * | 2017-12-27 | 2020-09-04 | 北京航空航天大学 | 一种基于Keystone和时频变换的地球同步轨道SAR运动目标成像处理装置 |
CN109633622B (zh) * | 2019-01-24 | 2022-11-29 | 中国人民解放军海军航空大学 | 基于时间反转二阶Keystone变换的雷达机动目标跨单元积累方法 |
CN109541568B (zh) * | 2019-01-24 | 2022-08-05 | 中国人民解放军海军航空大学 | 一种雷达机动目标跨距离和多普勒单元快速相参积累检测方法 |
CN111638516B (zh) * | 2019-03-01 | 2022-02-18 | 中国科学院电子学研究所 | 基于双频共轭处理技术的太赫兹频段sar运动补偿算法 |
CN110471060B (zh) * | 2019-07-08 | 2021-03-26 | 中国空间技术研究院 | 一种高效的雷达动目标调频率估计方法 |
CN112099021A (zh) * | 2020-08-14 | 2020-12-18 | 杭州电子科技大学 | 一种合成孔径声呐的点目标成像方法 |
CN112946642A (zh) * | 2021-01-27 | 2021-06-11 | 北京理工大学重庆创新中心 | 一种多通道uwb sar运动目标二维速度快速估计方法 |
CN113219458B (zh) * | 2021-05-26 | 2022-05-03 | 电子科技大学 | 双基地合成孔径雷达盲定位方法 |
Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778681A (zh) * | 2012-07-19 | 2012-11-14 | 电子科技大学 | 固定发射站双基地前视合成孔径雷达成像方法 |
CN102901966A (zh) * | 2012-10-22 | 2013-01-30 | 西安电子科技大学 | 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法 |
CN103983971A (zh) * | 2014-04-10 | 2014-08-13 | 西安电子科技大学 | 高效的星载sar的地面运动目标成像方法 |
CN104251990A (zh) * | 2014-09-15 | 2014-12-31 | 电子科技大学 | 合成孔径雷达自聚焦方法 |
CN104330795A (zh) * | 2014-10-20 | 2015-02-04 | 北京理工大学 | 一种基于Keystone变换的地基合成孔径雷达快速成像方法 |
CN105116395A (zh) * | 2015-07-02 | 2015-12-02 | 北京理工大学 | 一种天基双基地雷达的动目标长时间相参积累方法 |
Family Cites Families (1)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US8704703B2 (en) * | 2011-08-16 | 2014-04-22 | The Mitre Corporation | Pseudosynchronizer for an unsynchronized bistatic radar system |
-
2018
- 2018-01-29 CN CN201810083374.2A patent/CN108427115B/zh active Active
Patent Citations (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102778681A (zh) * | 2012-07-19 | 2012-11-14 | 电子科技大学 | 固定发射站双基地前视合成孔径雷达成像方法 |
CN102901966A (zh) * | 2012-10-22 | 2013-01-30 | 西安电子科技大学 | 基于解斜和楔石处理的动目标瞬时距离多普勒成像方法 |
CN103983971A (zh) * | 2014-04-10 | 2014-08-13 | 西安电子科技大学 | 高效的星载sar的地面运动目标成像方法 |
CN104251990A (zh) * | 2014-09-15 | 2014-12-31 | 电子科技大学 | 合成孔径雷达自聚焦方法 |
CN104330795A (zh) * | 2014-10-20 | 2015-02-04 | 北京理工大学 | 一种基于Keystone变换的地基合成孔径雷达快速成像方法 |
CN105116395A (zh) * | 2015-07-02 | 2015-12-02 | 北京理工大学 | 一种天基双基地雷达的动目标长时间相参积累方法 |
Non-Patent Citations (1)
Title |
---|
高速运动目标检测算法研究;于小龙;《中国优秀硕士学位论文全文数据库信息科技辑》;20140715;正文第46-48页 * |
Also Published As
Publication number | Publication date |
---|---|
CN108427115A (zh) | 2018-08-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN108051809B (zh) | 基于Radon变换的运动目标成像方法、装置及电子设备 | |
US7969345B2 (en) | Fast implementation of a maximum likelihood algorithm for the estimation of target motion parameters | |
CN103760558B (zh) | 一种太赫兹雷达isar成像方法 | |
CN108459321B (zh) | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 | |
CN102749621B (zh) | 一种双基地合成孔径雷达频域成像方法 | |
CN104749570B (zh) | 一种移不变机载双基合成孔径雷达目标定位方法 | |
CN111142105A (zh) | 复杂运动目标isar成像方法 | |
CN108226891B (zh) | 一种扫描雷达回波计算方法 | |
CN109270528B (zh) | 基于全解析距离模型的一站固定式双站sar成像方法 | |
CN110596701B (zh) | 基于二次椭圆模型的非平飞双站sar频域fenlcs成像方法 | |
CN103983974A (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN104833972A (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN114384520B (zh) | 一种机动平台对海面船只精细化雷达成像方法 | |
CN106054188A (zh) | 无人机合成孔径雷达成像的图像偏移自聚焦方法 | |
CN108008387B (zh) | 一种机载阵列天线下视三维成像方法 | |
CN114545411A (zh) | 一种基于工程实现的极坐标格式多模高分辨sar成像方法 | |
CN107153191B (zh) | 一种针对隐形飞机的双基isar成像检测方法 | |
CN111856466A (zh) | 一种高效的复杂运动目标isar平动补偿方法 | |
CN113671492A (zh) | 一种面向机动平台前视成像的samp重构方法 | |
CN106707278B (zh) | 一种基于稀疏表示的多普勒波束锐化成像方法及装置 | |
CN111060909A (zh) | 一种机载雷达斜前视超分辨成像方法 | |
CN109188436B (zh) | 适用于任意平台轨迹的高效双基sar回波生成方法 | |
CN110879391A (zh) | 基于电磁仿真和弹载回波仿真的雷达图像数据集制作方法 | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 |
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 |