CN107656274B - Sar动目标回波徙动校正与方位能量积累方法 - Google Patents
Sar动目标回波徙动校正与方位能量积累方法 Download PDFInfo
- Publication number
- CN107656274B CN107656274B CN201710851093.2A CN201710851093A CN107656274B CN 107656274 B CN107656274 B CN 107656274B CN 201710851093 A CN201710851093 A CN 201710851093A CN 107656274 B CN107656274 B CN 107656274B
- Authority
- CN
- China
- Prior art keywords
- moving target
- echo
- sar
- migration
- parameter
- 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
- G01S7/00—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00
- G01S7/02—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00
- G01S7/41—Details of systems according to groups G01S13/00, G01S15/00, G01S17/00 of systems according to group G01S13/00 using analysis of echo signal for target characterisation; Target signature; Target cross-section
Landscapes
- Engineering & Computer Science (AREA)
- Remote Sensing (AREA)
- Radar, Positioning & Navigation (AREA)
- Physics & Mathematics (AREA)
- Computer Networks & Wireless Communication (AREA)
- General Physics & Mathematics (AREA)
- Electromagnetism (AREA)
- Radar Systems Or Details Thereof (AREA)
- Geophysics And Detection Of Objects (AREA)
Abstract
本发明公开一种SAR动目标回波徙动校正与方位能量积累方法,本发明的方法利用改进的随机Hough变换,对具有噪声的距离脉冲压缩后回波曲线轨迹参数进行估计,得到曲线参考参数值,构造信号补偿方程,对原始回波信号进行相位补偿,从而解决了传统机载SAR不能同时具有一阶和二阶徙动量的运动目标回波校正问题,并且使运动目标回波能量在参数域上积累,更有利于对运动目标检测。本发明的方法与现有传统机载SAR雷达对地面动目标回波轨迹徙动校正方法相比,操作更简单,并且同时能校正一阶和二阶徙动。本发明的方法可以应用于机载SAR动目标检测和成像领域。
Description
技术领域
本发明属于雷达技术领域,具体涉及机载合成孔径雷达对地面动目标检测和成像技 术。
背景技术
合成孔径雷达(Synthetic Aperture Radar,SAR)是一种全天候、高分辨率微波遥感成像雷 达,它利用雷达天线和目标区域间的相对运动来获得空间的高分辨。目前在地形测绘,灾 区探测等领域得到广泛应用。
地面运动目标回波包含一阶和二阶徙动量,传统机载SAR成像算法对运动目标成像会 造成成像结果散焦、移位等情况。由于运动目标运动信息未知,传统成像算法对并不能对 动目标回波轨迹进行校正。
在文献:“Perry R P,Dipietro R C,Fante R.SAR imaging of moving targets[J].IEEE Transactions on Aerospace&Electronic Systems,1999,35(1):188-200”中,提出了一阶 Keystone变换,该方法在地面运动目标的运动信息未知的情况下进行一阶徙动量校正,但 是这种方法并不能对二阶徙动量进行校正,运动目标回波轨迹二阶徙动量依然会影响能量 积累和成像效果。
在文献:“Zhou F,Wu R,Xing M,et al.Approach for single channel SARground moving target imaging and motion parameter estimation.Iet Radar Sonar&Navigation,2007, 1(1):59-66”中,提出了二阶Keystone变换校正动目标二阶徙动量,该方法可以消除二阶徙 动量对能量积累和成像结果的影响,但是这种方法不能校正运动目标回波的一阶徙动量。
发明内容
本发明的目的是针对现有技术存在的缺陷,提出一种SAR动目标回波一阶和二阶徙动 校正和方位能量积累方法,以解决机载SAR不能同时对运动目标回波的一阶和二阶徙动量 校正的问题。
本发明的技术方案为:一种SAR动目标回波徙动校正与方位能量积累方法,具体包括 如下步骤:
步骤一:建立机载单基地正侧视SAR空间几何结构,并完成初始化;
在直角坐标系中,O为坐标原点,设(X0,0,0)为地面目标初始位置,其中,X0、0、0分别为地面运动目标的X轴,Y轴和Z轴坐标;
SAR飞行平台的位置坐标为(0,0,H0),其中,0、0和H0分别为SAR飞行平台的X轴,Y轴和Z轴坐标;
步骤二:根据几何模型获取运动目标到飞行平台的瞬时距离历史RM(t),可表示为:
其中,t为方位向时间变量;
步骤三:获取地面运动目标回波S(τ,t),并对回波进行下变频,具体公式为:
其中,rect(·)为距离向和方位向矩形窗,Tsar为合成孔径时间,τ为距离向时间变量, Δτ=2RM(t)/c为线性调频信号的时延,Tτ为发射信号时宽,kτ为距离向调频率,c为光速, fc为载波频率,为虚数单位。
步骤四:将地面运动目标回波进行距离向匹配滤波,滤波后回波信号相位可表示为:
其中,fτ为距离频率变量,θ为SAR天线入射角。
步骤五:构造地面运动目标回波信号初始补偿方程,可表示为:
y=c+bt+at2
其中,c为常数项系数,b一阶项系数,a二阶项系数。
步骤六:对脉冲压缩后的回波进行边缘检测,构造一个二维边缘点集合M,设定一个 幅度阈值N1,当边缘点幅度小于N1时将该点幅度置零;
步骤七:定义一个参数累加器A(a,b,c,n),用于存储参数值a,b,c和对应参数累加次数n,并将数组初始化;
步骤八:抽取二维边缘点集合中随机若干个不为零的数据点,计算出一组曲线参数(a1,b1,c1),参数结果存入参数累加器中,并对该组参数进行计数;
步骤九:在二维边缘点集合中去除抽取的若干个点;
步骤十:循环执行步骤八——步骤九,当边缘点集合中剩余的不为零数据点的个数不 足时停止循环;
步骤十一:设定一个参数累加次数阈值N2,累加次数n超过了阈值N2时,则将对应参数进行验证,用边缘点集合的所有数据点进行验证;
步骤十二:设定数据点数阈值N3,符合参数的数据点大于阈值N3时,则将估计出的参 数设定为原信号的徙动参数;
步骤十三:利用估计的徙动参数构造补偿方程,将补偿方程与原信号回波相乘,实现 运动目标回波轨迹徙动校正;
步骤十四:将估计出运动目标回波轨迹上数据点的值进行相干叠加,在曲线参数域可 以实现能量积累;
步骤十五:对目标回波信号经过徙动校正,获取校正后回波信号相位,徙动完全校正 的情况下,具体为:
从上式中可以看出,在机载单基地正侧视SAR空间几何构型下,距离向频率fτ和方位 项时间t的一阶和二阶耦合项同时被去除,剩余项仅为距离门、多普勒质心和多普勒调频率 项。因此可得,通过补偿运动目标回波信号,校正了SAR地面运动目标的回波一阶和二阶 徙动,去除了距离向频率fτ和方位项时间t的相互作用关系。
本发明的有益效果:本发明的方法利用改进的随机Hough变换,对具有噪声的距离脉 冲压缩后回波曲线轨迹参数进行估计,得到曲线参考参数值,构造信号补偿方程,对原始 回波信号进行相位补偿,从而解决了传统机载SAR不能同时具有一阶和二阶徙动量的运动 目标回波校正问题,并且使运动目标回波能量在参数域上积累,更有利于对运动目标检测。 本发明的方法与现有传统机载SAR雷达对地面动目标回波轨迹徙动校正方法相比,操作更 简单,并且同时能校正一阶和二阶徙动。本发明的方法可以应用于机载SAR动目标检测和 成像领域。
附图说明
图1是本发明提供的流程框图。
图2是本发明提供的机载SAR和地面运动目标空间几何结构。
图3是回波信噪比为0dB时回波轨迹校正前图像
图4是回波信噪比为0dB时回波轨迹校正后图像
图5是参数域能量积累图像
具体实施方式
本发明主要采用仿真实验的方法进行验证,所有步骤、结论都在Matlab2015上验证正 确。本发明的流程示意图如图1所示,具体包括如下步骤:
步骤一:建立机载单基地正侧视SAR空间几何结构,并完成初始化;
如图2所示,在直角坐标系中,设地面目标初始位置为(4km,0km,0km)。SAR平台的位 置坐标为(0km,0km,5km)。飞行平台的飞行速度为340m/s,地面运动目标的方位向速度和加 速度分别为10m/s和5m/s2,地面运动目标的距离向速度和加速度为17m/s和5m/s2,飞行平台到运动目标的最短斜距为6.4km,具体仿真参数如表1所示:
表1
参数 | 符号 | 数值 |
飞行平台位置坐标 | (0,0,H<sub>0</sub>) | 0km,0km,5km |
运动目标位置坐标 | (X<sub>0</sub>,0,0) | 4km,0km,0km |
飞行平台飞行速度 | V | 340m/s |
运动目标方位向速度 | V<sub>y</sub> | 10m/s |
运动目标方位向加速度 | a<sub>y</sub> | 5m/s<sup>2</sup> |
运动目标距离向速度 | V<sub>x</sub> | 17m/s |
运动目标距离向加速度 | a<sub>x</sub> | 5m/s<sup>2</sup> |
脉冲重复频率 | PRF | 800Hz |
发射信号时宽 | T<sub>r</sub> | 1.5μs |
发射信号带宽 | B<sub>r</sub> | 150MHz |
载波频率 | f<sub>c</sub> | 3GHz |
步骤二:根据几何模型获取运动目标到飞行平台的瞬时距离历史RM(t),可表示为:
步骤三:获取地面运动目标回波S(τ,t),并对回波进行下变频,具体公式为:
步骤四:将地面运动目标回波进行距离向匹配滤波,滤波后回波信号相位可表示为:
步骤五:构造地面运动目标回波信号初始补偿方程,可表示为:
y=c+bt+at2
步骤六:对脉冲压缩后的回波进行边缘检测,构造一个二维边缘点集合M,设定一个 幅度阈值N1,当边缘点幅度小于N1时将该点幅度置零。
步骤七:定义一个参数累加器A(a,b,c,n),用于存储参数值a,b,c和对应参数累加次数n,并将数组初始化。
步骤八:抽取二维边缘点集合中随机若干个不为零的数据点,计算出一组曲线参数 (a1,b1,c1),参数结果存入参数累加器中,并对该组参数进行计数,这里具体可以选取5个 不为零的数据点。
步骤九:在二维边缘点集合中去除抽取的5个点。
步骤十:循环执行步骤八——步骤九,当边缘点集合中剩余的不为零数据点不足5个 时停止循环。
步骤十一:设定一个参数累加次数阈值N2,累加次数n超过了阈值N2时,则将对应参数进行验证,用边缘点集合的所有数据点进行验证。
步骤十二:设定数据点数阈值N3,符合参数的数据点大于阈值N3时,把估计出的参数 设定为原信号的徙动参数。
步骤十三:利用估计的徙动参数构造补偿方程,将补偿方程与原信号回波相乘,实现 运动目标回波轨迹徙动校正。
步骤十四:估计出运动目标回波轨迹上数据点的值进行相干积累,在曲线参数域实现 能量积累。
步骤十五:对目标回波信号经过徙动校正,获取校正后回波信号相位,徙动完全校正 的情况下,其表达式为:
经过上述步骤处理,即可完成对机载SAR对运动目标回波的徙动校正,回波信噪比在 0dB时的校正后结果如图4所示,对比图3和图4可以看出,运动目标弯曲倾斜的徙动轨迹得到有效校正。如图5所示,信号能量能够聚集在参数域。
Claims (1)
1.一种SAR动目标回波徙动校正与方位能量积累方法,具体包括如下步骤:
步骤一:建立机载单基地正侧视SAR空间几何结构,并完成初始化;
在直角坐标系中,O为坐标原点,设(X0,0,0)为地面目标初始位置,其中,X0、0、0分别为地面运动目标的X轴,Y轴和Z轴坐标;
SAR飞行平台的位置坐标为(0,0,H0),其中,0、0和H0分别为SAR飞行平台的X轴,Y轴和Z轴坐标;
步骤二:根据几何模型获取运动目标到飞行平台的瞬时距离历史RM(t),可表示为:
其中,t为方位向时间变量;
步骤三:获取地面运动目标回波S(τ,t),并对回波进行下变频,具体公式为:
其中,rect(·)为距离向和方位向矩形窗,Tsar为合成孔径时间,τ为距离向时间变量,Δτ=2RM(t)/c为线性调频信号的时延,Tτ为发射信号时宽,kτ为距离向调频率,c为光速,fc为载波频率,为虚数单位;
步骤四:将地面运动目标回波进行距离向匹配滤波,滤波后回波信号相位可表示为:
其中,fτ为距离频率变量,θ为SAR天线入射角;
步骤五:构造地面运动目标回波信号初始补偿方程,可表示为:
y=c+bt+at2
其中,c为常数项系数,b一阶项系数,a二阶项系数;
步骤六:对脉冲压缩后的回波进行边缘检测,构造一个二维边缘点集合M,设定一个幅度阈值N1,当边缘点幅度小于N1时将该点幅度置零;
步骤七:定义一个参数累加器A(a,b,c,n),用于存储参数值a,b,c和对应参数累加次数n,并将数组初始化;
步骤八:抽取二维边缘点集合中随机若干个不为零的数据点,计算出一组曲线参数(a1,b1,c1),参数结果存入参数累加器中,并对该组参数进行计数;
步骤九:在二维边缘点集合中去除抽取的若干个点;
步骤十:循环执行步骤八——步骤九,当边缘点集合中剩余的不为零数据点的个数不足时停止循环;
步骤十一:设定一个参数累加次数阈值N2,累加次数n超过了阈值N2时,则将对应参数进行验证,用边缘点集合的所有数据点进行验证;
步骤十二:设定数据点数阈值N3,符合参数的数据点大于阈值N3时,则将估计出的参数设定为原信号的徙动参数;
步骤十三:利用估计的徙动参数构造补偿方程,将补偿方程与原信号回波相乘,实现运动目标回波轨迹徙动校正;
步骤十四:将估计出运动目标回波轨迹上数据点的值进行相干叠加,在曲线参数域可以实现能量积累;
步骤十五:对目标回波信号经过徙动校正,获取校正后回波信号相位,徙动完全校正的情况下,具体为:
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851093.2A CN107656274B (zh) | 2017-09-20 | 2017-09-20 | Sar动目标回波徙动校正与方位能量积累方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201710851093.2A CN107656274B (zh) | 2017-09-20 | 2017-09-20 | Sar动目标回波徙动校正与方位能量积累方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN107656274A CN107656274A (zh) | 2018-02-02 |
CN107656274B true CN107656274B (zh) | 2020-02-18 |
Family
ID=61129572
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201710851093.2A Active CN107656274B (zh) | 2017-09-20 | 2017-09-20 | Sar动目标回波徙动校正与方位能量积累方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN107656274B (zh) |
Families Citing this family (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN110379178B (zh) * | 2019-07-25 | 2021-11-02 | 电子科技大学 | 基于毫米波雷达成像的无人驾驶汽车智能泊车方法 |
CN110456351B (zh) * | 2019-08-29 | 2021-06-15 | 哈尔滨工业大学 | 基于时变幅值lfm信号参数估计的机动目标isar成像方法 |
CN112505647B (zh) * | 2020-09-25 | 2023-11-10 | 北京空间飞行器总体设计部 | 一种基于序贯子图像序列的动目标方位速度估计方法 |
Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749620A (zh) * | 2012-06-13 | 2012-10-24 | 北京航空航天大学 | 一种弹载/机载雷达单脉冲前视成像处理方法 |
KR101207220B1 (ko) * | 2012-10-17 | 2012-12-04 | 국방과학연구소 | 고해상도 sar 진폭 영상의 경계선 정밀 추출 방법 |
CN103412310A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 |
CN104360333A (zh) * | 2014-11-17 | 2015-02-18 | 西安电子科技大学 | 一种同时校正一阶和二阶距离徙动的相参积累检测方法 |
CN105093224A (zh) * | 2015-01-21 | 2015-11-25 | 电子科技大学 | 一种大斜视合成孔径雷达成像处理方法 |
CN105974410A (zh) * | 2016-04-28 | 2016-09-28 | 西安电子科技大学 | 机载雷达的多舰船目标sar和isar混合成像方法 |
CN106970371A (zh) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | 一种基于Keystone和匹配滤波的目标检测方法 |
-
2017
- 2017-09-20 CN CN201710851093.2A patent/CN107656274B/zh active Active
Patent Citations (7)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN102749620A (zh) * | 2012-06-13 | 2012-10-24 | 北京航空航天大学 | 一种弹载/机载雷达单脉冲前视成像处理方法 |
KR101207220B1 (ko) * | 2012-10-17 | 2012-12-04 | 국방과학연구소 | 고해상도 sar 진폭 영상의 경계선 정밀 추출 방법 |
CN103412310A (zh) * | 2013-08-26 | 2013-11-27 | 电子科技大学 | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 |
CN104360333A (zh) * | 2014-11-17 | 2015-02-18 | 西安电子科技大学 | 一种同时校正一阶和二阶距离徙动的相参积累检测方法 |
CN105093224A (zh) * | 2015-01-21 | 2015-11-25 | 电子科技大学 | 一种大斜视合成孔径雷达成像处理方法 |
CN105974410A (zh) * | 2016-04-28 | 2016-09-28 | 西安电子科技大学 | 机载雷达的多舰船目标sar和isar混合成像方法 |
CN106970371A (zh) * | 2017-04-28 | 2017-07-21 | 电子科技大学 | 一种基于Keystone和匹配滤波的目标检测方法 |
Non-Patent Citations (1)
Title |
---|
一种径向匀加速目标包络徙动补偿新方法;蒋千 等;《雷达科学与技术》;20130228;第11卷(第1期);第45-50页 * |
Also Published As
Publication number | Publication date |
---|---|
CN107656274A (zh) | 2018-02-02 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN107966688B (zh) | 基于相位干涉技术的宽带雷达目标速度解模糊方法 | |
CN103744068B (zh) | 双通道调频连续波sar系统的动目标检测成像方法 | |
CN103278820B (zh) | 临近空间慢速平台sar动目标检测方法及成像方法 | |
CN103412310B (zh) | 双基地前视合成孔径雷达地面动目标检测方法与成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN107656274B (zh) | Sar动目标回波徙动校正与方位能量积累方法 | |
CN109324315B (zh) | 基于双层次块稀疏性的空时自适应处理雷达杂波抑制方法 | |
CN112444811A (zh) | 一种综合mimo雷达和isar的目标检测和成像方法 | |
CN103983974B (zh) | 双站调频连续波合成孔径雷达成像方法 | |
CN109633644A (zh) | 一种机动目标isar成像方法 | |
CN109581362B (zh) | 合成孔径雷达高度计在可变脉冲簇模式下的信号处理方法 | |
CN102288963A (zh) | 基于子孔径参数估计的双基地isar图像融合方法 | |
CN107976676B (zh) | 一种机载合成孔径雷达动目标成像方法 | |
CN108196241B (zh) | 一种基于Hough变换的高速动目标速度估计方法 | |
CN104833972A (zh) | 一种双基地调频连续波合成孔径雷达频率变标成像方法 | |
CN110208796B (zh) | 基于奇异值逆滤波的扫描雷达超分辨成像方法 | |
CN109655819B (zh) | 一种基于实孔径多普勒波束锐化的杂波抑制三维成像方法 | |
CN103064084A (zh) | 基于距离频域的解模糊方法 | |
CN103576148B (zh) | 模拟星载sar距离模糊噪声图像的方法 | |
CN104931965B (zh) | 一种固定发射站双基正侧视sar成像方法 | |
CN103245949B (zh) | 一种基于改进理想滤波器的sar方位模糊抑制方法 | |
CN110824473A (zh) | 一种基于子空间的高分辨宽测绘带sar-gmti杂波抑制方法 | |
CN104182942B (zh) | Sar图像方位模糊抑制方法 | |
CN103616669B (zh) | 二维频域中对非均匀场景进行通道均衡方法 | |
CN112230220B (zh) | 基于Radon变换的Deramp-STAP动目标检测方法 |
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 |