一种用于sar匹配的特征区域检测方法
【技术领域】
[0001] 本发明涉及图像处理领域，尤其涉及合成孔径雷达图像处理及分析方法。
【背景技术】
[0002] 以下对本发明的相关技术背景进行说明，但这些说明并不一定构成本发明的现有 技术。
[0003] 在飞行器SAR成像导航中，通常采用星载SAR图像作为基准图或参考图。而哪处 区域的基准图可满足景象匹配的要求，是导航需要考虑的重要因素。不同的匹配算法对基 准图或参考图有不同的要求。目前SAR匹配方法多采用基于特征的匹配算法。因此，景象 区域应选择尽量包含弱散射地物、强散射地物共同构成的特征区域。如果选择的是单一纹 理、自相似性高的景象区域，如水平面或平坦的沙漠地带，都难以得到较好的匹配性能。
[0004] 为提取得到强、弱散射特性的区域，通常采用阈值分割法，即用一个或几个阈值将 图像的灰度直方图分成几个区间。然而对于灰度等级较少的SAR图像，不同区域的强度值 有较大的重叠，普通的阈值分割法难以获得理想的效果。因此，目前广泛应用于SAR图像的 阈值分割方法是基于恒虚警的方法。该方法的核心思想是通过分析SAR图像中像素强度的 统计特性，从而估计出目标分割的阈值，再对图像进行分割，在恒虚警率的约束条件下，做 到不丢失目标和尽可能地滤除杂波。该类方法的优点是分割速度快，缺点是需要对待分割 图像的阈值进行估计，因此阈值估计的效果将会极大地影响最终适配区选取的结果。陈石 平等在《高分辨率SAR图像CFAR分割的改进方法》文中对基于恒虚警的分割方法进行了改 进研究。
【发明内容】
[0005] 本发明的目的在于提出一种用于SAR匹配的特征区域检测方法，能够避免SAR图 像阈值分割难的问题，尤其适用于斑点噪声较大、地物散射特性值间有较大重叠的SAR影 像。
[0006] 根据本发明的用于SAR匹配的特征区域检测方法，包括：
[0007] S1、选取待检测图像中的任一点为观测点，获取以所述观测点为中心的初始窗口 对应区域的标准差；
[0008] S2、逐次递增所述初始窗口的尺寸，获取每个递增窗口对应区域的标准差；
[0009] S3、对初始窗口对应区域的标准差以及每个递增窗口对应区域的标准差进行显著 性分析，当显著性大于预设阈值时，判定所述待检测图像适于匹配。
[0010] 优选地，步骤S1之前进一步包括：
[0011] 根据所述待检测图像的分辨率以及所述待检测图像对应的地物的尺寸，确定所述 初始窗口的尺寸。
[0012] 优选地，所述初始窗口或所述递增窗口的形状为：由直线和/或曲线构成的封闭 图形。
[0013] 优选地，每个递增窗口的形状、以及每个递增窗口与初始窗口的形状相同或不同。
[0014] 优选地，步骤S2之前进一步包括：
[0015] 根据所述待检测图像的分辨率、所述待检测图像对应的地物的尺寸以及所述初始 窗口的尺寸，确定递增步长。
[0016] 优选地，每个递增窗口对应的递增步长相同或不同。
[0017] 优选地，每个递增窗口四周不同边缘点的递增步长相同或不同。
[0018] 优选地，步骤S2包括：
[0019] 按照递增步长，沿着所述观测点与初始窗口四周每个边缘点之间连线的方向递 增，得到第1个递增窗口，并获取第1个递增窗口对应区域的标准差；
[0020] 按照递增步长，沿着所述观测点与第η个递增窗口四周每个边缘点之间连线的方 向递增，得到第η+1个递增窗口，并获取第η+1个递增窗口对应区域的标准差；其中，η为正 整数。
[0021] 优选地，步骤S2包括：
[0022] 按照递增步长，扩大所述初始窗口的宽度和高度，得到第1个递增窗口的尺寸；
[0023] 以所述观测点为中心确定第1个递增窗口对应的区域，并获取第1个递增窗口对 应区域的标准差；
[0024] 按照递增步长，扩大第η个递增窗口的宽度和高度，得到第η+1个递增窗口的尺 寸；
[0025] 以所述观测点为中心确定第η+1个递增窗口对应的区域，并获取第η+1个递增窗 口对应区域的标准差；其中，η为正整数。
[0026] 优选地，所述递增窗口的尺寸满足公式1 :
[0027]
[0028] 式中，Wn为第η个递增窗口的宽度，单位为像素；Wni为第η-1个递增窗口的宽度， 单位为像素；Ηη为第η个递增窗口的高度，单位为像素；Hni为第η-1个递增窗口的高度，单 位为像素；AW为每次递增时的宽度递增步长；Ah为每次递增时的高度递增步长；η正整 数，η= 1时，W。为初始窗口的宽度，Η。为初始窗口的高度。
[0029] 根据本发明的用于SAR匹配的特征区域检测方法，包括：选取待检测图像中的任 一点为观测点，获取以所述观测点为中心的初始窗口对应区域的标准差；逐次递增所述初 始窗口的尺寸，获取每个递增窗口对应区域的标准差；对初始窗口对应区域的标准差以及 每个递增窗口对应区域的标准差进行显著性分析，当显著性大于预设阈值时，判定所述待 检测图像适于匹配。本发明通过获取初始窗口对应区域以及每个递增窗口对应区域的标准 差，能够考察待检测图像局部区域的散射特性变化情况，避开基于阈值分割方法面临的对 整幅待检测图像的阈值设定难、分割难的问题；通过标准差检测特征区域的匹配性，使得检 测方法不受观测点散射特性强或弱的影响。本发明的特征区域检测方法适用范围广，尤其 适用于斑点噪声较大、地物散射特性值间有较大重叠的SAR影像。
【附图说明】
[0030] 通过以下参照附图而提供的【具体实施方式】部分，本发明的特征和优点将变得更加 容易理解，在附图中：
[0031] 图1是示出根据本发明的用于SAR匹配的特征区域检测方法的流程图；
[0032] 图2是根据本发明实施例的待检测图像的示意图；
[0033] 图3是根据本发明实施例的待检测图像的监测结果示意图。
【具体实施方式】
[0034] 下面参照附图对本发明的示例性实施方式进行详细描述。对示例性实施方式的描 述仅是出于示范目的，而绝不是对本发明及其应用或用法的限制。
[0035] 本发明中，SAR是指合成孔径雷达，即SyntheticApertureRadar;本发明中的分 辨率是指每个待检测图像的像素对应的地物的面积。
[0036] 本发明通过考察局部区域的散射特性变化情况检测特征区域的匹配性，避开了待 检测图像阈值分割难的问题；通过标准差检测特征区域的匹配性，使得检测方法不受观测 点散射特性强或弱的影响。
[0037] 根据本发明的用于SAR匹配的特征区域检测方法起始于步骤S1，参见图1。步骤 S1选取待检测图像中的任一点为观测点，获取以所述观测点为中心的初始窗口对应区域的 标准差。标准差是相对量，以标准差作为检测待检测图像图像特征值的指标，使得检测方法 不依赖于特征区域的绝对图像特征值，而敏感于图像特征值的变化量。
[0038] 目前广泛应用于SAR图像的恒虚警阈值分割方法需要对待检测图像的阈值进行 估计，难以分割得到精确的强、弱散射区，因而影响了待检测图像特征区域的检测效果。为 了避免待检测图像阈值分割难的问题，本发明通过初始窗口以及多个递增窗口分别考察待 检测图像不同区域散射特性变化情况，从而检测特征区域的匹配性。
[0039] 无论观测点具备强或弱的散射特性，随着检测窗口尺寸的增大，若检测得到某个 递增窗口的标准差有较大变化，则该递增窗口对应区域必然出现了散射特性差异较大的不 同地物。因此，本发明的特征区域检测方法能够适用于不同波段、不同时相散射特性差异较 大的SAR遥感影像。
[0040] 根据本发明的优选实施例，初始窗口的形状为由直线和/或曲线构成的封闭图 形，比如，根据实时应用场景的不同，初始窗口的形状可以是椭圆形、圆形、等边多边形或者 其它不规则形状。
[0041] 初始窗口的尺寸大小直接影响根据本发明的特征区域检测方法的准确性，因此根 据本发明的优选实施例，步骤S1之前进一步包括：根据待检测图像的分辨率以及待检测图 像对应的地物的尺寸，确定初始窗口的尺寸。一方面，待检测图像的分辨率越高，待检测图 像中每个像素点对应的地物的面积越大，为了获得较好的检测效果，初始窗口的尺寸越小。 比如，对于5m分辨率的待检测图像，每个像素点对应的地物的面积为5mX5m，初始窗口的 尺寸越小，待检测图像中包含的初始窗口和递增窗口的数量越多，用于计算标准差的样本 数量越多，因此检测结果的准确性也就越高。另一方面，在保证用于计算标准差的样本数量 足够满足对检测结果准确性要求的情况下，待检测图像对应的地物的尺寸越大，用于计算 标准差的样本数量越多，计算过程越繁琐，为了适当降低计算过程的复杂程度，可以增大初 始窗口的尺寸。
[0042] 观测点可以是待检测图像的几何中心，也可以位于待检测图像