CN109313261A - 用于干涉测量分析的sar成像方法 - Google Patents

用于干涉测量分析的sar成像方法 Download PDF

Info

Publication number
CN109313261A
CN109313261A CN201780037690.XA CN201780037690A CN109313261A CN 109313261 A CN109313261 A CN 109313261A CN 201780037690 A CN201780037690 A CN 201780037690A CN 109313261 A CN109313261 A CN 109313261A
Authority
CN
China
Prior art keywords
sar
acquisition
doppler
radiancy
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.)
Pending
Application number
CN201780037690.XA
Other languages
English (en)
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.)
Thales Alenia Space Italia SpA
Original Assignee
Thales Alenia Space Italia SpA
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 Thales Alenia Space Italia SpA filed Critical Thales Alenia Space Italia SpA
Publication of CN109313261A publication Critical patent/CN109313261A/zh
Pending legal-status Critical Current

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
    • 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/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G01S7/4021—Means for monitoring or calibrating of parts of a radar system of receivers
    • 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/9023—SAR image post-processing techniques combined with interferometric 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/904—SAR modes
    • G01S13/9054—Stripmap 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/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/04—Display arrangements
    • 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/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G01S7/4026—Antenna boresight
    • 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/40—Means for monitoring or calibrating
    • G01S7/4004—Means for monitoring or calibrating of parts of a radar system
    • G01S7/4026—Antenna boresight
    • G01S7/403—Antenna boresight in azimuth, i.e. in the horizontal plane

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

本发明涉及一种用于干涉分析的SAR成像方法(40)，包括：接收与通过一个或多个合成孔径雷达执行的对地球表面的同一区域的两次或更多次SAR采集相关的原始SAR数据；以及处理所述原始SAR数据以生成SAR图像。所述方法的其特征在于，对于每次SAR采集，基于不同的两组处理参数处理与所述SAR采集相关的相应原始SAR数据，所述不同的两组处理参数包括：第一组，所述第一组对于所有SAR采集而言都相同，并且包括基于与所有SAR采集相关的物理多普勒参数计算的聚焦多普勒参数；以及相应的第二组，所述相应的第二组包括与所述SAR采集相关并且基于与所述SAR采集相关的相应物理多普勒参数来计算的相应的辐射度均衡多普勒参数。具体地，处理包括：基于所述聚焦多普勒参数聚焦与所有SAR采集相关的原始SAR数据；以及，对于每次SAR采集，基于相应的辐射度均衡多普勒参数将相应的辐射度均衡应用于相应的SAR数据，以在不降低方位分辨率并且不引入辐射度失真的情况下，补偿在用于执行SAR采集的所述合成孔径雷达的指向方面的可能差异。

Description

用于干涉测量分析的SAR成像方法
技术领域
本发明总体涉及借助于合成孔径雷达(SAR)的遥感，具体涉及基于SAR的干涉测量。
更具体地，本发明涉及一种用于生成用于干涉测量分析的SAR图像的技术。本发明在处理与以突发模式执行的两次或更多次SAR采集相关的原始SAR数据方面发现了虽然是非排它性的但是有利的应用。
背景技术
如已知的，SAR系统领域最广泛的应用之一是干涉测量，其不同的解决方案包括估计受关注区域的数字高程模型(DEM)，估计目标的高度变化(差分干涉测量)以及检测两个场景中的变化(相干变化检测)。
用于干涉测量的SAR采集的主要模式是所谓的条带图，因为采集特性(具体地，方位角角度参数)对于沿方位角布置的目标基本上不变。虽然稍微复杂一些，但其它操作模式(具体地Spotlight、ScanSAR、PingPong和TOPS)也可用于干涉测量应用。
图1中示出了用于生成地球表面SAR图像的典型参考几何结构。关于这一点，希望强调的是，在图1中(以及将在下文呈现和描述的以下图中)，为了方便和简单的说明和描述，地球的表面(并且将被)示出为“平坦的”，而不失一般性。
具体地，图1示意性地示出了合成孔径雷达10(为了简化描述，在下文中称为SAR传感器)，其在假定为基本恒定的高度h(相对于地球表面)的飞行方向d上移动。如已知的，SAR传感器10的高度h是沿着通过所述SAR传感器10的最低点方向z测量的(具体地它通过SAR传感器10的天线的相位中心)并且与地球表面和飞行方向d正交。方便地，SAR传感器10由诸如例如飞机、无人驾驶飞行器(UAV)或卫星之类的航空/太空平台(为了简化图示未在图1中示出)在飞行/轨道中运输。飞行方向d的地面轨迹标识了与所述飞行方向d平行且与最低点方向z正交的方位角方向x，而与最低点方向z和方位角方向x正交的交叉轨道方向y与方位角方向x一起标识了与地球表面相切的x-y平面。在使用中，借助于适当的天线(为了简化图示未在图1中示出)，SAR传感器10传输雷达脉冲并在采集方向sr上接收相关的后向散射信号，采集方向sr标识倾斜范围并且形成与最低点方向z的仰角θ和与飞行方向d的斜视角(或等效地，方位角方向x)，在图1所示的SAR采集几何结构中，其等于(大约)90°(所谓的“几乎零多普勒”几何结构)。
具体地，在图1中所示的条带图类型的SAR采集几何结构中，SAR传感器10利用雷达脉冲照射地球表面的条带(称为刈幅)，然后从刈幅接收相关的后向散射信号，刈幅主要平行于方位角方向x延伸并且沿着交叉轨道方向y具有给定的宽度W。为了更加清晰，图2示出了x-y平面中条带图模式下SAR采集的几何结构，其中可以观察到斜视角是如何全部相同的；具体地，在图2所示的示例中，斜视角全是直角(“几乎零多普勒”几何结构)。
此外，图3从辐射度角度示意性地示出了如何通过SAR传感器10的天线射束(由附图标记12表示)以条带图模式照射位于地球表面上的目标(由附图标记11表示)。具体地，如图3中所示(其中参考先前引入的xyz笛卡尔参考系统的x-z平面，已经在图1和图2中使用)，目标11由整个天线射束12照射。更普遍地，沿着方位角方向x的图像的每个点用相同的射束部分照射，即它接收相同的光焦度；因此，在条带图模式中不存在所谓的“扇形”效应。
另外，图3还示出了应用/适用于与目标11相关的信号的辐射度校正(由附图标记13表示)(即使在条带图模式中并非严格必需该辐射度校正的情况下，因为如刚才所解释的，沿着方位角方向x用相同射束部分照射所有点。
如已知的，条带图模式下SAR采集的方位角分辨率是角度孔径(或角度差-Δ角度)的函数，SAR传感器通过该角度孔径观察目标；或者，以等效的方式，方位角分辨率也可以看作是与SAR传感器的速度相关联的时间差(Δ时间)的函数，利用该函数来观察目标。具体地，方位角分辨率可以由以下等式表示(对于进一步的细节，可以参考例如G.Franceschetti和R.Lanari的“Synthetic Aperture RADAR Processing”，CRC Press，1999年3月)：
其中resStrip指示条带图方位角分辨率，λ指示SAR传感器使用的波长，并且delta_angle指示SAR传感器观察目标的角度孔径(或角度差-Δ角度)。
假设角度为天线的3dB孔径(单向)(＝0.886λ/L，其中L指示沿着SAR传感器天线的方位角方向的物理或等效长度)，可以获得传统上与条带图模式的方位角分辨率相关联的约束条件，其等于L/2(更多细节请再次参考“Synthetic Aperture RADAR Processing”)。
目前，使用非常宽的天线射束来改善方位角分辨率，这些通过使用小尺寸或欠照射的天线或具有幅度和/或相位调制以便减小等效尺寸，或者通过使用所谓的聚光灯模式来实现，聚光灯模式的采集逻辑在图4中示意性地示出。
具体地，如图4中所示(其中参考已经在图1至图3中使用的先前引入的xyz笛卡尔参考系统的x-y平面)，聚光灯模式下的SAR采集逻辑设想在SAR传感器10的飞行运动期间使用连续的或准连续的天线射束转向(通过动态调整斜视角的值)，以便始终用雷达脉冲照射相同的受关注区域，然后从其中接收相关的后向散射信号，以这种方式增加SAR传感器10在所述受关注区域上的持续时间，从而改善方位角分辨率。
用于改善方位角分辨率的上述两种方法都存在一些缺点。具体地，使用非常宽的天线射束需要使用非常高的传输功率，而聚光灯模式引入了对刈幅的方位角长度的限制。
国际专利申请WO 2014/122624 A1提供了对上述缺点的一种解决方案，其涉及一种用于生成SAR图像的方法，该方法包括借助于合成孔径雷达在地球表面区域的条带图模式中执行N次SAR采集，该合成孔径雷达通过天线或卫星平台传输并且包括单个非分区天线和耦合到所述单个非分区天线的单个接收机，其中N是大于1的整数。
具体地，根据WO 2014/122624 A1，以相对于合成孔径雷达的最低点的相同的预定仰角执行条带图模式中的所有N次SAR采集，使得在条带图模式中的所述N次SAR采集都考虑到地球表面的相同特定刈幅。相反，使用相对于合成孔径雷达的飞行方向的相应斜视角执行条带图模式中的每次SAR采集，所述相应的斜视角不同于用于在条带图模式中执行其它N-1次SAR采集的斜视角。
此外，根据WO 2014/122624 A1，在条带图模式中执行的每次SAR采集包括相应的雷达发射和接收操作：
·在条带模式下执行的其它N-1次SAR采集的单个或一组雷达发射和接收操作的情况下是单独或成组交时的；以及
·包括相应采集方向上的相应雷达射束的发射和接收，所述采集方向由预定仰角和用于条带图模式中的所述SAR采集的相应斜视角限定，从而导致所述相应采集方向彼此平行并且与在条带图模式下执行的其它N-1次SAR采集的采集方向不平行。
此外，根据WO 2014/122624 A1，在两个紧接着的连续时刻中执行的并且与条带图模式中的两次不同SAR采集相关的两个雷达传输和接收操作中传输和接收的雷达射束沿着方位角是连续的，以这种方式相对于经由在条带映射模式中执行的N次SAR采集中的任何一次可获得的积分时间，增加积分时间。
最后，根据WO 2014/122624 A1的方法还包括基于在条带图模式中执行的所有N次SAR采集来生成特定刈幅的区域的SAR成像，所述SAR图像具有相对于标称条带图方位分辨率增强高达N因子的方位角分辨率，其等于沿合成孔径雷达的单个非分区天线的方位角方向的物理或等效长度的一半。
如已知的，为了最大限度地利用SAR系统的能力，除了条带图和聚光灯模式之外，多年来已经提出了各种其它技术来使用SAR传感器，包括所谓的突发模式。
具体地，突发模式基于时分的概念，即SAR采集被分成单独的时间块(所谓的“突发”)，其恰当地聚合，使得能够观察到：
·若干个子刈幅，以这种方式增加观察到的表面(ScanSAR模式)；或者
·使用不同电磁特性的相同区域，例如，使用不同的极化(PingPong模式)。
主突发模式是ScanSAR模式，它将时间分成同步突发(即以常规速率)。在每个突发中获取天线的标称脉冲重复频率(PRF)的场景部分，其方式使得方位角频谱被正确地采样，但是在更短的时间内(因此更小的Δ角度)。该时分使得能够在随后的突发中沿其它方向切换射束，以便增加范围内的刈幅。实际上，“方位角时间”被分成Nb个突发，并且在每个突发中获得不同的子刈幅。该模式的最大禁忌症是通过减小Δ角度，方位角分辨率降低。
具体地，为了确保采集在方位角上没有间隙(即为了确保连续的方位角覆盖，或者没有“空洞”)，最大采集时间和分辨率必须满足以下要求关系：
其中TScanSar指示最大ScanSAR采集时间，TStrip指示最大条带图采集时间，resScanSar指示ScanSAR分辨率，resstrip指示条带图参考分辨率(即L/2)，并且如前所述，Nb表示每个采集周期细分的突发数量，即观察到的子刈幅的数目(对于更多细节，请再次参考“SyntheticAperture RADAR Processing”，第225至229页)。
关于这方面，图5示意性地示出了ScanSAR模式中的SAR采集逻辑的示例。具体地，如图5中所示(其中参考已在图1至图4中使用的先前引入的xyz笛卡尔参考系统的x-z平面)，SAR传感器(为简单起见未在图5中示出)在沿着飞行方向d的运动期间，连续地执行一系列采集周期Ci(其中i＝1，2，3，...)，其中每个采集周期被分成三个突发(即，在图5的示例中的Nb＝3)。详细地，在每个采集周期Ci的第一突发中，SAR传感器获取(或观察)第一子刈幅，在每个采集周期Ci的第二突发中，SAR传感器获取(或观察)第二子刈幅，并且在每个采集周期Ci的第三脉冲中，SAR传感器获取(或观察)第三子刈幅。显然，必须执行相同子刈幅的采集，以确保沿方位角方向x的采集连续性。
因此，在ScanSAR模式中，相对于条带映射模式中发生的情况，每个目标被照亮的时间更短。此外，总是关于条带图模式，在ScanSAR模式中，每个目标由天线射束的不同部分照射。关于这一点，图6从辐射度角度示意性地示出了三个目标(分别指示为P1、P2和P3)如何被照亮，该三个目标具有不同的方位角位置并位于图5中的ScanSAR采集的示例中的第二采集周期C2的第一突发中观察到的第一子刈幅的部分中。详细地，在图6中，三个目标P1、P2和P3中的每一个的相应照射天线射束分别由附图标记14、15和16指示。如图6中所示，每个成像突发的各个点由不同的射束部分照射(即，照射光焦度沿方位角变化)。
因此，ScanSAR模式中存在扇形。实际上，由于沿着方位角由天线射束的不同部分照射表面的每个部分，因此必须在处理阶段中均衡天线模式以便获得正确的辐射度值。为了执行这种均衡，有必要知道天线的转向；如果天线转向的知识/估计出现错误，则会产生称为扇形效应的辐射度伪影。
关于这方面，图7再次示出了图5中的ScanSAR采集的示例，以及表示第一采集周期C1和第二采集周期C2的第一突发和第二突发中的残余辐射度误差的曲线图；具体地，在图7中，所述突发与之间的辐射度跳跃导致扇形非常明显。
另外，图8再次示出了天线射束14、15和16，其照射位于图5中示出的ScanSAR采集的示例中的第二采集周期C2的第一突发中观察到的第一子刈幅的部分中的目标P1、P2和P3(图6中所示)。图8还示出了相应的辐射度校正(对于三个目标P1、P2和P3中的每一个，相应的辐射度校正分别由附图标记17、18和19指示)，以及相应的辐射度误差。
关于ScanSAR模式的进一步细节，可以参考例如Andrea Monti-Guarnieri和Pietro Guccione的“Optimal“focus”for low resolution ScanSAR”，IEEE Transactionson Geoscience and Remote Sensing，2001年3月1日第39卷第3期，XP011021705，ISSN：0196-2892。本文讨论了低分辨率ScanSAR数据的聚焦，用于检测幅度图像和干涉测量应用。具体地，如本文摘要所述，当ScanSAR回波的所谓方位角时间带宽乘积很大时，常规ScanSAR聚焦技术可以非常有效地补偿方位角天线模式(例如没有扇形)，但当突发持续时间非常短时(在这种情况下，它被减少到输出的无效加权)则无法执行补偿。结果是聚焦脉冲响应的方位角变化失真，在多视图平均值(不可用于干涉测量应用)中部分地补偿失真，其代价是处理的多普勒带宽的减少。
相反，上述“Optimal“focus”for low resolution ScanSAR”文章提出了一种方法，其中使用一组短内核，每个短内核适合于聚焦在特定方位角值。这些内核经过优化，可以用最小均方误差重建反射率。伪反转收敛于用于广泛和短暂突发的“常规”聚焦。可以适当地调整这些方位角变化内核以满足分辨率/旁瓣权衡中的约束条件，并且能够更好地利用可用的多普勒带宽。
换句话说，上述“Optimal“focus”for low resolution ScanSAR”文章描述了一种用于在具有短方位角时间带宽乘积的ScanSAR模式下处理采集而不增加脉冲响应失真的方法，证明了算法的有效性也用于干涉测量应用。关于这方面，重要的是要注意，在所述“Optimal“focus”for low resolution ScanSAR”文章中，假设将始终使用零多普勒几何结构，其中使用包括相应的多普勒速率的单组处理参数来聚焦每个采集，而对于凭借所使用的几何结构(即，零多普勒几何结构)的所有采集，假设多普勒中心为零。
突发模式还包括TOPS模式，其使用变化的方位角天线转向(对于进一步的细节，可以参考例如F.De Zan和A.Monti Guarnieri的“TOPSAR：Terrain Observation byProgressive Scans”，IEEE Transactions on Geoscience and Remote Sensing，第44卷第9期，2006年9月，第2352至2360页，以及Adriano Meta等人的“TOPS Imaging WithTerraSAR-X：Mode Design and Performance Analysis”，IEEE Transactions onGeoscience and Remote Sensing，第48卷第2期，2010年2月1日，第759至769页，XP011296215，ISSN：0196-2892)。具体地，在TOPS模式中，扫描基本上在与聚光灯模式相反的方向上执行。也就是说，在每个突发中执行方位角控制，使得SAR传感器在每个突发的开始时“向后看”并且在突发的结束时向前指向。与常规ScanSAR模式不同，这确保了每个目标都被整个天线模式照亮，从而允许更好的辐射度均衡，这样使得扇形效应更不危急。
关于这方面，图9示意性地示出了TOPS模式中的SAR采集逻辑的示例。具体地，如图9中所示(其中参考已经在图1至图7中使用的先前引入的xyz笛卡尔参考系统的x-z平面)，SAR传感器(为了简化图示，未在图9中示出)在沿着飞行方向d的运动期间，连续地执行一系列采集周期Ci(其中i＝1，2，3，...)，每个采集周期被分成三个突发)。详细地，在每个采集周期Ci的第一突发中，SAR传感器获取(或观察)第一子刈幅，在每个采集周期Ci的第二突发中，SAR传感器获取(或观察)第二子刈幅，并且在每个采集周期Ci的第三脉冲中，SAR传感器获取(或观察)第三子刈幅。显然，在TOPS的情况下，必须执行相同子刈幅的采集，以确保沿方位角方向x的采集连续性。
如上所述，在每个突发中执行方位角转向，使得每个目标都被整个天线模式照射，以这种方式使得扇形效应更不危急。关于这方面：
·图10从辐射度角度示意性地示出了当在方位角上定位在图9的TOPS采集示例中的第一采集周期的C1第一突发中观察到的第一子刈幅的部分的中心时如何照射第一目标(由附图标记20指示)；以及
·图11从辐射度角度示意性地示出了当在方位角上布置在图9的TOPS采集示例中的第一采集周期C1的第一突发中观察到的第一子条带部分的横向位置中时如何照射第二目标(由附图标记21指示)。
正如规范的ScanSAR模式中所发生的那样，在TOPS模式中，传感器也始终以小于条带图模式典型角的Δ角度(或Δ时间)看到目标。因此，TOPS模式也允许像ScanSAR一样(即使具有关于规范的ScanSAR的不同等式/约束条件)，以相对于条带图模式降低方位角分辨率为代价来扩展范围内的刈幅。
如已知的，在文献中普遍需要扫描具有共同方位角角度的区域，其中经常参考所谓的方位角频谱重叠。具体地，在文献中推测存在物理约束条件，其中方位角重叠的缺失导致SAR图像之间的完全不相干。因此，这种不相干将使得任何干涉测量都不可能。
用于干涉测量应用的基于“几乎零多普勒”几何结构(即，基于传统的方位角频谱重叠约束条件)的SAR图像采集技术的示例在以下提供：
·Marcus 等人的“Study on Concepts for RadarInterferometry from Satellites for Ocean(and Land)Applications(KoRIOLIs)”中的“第5节：技术问题”，最终报告，2002年4月1日，第1至34页，其中明确断言各种采集的多普勒频谱重叠(因此，使用相同的斜视角)是实现数据一致性的必要先决条件，因此能够执行干涉测量处理(关于这方面，例如，请参阅“Study on Concepts for Radar Interferometryfrom Satellites for Ocean(and Land)Applications(KoRIOLIs)”中的“第5节：技术问题”第5至15页所述的内容)；以及
·A.Ferretti等人的“InSAR Principles：guidelines for SAR InterferometryProcessing and Interpretation”，欧洲航天局(ESA)出版，TM-19，2007年2月，其中再次明确断言只有存在频谱重叠才有可能干涉测量(关于这方面，例如请参考“InSARPrinciples：guidelines for SAR Interferometry Processing and Interpretation”第C-103页所述的内容)。
在许多情况下，对于干涉测量目的，通过强制每个目标以相同的方位角角度聚焦(以这种方式，不幸地降低性能，尤其是分辨率)来消除不常见的两次或更多次SAR采集的频谱部分。这种类型的处理也称为公共带通滤波，并且具有仅允许极少量数据用于干涉测量应用的大缺点。显然，在频谱重叠最小或为零的情况下，即在突发模式下的SAR采集的典型情况下(例如ScanSAR和PingPong)，公共频带滤波是不可用的。
具体地，过去还提出了各种技术来管理频谱重叠，并因此提出在突发模式中获得的SAR采集的干涉测量处理，例如：
·A.Monti Guarnieri和C.Prati的“ScanSAR Focusing and Interferometry”，IEEE Transactions on Geoscience and Remote Sensing，第34卷第4期，1996年7月，第1029至1038页，其中仅设想了具有非零方位角重叠和仅对公共方位角频谱部分进行滤波的数据的专用。
·J.Holzner和R.Bamler的“Burst-Mode and ScanSAR Interferometry”IEEETransactions on Geoscience and Remote Sensing，第40卷第9期，2002年9月，第1917至1934页，其中方位角重叠的重要性在基于ScanSAR模式的SAR采集的干涉测量应用中得到了强调，并且提出了通过突发的同步滤波和通过使用相同的中心处理原始数据，但不对由天线模式引起的辐射度畸变(扇形)进行任何补偿，以这种方式引入去相关现象；
·R.Grandin的“Interferometric Processing of SLC Sentinel-1 TOPSData”，欧洲航天局会议记录，SP-371，FRINGE 2015，2015年3月23日至27日，意大利弗拉斯卡蒂(罗马)，其中强调了具有相同多普勒频谱的采集重要性，并提出了一些基于间接估计现象的补偿技术；具体地，这些技术需要高容量的图像配准，因为即使配准中的最小误差也会导致最终干涉测量产品中的显著误差(事实上，对于所提出的示例，选择具有最小差值的多普勒中心的采集，即现象学影响有限的采集)；以及
·P.Prats、R.Scheiber和F.De Zan的“Interferometric TOPS ChainDescription”，SEOM-INSARAP：Sentinel-1 InSAR Performance Study with TOPS Data，2014年12月10日的INSARAP研讨会，其中强调了在中心差值的情况下，数据过滤是必要的，因此，信息会丢失；即使没有强调，这也意味着在中心存在显著差值的情况下不可能使用数据。
最后，在US 9,019,144 B2中以非常彻底的方式处理与方位角频谱重叠相关的约束条件，其中证明了实际上可以消除该约束条件，因为在文献中观察到的现象是由于数学关系而不是物理特性引起的。具体地，US 9,019，144 B2涉及一种用于获取用于干涉测量处理的SAR图像的方法，其包括经由在飞行中传输的一个或多个SAR传感器获取具有采集几何结构的相同区域的SAR图像，诸如以允许所述SAR图像的干涉测量处理。详细地，根据US 9,019，144 B2的方法的特征在于采集几何结构，其中该区域的每个SAR图像在相应的采集方向上获取，该采集方向限定相对于飞行方向的相应的斜视角，并且其中斜视角用于确定除零之外的平均斜视角。此外，SAR图像由单个SAR传感器获取，该SAR传感器在飞行中由航空/卫星平台进行传输，其使用单个天线并在航空/卫星平台的单次通过中获取SAR图像。
发明内容
本发明的第一个目的是提供一种技术，能够至少部分地减轻当前出于干涉测量目应用于突发模式SAR采集的技术的上述技术缺陷。
本发明的第二个目的是提供一种SAR处理技术，能够确保频谱重叠，即使频谱重叠在物理上不存在时(即当传感器以不同的斜视角观察场景时)。
本发明实现了这些和其它目的，只要本发明涉及如所附权利要求中所限定的用于干涉测量分析的SAR成像方法。
具体地，根据本发明的SAR成像方法包括：
·接收与通过一个或多个合成孔径雷达执行的对地球表面的同一区域的两次或更多次SAR采集相关的原始SAR数据；以及
·处理原始SAR数据以生成SAR图像。
根据本发明的方法的特征在于，对于每次SAR采集，基于不同的两组处理参数处理与所述SAR采集相关的相应原始SAR数据，所述不同的两组处理参数包括：
·第一组，所述第一组对于所有SAR采集而言都相同，并且包括基于与所有SAR采集相关的物理多普勒参数计算的聚焦多普勒参数；以及
·相应的第二组，所述相应的第二组包括与所述SAR采集相关并且基于与所述SAR采集相关的相应物理多普勒参数来计算的相应的辐射度均衡多普勒参数。
具体地，处理包括：
·基于所述聚焦多普勒参数聚焦与所有SAR采集相关的原始SAR数据；以及
·对于每次SAR采集，基于相应的辐射度均衡多普勒参数将相应的辐射度均衡应用于相应的SAR数据，以在不降低方位角分辨率并且不引入辐射度失真的情况下，补偿在用于执行SAR采集的合成孔径雷达的指向方面的可能差异。
方便地，所述SAR成像方法还包括：
·对于每次SAR采集，估计对用于执行所述SAR采集的合成孔径雷达的指向加以指示的相应的物理多普勒参数；
·对于每次SAR采集，基于估计的相应物理多普勒参数来计算相应的辐射度均衡多普勒参数；以及
·基于针对所有SAR采集而估计的物理多普勒参数来计算聚焦多普勒参数。
附图说明
为了更好地理解本发明，现在将参考附图(未按比例)描述通过非限制性示例提供的一些优选实施例，其中：
·图1、图2和图3示意性地示出了条带图模式下的SAR采集的示例；
·图4示意性地示出了聚光灯模式下的SAR采集示例；
·图5、图6、图7和图8示意性地示出了ScanSAR模式下的SAR采集的示例；
·图9、图10和图11示意性地示出了TOPS模式下的SAR采集的示例；
·图12示意性地示出了在条带图模式下的两次SAR采集，其中不存在频谱重叠；
·图13示意性地示出了在ScanSAR模式下的两次SAR采集，其中不存在频谱重叠；以及
·图14示意性地示出了根据本发明优选实施例的SAR成像方法。
具体实施方式
提供以下描述以使本领域的专家能够具体实施和使用本发明。对所示实施例的各种修改对于专家来说将是显而易见的，并且在不脱离本发明的保护范围的情况下，本文描述的一般原理可以应用于其它实施例和应用。
因此，本发明并不旨在仅限于本文描述和示出的实施例，而是与符合本文公开的和所附权利要求中限定的原理和特征的最宽范围相一致。
本发明涉及一种SAR成像方法，包括：
·接收与通过一个或多个合成孔径雷达执行的对地球表面的同一区域的两次或更多次SAR采集相关的原始SAR数据(优选地以突发模式执行SAR采集；甚至更优选地，ScanSAR、PingPong或TOPS类型的采集)；以及
·处理原始SAR数据以生成SAR图像。
具体地，处理包括：
·基于聚焦多普勒参数聚焦与两次或更多次SAR采集相关的原始SAR数据，所述聚焦多普勒参数是基于与所述两次或更多次SAR采集相关的物理多普勒参数来计算的；以及
·对于每次SAR采集，基于相应的辐射度均衡多普勒参数将相应的辐射度均衡应用于相应的SAR数据，所述辐射度均衡多普勒参数是基于与所述SAR采集相关的相应物理多普勒参数来计算的。
方便地，所述SAR成像方法还包括：
·对于每次SAR采集，估计指示用于执行所述SAR采集的合成孔径雷达的指向的相应物理多普勒参数(优选地，相应物理多普勒中心和/或相应物理多普勒速率)；
·对于每次SAR采集，基于估计的相应物理多普勒参数来计算相应的辐射度均衡多普勒参数(优选地，相应的辐射度均衡多普勒中心和/或相应的辐射度均衡多普勒速率)；以及
·基于针对两次或更多次SAR采集估计的物理多普勒参数来计算聚焦多普勒参数(优选地，聚焦多普勒中心和/或聚焦多普勒速率)。
更方便地，所述SAR成像方法还包括：对于每次SAR采集，接收指示所述SAR采集的相应空间和时间特性的相应数据；指示每次SAR采集的相应空间和时间特性的相应数据优选地包括指示用于执行所述SAR采集的合成孔径雷达的指向和位置以及指示对应的采集时间的数据。
另外，对于每次SAR采集，基于指示所述SAR采集的相应空间和时间特性的相应数据来方便地估计相应物理多普勒参数。
优选地，两次或更多次SAR采集是以突发模式执行的SAR采集，并且所述SAR成像方法还包括：
·对于每次SAR采集，基于指示所述SAR采集的相应空间和时间特征的相应数据确定每个相应突发的位置；以及
·根据每次SAR采集的突发位置确定突发的参考位置；其中，还基于突发的参考位置处理原始SAR数据。
更优选地，所述SAR成像方法还包括：
·配准生成的SAR图像；
·基于配准的SAR图像检测和估计突发位置误差；
·基于估计的突发位置误差校正突发的参考位置；以及
·基于突发的校正参考位置再次处理原始SAR数据。
方便地，所述SAR成像方法还包括：
·基于生成的SAR图像检测和估计辐射度均衡多普勒参数和/或聚焦多普勒参数中的误差；
·基于估计的误差更新辐射度均衡多普勒参数和/或聚焦多普勒参数；以及
·基于更新的辐射度均衡多普勒参数和/或聚焦多普勒参数再次处理原始SAR数据。
本发明源于申请人的见解，根据该见解，文献中观察到的现象实际上是由于数学关系而不是物理特性引起的，因此即使不存在物理重叠，也可以处理图像以确保频谱重叠。由于本发明，因此可以避免性能局限性并扩展干涉测量分析，这甚至对突发模式下的采集设想数据一致性。
具体地，与文献中的规定不同，根据本发明，通过将用于辐射度补偿的多普勒参数(方便地为多普勒中心和多普勒速率)与用于聚焦的那些分离，可以从辐射度的角度以最大频谱重叠(即相干图像)获得正确的图像。
实际上，在突发模式中，具体地在TOPS模式中，中心根据目标的方位角位置而快速变化(多普勒速率)；因此，如果在处理干涉测量采集中使用不同的多普勒速率，例如与各次采集中的天线的标称转向变化相关联，则使用不同的并因此不相干的角度处理各个目标。
在数学项中，将一般时刻沿着相对于卫星的零多普勒的方位角表示为贝塔(β)，在去除与方位角聚焦相关的二次方之后，从与第m个图像相关的等于德尔塔(Δ)的半孔径的积分得到的一般贝塔响应(β)可写为：
其中Ai、θi、分别表示一般第i个目标的幅度、位置和相位，Ga是天线的双向增益，Fw表示白化滤波器(或天线方向图的辐射度均衡滤波器)，以及
·DC表示与获取通用目标的斜视角相关的多普勒中心，
·表示用于与数据的聚焦方向相关的处理的多普勒中心，并且
·表示在处理期间用于补偿由天线增益引起的辐射度变化的多普勒中心。
当前的处理算法使用仅一个处理中心用于聚焦和辐射度均衡(即白化滤波器)，并且基于该处理中心必须与物理采集的处理中心重合的假设；因此，在这个假设并且还假设优选地知道天线方向图以及假设优选地地估计物理多普勒中心的情况下，发现
Fw＝Ga
这得出
考虑用不同多普勒中心获取的目标相位差Δint进行第二干涉采集，得出：
因此，两次采集之间的相位差由下式给出：
项
导致数据不一致。
为了更好地理解刚刚描述的内容，图12(其中参考了已在图1至7和图9至图11中使用的先前引入的xyz笛卡尔参考系统的x-z平面)示意性地示出了在相同目标22的条带图模式中的两次SAR采集，其中不存在频谱重叠。具体地，第一SAR采集由沿第一飞行方向d1移动的第一卫星23执行，而第二SAR采集由沿与第一飞行方向d1平行的第二飞行方向d2移动的第二卫星24执行。详细地，在图12所示的示例中，出现Δint＝2Δ。
在ScanSAR、TOPS和PingPong突发模式中，取决于目标的位置，多普勒参数存在显著变化，可以通过插入处理多普勒速率实际采集多普勒速率(DR)和在处理中用于白化滤波器的多普勒速率的概念来修改等式(1)。将多普勒速率表示为方位角角度与标称采集角度贝塔(β)之间的比率，得出：
在这种情况下，项(6)变为：
根据本发明的优选实施例，使用两个多普勒中心和两个多普勒速率，其中：
·第一多普勒中心和第一多普勒速率用于聚焦(具体地，用于确定聚焦角度)，以及
·第二多普勒中心和第二多普勒速率取而代之用于白化滤波器，即用于均衡辐射模式。
通过使用相同的处理多普勒中心和相同的处理多普勒速率来聚焦两个图像，项(6)(或项(8))变为空并且两个图像保持一致性。
如前所述，目前采用公共频带滤波以克服频谱重叠问题，即对图像/SAR采集滤波，仅留下共同角度。该技术仅在采集具有显著的频谱重叠时才起作用并且在任何情况下对性能具有严重影响。在突发模式中(例如，在ScanSAR模式中)，角度孔径非常小，以致于SAR传感器的最小指向误差可以导致完全没有频谱重叠，从而有效地防止了干涉测量应用的使用。
关于这方面，图13(其中参考已经在图1至图7和图9至图12中使用的先前引入的xyz笛卡尔参考系统的x-z平面)示意性地示出了相同目标25的ScanSAR模式下的两次SAR采集，其中在不存在频谱重叠。具体地，第一ScanSAR采集由沿第一飞行方向d1移动的第一卫星26执行，而第二ScanSAR采集由沿与第一飞行方向d1平行的第二飞行方向d2移动的第二卫星27执行。此外，所考虑的目标25在方位角上布置在第一突发中观察到的子刈幅的边缘处。
具体地，如图13中所示，利用当前的处理技术，由于天线的不同转向(并且因此不同的多普勒中心)，在不同的干涉测量采集中，在非同源突发中，因此在完全不同的方位角角度情况下，方位角边缘处的区域可以被聚焦。这将导致无法将这两次ScanSAR采集用于干涉测量分析。
相反，根据本发明的优选实施例，基于一对相应多普勒参数(方便地，多普勒中心和多普勒速率)的相应辐射度均衡应用于两次ScanSAR采集中的每一个，该对相应多普勒参数是基于ScanSAR采集的相应物理多普勒参数(即，基于相应实际采集几何结构)计算的。此外，总是根据本发明的所述优选实施例，两次ScanSAR采集都是基于同一对多普勒参数(方便地，多普勒中心和多普勒速率)被聚焦的，所述同一对多普勒参数是基于这两次ScanSAR采集的物理多普勒参数来计算的。
在TOPS模式中采集的情况下，多普勒中心快速变化；因此，如果不使用公共的多普勒频率，则将以不同的角度处理目标，从而导致不连贯。
本发明有效地克服了这种限制，并且甚至使得能够通过恢复从处理角度孔径到天线的角度孔径的约束来使用在突发模式中获取的数据用于干涉测量目的，这些参数通常不重合(具体地，对于突发模式或分辨率降低的模式)。
重要的是强调本发明解决了原点处频谱重叠的问题，因此，与根据US 9,019,144B2的方法不同，根据本发明的技术：
·用于任何平均斜视角，只要各种图像之间的最大斜视差值小于天线的角度孔径；以及
·使得能够针对各个采集以相同角度“整合”每个目标，即使这些采集是以不同的斜视角执行的，而无需设想进一步的数据处理以确保一致性，诸如求平均操作，因此没有关联的分辨率损耗。
为了更好地理解本发明，图14示出了根据本发明优选实施例的用于生成SAR图像的方法的流程图(由附图标记40整体指示)。
具体地，方法40包括：
·接收与由一个或多个卫星通过一个或多个SAR传感器优选地在突发模式下执行的地球表面的同一区域的两次或更多次SAR采集相关的输入数据(宏块41)；对于每次SAR采集，所述输入数据方便地包括相应的原始SAR数据(块411)(即，由执行采集的相应SAR传感器产生的数据)、指示SAR采集的相应特性(方便地，相应SAR传感器的指向，所使用的相应天线射束的特性等)的相应姿态数据(块412)以及相应的轨道数据和采集时间(块413)，所述相应的轨道数据与执行SAR采集的相应卫星相关并且所述相应采集时间指示SAR采集的相应时间特性；
·基于输入数据估计SAR采集的物理参数(宏块42)；
·基于物理参数确定处理参数(宏块43)；
·基于物理和处理参数处理(宏块44)输入数据；以及
·提供SAR图像(块451)作为输出(宏块45)，准备好通过一种或多种干涉测量分析技术进行处理(块50)，例如用于估计DEM(数字高程模型)，或用于差分干涉测量和/或相干变化检测。
详细地，关于估计物理多普勒参数(宏块42)，从每单个SAR采集的原始SAR数据(块411)和姿态数据(块412)开始，获得指示与单个SAR采集相关的所有方位角转向的物理多普勒参数(多普勒中心和/或多普勒速率)的估计(块421)。此外，基于每单个SAR采集的轨道数据和相关采集时间(块413)，确定每单个SAR采集的每单个突发的空间位置(或以等效方式，时间位置)(块422)。
取而代之地，关于确定处理参数(宏块43)，基于物理多普勒参数(多普勒中心和多普勒速率)，具体地基于与单个SAR采集相关的方位角转向，计算处理多普勒参数(块431)(多普勒中心和多普勒速率)，包括：
·对于每次SAR采集，基于与所述SAR采集相关的相应物理多普勒参数，具体地基于在执行所述SAR采集时SAR传感器的实际指向，计算相应的辐射度均衡多普勒参数；以及
·聚焦基于与所有SAR采集相关的物理多普勒参数计算的多普勒参数。
方便地，确定处理参数(宏块43)还包括确定等效主控(事实上，如果每次采集的多普勒中心值是已知的，则可以定义参考值，例如平均值)。
此外，突发的参考空间/时间位置(块432)(例如，最大化空间/时间重叠的突发的参考空间/时间位置)是基于每次SAR采集的每单个突发的空间/时间位置来确定的。
取而代之地，关于处理(宏块44)，基于物理多普勒参数、处理多普勒参数和突发的参考空间/时间位置来处理与两次或更多次SAR采集相关的原始SAR数据(块411)，以便形成SAR图像(块451)。具体地，处理(宏块44)包括(块441)：
·基于聚焦多普勒参数聚焦与所有SAR采集相关的原始SAR数据，以这种方式获得聚焦在雷达域(倾斜范围-方位角)中的相应图像(块451)；以及
·对于每个聚焦的SAR采集，基于与所述SAR采集相关的相应的辐射度均衡多普勒参数，将相应的辐射度均衡应用于相应SAR数据。
通过该处理(宏块44)，可以获得具有最小辐射度伪影的最终产品，并且使得每个目标变成针对各个采集以相同角度被“整合”。
此外，在聚焦和辐射度均衡之后(块441)，可以方便地校正物理和/或处理多普勒参数的估计(块442)，以便随后用于执行新的聚焦和辐射度均衡(块441)。具体地，可以通过分析SAR图像中的任何伪影残余来校正物理和/或处理多普勒参数。
此外，SAR图像是配准的，并且在方位角配准之后，可以估计任何突发位置误差(例如，由于SAR传感器的标称位置的估计中的误差而引起的)，并且如果必要时，校正突发的参考空间/时间位置(块452和432)，然后重新使用它们以执行新的聚焦和辐射度均衡(块441)。
最后，获得与干涉测量处理的标准算法兼容的最终产品(块50)。
方法40优选地通过适当配置和编程以执行所述方法40的上述步骤的电子处理装置来实现。方便地，所述电子处理装置被配置成直接从SAR遥感系统和/或数据库接收输入数据(宏块41)。所述电子处理装置可以被配置成将所生成的SAR图像(块451)提供给处理系统以进行干涉测量分析，或者可以适当地编程以基于所述SAR图像执行一个或干涉测量处理操作(块50)(例如，用于计算一个或多个高度，或者更一般地，用于估计DEM，用于差分干涉测量和/或变化检测，或计算干涉图、一致性图等的目的)。
显然，虽然关于卫星SAR采集进行了描述，但是方法40可以与空间SAR采集(例如，由卫星或航天器执行)和航空SAR采集(例如，由飞机、UAV、直升机等执行)一起使用。
以相同的方式，虽然已经关于突发模式中的SAR采集描述了该方法，但是方法40经过必要的修改显然还可用于条带图和聚光灯模式中的SAR采集。
从前面的描述中可以立即明显看出本发明的技术优点。
具体地，重要的是强调如下事实：如前所述，本发明有效地克服了当前已知技术的上述限制，因为即使频谱重叠物理上不存在，本发明也能够确保频谱重叠。
同样重要的是要注意，由于方位角天线方向图，没有已知技术使用不同的多普勒参数进行聚焦和辐射度补偿。
具体地，上述文章“Burst-Mode and ScanSAR Interferometry”提出了如在本发明中用于处理原始数据的相同中心的用途，但是，与本发明不同，它放弃以补偿由于天线模式引起的辐射度畸变(扇形)，引入去相关现象，并因此导致图像之间的差异，这降低了最终干涉测量产品的性能。此外，重要的是要强调“Burst-Mode and ScanSAR Interferometry”没有提及使用相同的多普勒速率，专门处理ScanSAR模式，并且绝不假定使用双参数进行聚焦和辐射度补偿(未实现)。
如前所述，“Optimal“focusing”for low resolution ScanSAR”也教示使用单组处理参数，而使用两组不同的处理参数(一组用于聚焦，一组用于辐射度均衡)的可能性未设想也未在所述文章中进行假定。
具体地，在“Optimal“focusing”for low resolution ScanSAR”的几个部分中，作者讨论了每次采集的单组聚焦参数，而没有引入技术来补偿各种干涉测量采集之间方位角控制的可能差值。
详细地，在“Optimal“focus”for low resolution ScanSAR”中，多普勒中心被认为是具有单个值的每次采集的特定单个参数(具体地，假设多普勒中心等于零，因为在文章中假设了“零多普勒”类型的采集几何结构)。
这同样适用于多普勒速率，其在“Optimal“focusing”for low resolutionScanSAR”被认为是单个参数。
该文章没有提供与使用不同方位角采集方向(即使用不同的斜视角)执行的SAR采集相关的问题的任何解决方案。
具体地，在“Optimal“focuS”for low resolution ScanSAR”第487页的“PhasePreserving Processing”的第V-A段中指出：“最终干涉测量质量的重要因素是处理器相位失真。在低TBPaz[方位角时间带宽乘积]ScanSAR干涉测量中，通过在两个配准的聚焦图像中的不同FIRF[聚焦脉冲响应函数]形状引入去相关源。当两个AAP[方位角天线模式]被移位时或者当两个图像中的一个从SAR聚焦图像合成时可能是这种情况...”。
“Optimal“focusing”for low resolution ScanSAR”的这一段V-A是文章中讨论干涉测量目的的不同采集之间关系的唯一部分。由于使用不同的方位角天线模式引起的去相关在“Optimal“focusing”for low resolution ScanSAR”中作为约束条件被提及，而没有提供任何解决方案来补偿这种去相关。事实上，在“Optimal“focusing”for lowresolution ScanSAR”的同一段(即VA)中，作者陈述：“在这种情况下，MMSE[最小均方误差]解决方案提供了与余弦窗参考相比的去相关性，然而MMSE表现更好，因为它探索更大的带宽。”。
总之，“Optimal“focus”for low resolution ScanSAR”中提出的技术不会引入额外的相位噪声，但是，与本发明不同，它也没有提供任何解决方案来补偿使用不同方位角采集方向(即使用不同的斜视角)执行的干涉测量采集的相位失真。
作为刚才解释的内容的补充证明，重要的是要注意“Optimal“focusing”for lowresolution ScanSAR”第482页的表II列出了用于ScanSAR聚焦的单组参数的值(其中指示单个多普勒速率，而没有指示多普勒中心，因为它被假设为零)。以这种方式，在“Optimal“focus”for low resolution ScanSAR”中使用的处理参数(具体地，聚焦参数)可以针对每次干涉测量采集而改变(在方位角采集几何形状改变的情况下)。“Optimal“focus”for lowresolution ScanSAR”中使用的这些聚焦参数实质上对应于每次采集的上述特定物理多普勒参数，根据本发明，这些参数用于计算每次采集的上述双组处理参数，即：
·上述聚焦多普勒参数对于所有SAR采集是相同的，并且基于与所有SAR采集相关的物理多普勒参数来计算；以及
·上述相应的辐射度均衡多普勒参数与所述SAR采集相关，并且基于与所述SAR采集相关的相应物理多普勒参数来计算。
此外，还应该注意的是，“Interferometric Processing of SLC Sentinel-1TOPSData”、“Interferometric TOPS Chain Description”，甚至“TOPS Imaging WithTerraSAR-X：Mode Design and Performance Analysis”都假定使用用于聚焦和辐射度补偿的双重参数。
总之，不同干涉测量采集之间的多普勒频移导致去相关。如前所述，为了处理这种去相关，传统文献教导使用所谓的公共频带滤波，但是该技术导致方位角分辨率降低，在突发模式下执行的采集不能处理该降低(例如在ScanSAR模式中)。
相反，本发明教导使用不同的两组处理参数来处理每次采集；具体地，相同的一组聚焦参数用于聚焦与要处理的所有SAR采集相关的原始SAR数据，而对于每次SAR采集，在不降低方位角分辨率并且不引入辐射度失真(即，不引入辐射度性能的降低)的情况下，还基于与所述SAR采集相关的相应组辐射度均衡参数执行相应辐射度均衡以补偿在用于执行各种SAR采集的SAR传感器的可能的指向方面的差异。
换句话说，通过解决由于原点处的不同采集几何结构引起的问题，本发明消除了由于各种采集之间的多普勒速率的不同姿态或差值而导致的问题，从而即使在物理上不存在时也保证了频谱重叠并且克服了目前已知技术的上述限制。
同样重要的是要注意，本发明的另外技术优点源于在聚焦和辐射度均衡(图14中的块441)之后，例如通过分析SAR图像中可能的伪影残余物来校正(图14中的块442)物理和/或处理多普勒参数(即聚焦和辐射度均衡参数)的估计的可能性。以这种方式，可以迭代地执行处理(即，聚焦和辐射度均衡-图14中的块441)，逐步细化物理和/或处理多普勒参数的估计(即，聚焦和辐射度均衡参数-图14中的块442)，因此，逐步改善最终产品的质量(即生成的SAR图像-图14中的块451)。
此外，如前所述，在SAR图像的方位角配准之后，根据本发明，还可以估计可能的突发位置误差(例如，由于SAR传感器的标称位置的估计中的误差引起的)并且，如果必要，校正突发的参考空间/时间位置(图14中的块452和432)，然后使用它们来执行新的聚焦和辐射度均衡(块441)。
综上所述，应清楚，在不脱离所附权利要求限定的本发明的范围的情况下，可以对本发明进行各种修改。

Claims (11)

1.一种用于干涉分析的SAR成像方法(40)，包括：
接收与通过一个或多个合成孔径雷达执行的对地球表面的同一区域的两次或更多次SAR采集相关的原始SAR数据；以及
处理所述原始SAR数据以生成SAR图像；
其特征在于，对于每次SAR采集，基于不同的两组处理参数处理与所述SAR采集相关的相应原始SAR数据，所述不同的两组处理参数包括：
第一组，所述第一组对于所有SAR采集而言都相同，并且包括基于与所有SAR采集相关的物理多普勒参数计算的聚焦多普勒参数；以及
相应的第二组，所述相应的第二组包括与所述SAR采集相关并且基于与所述SAR采集相关的相应物理多普勒参数来计算的相应的辐射度均衡多普勒参数；
并且其中处理包括：
基于所述聚焦多普勒参数聚焦与所有SAR采集相关的原始SAR数据；以及，
对于每次SAR采集，基于相应的辐射度均衡多普勒参数将相应的辐射度均衡应用于相应的SAR数据，以在不降低方位分辨率并且不引入辐射度失真的情况下，补偿在用于执行SAR采集的所述合成孔径雷达的指向方面的可能差异。
2.根据权利要求1所述的方法，还包括：
对于每次SAR采集，估计对用于执行SAR采集的所述合成孔径雷达的指向加以指示的相应的物理多普勒参数；
对于每次SAR采集，基于所估计的相应的物理多普勒参数计算相应的辐射度均衡多普勒参数；以及
基于针对所有SAR采集而估计的物理多普勒参数来计算聚焦多普勒参数。
3.根据权利要求2所述的方法，还包括：对于每次SAR采集，接收指示所述SAR采集的相应空间和时间特性的相应数据，并且所述相应数据包括指示用于执行所述SAR采集的合成孔径雷达的指向和位置的数据以及指示对应的采集时间的数据；
并且其中，对于每次SAR采集，基于指示所述SAR采集的相应空间和时间特性的所述相应数据来估计相应的物理多普勒参数。
4.根据权利要求3所述的方法，其中所述两次或更多次SAR采集是以突发模式执行的SAR采集。
5.根据权利要求4所述的方法，其中所述两次或更多次SAR采集是ScanSAR、PingPong或TOPS类型的采集。
6.根据权利要求4或5所述的方法，还包括：
对于每次SAR采集，基于指示所述SAR采集的相应空间和时间特性的相应数据确定每个相应突发的位置；以及
基于每次SAR采集的突发的位置确定突发的参考位置；
其中还基于突发的参考位置处理原始SAR数据。
7.根据权利要求6所述的方法，还包括：
配准所生成的SAR图像；
基于配准的SAR图像检测和估计突发位置误差；
基于所估计的突发位置误差校正突发的参考位置；以及
基于所校正的突发的参考位置再次处理原始SAR数据。
8.根据权利要求2至7中任一项所述的方法，其中：
对于每次SAR采集，所述相应的物理多普勒参数包括对用于执行所述SAR采集的合成孔径雷达的指向加以指示的相应物理多普勒中心和/或相应物理多普勒速率；
对于每次SAR采集，所述相应的辐射度均衡多普勒参数包括基于所述相应物理多普勒中心和/或所述相应物理多普勒速率计算的相应的辐射度均衡多普勒中心和/或相应的辐射度均衡多普勒速率；以及
所述聚焦多普勒参数包括基于与所述两次或更多次SAR采集相关的物理多普勒中心和/或物理多普勒速率计算的聚焦多普勒中心和/或聚焦多普勒速率。
9.根据前述任一项权利要求所述的方法，还包括：
基于所生成的SAR图像检测和估计所述辐射度均衡多普勒参数和/或所述聚焦多普勒参数中的误差；
基于所估计的误差更新所述辐射度均衡多普勒参数和/或所述聚焦多普勒参数；以及
基于所更新的辐射度均衡多普勒参数和/或聚焦多普勒参数再次处理原始SAR数据。
10.一种电子处理装置，被配置成执行根据前述任一项权利要求所述的SAR成像方法(40)。
11.一种计算机程序产品，包括软件代码部分，所述软件代码部分能够由电子处理装置执行，并且在被执行时使所述电子处理装置实现根据权利要求1至9中任一项所述的SAR成像方法(40)。
CN201780037690.XA 2016-06-21 2017-06-21 用于干涉测量分析的sar成像方法 Pending CN109313261A (zh)

Applications Claiming Priority (3)

Application Number Priority Date Filing Date Title
IT102016000064085 2016-06-21
ITUA2016A004568A ITUA20164568A1 (it) 2016-06-21 2016-06-21 Metodo di formazione di immagini sar per analisi interferometriche
PCT/IB2017/053689 WO2017221166A1 (en) 2016-06-21 2017-06-21 Sar imaging method for interferometric analyses

Publications (1)

Publication Number Publication Date
CN109313261A true CN109313261A (zh) 2019-02-05

Family

ID=57750371

Family Applications (1)

Application Number Title Priority Date Filing Date
CN201780037690.XA Pending CN109313261A (zh) 2016-06-21 2017-06-21 用于干涉测量分析的sar成像方法

Country Status (11)

Country Link
US (1) US11209540B2 (zh)
EP (1) EP3472643B1 (zh)
JP (1) JP2019519784A (zh)
KR (1) KR102391415B1 (zh)
CN (1) CN109313261A (zh)
AR (1) AR108852A1 (zh)
CA (1) CA3027670A1 (zh)
IL (1) IL263708B (zh)
IT (1) ITUA20164568A1 (zh)
RU (1) RU2018143986A (zh)
WO (1) WO2017221166A1 (zh)

Cited By (3)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
CN113129318A (zh) * 2021-04-25 2021-07-16 水利部信息中心 一种利用sar影像计算蓄滞洪区蓄水量的方法
CN113805175A (zh) * 2021-09-17 2021-12-17 北京工商大学 一种条带模式sar成像重聚焦方法
TWI807528B (zh) * 2021-01-12 2023-07-01 美商雷森公司 用於檢查水下載具的方法、設備及物品

Families Citing this family (10)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US11175398B2 (en) * 2016-12-21 2021-11-16 The Boeing Company Method and apparatus for multiple raw sensor image enhancement through georegistration
US10802135B2 (en) * 2016-12-21 2020-10-13 The Boeing Company Method and apparatus for raw sensor image enhancement through georegistration
US10989802B2 (en) * 2017-10-12 2021-04-27 Honeywell International Inc. Altimeter with high-resolution radar
US11933882B2 (en) * 2019-08-20 2024-03-19 Institut National D'optique Method and system for detection and synthetic aperture imaging of a target
KR20210062453A (ko) 2019-11-21 2021-05-31 삼성전자주식회사 운동 정보 결정 방법 및 장치
CN113993199B (zh) * 2020-07-27 2024-05-14 北京小米移动软件有限公司 天线辐射功率的调节方法及装置、电子设备
RU2748760C2 (ru) * 2020-11-09 2021-05-31 Виктор Андреевич Кузнецов Способ получения трехмерного радиолокационного изображения земной поверхности в двухпроходном интерферометрическом режиме съемки с беспилотного летательного аппарата
CN114217278B (zh) * 2021-12-14 2023-03-24 齐鲁空天信息研究院 一种基于动态目标的相控阵雷达幅相校准方法及系统
WO2023233563A1 (ja) * 2022-06-01 2023-12-07 三菱電機株式会社 飛翔体に搭載されるセンサシステム、及びセンサシステム用のデータ処理システム
CN116859344B (zh) * 2023-08-28 2023-11-03 中国电子科技集团公司第十四研究所 一种面向相干最优的能谱自适应分布式InSAR空间同步方法

Family Cites Families (2)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
ITTO20130108A1 (it) 2013-02-08 2014-08-09 Thales Alenia Space Italia S P A C On Unico Socio Innovativo metodo per generare immagini sar in modalita' stripmap
US10481257B2 (en) * 2014-10-30 2019-11-19 Mitsubishi Electric Corporation Synthetic aperture radar apparatus

Cited By (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
TWI807528B (zh) * 2021-01-12 2023-07-01 美商雷森公司 用於檢查水下載具的方法、設備及物品
US11753126B2 (en) 2021-01-12 2023-09-12 Raytheon Company Underwater vehicle inspection
CN113129318A (zh) * 2021-04-25 2021-07-16 水利部信息中心 一种利用sar影像计算蓄滞洪区蓄水量的方法
CN113805175A (zh) * 2021-09-17 2021-12-17 北京工商大学 一种条带模式sar成像重聚焦方法
CN113805175B (zh) * 2021-09-17 2022-04-19 北京工商大学 一种条带模式sar成像重聚焦方法

Also Published As

Publication number Publication date
IL263708A (en) 2019-01-31
US20190196005A1 (en) 2019-06-27
WO2017221166A1 (en) 2017-12-28
US11209540B2 (en) 2021-12-28
CA3027670A1 (en) 2017-12-28
AR108852A1 (es) 2018-10-03
EP3472643B1 (en) 2021-12-08
EP3472643A1 (en) 2019-04-24
IL263708B (en) 2021-12-01
KR20190035624A (ko) 2019-04-03
KR102391415B1 (ko) 2022-04-27
JP2019519784A (ja) 2019-07-11
ITUA20164568A1 (it) 2017-12-21
RU2018143986A (ru) 2020-07-23

Similar Documents

Publication Publication Date Title
CN109313261A (zh) 用于干涉测量分析的sar成像方法
Jao Theory of synthetic aperture radar imaging of a moving target
Duersch Backprojection for synthetic aperture radar
USH1720H (en) Time frequency processor for radar imaging of moving targets
Guarnieri Adaptive removal of azimuth ambiguities in SAR images
CN104076360B (zh) 基于压缩感知的二维sar稀疏目标成像方法
CN108459321A (zh) 基于距离-方位圆模型的大斜视高分辨率sar成像方法
EP1707987B1 (en) Method for mapping a target scene using scanning radar
US20080297405A1 (en) Synthetic Aperture focusing techniques
EP2817655B1 (en) Systems and methods for image sharpening
CN103901416A (zh) 一种基于稳健主成分分析法的多通道杂波抑制方法
CN113589287A (zh) 合成孔径雷达稀疏成像方法、装置、电子设备及存储介质
CN105022060A (zh) 针对快速空天目标的步进isar成像方法
Marston et al. Spatially variant autofocus for circular synthetic aperture sonar
CN113608218A (zh) 一种基于后向投影原理的频域干涉相位稀疏重构方法
CN113484859A (zh) 一种基于融合技术的二维超分辨雷达成像方法
CN102680956A (zh) 星载合成孔径雷达方位向第一模糊区回波信号能量抑制系统
CN110045374B (zh) 一种基于Chirp Scaling的多角度超长聚束SAR超分辨率成像算法
Sowjanya et al. Performance comparison of DAS and BP algorithms for Synthetic Aperture Radar Imaging
DE102004045273B4 (de) Verfahren zur Bilderzeugung bei einem Synthetischen Apertur Radar
Tema et al. Factors to Consider for Radiometric Calibration of Airborne FMCW SAR imagery
Kubica et al. Passive SAR imaging using the ASAR instrument of ENVISAT as transmitter of opportunity
Wang Pulse Regime Airborne Millimeter Wave Synthetic Aperture Radar
Zhang et al. Very high-resolution bistatic SAR imaging with TerraSAR-X as the illuminator in ST mode
Richter Coherent Azimuth Ambiguity Removal from Along-Track Interferometric Synthetic Aperture Radar Data for the Harmony Mission

Legal Events

Date Code Title Description
PB01 Publication
PB01 Publication
WD01 Invention patent application deemed withdrawn after publication

Application publication date: 20190205

WD01 Invention patent application deemed withdrawn after publication