CN112034459A - 基于天线方向图补偿的直线合成孔径雷达成像方法和系统 - Google Patents

基于天线方向图补偿的直线合成孔径雷达成像方法和系统 Download PDF

Info

Publication number
CN112034459A
CN112034459A CN202010816258.4A CN202010816258A CN112034459A CN 112034459 A CN112034459 A CN 112034459A CN 202010816258 A CN202010816258 A CN 202010816258A CN 112034459 A CN112034459 A CN 112034459A
Authority
CN
China
Prior art keywords
synthetic aperture
linear synthetic
imaging
antenna
radar
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
Application number
CN202010816258.4A
Other languages
English (en)
Other versions
CN112034459B (zh
Inventor
王鹏宇
宋千
Current Assignee (The listed assignees may be inaccurate. Google has not performed a legal analysis and makes no representation or warranty as to the accuracy of the list.)
Hunan Jianwei Intelligent Technology Co ltd
Original Assignee
Individual
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by Individual filed Critical Individual
Priority to CN202010816258.4A priority Critical patent/CN112034459B/zh
Publication of CN112034459A publication Critical patent/CN112034459A/zh
Application granted granted Critical
Publication of CN112034459B publication Critical patent/CN112034459B/zh
Active legal-status Critical Current
Anticipated expiration legal-status Critical

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/904—SAR modes
    • G01S13/9056—Scan SAR mode
    • 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)

Abstract

本发明涉及基于天线方向图补偿的直线合成孔径雷达成像方法和系统。一种基于天线方向图补偿的直线合成孔径雷达成像方法，包括步骤：S1、对直线合成孔径雷达天线方向图进行精确测量；S2、形成若干成像像素点；划定为多个直线合成孔径；S3、计算每个所述成像像素点与所述直线合成孔径之间的距离延迟、方位角与俯仰角；S4、多个所述相位补偿因子形成相位补偿矩阵；S5、得到所述一维距离像；S6、得到所述成像像素点的补偿图像；S7、得到所述直线合成孔径的所述子图像；S8、所述子图像进行相干叠加得到最终成像。本发明通过对所述成像像素点进行补偿，将由天线相位方向图引入的相位不一致性误差进行消除，提高了直线孔径SAR大积累角条件下的方位分辨率。

Description

基于天线方向图补偿的直线合成孔径雷达成像方法和系统
技术领域
本发明涉及雷达成像领域，尤其涉及基于天线方向图补偿的直线合成孔径雷达成像方法和系统。
背景技术
雷达(Radio Detection and Ranging,Radar)通过主动辐射电磁波来获取目标的距离、方位和速度等信息。为获取更多的观测信息和目标特性，实孔径雷达成像技术应运而生，其原理是：对发射的宽带信号进行脉冲压缩分辨距离维远近分布的点目标；通过波束形成分辨平行于天线方向分布的点目标。在实孔径雷达成像中，方位分辨力由雷达波束宽度决定的。孔径越长，分辨率越高；而距离越远，分辨率越低。对于远距离作用雷达，为了获取足够高的方位分辨率，天线孔径需要达到足够的长度。
1951年6月，美国Goodyear Aerospace公司的C.Wiley提出利用多普勒波束锐化技术(Doppler Beam Sharpening,DBS)来提高方位分辨率，为雷达信号对地表进行二维高分辨成像提供依据。DBS利用空间采样构建了一个等效的长天线孔径，突破了实孔径技术对方位分辨率的限制，标志着合成孔径雷达(Synthetic Aperture Radar,SAR)成像技术的诞生。SAR利用雷达与目标之间的相对运动获取不同时刻、不同空间位置的观测信息，通过相干处理实现方位高分辨成像。与光学成像相比，SAR成像具有全天时、全天候和远距探测等优点。
根据工作平台的不同，SAR可以分为机载SAR、星载SAR、弹载SAR、车载SAR和球载SAR等。同时SAR系统根据具体工作模式的不同又分为条带SAR、聚束SAR、循序扫描地形观测(Terrain Observation by Progressive Scans,TOPS)SAR以及滑动聚束SAR，并且每种工作模式又可以衍生出正侧视和斜视两种几何构型。其中，正侧视条带SAR是最常见的工作模式，在该模式下，雷达安装于搭载平台，天线指向角度固定，当雷达沿直线运动时，天线波束均匀扫过地面，形成带状观测区域。
在正侧视条件下，雷达波束中心线与测绘带垂直，设雷达在方位向波束角宽度为θ，合成孔径的长度为L。建立雷达-目标坐标系，以雷达坐标为(X,0)、目标P的坐标为P(X0,Y0)，雷达沿X轴做匀速运动，速度为v，设探测区域内有一点目标P(X0,Y0)，雷达到目标的距离为R(t)，最短径向距离设为R0，则时刻t雷达到目标的瞬时距离为：
Figure BDA0002632823000000011
其中X＝vt，在理想运动条件下，且当目标点和雷达满足R0＞＞(X0-X)时，用菲涅尔近似对式(1)进行简化得：
Figure BDA0002632823000000021
假设雷达发射连续正弦波s(t)＝Aexp(jωct)，则理想点目标回波信号r(t)为：
r(t)＝Aexp[jωc(t-τ)] (3)
式中，A为天线幅度增益，ωc为发射信号的载频，τ＝2R(t)/c为目标回波信号延迟，c为光速。对式(3)进行相干检波并幅度归一化得：
r(t)＝exp[-jωcτ] (4)
在SAR孔径合成过程中，雷达与目标发生相对运动，即回波存在多普勒频移，对式(4)的相位求导数可得目标方位向回波的多普勒频率为：
Figure BDA0002632823000000022
从式子(5)可以看出回波信号在方位向上为线性调频信号，对式(5)求导数可得多普勒调频率：
Figure BDA0002632823000000023
由于受天线相位方向图的影响，SAR方位积累过程中通常考虑3dB波束宽度的积累角，因此合成孔径长度、波束角以及目标径向垂直距离的关系近似满足L≈R0θ，则相干积累时间T≈L/v。由多普勒调频率可以得出多普勒带宽，即
Figure BDA0002632823000000024
由式(7)可进一步得到系统的方位向分辨率为：
Figure BDA0002632823000000025
由式(8)不难发现，当发射信号一定的条件下，系统方位分辨率与积累角成反比。积累角越大，合成孔径越长，方位分辨率就越高。与方位分辨率相比，SAR距离向分辨率ρr只与发射信号带宽B相关，即ρr＝c/2B。因此，SAR系统在系统参数和几何配置一定的条件下，增大积累角有利于系统成像分辨率的提高。随着积累角的逐渐增大，系统方位分辨率逐渐提高，当积累约束突破天线3dB波束限制时，由天线相位方向图引入的相位不一致性误差将变得无法忽略，就会造成图像分辨率降低的情况发生。
申请号为CN201210166779.5的专利文献公开了一种线性调频连续波雷达高度表的数字化处理方法，属于无线电高度表领域，包括步骤一：从雷达高度表的前端正交混频器得到差拍信号，进行A/D采样得到数字信号序列；步骤二：对数字信号序列进行傅里叶变换并计算谱线的幅度；步骤三：根据幅度谱谱线幅度进行判决并测量数字信号序列的差拍频率；步骤四：进行频谱细化，计算出锯齿波斜率修正量，对锯齿波电压信号的斜率进行修正；步骤五：根据锯齿波电压信号斜率计算当前高度表的高度值。但是依然没有解决上述问题。
因而现有的雷达成像技术还存在不足，有待改进和提高。
发明内容
鉴于上述现有技术的不足之处，本发明的目的在于提供基于天线方向图补偿的直线合成孔径雷达成像方法和系统，能够解决在大积累角条件下直线合成孔径雷达方位分辨率受天线相位方向图不一致性影响而下降的问题。
为了达到上述目的，本发明采取了以下技术方案：
一种基于天线方向图补偿的直线合成孔径雷达成像方法，包括步骤：
S1、对直线合成孔径雷达天线方向图进行精确测量，获得雷达天线相位特性矩阵；
S2、对雷达探测区域进行成像网络剖分，形成若干成像像素点；确定直线合成孔径雷达的所有直线合成孔径位置；
S3、获取所述直线合成孔径雷达在单一所述直线合成孔径上波束照射范围内的所有所述成像像素点，计算每个所述成像像素点与所述直线合成孔径之间的距离延迟、方位角与俯仰角；
S4、根据步骤S3中每个所述成像像素点的所述方位角和所述俯仰角结合所述雷达天线相位特性矩阵计算得到每个像素点的相位补偿因子；多个所述相位补偿因子形成相位补偿矩阵；
S5、所述直线合成孔径雷达对每个所述成像像素点的采样回波进行脉冲压缩处理得到所述一维距离像；
S6、根据所述距离延迟将所述一维距离像中的复数值投影至相应的所述成像像素点，同时对所述成像像素点进行多普勒相位和剩余视频相位的空变补偿，得到所述成像像素点的补偿图像；
S7、将得到所述补偿图像与所述相位补偿矩阵进行复共轭相乘，得到所述直线合成孔径的所述子图像；
S8、将步骤S2中确定的所有所述直线合成孔径均执行步骤S3-S7，将各个所述直线合成孔径得到的所述子图像进行相干叠加得到最终成像。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，所述步骤S5中，所述一维距离像的获取步骤为：
S51、所述直线合成孔径雷达对采样回波进行混频和幅度归一化后，得到差频信号；
S52、进而对所述差频信号进行傅立叶变换后得到采样回波频谱；
S53、再根据频率与目标距离的对应关系，得到所述采样回波的一维距离像。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，所述步骤S51中，所述差频信号的求解公式为：
所述采样回波信号表达式为：
Figure BDA0002632823000000041
其中，t为时刻；R为时刻t直线合成孔径雷达到目标的瞬时距离；c为光速；TP为信号脉冲；Kr表示调频斜率；α为方位角；β为俯仰角；ω(α,β)为天线相位因子，即在方位角α和俯仰角β方向上天线的相位值；
所述差频信号的表达式为：
Figure BDA0002632823000000042
其中，fc为发射信号的载频。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，所述步骤S52中，所述采样回波频谱的表达式为：
Figure BDA0002632823000000043
其中，f为频率。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，所述步骤S53中，所述频率与目标距离的对应关系表达是为：
f＝2rKr/c；
其中，f为频率；Kr表示调频斜率；c为光速；r为目标距离。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，所述一维距离像的表达式为：
Figure BDA0002632823000000051
其中，r为目标距离；Kr表示调频斜率；c为光速；ω(α,β)为天线相位因子；R为时刻t直线合成孔径雷达到目标的瞬时距离。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，在所述步骤S2中，所述直线合成孔径的数量根据方位采样间隔和雷达扫描探测范围确定。
优选的所述的基于天线方向图补偿的直线合成孔径雷达成像方法，在所述步骤S3中，在计算所述成像像素点的俯仰角时，叠加当前雷达天线探测俯视角进行计算。
一种使用所述的基于天线方向图补偿的直线合成孔径雷达成像方法的直线合成孔径雷达成像系统，直线合成孔径雷达装设在可移动装置上，所述可移动装置用于按照预定的直线合成孔径移动。
优选的所述的直线合成孔径雷达成像系统，所述可移动装置按照预定的直线合成孔径直线移动。
相较于现有技术，本发明提供的基于天线方向图补偿的直线合成孔径雷达成像方法和系统，具有以下有益效果：
本发明通过对成像区域进行剖分得到若干成像像素点，进而基于若干所述成像像素点进行补偿，实现对天线相位方向图的相位不一致性的补偿，将由天线相位方向图引入的相位不一致性误差进行消除，突破了3dB波束角的限制，提高了直线孔径SAR大积累角条件下的方位分辨率，改善了二维成像质量。
附图说明
图1是本发明提供的基于天线方向图补偿的直线合成孔径雷达成像方法的流程图；
图2是本发明提供的直线SAR测量几何示意图；
图3是本发明提供的成像仿真得到的天线相位方向图；
图4是本发明提供的点目标未叠加天线相位方向图二维成像结果；
图5是本发明提供的图4点目标二维图像的一维剖面图；
图6是本发明提供的点目标叠加天线相位方向图二维成像结果；
图7是本发明提供的图6点目标二维图像的一维剖面图；
图8是本发明提供的基于像素点的天线相位方向图补偿点目标二维成像结果；
图9是本发明提供的图8点目标二维图像的一维剖面图。
具体实施方式
为使本发明的目的、技术方案及效果更加清楚、明确，以下参照附图并举实施例对本发明进一步详细说明。应当理解，此处所描述的具体实施例仅用以解释本发明，并不用于限定本发明。
请一并参阅图1-图9，本发明提供一种基于天线方向图补偿的直线合成孔径雷达成像方法，包括步骤：
S1、对直线合成孔径雷达天线方向图进行精确测量，获得雷达天线相位特性矩阵；
S2、对雷达探测区域进行成像网络剖分，形成若干成像像素点；确定直线合成孔径雷达的所有直线合成孔径位置；具体的，所述雷达探测区域是雷达波束扫描过的区域，在本实施例中，使用直线合成孔径雷达进行探测，所述直线合成孔径雷达沿直线移动，则形成的是长方形的探测区域，将所述探测区域进行网络剖分对探测区域进行细化处理，进而能够基于对每个所述成像像素点的补偿，实现对探测区域进行精细化补偿，进而提高在大积累角条件下的方位分辨率；剖分的精度视现场的具体情况进行设定，本发明不做限定；直线合成孔径是所述直线合成孔径雷达在探测过程中形成的合成孔径，在本实施例中，所述直线合成孔径是每个直线合成孔径的中线位置；
S3、获取所述直线合成孔径雷达在单一所述直线合成孔径上波束照射范围内的所有所述成像像素点，计算每个所述成像像素点与所述直线合成孔径之间的距离延迟、方位角与俯仰角；所述距离延迟是所述直线合成孔径雷达发出的发射信号与相应接收信号之间的时延，即距离向接收时延；
S4、根据步骤S3中每个所述成像像素点的所述方位角和所述俯仰角结合所述雷达天线相位特性矩阵计算得到每个像素点的相位补偿因子；多个所述相位补偿因子形成相位补偿矩阵；
S5、所述直线合成孔径雷达对每个所述成像像素点的采样回波进行脉冲压缩处理得到所述一维距离像；
S6、根据所述距离延迟将所述一维距离像中的复数值投影至相应的所述成像像素点，同时对所述成像像素点进行多普勒相位和剩余视频相位的空变补偿，得到所述成像像素点的补偿图像；
S7、将得到所述补偿图像与所述相位补偿矩阵进行复共轭相乘，得到所述直线合成孔径的所述子图像；
S6、将步骤S2中确定的所有所述直线合成孔径均执行步骤S3-S5，将各个所述直线合成孔径得到的所述子图像进行相干叠加得到最终成像。
具体的，本发明提供的直径合成孔径雷达成像方法首选根据成像区域和方位孔径位置计算每个成像像素点对应当前孔径的方位角与俯仰角；然后依据天线相位方向图计算天线相位方向图相位补偿因子；最后在二维成像积累过程中对每个像素点进行天线相位方向图相位补偿。能够有效解决天线相位方向图相位不一致性而造成直线合成孔径雷达方位向聚焦不佳的问题，通过对像素点进行的天线相位方向图空变相位补偿，降低了天线相位方向图相位不一致性对雷达系统方位聚焦和目标RCS(Radar Cross section，雷达散射截面)散射特性影响，增大了直线合成孔径雷达成像方位积累角，突破了常规直线合成孔径雷达3dB波束宽度的成像制约，提高了方位分辨率，改善了目标成像质量。应当说明的是，所述积累角为雷达检测中的天线波束角，使用本发明提供的成像方法进行雷达成像，突破了所述积累角制约在3dB的限制，此时所述积累角即使扩大几倍，成像依然清晰，极大的扩大了雷达检测的方位分辨率。
作为优选方案，本实施例中，所述一维距离像的获取步骤包括：
S51、所述直线合成孔径雷达对采样回波进行混频和幅度归一化后，得到差频信号；其中，所述采样回波信号表达式为：
Figure BDA0002632823000000071
其中，t为时刻；R为时刻t直线合成孔径雷达到目标的瞬时距离；c为光速；TP为信号脉冲；Kr表示调频斜率；α为方位角；β为俯仰角；ω(α,β)为天线相位因子，即在方位角α和俯仰角β方向上天线的相位值；
所述差频信号的表达式为：
Figure BDA0002632823000000072
其中，fc为发射信号的载频；
S52、进而对所述差频信号进行傅立叶变换后得到采样回波频谱；其中，所述采样回波频谱的表达式为：
Figure BDA0002632823000000081
其中，f为频率；
S53、再根据频率与目标距离的对应关系f＝2rKr/c，得到所述采样回波的一维距离像，一维距离像的表达式为：
Figure BDA0002632823000000082
其中，r为目标距离。上述所有公式中，相同的项均表达相同的物理量，未重复标示。
具体的，下面结合附图并举例对本发明进行详细描述。图2给出了本发明直线SAR的测量几何示意图，表1给出了参与测试的直线孔径SAR系统的几何参数。
表1直线SAR系统、几何参数
参数名称 符号 参数值
载频 f<sub>c</sub> 10GHz
发射信号带宽 B 1GHz
信号脉宽 T<sub>p</sub> 800us
调频率 K<sub>r</sub> K<sub>r</sub>＝B/T<sub>p</sub>
积累角 θ 50°
光速 c 3×10<sup>8</sup>m/s
目标坐标 (X<sub>0</sub>,Y<sub>0</sub>,Z<sub>0</sub>) (0,200,-100)
请着重参阅图1，所述直线合成孔径雷达沿X轴正向(箭头方向)匀速运动，t时刻的坐标为(X,Y,Z)，直线合成孔径雷达测量几何如图1所示。
假设理想点目标P位于探测区域的中心，坐标为(X0,Y0,Z0)，所述直线合成孔径雷达发射调频连续波信号，收发天线共用同一个雷达天线(即所述调频连续波信号的发射点和采样回波的接收点为相同的位置，例如在图1中为坐标(X,Y,Z))，可以得到雷达天线的相位方向图的相位特性因子为exp[jω(α,β)]，点目标P的采样回波信号sr(t)可表示为：
Figure BDA0002632823000000083
此时，根据所述采样回波信号和发射信号(即调频连续信号)经过混频和幅度归一化后得到的差频信号为：
Figure BDA0002632823000000084
其中，
Figure BDA0002632823000000091
为t时刻直径合成孔径雷达与目标P之间的瞬时径向距离，Kr表示调频斜率。在一个周期内，R为一常数；所述差频信号表达式中，
Figure BDA0002632823000000092
为第一项，表示所述瞬时径向距离所对应的相位；
Figure BDA0002632823000000093
为第二项，表示回波的多普勒效应，这是进行方位向脉压所必须处理的；
Figure BDA0002632823000000094
为第三项，是解线性调频方法所特有，称为剩余视频相位，一般为常数：
Figure BDA0002632823000000095
当第三项为常数时，需要补偿的相位项可表示为：
Figure BDA0002632823000000096
对式(2)进行傅立叶变换得采样回波频谱：
Figure BDA0002632823000000097
根据频率与目标距离的对应关系f＝2rKr/c，可以得到采样回波一维距离像的表达式为：
Figure BDA0002632823000000098
从式(5)可以看出：当目标与直线合成孔径雷达相对固定时，雷达天线相位方向图相位特性因子均可视为常数。由于雷达天线的相位方向图相位特性因子是关于方位角α和俯仰角β的空变函数，因此需要在成像过程中逐像素点进行补偿。
以下结合图2-图9，使用后向投影(Back Projection,BP)算法为例对基于天线方向图补偿的直线合成孔径雷达成像方法进行详细说明：
1)根据表1中的几何参数进行直线合成孔径雷达进行设定，设定后的直线合成孔径雷达的天线方向图如图3所示，所述直线合成孔径雷达(X,Y,Z)在t时刻与点目标P(X0,Y0,Z0)之间的位置关系如图2所示，计算此时点目标P在对应所述直线合成孔径雷达在此位置形成的孔径的方位角α和俯仰角β，在仿真回波上叠加雷达天线相位方向图的相位特性因子；
2)对所有直线合成孔径的采样回波均进行脉冲压缩处理，即均执行步骤S511-S513得到所述直线合成孔径雷达位于所有直线合成孔径的一维距离像；此处应当说明，所述直线合成孔径的数量根据方位采样间隔和雷达扫描探测范围进行确定，数量设为N；
3)对雷达探测区域进行成像网络剖分，得到若干成像像素点；
4)计算成像网格中各个所述成像像素点到第i个所述直线合成孔径的距离延迟、方位角和俯仰角；具体的，每个所述直线合成孔径均对应多个所述成像像素点；
5)根据步骤4计算的距离延迟将第i(1≤i≤N)个方位孔径对应的一维距离像的复数值投影至成像网格中对应像素点，同时对各像素点进行多普勒相位和剩余视频相位的空变补偿，获取第i个补偿图像；
6)根据步骤4计算的方位角和俯仰角，结合所述雷达天线相位特性举证(由所述天线方位图得到)得到每个所述成像像素点的相位补偿因子，进而整合成所述相位补偿矩阵，并对步骤5)获取的第i个所述补偿图像进行复共轭相乘，完成天线相位方向图相位补偿得到第i个所述直线合成孔径的子图像；
7)重复步骤4)至步骤6)，直至遍历N个方位孔径，然后将N个子图像进行相干叠加完成最终成像处理。
综上，设天线相位方向图相位特性因子为exp[jω(α,β)]，其中α和β分别对应点目标P在当前方位孔径下的方位角与俯仰角，那么理想点目标回波信号r(t)则变为：
r(t)＝exp[jω(α,β)]exp[jωc(t-τ)]；
由于采样回波叠加的天线相位方向图相位特性因子是关于方位角α和俯仰角β的空变函数，而非关于时间t的时变函数，因此天线相位方向图相位补偿可考虑在二维时域成像过程中进行补偿。不难发现，通过对天线相位方向图相位一致性的补偿，能够增大方位积累角，改善SAR方位分辨率。
图3为本发明示例仿真所使用的仿真天线相位方向图。图4-图9给出了利用本发明实现的点目标仿真成像结果。图4和图5分别为无天线相位方向图影响的点目标二维成像结果和一维剖面图，即结果目标图。图6和图7分别为增加天线相位方向图影响的点目标二维成像结果和一维剖面图，即本领域传统的简单叠加天线相位方向图的结果，从结果可以看出，天线相位方向图的相位不一致性影响了目标的方位向聚焦。图8和图9分别为利用本发明得到的点目标二维成像结果和一维剖面图。可以看出，本发明能够有效补偿天线相位方向图的相位不一致性，实现目标的精确聚焦。
表2给出了三种仿真条件下的成像分辨率分析结果，对比可知，天线相位方向图相位不一致性会造成二维分辨率的降低，而通过本发明则可以获取与理论分辨率基本一致的成像效果。
表2成像分辨率对比表
Figure BDA0002632823000000101
相应的，本发明还提供一种使用所述的基于天线方向图补偿的直线合成孔径雷达成像方法的直线合成孔径雷达成像系统，直线合成孔径雷达装设在可移动装置上，所述可移动装置用于按照预定的直线合成孔径移动。
作为优选方案，本实施例中，所述可移动装置按照预定的直线合成孔径直线移动。
作为优选方案，本实施例中，所述步骤S2中，多个预定的直线合成孔径按照顺序呈一条直线排列。
作为优选方案，本实施例中，在所述步骤S2中，所述直线合成孔径的数量根据方位采样间隔和雷达扫描探测范围确定。
作为优选方案，本实施例中，在所述步骤S3中，在计算所述成像像素点的俯仰角时，叠加当前雷达天线探测俯视角进行计算。
具体的，在所述成像系统实施过程中，在应用所述基于天线方向图补偿的直线合成孔径雷达成像方法的流程为成像过程中实时补偿，具体请参阅上述实施例。
可以理解的是，对本领域普通技术人员来说，可以根据本发明的技术方案及其发明构思加以等同替换或改变，而所有这些改变或替换都应属于本发明所附的权利要求的保护范围。

Claims (10)

1.一种基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，包括步骤：
S1、对直线合成孔径雷达天线方向图进行精确测量，获得雷达天线相位特性矩阵；
S2、对雷达探测区域进行成像网络剖分，形成若干成像像素点；确定直线合成孔径雷达的所有直线合成孔径位置；
S3、获取所述直线合成孔径雷达在单一所述直线合成孔径上波束照射范围内的所有所述成像像素点，计算每个所述成像像素点与所述直线合成孔径之间的距离延迟、方位角与俯仰角；
S4、根据步骤S3中每个所述成像像素点的所述方位角和所述俯仰角结合所述雷达天线相位特性矩阵计算得到每个像素点的相位补偿因子；多个所述相位补偿因子形成相位补偿矩阵；
S5、所述直线合成孔径雷达对每个所述成像像素点的采样回波进行脉冲压缩处理得到所述一维距离像；
S6、根据所述距离延迟将所述一维距离像中的复数值投影至相应的所述成像像素点，同时对所述成像像素点进行多普勒相位和剩余视频相位的空变补偿，得到所述成像像素点的补偿图像；
S7、将得到所述补偿图像与所述相位补偿矩阵进行复共轭相乘，得到所述直线合成孔径的所述子图像；
S8、将步骤S2中确定的所有所述直线合成孔径均执行步骤S3-S7，将各个所述直线合成孔径得到的所述子图像进行相干叠加得到最终成像。
2.根据权利要求1所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，所述步骤S5中，所述一维距离像的获取步骤为：
S51、所述直线合成孔径雷达对采样回波进行混频和幅度归一化后，得到差频信号；
S52、进而对所述差频信号进行傅立叶变换后得到采样回波频谱；
S53、再根据频率与目标距离的对应关系，得到所述采样回波的一维距离像。
3.根据权利要求2所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，所述步骤S51中，所述差频信号的求解公式为：
所述采样回波信号表达式为：
Figure FDA0002632822990000021
其中，t为时刻；R为时刻t直线合成孔径雷达到目标的瞬时距离；c为光速；TP为信号脉冲；Kr表示调频斜率；α为方位角；β为俯仰角；ω(α,β)为天线相位因子，即在方位角α和俯仰角β方向上天线的相位值；
所述差频信号的表达式为：
Figure FDA0002632822990000022
其中，fc为发射信号的载频。
4.根据权利要求2所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，所述步骤S52中，所述采样回波频谱的表达式为：
Figure FDA0002632822990000023
其中，f为傅立叶频率。
5.根据权利要求2所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，所述步骤S53中，所述频率与目标距离的对应关系表达是为：
f＝2rKr/c；
其中，f为频率；Kr表示调频斜率；c为光速；r为目标距离。
6.根据权利要求5所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，所述一维距离像的表达式为：
Figure FDA0002632822990000024
其中，r为目标距离；Kr表示调频斜率；c为光速；ω(α,β)为天线相位因子；R为时刻t直线合成孔径雷达到目标的瞬时距离。
7.根据权利要求1所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，在所述步骤S2中，所述直线合成孔径的数量根据方位采样间隔和雷达扫描探测范围确定。
8.根据权利要求1所述的基于天线方向图补偿的直线合成孔径雷达成像方法，其特征在于，在所述步骤S3中，在计算所述成像像素点的俯仰角时，叠加当前雷达天线探测俯视角进行计算。
9.一种使用权利要求1-8任一所述的基于天线方向图补偿的直线合成孔径雷达成像方法的直线合成孔径雷达成像系统，其特征在于，直线合成孔径雷达装设在可移动装置上，所述可移动装置用于按照预定的直线合成孔径移动。
10.根据权利要求9所述的直线合成孔径雷达成像系统，其特征在于，所述可移动装置按照预定的直线合成孔径直线移动。
CN202010816258.4A 2020-08-14 2020-08-14 基于天线方向图补偿的直线合成孔径雷达成像方法和系统 Active CN112034459B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202010816258.4A CN112034459B (zh) 2020-08-14 2020-08-14 基于天线方向图补偿的直线合成孔径雷达成像方法和系统

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202010816258.4A CN112034459B (zh) 2020-08-14 2020-08-14 基于天线方向图补偿的直线合成孔径雷达成像方法和系统

Publications (2)

Publication Number Publication Date
CN112034459A true CN112034459A (zh) 2020-12-04
CN112034459B CN112034459B (zh) 2024-05-28

Family

ID=73578212

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202010816258.4A Active CN112034459B (zh) 2020-08-14 2020-08-14 基于天线方向图补偿的直线合成孔径雷达成像方法和系统

Country Status (1)

Country Link
CN (1) CN112034459B (zh)

Cited By (4)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112799064A (zh) * 2020-12-30 2021-05-14 内蒙古工业大学 柱面孔径非线性渐进式相位迭代成像的方法及装置
CN113156436A (zh) * 2021-03-26 2021-07-23 中国科学院空天信息创新研究院 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备
CN115508803A (zh) * 2022-11-23 2022-12-23 深圳市中科海信科技有限公司 基于dsp数字信号处理板的波束锐化处理方法及系统
CN116520323A (zh) * 2023-06-29 2023-08-01 中国科学院空天信息创新研究院 月基合成孔径雷达的对地观测方法及装置

Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4134113A (en) * 1977-04-18 1979-01-09 Westinghouse Electric Corporation Monopulse motion compensation for a synthetic aperture radar
JP2001242244A (ja) * 2000-02-25 2001-09-07 Mitsubishi Electric Corp 合成開口レーダ装置及びその像再生方法
US20070001896A1 (en) * 2003-04-25 2007-01-04 Hans Hellsten Rfi suppression in sar
CN101975948A (zh) * 2010-10-28 2011-02-16 电子科技大学 一种遥感卫星照射源前视合成孔径雷达成像方法
KR20110113926A (ko) * 2010-04-12 2011-10-19 한국과학기술원 연속형 광역 고해상도 영상 구현을 위한 합성 개구면 레이더 시스템 및 그 방법
CN103454633A (zh) * 2013-07-12 2013-12-18 电子科技大学 一种基于后向投影算法的干涉sar动基线处理方法
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
CN108387900A (zh) * 2018-05-09 2018-08-10 西安电子科技大学 一种直升机载旋转式合成孔径雷达的振动误差补偿方法
CN108872986A (zh) * 2018-04-19 2018-11-23 西安电子科技大学 一种方位角波数均匀化处理的极坐标sar运动补偿成像方法
CN110554385A (zh) * 2019-07-02 2019-12-10 中国航空工业集团公司雷华电子技术研究所 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统

Patent Citations (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US4134113A (en) * 1977-04-18 1979-01-09 Westinghouse Electric Corporation Monopulse motion compensation for a synthetic aperture radar
JP2001242244A (ja) * 2000-02-25 2001-09-07 Mitsubishi Electric Corp 合成開口レーダ装置及びその像再生方法
US20070001896A1 (en) * 2003-04-25 2007-01-04 Hans Hellsten Rfi suppression in sar
KR20110113926A (ko) * 2010-04-12 2011-10-19 한국과학기술원 연속형 광역 고해상도 영상 구현을 위한 합성 개구면 레이더 시스템 및 그 방법
CN101975948A (zh) * 2010-10-28 2011-02-16 电子科技大学 一种遥感卫星照射源前视合成孔径雷达成像方法
CN103454633A (zh) * 2013-07-12 2013-12-18 电子科技大学 一种基于后向投影算法的干涉sar动基线处理方法
CN103901429A (zh) * 2014-04-09 2014-07-02 西安电子科技大学 基于稀疏孔径的机动目标逆合成孔径雷达成像方法
CN108872986A (zh) * 2018-04-19 2018-11-23 西安电子科技大学 一种方位角波数均匀化处理的极坐标sar运动补偿成像方法
CN108387900A (zh) * 2018-05-09 2018-08-10 西安电子科技大学 一种直升机载旋转式合成孔径雷达的振动误差补偿方法
CN110554385A (zh) * 2019-07-02 2019-12-10 中国航空工业集团公司雷华电子技术研究所 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
吴星辉 等: "地基合成孔径雷达的发展现状及应用", 武汉大学学报, vol. 44, no. 7, pages 1073 - 1081 *
张群英 等: "太赫兹合成孔径雷达成像运动补偿算法", 电子与信息学报, vol. 89, no. 1, pages 179 - 187 *

Cited By (7)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN112799064A (zh) * 2020-12-30 2021-05-14 内蒙古工业大学 柱面孔径非线性渐进式相位迭代成像的方法及装置
CN113156436A (zh) * 2021-03-26 2021-07-23 中国科学院空天信息创新研究院 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备
CN113156436B (zh) * 2021-03-26 2023-04-25 中国科学院空天信息创新研究院 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备
CN115508803A (zh) * 2022-11-23 2022-12-23 深圳市中科海信科技有限公司 基于dsp数字信号处理板的波束锐化处理方法及系统
CN115508803B (zh) * 2022-11-23 2023-02-03 深圳市中科海信科技有限公司 基于dsp数字信号处理板的波束锐化处理方法及系统
CN116520323A (zh) * 2023-06-29 2023-08-01 中国科学院空天信息创新研究院 月基合成孔径雷达的对地观测方法及装置
CN116520323B (zh) * 2023-06-29 2024-04-09 中国科学院空天信息创新研究院 月基合成孔径雷达的对地观测方法及装置

Also Published As

Publication number Publication date
CN112034459B (zh) 2024-05-28

Similar Documents

Publication Publication Date Title
CN112034459B (zh) 基于天线方向图补偿的直线合成孔径雷达成像方法和系统
CN107741586B (zh) 基于DBF-TOPS加权的星载Ka InSAR信号处理方法
JPH045157B2 (zh)
CN108490443B (zh) 基于解析解及NUFFT的多子阵合成孔径声纳ωk成像算法
CN110146889B (zh) 基于最优回波子区域选择的大转角isar成像方法
CN102608597A (zh) 基于非完全数据解卷积的实孔径前视成像方法
Fornaro et al. Role of processing geometry in SAR raw data focusing
CN113238225B (zh) 一种用于动目标检测的雷达
CN114545411B (zh) 一种基于工程实现的极坐标格式多模高分辨sar成像方法
Goh et al. Bistatic SAR experiment with the Ingara imaging radar
CN112305510A (zh) 一种基于dem匹配的合成孔径雷达影像几何定标方法
CN108226925A (zh) 一种适用于弹载大前斜视时变参数sar的重叠子孔径算法
CN111257878A (zh) 基于俯仰维频率脉内扫描高分宽幅sar的波形设计方法
CN116794612A (zh) 基于线性变标的合成孔径雷达距离-孔径空变运动补偿方法
CN108732555B (zh) 一种自动驾驶阵列微波成像运动补偿的方法
CN111880179A (zh) 一种弹载弧线俯冲大斜视tops sar的成像方法
CN107607950B (zh) 基于DBF-TOPS加权的星载Ka InSAR波束形成方法
RU2614041C1 (ru) Способ формирования изображения земной поверхности в радиолокационной станции с синтезированием апертуры антенны
CN116148891A (zh) 基于顺轨干涉的星地双基海面二维流场测量方法
Raj et al. Velocity-ISAR: On the application of ISAR techniques to multichannel SAR imaging
RU2707556C1 (ru) Способ определения высоты рельефа местности радиолокатором с синтезированной апертурой антенны
Amin et al. Mixing Matrix Calibration for BSS Range Ambiguity Suppression in Multichannel SAR Systems
Freeman et al. An application of the monopulse principle to determining elevation angles in SAR images
Chen et al. Very High-Resolution Synthetic Aperture Radar Systems and Imaging: A Review
RU2626012C1 (ru) Способ картографирования земной поверхности бортовой радиолокационной станцией

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
TA01 Transfer of patent application right
TA01 Transfer of patent application right

Effective date of registration: 20240410

Address after: 410000, Block 1-10, B1E1, Beichen Fenghuang Tianjie Garden, No. 68 Qinglan Road, Xinhe Street, Kaifu District, Changsha City, Hunan Province, China, along with an underground garage and connecting platform 8005

Applicant after: Hunan Jianwei Intelligent Technology Co.,Ltd.

Country or region after: China

Address before: Room 307, building 4, Hunan scientific and technological achievements transformation center, Longping high tech park, Furong district, Changsha City, Hunan Province

Applicant before: Song Qian

Country or region before: China

Applicant before: Wang Pengyu

GR01 Patent grant