机载干涉合成孔径雷达影像的几何校正方法

摘要

本公开提供一种机载干涉合成孔径雷达影像的几何校正方法，包括：识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域；对所述非建筑物区域中的各像素进行几何校正；根据所述建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块，对于每一子块，计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数；根据所述仿射变换系数对子块内的各像素进行几何校正，得到高斯坐标下校正后的幅度图像；对所述幅度图像进行插值处理。

权利要求

1.一种机载干涉合成孔径雷达影像的几何校正方法，包括：

识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域；

对所述非建筑物区域中的各像素进行几何校正，包括：获取所述非建筑物区域中的各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔；根据各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔计算各像素的地面坐标；计算各像素的地面坐标在高斯网格中对应的坐标；根据各像素的坐标将各像素的图像幅度插入所述高斯网格中对应的位置；

根据所述建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块，对于每一子块，计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数；

根据所述仿射变换系数对子块内的各像素进行几何校正，得到高斯坐标下校正后的幅度图像，包括：根据所述仿射变换系数计算子块内各像素在高斯网格中对应的坐标，根据各像素的坐标将各像素的图像幅度插入所述高斯网格中对应的位置，得到高斯坐标下校正后的幅度图像；

对所述高斯坐标下校正后的幅度图像进行双线性插值消除采样不均匀的现象。

2.根据权利要求1所述的几何校正方法，其中，所述识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域包括：计算所述机载干涉合成孔径雷达影像的变差系数；

基于所述变差系数，采用最大类间方差法对所述机载干涉合成孔径雷达影像中的像素进行分类，识别所述建筑物区域和非建筑物区域。

3.根据权利要求1所述的几何校正方法，其中，所述根据各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔计算各像素的地面坐标包括：根据

计算各像素的地面坐标，其中，(XS，YS，ZS)为像素成像时的载机位置，(VX，VY，VZ)为像素成像时的载机速度，R0为图像近距，ρr为距离向像素间隔，fd为多普勒中心频率，λ为波长，(XG，YG)为像素的地面坐标，ZG像素的高程值，x为非建筑物区域中像素的横坐标。

4.根据权利要求1所述的几何校正方法，其中，所述计算各像素的地面坐标在高斯网格中对应的坐标包括：根据坐标转换关系将各像素的地面坐标转化为高斯坐标系下的坐标；

根据公式：

计算像素的地面坐标在高斯网格中对应的坐标，其中，(Xgauss，Ygauss)为各像素的地面坐标对应的高斯坐标系下的坐标，(Xgausss_lt，Ygauss_lt)为高斯网格一顶角的坐标，resX为高斯网格中X方向的采样间隔，resY为高斯网格中Y方向的采样间隔，(i，j)为高斯网格下的坐标。

5.根据权利要求1所述的几何校正方法，其中，所述根据所述建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块包括：根据各个建筑物的外接矩形将所述建筑物区域划分为至少一个子块，其中，每一个子块的四个角点落在所述非建筑物区域。

6.根据权利要求5所述的几何校正方法，其中，所述计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数包括：对于每一子块，计算子块的四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标；

建立高斯坐标系下子块与斜距坐标系下子块之间的仿射变换关系；

将所述四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标代入所述仿射变换关系求解所述仿射变换系数。

7.根据权利要求6所述的几何校正方法，其中，所述仿射变换关系为：其中，(x，y)为子块的角点在斜距坐标系下的坐标，(i，j)为子块的角点在高斯网格中的坐标，a0，a1，a2，a3，b0，b1，b2，b3为仿射变换系数。

说明书

机载干涉合成孔径雷达影像的几何校正方法

技术领域

[0001] 本公开涉及电子信息技术雷达技术领域，特别涉及一种机载干涉合成孔径雷达影像的几何校正方法。

背景技术

[0002] 合成孔径雷达(Synthetic Aperture Radar，SAR)影像几何校正的目的是修正SAR原始图像的几何变形，生成符合某种地图投影或图形表达要求的新图像。目前，机载SAR系统具备高精度的GPS定位系统，因此，几何校正过程需要重点考虑地形起伏的影响。

[0003] 干涉合成孔径雷达(Interferometric SAR，InSAR)利用合成孔径雷达(SAR)两个通道的干涉相位信息提取地表的高程信息，能够将SAR的测量拓展到三维空间。因而，机载干涉SAR能利用自身获得的高精度DEM进行影像的几何校正，对消除地形起伏对成像的影响具有优势。

[0004] 机载干涉SAR几何校正的通常的方法是根据已知的DEM，逐个像素联立距离多普勒方程求解其平面定位。这一方法存在的问题是：一是计算效率低，对大面积影像，逐个像素计算耗费时间长，另外，干涉SAR获得的DEM本身存在一定的噪声，在建筑物区域，受叠掩、阴影几何畸变影响，相位噪声大、DEM数据噪声较大，由此会造成平面定位的偏差，尤其在建筑物的边缘，定位偏差在几何校正后的影像上就会表现为边缘的扭曲，影响目视解译。

发明内容

[0005] 有鉴于此，本公开提供一种机载干涉合成孔径雷达影像的几何校正方法，包括：识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域；对所述非建筑物区域中的各像素进行几何校正；根据所述建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块，对于每一子块，计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数；根据所述仿射变换系数对子块内的各像素进行几何校正，得到高斯坐标下校正后的幅度图像；对所述幅度图像进行插值处理。

[0006] 根据本公开的实施例，其中，所述识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域包括：计算所述机载干涉合成孔径雷达影像的变差系数；基于所述变差系数，采用最大类间方差法对所述机载干涉合成孔径雷达影像中的像素进行分类，识别所述建筑物区域和非建筑物区域。

[0007] 根据本公开的实施例，其中，所述对所述非建筑物区域中的各像素进行几何校正包括：获取所述非建筑物区域中的各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔；根据各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔计算各像素的地面坐标；计算各像素的地面坐标在高斯网格中对应的坐标；根据各像素的坐标将各像素的图像幅度插入所述高斯网格中对应的位置。

[0008] 根据本公开的实施例，其中，所述根据各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔计算各像素的地面坐标包括：根据

[0009]

[0010] 计算各像素的地面坐标，其中，(XS，YS，ZS)为像素成像时的载机位置，(VX，VY，VZ)为像素成像时的载机速度，R0为图像近距，ρr为距离向像素间隔，fd为多普勒中心频率，λ为波长，(XG，YG)为像素的地面坐标，ZG像素的高程值。

[0011] 根据本公开的实施例，其中，所述计算各像素的地面坐标在高斯网格中对应的坐标包括：根据坐标转换关系将各像素的地面坐标转化为高斯坐标系下的坐标；根据公式：

[0012]

[0013] 计算像素的地面坐标在高斯网格中对应的坐标，其中，(Xgauss，Ygauss)为各像素的地面坐标对应的高斯坐标系下的坐标，(Xgauss_lt，Ygauss_lt)为高斯网格一顶角的坐标，resX为高斯网格中X方向的采样间隔，resY为高斯网格中Y方向的采样间隔，(i，j)为高斯网格下的坐标。

[0014] 根据本公开的实施例，其中，所述根据所述建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块包括：根据各个建筑物的外接矩形将所述建筑物区域划分为至少一个子块，其中，每一个子块的四个角点落在所述非建筑物区域。

[0015] 根据本公开的实施例，其中，所述计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数包括：对于每一子块，计算子块的四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标；建立高斯坐标系下子块与斜距坐标系下子块之间的仿射变换关系；将所述四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标代入所述仿射变换关系求解所述仿射变换系数。

[0016] 根据本公开的实施例，其中，所述仿射变换关系为：

[0017]

[0018] 其中，(x，y)为子块的角点在斜距坐标系下的坐标，(i，j)为子块的角点在高斯网格中的坐标。

[0019] 根据本公开的实施例，其中，所述根据所述仿射变换系数对子块内的各像素进行几何校正包括：根据所述仿射变换系数计算子块内各像素在高斯网格中对应的坐标；根据各像素的坐标将各像素的图像幅度插入所述高斯网格中对应的位置。

[0020] 根据本公开的实施例，其中，采用双线性插值法对所述幅度图像进行插值处理。

[0021] 本公开实施例提供的几何校正方法针对机载干涉SAR影像几何校正存在的难点，通过一定的策略判断出影像中建筑物等DEM噪声较大的区域，对于非建筑物区域，采用像素级的几何校正方法，保证几何校正的精度；对于建筑物区域，则通过块校正的方法，避免DEM噪声导致的校正质量下降，同时能够提高校正的速度，从而克服了像素级的几何校正方法效率低且在建筑物等区域校正产品质量较差的缺点，提高了机载干涉SAR影像几何校正的效率和校正产品的可判读性。

附图说明

[0022] 通过以下参照附图对本公开实施例的描述，本公开的上述以及其他目的、特征和优点将更为清楚，在附图中：

[0023] 图1示意性示出了根据本公开实施例提供的机载干涉合成孔径雷达影像的几何校正方法的总体流程图。

[0024] 图2示意性示出了根据本公开实施例提供的具体像素级校正和块校正结合校正的流程图

具体实施方式

[0025] 为使本公开的目的、技术方案和优点更加清楚明白，以下结合具体实施例，并参照附图，对本公开进一步详细说明。显然，所描述的实施例是本公开一部分实施例，而不是全部的实施例。基于本公开中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本公开保护的范围。

[0026] 在此使用的术语仅仅是为了描述具体实施例，而并非意在限制本公开。在此使用的术语“包括”、“包含”等表明了所述特征、步骤、操作和/或部件的存在，但是并不排除存在或添加一个或多个其他特征、步骤、操作或部件。

[0027] 在本公开中，除非另有明确的规定和限定，术语“安装”、“相连”“连接”、“固定”等术语应做广义理解，例如，可以是固定连接，也可以是可拆卸连接，或成一体；可以是机械连接，也可以是电连接或可以互相通讯；可以是直接连接，也可以通过中间媒介间接相连，可以是两个元件内部的连通或两个元件的相互作用关系。对于本领域的普通技术人员而言，可以根据具体情况理解上述术语在本公开中的具体含义。

[0028] 在本公开的描述中，需要理解的是，术语“纵向”、“长度”、“周向”、“前”、“后”、“左”、“右”、“顶”、“底”、“内”、“外”等指示的方位或位置关系为基于附图所示的方位或位置关系，仅是为了便于描述本公开和简化描述，而不是指示或暗示所指的子系统或元件必须具有特定的方位、以特定的方位构造和操作，因此不能理解为对本公开的限制。

[0029] 贯穿附图，相同的元素由相同或相近的附图标记来表示。可能导致本公开的理解造成混淆时，将省略常规结构或构造。并且图中各部件的形状、尺寸、位置关系不反映真实大小、比例和实际位置关系。另外，在权利要求中，不应将位于括号之间的任何参考符号构造成对权利要求的限制。

[0030] 类似地，为了精简本公开并帮助理解各个公开方面中的一个或多个，在上面对本公开示例性实施例的描述中，本公开的各个特征有时被一起分到单个实施例、图或者对其描述中。参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本公开的至少一个实施例或示例中。本说明书中，对上述术语的示意性表述不一定指的是相同的实施例或示例。而且，描述的具体特征、结构、材料或者特点可以在任何的一个或者多个实施例或示例中以合适的方式结合。

[0031] 此外，术语“第一”、“第二”仅用于描述目的，而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。因此，限定有“第一”、“第二”的特征可以明示或者隐含地包括一个或者更多个该特征。在本公开的描述中，“多个”的含义是至少两个，例如两个、三个等，除非另有明确具体的限定。

[0032] 本公开的目的在于提供一种像素级校正与块校正相结合的机载干涉SAR影像几何校正方法，首先通过一定的策略判断出影像中建筑物等数字高程模型(DEM)噪声较大的区域，然后对于非建筑物区域，采用常规的像素级的几何校正方法，保证几何校正的精度；对于建筑物区域，则通过块校正的方法，避免DEM噪声导致的校正质量下降，同时能够提高校正的速度。从而提高了机载干涉SAR影像几何校正的效率和校正产品的可判读性。

[0033] 图1示意性示出了根据本公开实施例提供的机载干涉合成孔径雷达影像的几何校正方法的总体流程图。图2示意性示出了根据本公开实施例提供的具体像素级校正和块校正结合校正的流程图。

[0034] 如图1和图2所示，该方法例如可以包括操作S101～操作S106。

[0035] 在操作S101，识别机载干涉合成孔径雷达影像中的建筑物区域和非建筑物区域。

[0036] 在本公开实施例中，可以通过计述机载干涉合成孔径雷达影像的变差系数，基于所述变差系数，采用最大类间方差法对机载干涉合成孔径雷达影像中的像素进行分类，识别建筑物区域和非建筑物区域。

[0037] 示例性的，可以通过：

[0038]

[0039] 计算图像的变差系数，进而得到图像的变差系数图。然后，采取最大类间方差法(OSTU)自动选取阈值对其进行分类，得到建筑物区域和非建筑物区域。

[0040] 在操作S102，对非建筑物区域中的各像素进行几何校正。

[0041] 在本公开实施例中，对非建筑物区域中的每一个像素(x，y)进行几何校正可以为：获取非建筑物区域中的各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔。根据各像素的高程值、各像素成像时的载机位置和载机速度、图像近距、距离向像素间隔计算各像素的地面坐标。计算各像素的地面坐标在高斯网格中对应的坐标。根据各像素的坐标将各像素的图像幅度插入高斯网格中对应的位置。

[0042] 具体的，首先，可以读取利用机载干涉SAR数据生成的DEM中该像素的高程值ZG。

[0043] 然后，根据

[0044]

[0045] 计算各像素的地面坐标，其中，(XS，YS，ZS)为像素成像时的载机位置，(VX，VY，VZ)为像素成像时的载机速度，R0为图像近距，ρr为距离向像素间隔，fd为多普勒中心频率，λ为波长，(XG，YG)为像素的地面坐标，ZG像素的高程值。由此可见，上述方程中除(XG，YG)之外，其他参数均为已知的，因此，解该方程即可以得出像素(x，y)的地面坐标(XG，YG)。

[0046] 接下来，根据坐标转换关系将该像素的地面坐标(XG，YG，ZG)转换到高斯坐标系下的坐标(Xgauss，Ygauss，Zgauss)，计算出该坐标在高斯网格中的位置为：

[0047]

[0048] 计算像素的地面坐标在高斯网格中对应的坐标，其中，(Xgauss，Ygauss)为各像素的地面坐标对应的高斯坐标系下的坐标，(Xgauss_lt，Ygauss_lt)为高斯网格一顶角的坐标，resX为高斯网格中X方向的采样间隔，resY为高斯网格中Y方向的采样间隔，(i，j)为高斯网格下的坐标。

[0049] 最后，根据非建筑物区域中各像素的坐标将各像素的图像幅度插入所述高斯网格中对应的位置，完成非建筑物区域中各像素的几何校正。

[0050] 在操作S103，根据建筑物区域中各建筑物的大小将所述建筑物区域划分为至少一个子块。

[0051] 在本公开实施例中，针对建筑物区域，采用块校正的方法进行几何校正。可以根据每一个建筑物区域的大小，设定子块的尺寸和位置。子块一般设定为建筑物区域的外接矩形，例如200×200，如果建筑物区域范围较大，可以划分为多个子块。子块的四个角点需落在非建筑物区域。

[0052] 在操作S104，对于每一子块，计算高斯坐标系下子块与斜距坐标系下子块之间的仿射变换系数。

[0053] 在本公开实施例中，对于每一子块，计算子块的四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标。建立高斯坐标系下子块与斜距坐标系下子块之间的仿射变换关系。将四个角点在高斯网格中的坐标以及在斜距坐标系下的坐标代入仿射变换关系求解所述仿射变换系数。

[0054] 具体的，对划分出的每一个子块，对该块的四个角点，分别按操作S102计算出四个角点在高斯网格中的位置(i1，j1)，(i2，j2)，(i3，j3)，(i4，j4)。

[0055] 对划分出的每一个子块，建立高斯坐标系下子块与斜距坐标系下子块之间的仿射变换关系可以为：

[0056]

[0057] 其中，(x，y)为子块的角点在斜距坐标系下的坐标，(i，j)为子块的角点在高斯网格中的坐标。

[0058] 将四个角点的斜距图像坐标(x1，y1)，(x2，y2)，(x3，y3)，(x4，y4)和高斯网格坐标(i1，j1)，(i2，j2)，(i3，j3)，(i4，j4)代入上式，构成如下线性方程组：

[0059]

[0060] 对划分出的每一个子块，求解对应的线性方程组得出仿射变换系数a0，a1，a2，a3，b0，b1，b2，b3。

[0061] 在操作S105，根据仿射变换系数对子块内的各像素进行几何校正，得到高斯坐标下校正后的幅度图像。

[0062] 对划分出的每一个子块，对块内的每一个像素(x，y)，利用计算出的仿射变换系数，计算几何校正后的图像高斯坐标(i，j)。根据各像素的坐标将各像素的图像幅度插入高斯网格中对应的位置，完整建筑区的块校正。

[0063] 在操作S106，对幅度图像进行插值处理。

[0064] 在本公开实施例中，由于干涉SAR侧视成像的特点，干涉SAR影像是在斜距方向进行等间隔采样，因此几何校正时将侧视方向转换成正射方向，会造成采样不均匀，需要对校正后的幅度图像进行插值处理。示例性的，可以采用双线性插值法进行插值。

[0065] 综上所述，本公开基于像素级校正与块校正相结合的方法对机载干涉SAR影像进行几何校正。首先利用变差系数图的最大类间方差分类区分出影像中的建筑物区域和非建筑物区域。对非建筑物区域进行像素级的几何校正；而对建筑物区域，则将影像划分为子块，根据四个角点的平面定位建立该子块内部的仿射变换关系，子块内部根据仿射变换关系进行几何校正。这样可以克服像素级的几何校正方法效率低且在建筑物等区域校正产品质量较差的缺点，提高了机载干涉SAR影像几何校正的效率和校正产品的可判读性。

[0066] 以上所述的具体实施例，对本公开的目的、技术方案和有益效果进行了进一步详细说明，所应理解的是，以上所述仅为本公开的具体实施例而已，并不用于限制本公开，凡在本公开的精神和原则之内，所做的任何修改、等同替换、改进等，均应包含在本公开的保护范围之内。