211173157_旋转不变描述和差异化约束的机载SAR图像连接点提取

数字技术与应用
第 41 卷114
其中，W i 是第i 个径向区间的描述子向量，V i ,j 是区间(i ,j)的描述子向量，U i ,j ,k 是V i ,j 中第k 个描述子数值。

描述子的维度为(1+12×2)×6=150。

为了增强描述子的
辐射不变性，需要将描述子向量归一化为单位长度。

图2 密集GLOH 格网结构示意图
Fig.2 Schematic diagram of dense GLOH grid structure
上述步骤描述了如何构造特征向量，但没有考虑图像之间的旋转变化。

为了消除图像旋转差异的影响，设计了一种旋转不变性增强方法，该方案为参考图像中每个特征点计算一个描述符向量，为待匹配中每个特征点计算12个不同方向的描述子向量。

下面详细介绍生成第s 个方向描述子向量的过程。

第s 个方向描述子向量基于图像间旋转差异为360s /12=30s 度的假设计算。

由于梯度值不随图像旋转而变化，可以通过轮转外侧两个径向区间的描述子向量快速获取不同方向的描述子向量。

第s 个方向的外侧两个径向区间描述子向量W i s 如式（3）所示：
,1,12,1
,12,1,,,,0,,,,,,1,,11s s
i i s
i s s s s
i s i i i s V V s W V V V V s + =  = = （3）第s 个方向描述子向量D s 如式（4）所示：
12
3
,,s
s
s D W W W = （
4）总之，上述方案仅在密集GLOH 格网中捕获一次描述符向量，可以通过轮转描述子向量快速生成12个方向的描述子向量。

1.2描述子方向估计
参考图像获取一个方向描述子且待匹配图像获取12个方向描述子后，通过NNDR 方法选择对应匹配。

对于一个匹配对，由于待匹配图像中每个特征点有12个不同方向的描述子向量，因此该匹配对能够得到12个距离比值。

距离比较小的匹配对更可能是正确匹配，因此可从12个距离比值中选择最小值作为该匹配对的最终距离比值，最终距离比值对应一个取值为1,2,…,12的方向索引。

选择最终距离比值最小的前300个匹配对构成样本
集合，统计这300个样本方向索引的直方图，由于该样本集合的匹配正确率较高，最大匹配对数量对应的方向索引值可认为是描述子的正确方向，从而完成描述子方向估计。

确定正确方向后，对于待匹配图像中的每个特征点，只保留正确方向的描述子，剔除其他方向的描述。

1.3差异化约束误匹配剔除
经过描述子方向估计后，每个特征点都只对应一个描述子向量，可以采用NNDR 匹配策略获得初始匹配对，该匹配点对通常包含一定数量的误匹配，需要剔除误匹配以精化匹配结果。

SAR-SIFT 算法中采用改进的RANSAC 算法，如FSC 算法或AC-RANSAC 算法，
在少量迭代的情况下，从包含大量误匹配的结果中稳健地提取正确的匹配点对。

这些算法均采用严格的几何变换模型，通常是仿射变换模型建立匹配点对几何关系，即有如式（5）、式（6）所示：
201121x a a x a y =++ （5）
201121y b b x b y =++ （6）
其中，式(5)和式(6)分别为距离向和方位向坐标变换公式。

x 1和x 2为两幅SAR 图像距离向坐标，y 1和y 2为两幅SAR 图像方位向坐标。

a i 和b i （i =0,1,2）为变换模型参数。

正确匹配点对的判定准则如式（7）所示
：
d T =<
（7）其中，T 为判定阈值，通常为一个较小的值，u 和v
分别为距离向和方位向残差，011212u a a x a y x =
++−，011212v b b x b y y =++−。

由于斜距SAR 图像距离向几何畸变引起的几何差异较大，严格几何模型会导致匹配精化结果陷入局部最优，大量正确匹配点对漏检。

为改善这一问题，本文在距离向和方位向采用差异化约束阈值，如式（8）所示：
1
2u T v T < < （8）其中，T 1为距离向阈值，可设为一个较大的值，T 1=100；T 2为方位向阈值，可设为一个较小的值。

本文在FSC 算法中采用差异化约束进行误匹配剔除，由于距离向阈值较大，在FSC 算法剔除误匹配后需要基于计算距离向残差，采用3σ法则进一步剔除误匹配点。

误匹配剔除的具体步骤为：
（1）基于NNDR 匹配策略提取两个匹配点对集合，分别为样本匹配集和总匹配集。

样本匹配集对应选择距离比最小的前300个匹配对，正确匹配点对占比较高；总匹配集对应的距离比率阈值为1，包含了大量的正确匹
2023年第 4 期
115
张昊：旋转不变描述和差异化约束的机载SAR 图像连接点提取配点对。

（2）从样本匹配集中随机选取一个样本（包含3个匹配点对），利用该样本初始化两幅图像间的仿射变换模型。

（2）对于总匹配集中的每个匹配点对，分别计算方位向和距离向残差，距离向残差小于100像素、方位向残差小于1.5像素的匹配点对认为是正确匹配点对，所有正确匹配点对构成初始正确匹配集。

（4）重新选取样本，重复步骤（2）～（3）K 次，获取匹配点对数量最大的初始正确匹配集，以及该集合中每个正确匹配对应的距离向残差。

（5）计算距离向残差的均值和方差，将残差不满足3σ法则的匹配点对从初始正确匹配集中剔除。

（6）重复步骤（5），直到初始正确匹配集不再变化，输出得到最终的正确匹配集。

2实验
为了验证本文方法的正确性和有效性，采用机载SAR 系统数据进行连接点匹配实验。

实验中，SAR 图像对是在原始SAR 图像重叠区域截取小图像块，图像块大小为1200×1200像素，共获取3对SAR 图像，分别覆盖居民地、农田和山地3种典型场景，如图3所示，这类图像用于验证算法精匹配机载SAR
图像块的能力。

(a) SAR 图像对1 (b) SAR 图像对2 (c) SAR 图像对3
图3 SAR 图像对Fig.3 SAR image pairs
实验评价指标包括匹配点对数量和计算耗时。

计算耗时是评估算法优劣的重要指标，计算耗时在Intel i7-8750H 处理器、32GB 内存和英伟达 Quadro P1000显卡的笔记本电脑、采用MATLABR2018b 软件统计。

实验采用3种连接点匹配算法，包括SAR-SIFT 算法、SAR-SIFT-RI 算法和SAR-SIFT-RIDC 算法，其中SAR-SIFT 算法采用FSC 算法进行误匹配点对剔除，SAR-SIFT-RI 算法构造旋转不变描述子并采用FSC 算
法进行误匹配点对剔除，SAR-SIFT-RIDC 算法构造旋转不变描述子并采用差异化约束的FSC 算法进行误匹配点对剔除。

SAR 图像的连接点匹配结果如表1和图4所示。

SAR 图像的连接点匹配结果表明本文算法能够更加有效精匹配机载SAR 图像块。

表1 3对SAR 图像连接点提取结果
Tab.1 Extraction results of 3 pairs of SAR image connection points
方法
SAR 图像对1SAR 图像对2SAR 图像对3
匹配点数量耗时/秒匹配点数量耗时/秒匹配点数量耗时/
秒SAR-SIFT 7236.17333.44534.4SAR-SIFT-RI 8028.47926.64728.1SAR-SIFT-RI-DC
287
28.2
295
26.7
243
29.1
SAR-SIFT SAR-SIFT-RI SAR-SIFT-RIDC
SAR 图像对1 SAR 图像对2 SAR 图像对3
图4 SAR 图像连接点提取结果
Fig.4 Extraction results of SAR image connection points
3结语
本文设计了一种利用旋转不变描述和差异化约束的机载SAR 图像连接点提取方法。

通过计算旋转不变梯度并在密集GLOH 网格中捕获梯度直方图，及在FSC 算法中引入方位向强约束、距离向弱约束的差异化约束模型，在误匹配剔除过程中有效避免正确匹配点对的漏检。

采用国产机载高分辨率毫米波SAR 数据进行连接点提取实验，实验表明旋转不变描述能够提高连接点提取效率，差异化约束能够增加提取的连接点数量。

引用
[1] 靳国旺,徐青,张红敏.合成孔径雷达干涉测量[M].北京:国防工业出版社,2014.
[2] 靳国旺,张红敏,徐青.雷达摄影测量[M].北京:测绘出版社, 2015.
[3] MIKOLAJCZYK K,SCHMID C.A Performance Evaluation of Local Descriptors[J].IEEE Transactions on Pattern Analysis and Machine Intelligence,2005,27(10):1615-1630.。