CN113156436A - 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备 - Google Patents

圆迹合成孔径雷达自聚焦成像方法、系统及电子设备 Download PDF

Info

Publication number
CN113156436A
CN113156436A CN202110327795.7A CN202110327795A CN113156436A CN 113156436 A CN113156436 A CN 113156436A CN 202110327795 A CN202110327795 A CN 202110327795A CN 113156436 A CN113156436 A CN 113156436A
Authority
CN
China
Prior art keywords
image
projection matrix
synthetic aperture
imaging
phase error
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
CN202110327795.7A
Other languages
English (en)
Other versions
CN113156436B (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.)
Aerospace Information Research Institute of CAS
Original Assignee
Aerospace Information Research Institute of CAS
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 Aerospace Information Research Institute of CAS filed Critical Aerospace Information Research Institute of CAS
Priority to CN202110327795.7A priority Critical patent/CN113156436B/zh
Publication of CN113156436A publication Critical patent/CN113156436A/zh
Application granted granted Critical
Publication of CN113156436B publication Critical patent/CN113156436B/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/9052—Spotlight 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
    • 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

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

本发明的实施例提供了一种圆迹合成孔径雷达自聚焦成像方法、系统及电子设备。该方法包括：将成像区域划分为多个矩形区域；对多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；对投影矩阵进行图像偏移量补偿；对图像偏移量补偿后的投影矩阵进行相位误差估计和相位误差补偿；并且基于相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。本发明通过同时补偿包络误差和相位误差，提高了SAR成像的精度，同时通过对成像区域分块处理，减少了相位误差空变性对SAR成像的影响并且降低了SAR成像对存储空间的要求。

Description

圆迹合成孔径雷达自聚焦成像方法、系统及电子设备
技术领域
本发明涉及雷达成像技术领域，尤其涉及一种圆迹合成孔径雷达自聚焦成像方法、系统及电子设备。
背景技术
合成孔径雷达(Synthetic Aperture Radar，SAR)是一种微波主动遥感设备，具有全天时全天候的工作特点，因此在军事和民用各领域里均有广泛的应用。在传统的SAR工作模式中，雷达平台均沿直线飞行且合成孔径长度十分有限，这就造成了叠掩、阴影、透视缩短等SAR图像普遍存在的问题。为了解决这些问题，提出了圆迹SAR工作模式。
一般使用后向投影算法(Back Projection，BP)对圆迹SAR数据进行成像处理，而BP算法对地形和雷达平台的运动轨迹的精度有很高的要求。目前机载GPS和惯导系统与开源DEM数据很难自聚焦满足高波段圆迹SAR成像要求，因此需要使用自聚焦算法成像。
现有的自聚焦BP成像算法中，只针对运动误差引起的相位误差进行了估计和补偿，而忽略了对包络误差的补偿，这会导致补偿不够精确，此外，相位误差空变性也会对自聚焦BP效果产生影响，这使得SAR图像聚焦效果不够理想。同时，现有SAR图像聚焦过程产生的开销较大，对计算机内存要求比较高。
发明内容
针对现有技术中的问题，本发明的实施例提供了一种圆迹合成孔径雷达自聚焦成像方法、系统及电子设备。
根据本发明的一方面，提供了一种圆迹合成孔径雷达自聚焦成像方法，包括：
将成像区域划分为多个矩形区域；
对多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；
对投影矩阵进行图像偏移量补偿；
对图像偏移量补偿后的投影矩阵进行相位误差估计和相位误差补偿；并且
基于相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。
进一步地，对投影矩阵进行图像偏移量补偿包括：
将合成孔径划分为多个子孔径，并获取多个子孔径中的每一个的图像；
选择图像中的一个作为参考图像，计算参考图像与图像中的其余图像之间的归一化互相关系数，并获取归一化互相关系数最大时的二维偏移量；并且
基于二维偏移量对投影矩阵进行图像偏移量补偿。
进一步地，将合成孔径划分为多个子孔径包括：
根据雷达信噪比和分辨率将合成孔径划分为多个子孔径。
进一步地，通过以下公式获取多个子孔径中的每一个的图像：
Figure BDA0002995275920000021
其中，I(i，j)为像素点(i，j)的像素值，sij为像素点(i，j)对应的距离向脉冲压缩后的回波信号，fc为雷达中心频率，c为光速，Rij为像素点(i，j)到雷达平台的距离。
进一步地，通过以下公式计算Rij：
Figure BDA0002995275920000022
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置。
进一步地，通过以下公式对投影矩阵进行图像偏移量补偿：
Figure BDA0002995275920000031
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置，Δx，Δy为归一化互相关系数最大时的二维偏移量。
进一步地，通过以下公式计算参考图像与图像中的其余图像之间的归一化互相关系数：
Figure BDA0002995275920000032
其中，Gh为参考图像与第h个子孔径的图像之间的归一化互相关系数，p为横向的像素点数量，q为纵向的像素点数量，i为像素点的横坐标位置，j为像素点的纵坐标位置，I为图像矩阵，mean()代表取均值操作。
根据本发明的另一方面，提供了一种圆迹合成孔径雷达自聚焦成像系统，包括：
成像区域划分模块，用于将成像区域划分为多个矩形区域；
投影矩阵计算模块，用于对多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；
包络误差补偿模块，用于对投影矩阵进行图像偏移量补偿；
相位误差补偿模块，用于对图像偏移量补偿后的投影矩阵进行相位误差估计和相位误差补偿；以及
成像区域重组模块，用于基于相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。
进一步地，包络误差补偿模块包括：
成像模块，用于将合成孔径划分为多个子孔径，并获取多个子孔径中的每一个的图像；
包络误差计算模块，用于选择图像中的一个作为参考图像，计算参考图像与图像中的其余图像之间的归一化互相关系数，并获取归一化互相关系数最大时的二维偏移量；以及
包络误差补偿执行模块，用于基于二维偏移量对投影矩阵进行图像偏移量补偿。
根据本发明的另一方面，提供了一种电子设备，包括存储器、处理器及存储在存储器上并可在处理器上运行的计算机程序，处理器执行计算机程序时实现上述任一项圆迹合成孔径雷达自聚焦成像方法的步骤。
由以上技术方案可知，本发明的技术方案通过同时补偿包络误差和相位误差，提高了SAR成像的精度，同时通过对成像区域分块处理，减少了相位误差空变性对SAR成像的影响并且降低了SAR成像对存储空间的要求。
附图说明
为了更清楚地说明本发明的实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。
图1为本发明的实施例提供的SAR成像示意图；
图2为本发明的实施例提供的圆迹合成孔径雷达自聚焦成像方法的流程图；
图3为图2所示的包络误差补偿过程的示意图；
图4为本发明的实施例提供的圆迹合成孔径雷达自聚焦成像系统的示意图；以及
图5为本发明的实施例提供的电子设备的示意图。
具体实施方式
为使本发明的实施例的目的、技术方案和优点更加清楚，下面将结合本发明的实施例中的附图，对本发明的实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。
图1示出了本发明的实施例提供的SAR成像示意图，圆迹SAR以一个小天线作为单个辐射单元，在不同位置上接收同一地物的回波信号并进行相关解调压缩处理。一个小天线(例如，孔径1、孔径2、孔径3、......)通过“运动”方式就合成一个等效“大天线”，这样SAR就可以安装在卫星平台上而可以获取较高分辨率的SAR图像。
图2示出了本发明的实施例提供的圆迹合成孔径雷达自聚焦成像方法的流程图。
参考图2，根据本发明的实施例的圆迹合成孔径雷达自聚焦成像方法包括以下步骤：
S201：将成像区域划分为多个矩形区域；
S203：对所述多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；
S205：对所述投影矩阵进行图像偏移量补偿；
S207：对图像偏移量补偿后的所述投影矩阵进行相位误差估计和相位误差补偿；并且
S209：基于相位误差补偿后的所述投影矩阵，拼接得到所述成像区域的全场景图像。
在本实施例中，需要说明的是，可以根据计算机的存储性能来对成像区域进行划分。通常，在存储性能满足需求的情况下，划分后的区域的面积越大，相位估计的精度越高，但是，面积太大也可能会导致相位误差空变性影响更大。为了保证相位估计的精确性，可以选择C波段雷达每块矩形区域边长小于500米，X波段雷达每块区域边长小于300米，Ku波段边长小于100米，Ka波段小于50米。
为了便于说明，下面假设将成像区域划分为M*N个矩形区域。
在S203中，在通过S201获得M*N个矩形区域之后，对M*N个矩形区域中的每一个，针对每个脉冲计算投影矩阵。
在本实施例中，需要说明的是，可以通过传统的BP算法针对每个脉冲计算投影矩阵，然而，本发明不限于此。
在S205中，对投影矩阵进行图像偏移量补偿。
在本实施例中，需要说明的是，对投影矩阵进行图像偏移量补偿是指进行包络误差补偿。包络误差补偿的操作将结合图3进行详细描述。
在S207中，在通过S205对投影矩阵进行图像偏移量补偿之后，对投影矩阵进行相位误差估计和相位误差补偿。
在本实施例中，需要说明的是，可以通过传统的自聚焦BP算法实现相位误差估计和相位误差补偿，然而，本发明不限于此。
在S209中，基于通过S207进行相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。
下面将参考图3详细描述图2所示的包络误差补偿S205的具体步骤。如图3所示，对投影矩阵进行偏移量补偿可以包括以下步骤：
S2051：将合成孔径划分为多个子孔径，并获取多个子孔径中的每一个的图像；
S2053：选择所述图像中的一个作为参考图像，计算所述参考图像与所述图像中的其余图像之间的归一化互相关系数，并获取所述归一化互相关系数最大时的二维偏移量；并且
S2055：基于所述二维偏移量对所述投影矩阵进行图像偏移量补偿。
在本实施例中，需要说明的是，可以根据雷达信噪比、距离向分辨率、方位向分辨率等因素来对子孔径进行分割。为了保证子孔径分割之后能够进行后续的步骤并得到准确的偏移量，分割后的单幅子孔径中必须能看到场景信息。
在本实施例中，在分割获得多个子孔径之后，可以通过以下BP成像公式获取多个子孔径中的每一个的图像：
Figure BDA0002995275920000071
其中，I(i，j)为像素点(i，j)的像素值，sij为像素点(i，j)对应的距离向脉冲压缩后的回波信号，fc为雷达中心频率，c为光速，Rij为像素点(i，j)到雷达平台的距离。
在本实施例中，可以通过以下公式计算Rij：
Figure BDA0002995275920000072
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置。
为了便于说明，下面假设将合成孔径分为k个子孔径，其中k为大于1的自然数。
在S2053中，在通过S2051获得k个子孔径的图像之后，选择k个图像中的一个作为参考图像，计算参考图像与k个图像中的其余图像之间的归一化互相关系数，并获取归一化互相关系数最大时的二维偏移量。
在本实施例中，需要说明的是，可以通过以下公式计算参考图像与图像中的其余图像之间的归一化互相关系数：
Figure BDA0002995275920000081
其中，h为子孔径序号，h＝1、2、......、k，Gh为参考图像与第h个子孔径的图像之间的归一化互相关系数，p为横向的像素点数量，q为纵向的像素点数量，i为像素点的横坐标位置，j为像素点的纵坐标位置，I为图像矩阵，mean()代表取均值操作。
在本实施例中，需要说明的是，在计算出每个子孔径的图像与参考图像之间的归一化互相关系数之后，可以得到归一化相关系数最大时的二维偏移量。例如，当第h个子孔径的图像被计算为具有最大的归一化相关系数时，该二维偏移量可以被表示为Δxh,Δyh。。
为了便于说明，下面将假设Δxh,Δyh为归一化相关系数最大时的二维偏移量。
在S2055中，在通过S2053获得归一化相关系数最大时的二维偏移量之后，基于该二维偏移量对投影矩阵进行图像偏移量补偿。
在本实施例中，可以通过以下公式对投影矩阵进行图像偏移量补偿：
Figure BDA0002995275920000082
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置，Δx，Δy为归一化互相关系数最大时的二维偏移量，即Δxh,Δyh。
由以上可知，本实施例通过同时补偿包络误差和相位误差，可以提高SAR成像的精度，同时通过对成像区域分块处理，可以减少相位误差空变性对SAR成像的影响并且降低SAR成像对存储空间的要求。
图4示出了本发明的实施例的圆迹合成孔径雷达自聚焦成像系统的结构示意图。
参考图4，圆迹合成孔径雷达自聚焦成像系统400包括成像区域划分模块401、投影矩阵计算模块403、包络误差补偿模块405、相位误差补偿模块407和成像区域重组模块409。
成像区域划分模块401用于将成像区域划分为多个矩形区域。
在本实施例中，成像区域划分模块401可以根据计算机的存储性能来对成像区域进行划分。通常，在存储性能满足需求的情况下，划分后的区域的面积越大，相位估计的精度越高，但是，面积太大也可能会导致相位误差空变性影响更大。为了保证相位估计的精确性，可以选择C波段雷达每块矩形区域边长小于500米，X波段雷达每块区域边长小于300米，Ku波段边长小于100米，Ka波段小于50米。
投影矩阵计算模块403用于对多个矩形区域中的每一个，针对每个脉冲计算投影矩阵。
在本实施例中，投影矩阵计算模块403可以通过传统的BP算法针对每个脉冲计算投影矩阵，然而，本发明不限于此。
包络误差补偿模块405用于对投影矩阵进行图像偏移量补偿。
稍后将详细描述包络误差补偿模块405的详细配置。
相位误差补偿模块407用于对图像偏移量补偿后的投影矩阵进行相位误差估计和相位误差补偿。
在本实施例中，相位误差补偿模块407可以通过传统的自聚焦BP算法实现相位误差估计和相位误差补偿，然而，本发明不限于此。
成像区域重组模块409用于基于相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。
根据实施例，包络误差补偿模块405可以被配置为包括成像模块4051、包络误差计算模块4053和包络误差补偿执行模块4055。
成像模块4051用于将合成孔径划分为多个子孔径，并获取多个子孔径中的每一个的图像。
在本实施例中，成像模块4051可以根据雷达信噪比、距离向分辨率、方位向分辨率等因素来对子孔径进行分割。为了保证子孔径分割之后能够进行后续的步骤并得到准确的偏移量，分割后的单幅子孔径中必须能看到场景信息。
在本实施例中，在分割获得多个子孔径之后，成像模块4051可以通过以下BP成像公式获取多个子孔径中的每一个的图像：
Figure BDA0002995275920000101
其中，I(i，j)为像素点(i，j)的像素值，sij为像素点(i，j)对应的距离向脉冲压缩后的回波信号，fc为雷达中心频率，c为光速，Rij为像素点(i，j)到雷达平台的距离。
在本实施例中，成像模块4051可以通过以下公式计算Rij：
Figure BDA0002995275920000102
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置。
为了便于说明，下面假设成像模块4051将合成孔径分为k个子孔径，其中k为大于1的自然数。
包络误差计算模块4053用于选择图像中的一个作为参考图像，计算参考图像与图像中的其余图像之间的归一化互相关系数，并获取归一化互相关系数最大时的二维偏移量。
在本实施例中，包络误差计算模块4053可以选择k个图像中的一个作为参考图像，计算参考图像与k个图像中的其余图像之间的归一化互相关系数，并获取归一化互相关系数最大时的二维偏移量。
在本实施例中，包络误差计算模块4053可以通过以下公式计算参考图像与图像中的其余图像之间的归一化互相关系数：
Figure BDA0002995275920000111
其中，h为子孔径序号，h＝1、2、......、k，Gh为参考图像与第h个子孔径的图像之间的归一化互相关系数，p为横向的像素点数量，q为纵向的像素点数量，i为像素点的横坐标位置，j为像素点的纵坐标位置，I为图像矩阵，mean()代表取均值操作。
在本实施例中，在计算出每个子孔径的图像与参考图像之间的归一化互相关系数之后，可以得到归一化相关系数最大时的二维偏移量。例如，当第h个子孔径的图像被计算为具有最大的归一化相关系数时，该二维偏移量可以被表示为Δxh,Δyh。。
为了便于说明，下面将假设Δxh,Δyh为归一化相关系数最大时的二维偏移量。
包络误差补偿执行模块4055用于基于二维偏移量对投影矩阵进行图像偏移量补偿。
在本实施例中，包络误差补偿执行模块4055可以通过以下公式对投影矩阵进行图像偏移量补偿：
Figure BDA0002995275920000112
其中，xa，ya，za为雷达平台的空间坐标位置，xp，yp，zp为像素点(i，j)的空间坐标位置，Δx，Δy为归一化互相关系数最大时的二维偏移量，即Δxh,Δyh。
由以上可知，本实施例通过同时补偿包络误差和相位误差，可以提高SAR成像的精度，同时通过对成像区域分块处理，可以减少相位误差空变性对SAR成像的影响并且降低SAR成像对存储空间的要求。
在本实施例中，本发明的系统中的各个模块可以集成于一体，也可以分离部署。上述模块可以合并为一个模块，也可以进一步拆分成多个子模块。
图5示出了本发明的实施例提供的电子设备的示意图。
参考图5，电子设备500包括处理器501、存储器503、通信接口505和通信总线507，并且处理器501、存储器503、通信接口505通过通信总线507在彼此之间进行通信。
处理器501用于调用存储器503中存储的计算机程序，处理器501执行计算机程序时实现上述圆迹合成孔径雷达自聚焦成像方法的全部步骤，例如，处理器501执行计算机程序时实现下述过程：将成像区域划分为多个矩形区域；对多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；对投影矩阵进行图像偏移量补偿；对图像偏移量补偿后的投影矩阵进行相位误差估计和相位误差补偿；并且基于相位误差补偿后的投影矩阵，拼接得到成像区域的全场景图像。
可以理解的是，计算机程序可以执行的细化功能和扩展功能可参考上面实施例的描述。
此外，上述存储器中的逻辑指令可以通过软件功能单元的形式实现并作为独立的产品销售或使用时，可以存储在一个计算机可读取存储介质中。基于这样的理解，本发明的技术方案本质上或者说对现有技术做出贡献的部分或者该技术方案的部分可以以软件产品的形式体现出来，该计算机软件产品存储在一个存储介质中，包括若干指令以使得计算机设备(可以是个人计算机、服务器或者网络设备等)执行本发明各个实施例所述方法的全部或部分步骤，并且上述存储介质包括不但限于U盘、移动硬盘、只读存储器(Read-OnlyMemory，ROM)、随机存取存储器(Random Access Memory，RAM)、磁碟或者光盘等各种可以存储程序代码的介质。
以上所描述的装置的实施例仅仅是示意性的，其中所述作为分离部件说明的单元可以是或者也可以不是物理上分开的，作为单元显示的部件可以是或者也可以不是物理单元，即可以位于一个地方，或者也可以分布到多个网络单元上。可以根据实际的需要选择其中的部分或者全部模块来实现本发明的实施例方案的目的。本领域普通技术人员在不付出创造性的劳动的情况下，即可以理解并实施。
通过对以上实施方案的描述，本领域的技术人员可以清楚地了解到各实施方案可借助软件加必需的通用硬件平台的方式来实现，当然也可以通过硬件。基于这样的理解，上述技术方案本质上或者说对现有技术做出贡献的部分可以以软件产品的形式体现出来，该计算机软件产品可以存储在计算机可读存储介质中，如ROM/RAM、磁碟、光盘等，包括若干指令用以使得一台计算机设备(可以是个人计算机，服务器，或者网络设备等)执行各个实施例或者实施例的某些部分所述的圆迹合成孔径雷达自聚焦成像方法。
此外，在本发明中，参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。在本说明书中，对上述术语的示意性表述不必须针对的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外，在不相互矛盾的情况下，本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
最后应说明的是：以上实施例仅用以说明本发明的技术方案，而非对其限制；尽管参考前述实施例对本发明进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本发明各实施例技术方案的精神和范围。

Claims (10)

1.一种圆迹合成孔径雷达自聚焦成像方法，其特征在于，包括：
将成像区域划分为多个矩形区域；
对所述多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；
对所述投影矩阵进行图像偏移量补偿；
对图像偏移量补偿后的所述投影矩阵进行相位误差估计和相位误差补偿；并且
基于相位误差补偿后的所述投影矩阵，拼接得到所述成像区域的全场景图像。
2.根据权利要求1所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
所述对所述投影矩阵进行图像偏移量补偿包括：
将合成孔径划分为多个子孔径，并获取所述多个子孔径中的每一个的图像；
选择所述图像中的一个作为参考图像，计算所述参考图像与所述图像中的其余图像之间的归一化互相关系数，并获取所述归一化互相关系数最大时的二维偏移量；并且
基于所述二维偏移量对所述投影矩阵进行图像偏移量补偿。
3.根据权利要求2所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
所述将合成孔径划分为多个子孔径包括：
根据雷达信噪比和分辨率将所述合成孔径划分为所述多个子孔径。
4.根据权利要求3所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
通过以下公式获取所述多个子孔径中的每一个的图像：
Figure FDA0002995275910000021
其中，I(i，j)为像素点(i，j)的像素值，sij为所述像素点(i，j)对应的距离向脉冲压缩后的回波信号，fc为雷达中心频率，c为光速，Rij为所述像素点(i，j)到雷达平台的距离。
5.根据权利要求4所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
通过以下公式计算所述Rij：
Figure FDA0002995275910000022
其中，xa，ya，za为所述雷达平台的空间坐标位置，xp，yp，zp为所述像素点(i，j)的空间坐标位置。
6.根据权利要求5所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
通过以下公式对所述投影矩阵进行图像偏移量补偿：
Figure FDA0002995275910000023
其中，xa，ya，za为所述雷达平台的空间坐标位置，xp，yp，zp为所述像素点(i，j)的空间坐标位置，Δx，Δy为所述归一化互相关系数最大时的二维偏移量。
7.根据权利要求2所述的圆迹合成孔径雷达自聚焦成像方法，其特征在于，
通过以下公式计算所述参考图像与所述图像中的其余图像之间的归一化互相关系数：
Figure FDA0002995275910000024
其中，Gh为所述参考图像与第h个子孔径的图像之间的归一化互相关系数，p为横向的像素点数量，q为纵向的像素点数量，i为像素点的横坐标位置，j为像素点的纵坐标位置，I为图像矩阵，mean()代表取均值操作。
8.一种圆迹合成孔径雷达自聚焦成像系统，其特征在于，包括：
成像区域划分模块，用于将成像区域划分为多个矩形区域；
投影矩阵计算模块，用于对所述多个矩形区域中的每一个，针对每个脉冲计算投影矩阵；
包络误差补偿模块，用于对所述投影矩阵进行图像偏移量补偿；
相位误差补偿模块，用于对图像偏移量补偿后的所述投影矩阵进行相位误差估计和相位误差补偿；以及
成像区域重组模块，用于基于相位误差补偿后的所述投影矩阵，拼接得到所述成像区域的全场景图像。
9.根据权利要求8所述的圆迹合成孔径雷达自聚焦成像系统，其特征在于，
所述包络误差补偿模块包括：
成像模块，用于将合成孔径划分为多个子孔径，并获取所述多个子孔径中的每一个的图像；
包络误差计算模块，用于选择所述图像中的一个作为参考图像，计算所述参考图像与所述图像中的其余图像之间的归一化互相关系数，并获取所述归一化互相关系数最大时的二维偏移量；以及
包络误差补偿执行模块，用于基于所述二维偏移量对所述投影矩阵进行图像偏移量补偿。
10.一种电子设备，其特征在于，包括存储器、处理器及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现权利要求1至7中任一项所述圆迹合成孔径雷达自聚焦成像方法的步骤。
CN202110327795.7A 2021-03-26 2021-03-26 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备 Active CN113156436B (zh)

Priority Applications (1)

Application Number Priority Date Filing Date Title
CN202110327795.7A CN113156436B (zh) 2021-03-26 2021-03-26 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
CN202110327795.7A CN113156436B (zh) 2021-03-26 2021-03-26 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备

Publications (2)

Publication Number Publication Date
CN113156436A true CN113156436A (zh) 2021-07-23
CN113156436B CN113156436B (zh) 2023-04-25

Family

ID=76885064

Family Applications (1)

Application Number Title Priority Date Filing Date
CN202110327795.7A Active CN113156436B (zh) 2021-03-26 2021-03-26 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备

Country Status (1)

Country Link
CN (1) CN113156436B (zh)

Cited By (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114519778A (zh) * 2022-03-02 2022-05-20 中国科学院空天信息创新研究院 多角度sar数据的目标三维重构方法、装置、设备及介质
CN115185078A (zh) * 2022-05-30 2022-10-14 元潼(北京)技术有限公司 非相干孔径合成像差校正方法及装置

Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060109164A1 (en) * 2004-11-23 2006-05-25 Cho Kwang M Efficient autofocus method for swath SAR
CN104730520A (zh) * 2015-03-27 2015-06-24 电子科技大学 基于子孔径合成的圆周sar后向投影自聚焦方法
CN105974414A (zh) * 2016-06-24 2016-09-28 西安电子科技大学 基于二维自聚焦的高分辨聚束sar自聚焦成像方法
CN109917389A (zh) * 2019-04-16 2019-06-21 中国人民解放军国防科技大学 一种机载全息sar成像中的相位校正方法
CN110554385A (zh) * 2019-07-02 2019-12-10 中国航空工业集团公司雷华电子技术研究所 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统
CN111208513A (zh) * 2020-01-15 2020-05-29 西安电子科技大学 空间目标isar图像序列能量反向投影与三维重构方法
CN111880180A (zh) * 2020-07-03 2020-11-03 西安电子科技大学 一种高分辨运动舰船sar成像的自聚焦方法
CN112034459A (zh) * 2020-08-14 2020-12-04 宋千 基于天线方向图补偿的直线合成孔径雷达成像方法和系统

Patent Citations (8)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US20060109164A1 (en) * 2004-11-23 2006-05-25 Cho Kwang M Efficient autofocus method for swath SAR
CN104730520A (zh) * 2015-03-27 2015-06-24 电子科技大学 基于子孔径合成的圆周sar后向投影自聚焦方法
CN105974414A (zh) * 2016-06-24 2016-09-28 西安电子科技大学 基于二维自聚焦的高分辨聚束sar自聚焦成像方法
CN109917389A (zh) * 2019-04-16 2019-06-21 中国人民解放军国防科技大学 一种机载全息sar成像中的相位校正方法
CN110554385A (zh) * 2019-07-02 2019-12-10 中国航空工业集团公司雷华电子技术研究所 机动轨迹合成孔径雷达自聚焦成像方法、装置及雷达系统
CN111208513A (zh) * 2020-01-15 2020-05-29 西安电子科技大学 空间目标isar图像序列能量反向投影与三维重构方法
CN111880180A (zh) * 2020-07-03 2020-11-03 西安电子科技大学 一种高分辨运动舰船sar成像的自聚焦方法
CN112034459A (zh) * 2020-08-14 2020-12-04 宋千 基于天线方向图补偿的直线合成孔径雷达成像方法和系统

Non-Patent Citations (2)

* Cited by examiner, † Cited by third party
Title
DAIYIN ZHU ET AL: "Multi-Subaperture PGA for SAR Autofocusing", 《IEEE TRANSACTIONS ON AEROSPACE AND ELECTRONIC SYSTEMS》 *
朱岱寅等: "高分辨率极化合成孔径雷达成像研究进展", 《数据采集与处理》 *

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN114519778A (zh) * 2022-03-02 2022-05-20 中国科学院空天信息创新研究院 多角度sar数据的目标三维重构方法、装置、设备及介质
CN114519778B (zh) * 2022-03-02 2022-11-22 中国科学院空天信息创新研究院 多角度sar数据的目标三维重构方法、装置、设备及介质
CN115185078A (zh) * 2022-05-30 2022-10-14 元潼(北京)技术有限公司 非相干孔径合成像差校正方法及装置

Also Published As

Publication number Publication date
CN113156436B (zh) 2023-04-25

Similar Documents

Publication Publication Date Title
Ran et al. An autofocus algorithm for estimating residual trajectory deviations in synthetic aperture radar
CN111443339A (zh) 一种双基sar空变校正成像方法、装置、设备及存储介质
JP6249110B1 (ja) 合成開口レーダ装置
US8212714B1 (en) Using doppler radar images to estimate aircraft navigational heading error
CN113156436A (zh) 圆迹合成孔径雷达自聚焦成像方法、系统及电子设备
CN108090956B (zh) 一种几何校正方法及系统
CN106908792B (zh) 一种合成孔径雷达成像方法及装置
RU2661941C1 (ru) Способ построения радиолокационного изображения с помощью радиолокационной станции с синтезированной апертурой
CN111652931A (zh) 一种地理定位方法、装置、设备及计算机可读存储介质
Linnehan et al. Map-drift autofocus and scene stabilization for video-SAR
Roth et al. Geocoding of TerraSAR-X data
Bezvesilniy et al. Synthetic aperture radar systems for small aircrafts: Data processing approaches
KR102185307B1 (ko) Sar 영상의 객체 응답 초해상도화 방법 및 객체 응답 초해상도화 장치
CN116609781A (zh) 多星数据联合的北斗InSAR DEM误差补偿方法
Guo et al. Focusing multistatic geo sar with two stationary receivers using spectrum alignment and extrapolation
Li et al. Ground moving target detection and location based on SAR images for distributed spaceborne SAR
Klein et al. Pulse-to-pulse radar-aided positioning using multibeam autofocus
US20230152443A1 (en) Method for georeferencing of a digital elevation model
CN114187332A (zh) 一种雷达图像配准方法和系统
De Macedo et al. Recent LC-and X-band MetaSensing airborne SAR campaigns for emerging applications
Hu et al. A novel antenna phase center estimation method for synthetic aperture radar
Brotzer et al. Drone With Integrated Moving Baseline System and Time Domain Autofocus Algorithm for High Resolution SAR Images
CN117745779B (zh) 一种光学和sar共孔径一致性成像方法
Zhang et al. High-Precision Spotlight SAR Imaging with Joint Modified Fast Factorized Backprojection and Autofocus
CN115311363B (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