CN106291556B - 一种基于双曲线模型的大斜视sar成像算法 - Google Patents
一种基于双曲线模型的大斜视sar成像算法 Download PDFInfo
- Publication number
- CN106291556B CN106291556B CN201610766898.2A CN201610766898A CN106291556B CN 106291556 B CN106291556 B CN 106291556B CN 201610766898 A CN201610766898 A CN 201610766898A CN 106291556 B CN106291556 B CN 106291556B
- Authority
- CN
- China
- Prior art keywords
- signal
- distance
- orientation
- frequency
- indicate
- 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
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/9004—SAR image acquisition 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/9004—SAR image acquisition techniques
- G01S13/9005—SAR image acquisition techniques with optical processing of the SAR signals
-
- 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/9041—Squint mode
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)
- Signal Processing (AREA)
- Radar Systems Or Details Thereof (AREA)
Abstract
本发明涉及一种基于双曲线模型的大斜视SAR成像算法,采用更精确的雷达与目标之间的瞬时斜距的双曲线模型,且对距离压缩后的信号进行方位频率三次项的补偿,与时域校正线性距离走动的RD算法相比,本发明方法的成像目标结果不仅聚焦效果好、方位距离向旁瓣对称,而且距离向和方位向的中心目标峰值旁瓣比更高。
Description
技术领域
本发明属于雷达技术领域,具体涉及一种新的大斜视合成孔径雷达(SyntheticAperture Radar,SAR)成像算法。
背景技术
合成孔径雷达工作于大斜视模式下,发射信号与载机飞行方向不垂直,导致了回波信号存在严重的距离方位耦合,距离向调频率变化和距离徙动严重等问题,并且距离徙动呈现出大距离走动、小距离弯曲的特点。本算法利用时域校正距离走动方法补偿掉距离徙动中的线性距离走动分量,大大降低了距离向和方位向的耦合,并且对时域校正距离走动后所产生的方位频率三次相位项进行补偿,有效改善了方位向旁瓣不对称畸变的问题,可被用于大斜视SAR成像处理。
时域校正线性距离走动的距离-多普勒(Range-Doppler,RD)算法为:
步骤1 根据斜视SAR成像空间几何模型确定雷达与目标之间的瞬时斜距;
步骤2 依据雷达参数确定接收到的经解调后的回波基频信号;
步骤3 由步骤1中的雷达与目标之间的瞬时斜距泰勒展开式中确定线性距离走动项,并对其进行补偿以校正距离走动;
步骤4 在二维频域中队信号距离向匹配滤波和残余距离徙动校正,使得信号的距离方位项解耦合;
步骤5 对于距离徙动校正后二位解耦合的信号可以沿着方位向对信号进行方位向匹配滤波,从而得到聚焦后目标图像。
时域校正线性距离走动的距离-多普勒算法具有成像速度快、普适性强等优点。但是其在斜视角度增大的情况下聚焦后的目标会出现方位向旁瓣不对称畸变的情况,这就导致了当多目标成像时成像质量下降,不利于雷达成像系统对图像质量的要求。
发明内容
要解决的技术问题
为了克服成像目标聚焦后方位向旁瓣不对称的情况,本发明基于更精确的目标与雷达之间距离的双曲线模型,并详细分析了距离走动校正后二维频域中的信号,提出了一种利用更精确的数学模型分析二维频域中回波信号特点优化改善了大斜视SAR雷达成像系统性能指标的方法。
技术方案
一种基于双曲线模型的大斜视SAR成像算法,其特征在于步骤如下:
步骤1:根据距离等式采用双曲线模型,确定雷达与目标之间的瞬时斜距为
其中,R(t,R0)表示雷达与目标之间的瞬时斜距,t表示载机飞行的慢时间,R0表示载机航线与场景中心线之间的最短距离,V表示载机飞行的速度,θ0表示雷达发射信号的斜视角度;
步骤2:根据瞬时斜距和载机雷达工作的参数,确定经解调的二维时域中的回波信号
其中,τ是距离向时间,G(t)表示二维时域中回波信号的振幅,ma(t)表示方位向窗函数,mr(·)表示距离向窗函数,c表示光速,Kr表示发射信号的调频率,λ表示发射信号的波长;
步骤3:对回波信号进行距离向FFT,将信号变换到距离频域中,在距离频域中对信号乘以距离走动校正函数,所述的距离走动校正函数:
其中,fτ表示距离向频率,fc表示对发射信号进行调制的载频;
步骤4:对距离走动校正后的信号进行方位向FFT,将信号变换到二维频域中;对信号依次进行距离压缩、方位频率三次项补偿和残余距离徙动校正:
所述的距离压缩函数为:
其中,f为方位向频率,V′=V cosθ0;
对距离压缩后的信号进行方位频率三次项的补偿,这一步补偿了因距离走动校正而引起的压缩后目标方位向旁瓣的不对称畸变问题,所述的方位频率三次项补偿函数为:
所述的残余距离徙动校正函数为:
步骤5:对步骤4得到的信号进行距离向IFFT,将信号变换到距离多普勒域中,然后沿着方位向对信号进行方位向压缩;所述的方位向压缩函数为:
步骤6:对压缩后的信号再进行方位向IFFT,得到聚焦目标图像。
有益效果
本发明的有益效果是在时域校正线性距离走动的RD算法的基础上,通过详细分析大斜视下成像模型和雷达回波信号的特点,采用更精确的雷达与目标之间的瞬时斜距的双曲线模型,提出了一种新的大斜视SAR成像算法。与时域校正线性距离走动的RD算法相比,本发明方法的成像目标结果不仅聚焦效果好、方位距离向旁瓣对称,而且距离向和方位向的中心目标峰值旁瓣比更高。采用本法名的大斜视SAR雷达成像算法的雷达成像系统可以得到更加精确的目标图像。
附图说明
图1是大斜视机载SAR空间几何模型图。
H为载机飞行高度,Rb为航线距场景中心的距离,R0为波束相位中心的斜距,θ0为斜视角,R(t)为雷达与目标之间的瞬时斜距,Vt为载机飞行速度与时间。
图2是距离走动校正后的SAR空间几何等效示意模型图。
图3是本发明算法的流程图图。
图4是在斜视角为45°下时的3点目标成像效果图。图a分别是参考算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图。图b分别是本文算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图
图5是在斜视角为60°下时的3点目标成像效果图。图a分别是参考算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图。图b分别是本文算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图
图6是在斜视角为70°下时的3点目标成像效果图。图a分别是参考算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图。图b分别是本文算法的二维时域目标信号和中心目标等高线图、中心目标距离向截面图和方位向截面图。
图7中心目标峰值旁瓣比
具体实施方式
现结合实施例、附图对本发明作进一步描述:
利用以上所采取的距离等式更加精确的双曲线模拟雷达成像系统的雷达与目标之间的瞬时斜距,补偿二维频域中因距离走动校正产生的方位频率三次项对雷达回波数据进行处理,比较仿真系统成像性能。
步骤1:根据距离等式采用的双曲线模型,确定雷达与目标之间的瞬时斜距为
式(1)中R(t,R0)表示雷达与目标之间的瞬时斜距,t表示载机飞行的慢时间,R0表示载机航线与场景中心线之间的最短距离,V表示载机飞行的速度,θ0表示雷达发射信号的斜视角度。
对雷达与目标之间的瞬时斜距进行泰勒展开如下式所示
从式(3)中可以确定处距离徙动量中的线性距离走动分量为ΔR=-sinθ0Vt,距离弯曲分量为
步骤2:利用步骤1中的雷达与目标之间的瞬时斜距和载机雷达工作的参数调频率Kr和波长λ来确定经解调的二维时域中的回波信号
τ是距离向时间,G(t)表示二维时域中回波信号的振幅,ma(t)表示方位向窗函数,mr(·)表示距离向窗函数,c表示光速,Kr表示发射信号的调频率,λ表示发射信号的波长。
步骤3:将式(1)代入到式(2)中,并对步骤2中的二维时域中的回波信号进行距离向FFT,将信号变换到距离频域中。因此可在距离频域中对信号乘以函数Hw(fτ),对距离频域中的信号进行距离走动校正
其中,fτ表示距离向频率,fc表示对发射信号进行调制的载频;
步骤4:对步骤3在距离频域中进行距离走动校正后的信号对其做方位向FFT,将信号变换到二维频域中。对信号距离压缩、位频率三次项补偿和残余距离徙动校正。
首先对信号进行距离向压缩,压缩函数为
其中,f为方位向频率,V′=V cosθ0;
对距离压缩后的信号进行方位频率三次项的补偿,这一步补偿了因距离走动校正而引起的压缩后目标方位向旁瓣的不对称畸变问题,补偿函数为
通过对雷达与目标之间瞬时斜距的泰勒展开进行分析,在第三步中对距离徙动中的距离走动分量进行了补偿,为了对距离徙动进行完全的补偿,使信号的距离方位向解耦合,则必须要对信号中残余的距离徙动量(即距离弯曲量)进行校正,校正函数为
通过残余距离徙动校正后,信号的距离方位向解耦合,并将因距离走动校正在二维频域中引起的方位向频率三次项用共轭消除法进行补偿,使得信号的方位向压缩可以沿着每个距离门进行。
步骤5:对距离徙动校正后的二维解耦合的信号进行距离向IFFT,将信号变换到距离多普勒域中,然后沿着方位向对信号进行方位向压缩。方位向压缩函数为
步骤6:对压缩后的信号再进行方位向IFFT,得到聚焦目标图像。
图3为本发明所采用的基于时域去走动校正的RD算法流程图图,将接收到的回波信号进行距离向FFT,将信号转换到距离多普勒域中对其进行距离走动校正;将距离走动校正后的信号在二维频域进行距离向压缩、距离徙动校正和补偿因时域走动校正而引起的方位频率三次项的问题,最后在距离多普勒域中进行方位向压缩便可以得到目标的二维像。
在SAR工作于大斜视模式下时,分别选取参考算法和本发明中的算法不同斜视角下的回波数据进行方针处理,仿真结果参见图4、图5和图6,分别为45°、60°和70°的斜视角度下的结果对比。从结果中可以看出本发明中的算法在斜视角度变大的情况下仍能使聚焦后目标旁瓣对称且中心目标的峰值旁瓣比相较于参考算法有所提升,可以满足大斜视SAR的成像要求。
Claims (1)
1.一种基于双曲线模型的大斜视SAR成像算法,其特征在于步骤如下:
步骤1:根据距离等式采用双曲线模型,确定雷达与目标之间的瞬时斜距为
其中,R(t,R0)表示雷达与目标之间的瞬时斜距,t表示载机飞行的慢时间,R0表示载机航线与场景中心线之间的最短距离,V表示载机飞行的速度,θ0表示雷达发射信号的斜视角度;
步骤2:根据瞬时斜距和载机雷达工作的参数,确定经解调的二维时域中的回波信号
其中,τ是距离向时间,G(t)表示二维时域中回波信号的振幅,ma(t)表示方位向窗函数,mr(·)表示距离向窗函数,c表示光速,Kr表示发射信号的调频率,λ表示发射信号的波长;
步骤3:对回波信号进行距离向FFT,将信号变换到距离频域中,在距离频域中对信号乘以距离走动校正函数,所述的距离走动校正函数:
其中,fτ表示距离向频率,fc表示对发射信号进行调制的载频;
步骤4:对距离走动校正后的信号进行方位向FFT,将信号变换到二维频域中;对信号依次进行距离压缩、方位频率三次项补偿和残余距离徙动校正:
所述的距离压缩函数为:
其中,f为方位向频率,V′=Vcosθ0;
对距离压缩后的信号进行方位频率三次项的补偿,这一步补偿了因距离走动校正而引起的压缩后目标方位向旁瓣的不对称畸变问题,所述的方位频率三次项补偿函数为:
所述的残余距离徙动校正函数为:
步骤5:对步骤4得到的信号进行距离向IFFT,将信号变换到距离多普勒域中,然后沿着方位向对信号进行方位向压缩;所述的方位向压缩函数为:
步骤6:对方位向压缩后的信号再进行方位向IFFT,得到聚焦目标图像。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610766898.2A CN106291556B (zh) | 2016-08-29 | 2016-08-29 | 一种基于双曲线模型的大斜视sar成像算法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201610766898.2A CN106291556B (zh) | 2016-08-29 | 2016-08-29 | 一种基于双曲线模型的大斜视sar成像算法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN106291556A CN106291556A (zh) | 2017-01-04 |
CN106291556B true CN106291556B (zh) | 2018-10-19 |
Family
ID=57674706
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201610766898.2A Active CN106291556B (zh) | 2016-08-29 | 2016-08-29 | 一种基于双曲线模型的大斜视sar成像算法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN106291556B (zh) |
Families Citing this family (4)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN111208514B (zh) * | 2020-01-06 | 2023-05-05 | 华侨大学 | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 |
CN113329331B (zh) * | 2021-05-19 | 2022-04-29 | 天津大学 | 一种基于改进双曲线模型的rfid定位方法 |
CN115616576B (zh) * | 2022-12-16 | 2023-03-10 | 西安空间无线电技术研究所 | 大斜视成像下的sar天线方向图校正方法及装置 |
CN116299464B (zh) * | 2023-05-17 | 2023-07-21 | 西安电子科技大学 | 一种高速高机动宽幅扇扫sar成像方法 |
Family Cites Families (6)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
ITTO20110526A1 (it) * | 2011-06-15 | 2012-12-16 | Thales Alenia Space Italia S P A C On Unico Socio | Acquisizione di immagini sar per calcolare una quota o un modello digitale di elevazione tramite elaborazioni interferometriche |
CN102393518B (zh) * | 2011-10-27 | 2013-07-17 | 上海大学 | 一种适用于大斜视角的机载sar成像方法 |
CN102879784B (zh) * | 2012-10-25 | 2014-06-18 | 西安电子科技大学 | 四种模式的合成孔径雷达统一成像方法 |
CN104020472A (zh) * | 2014-05-26 | 2014-09-03 | 中国电子科技集团公司第三十八研究所 | 一种便于实时处理的方位ncs大斜视sar成像方法 |
CN104459693A (zh) * | 2014-12-01 | 2015-03-25 | 西安电子科技大学 | 基于gpu的弹载sar前斜视成像方法 |
CN105093224A (zh) * | 2015-01-21 | 2015-11-25 | 电子科技大学 | 一种大斜视合成孔径雷达成像处理方法 |
-
2016
- 2016-08-29 CN CN201610766898.2A patent/CN106291556B/zh active Active
Also Published As
Publication number | Publication date |
---|---|
CN106291556A (zh) | 2017-01-04 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
CN105259552B (zh) | 一种基于非线性调频信号的合成孔径雷达成像方法和装置 | |
CN106291556B (zh) | 一种基于双曲线模型的大斜视sar成像算法 | |
CN108459321B (zh) | 基于距离-方位圆模型的大斜视高分辨率sar成像方法 | |
CN102749620B (zh) | 一种弹载/机载雷达单脉冲前视成像处理方法 | |
CN108710111B (zh) | 一种机载双基前视sar方位相位的二维空变校正方法 | |
CN105759263B (zh) | 一种高分辨率大场景下的星载斜视sar雷达成像方法 | |
CN109799495B (zh) | 一种用于高保真阵列处理的宽带时延估计方法 | |
CN109143237B (zh) | 适用于任意平台轨迹的双基聚束sar的pfa波前弯曲校正方法 | |
CN109870686B (zh) | 一种基于改进斜距模型的曲线轨迹sar斜视成像方法 | |
CN103293520B (zh) | 基于等效距离方程的sar地面运动目标成像方法 | |
CN108427115B (zh) | 合成孔径雷达对运动目标参数的快速估计方法 | |
CN107607952B (zh) | 基于电磁涡旋波的三维合成孔径雷达成像方法 | |
CN110673144B (zh) | 一种基于时间变标的子孔径大斜视sar成像处理方法 | |
CN103869311A (zh) | 实波束扫描雷达超分辨成像方法 | |
CN105607055A (zh) | 一种基于天线方向图的机载雷达单脉冲前视成像方法 | |
CN103489176A (zh) | 一种对于严重几何畸变的sar图像进行同名点提取的方法 | |
CN106646466B (zh) | 一种基于主成分分析的加权后向投影算法的成像方法 | |
CN111722225B (zh) | 基于先验相位结构信息的双基sar两维自聚焦方法 | |
CN106054187A (zh) | 基于斜距模型下的大斜视sar曲线轨迹波数域成像方法 | |
CN113281727A (zh) | 一种基于水平线列阵的输出增强的波束形成方法及其系统 | |
CN111208514A (zh) | 一种曲线运动轨迹SAR的切比雪夫斜距模型和Chirp Scaling成像方法 | |
CN106772373B (zh) | 针对任意地面运动目标的sar成像方法 | |
CN108562898A (zh) | 一种前侧视sar的距离和方位两维空变自聚焦方法 | |
CN103792534A (zh) | 一种基于先验相位结构知识的sar两维自聚焦方法 | |
CN114994676A (zh) | 一种一站固定式双站低频超宽带sar运动目标成像方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
C06 | Publication | ||
PB01 | Publication | ||
C10 | Entry into substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant | ||
TR01 | Transfer of patent right | ||
TR01 | Transfer of patent right |
Effective date of registration: 20230626 Address after: 710075 Room 014, F2003, 20th Floor, Building 4-A, Xixian Financial Port, Fengdong New Energy Jinmao District, Xixian New District, Xi'an City, Shaanxi Province Patentee after: Shaanxi Debaoyun Technology Co.,Ltd. Address before: 710072 No. 127 Youyi West Road, Shaanxi, Xi'an Patentee before: Northwestern Polytechnical University |