面向可见光和SAR影像配准的特征点检测

面向可见光和SAR影像配准的特征点检测 Feature point detection for optical and SAR remote sensing images registration 1 first-author 王丽娜（1986-），女，黑龙江鹤岗人，博士，讲师，2021年于中国科学院长春光机所获得博士学位。现就职于长春理工大学机电工程学院，主要从事多源遥感图像配准等方面的研究。E-mail： wangln@cust.edu.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995942&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995950&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995946&type= 22.01333237 32.00399780 王丽娜 （1986-），女，黑龙江鹤岗人，博士，讲师，2021年于中国科学院长春光机所获得博士学位。现就职于长春理工大学机电工程学院，主要从事多源遥感图像配准等方面的研究。E-mail： wangln@cust.edu.ac.cn wangln@cust.edu.ac.cn 1 2 2 1 corresp E-mail： shiguangfeng@cust.edu.cn. E-mail： shiguangfeng@cust.edu.cn. 石广丰（1981-），男，吉林长春人，教授，博导，2010年于中国科学院长春光机所获得博士学位，现为长春理工大学机电工程学院教授，主要从事图像特征检测与识别方面的研究。E-mail： shiguangfeng@cust.edu.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995953&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995961&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995958&type= 22.01332855 32.00399780 石广丰 （1981-），男，吉林长春人，教授，博导，2010年于中国科学院长春光机所获得博士学位，现为长春理工大学机电工程学院教授，主要从事图像特征检测与识别方面的研究。E-mail： shiguangfeng@cust.edu.cn shiguangfeng@cust.edu.cn 长春理工大学 机电工程学院，吉林 长春 130022 长春理工大学 机电工程学院，吉林 长春 130022 College of Electro-Mechanical Engineering， Changchun University of Science and Technology，Changchun 130022， China College of Electro-Mechanical Engineering， Changchun University of Science and Technology，Changchun 130022， China 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 中国科学院 长春光学精密机械与物理研究所，吉林 长春 130033 Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China Changchun Institute of Optics， Fine Mechanics and Physics， Chinese Academy of Sciences， Changchun 130033， China 由于可见光和SAR影像间的非线性辐射差异和SAR斑点噪声的影响，现有算法在对SAR影像进行特征点提取时，难以保证特征点的重复率，导致匹配性能下降。针对上述问题，本文提出了一种基于相位一致性矩特征的分块Harris特征点提取算法。首先，对输入图像进行分块操作；其次，定义了相位一致性中间矩，并结合最大矩和最小矩对每个图像块构建相位一致性多矩图；最后，在相位一致性多矩图上设计了投票策略，选取了在多矩图上重复出现超过半数的特征点，即稳定且重复率高的特征点作为最终的特征点。本文采用仿真的可见光和SAR图像作为实验数据，选取了三种不同的特征点检测算法与本文算法进行对比，实验结果表明，该算法能够克服SAR斑点噪声的影响和影像间的非线性辐射差异，有效地提高了特征点的重复率。可见光和SAR图像的配准结果表明，匹配点数较其他三种测试算法分别提高了23、26和35对，均方根误差分别降低了12.6%、37.2%和40.8%，有效地提升了配准算法性能。 The influence of SAR speckle noise makes it difficult for the existing state-of-the art algorithms to guarantee the repeatability rate of feature points when extracting them from optical and SAR images owing to the nonlinear radiation differences between optical and SAR remote sensing images， which consequently reduce the matching performance. To address the above problems， a Harris feature point extraction algorithm based on phase congruency moment feature is proposed. Firstly， blocking strategy was used to divide the input image into several image blocks； secondly， phase congruency intermediate moments were defined； then， phase congruency multi-moment maps were calculated for each image block； and finally， a voting strategy was designed on the phase congruency multi-moment maps. The feature points that appeared more than half of the time on the multi-moment image were selected as the final feature points. In this study， the simulated optical and SAR images were used as experimental data， and three different feature point detection algorithms were selected for comparison with the proposed algorithm. Experimental results showed that the proposed algorithm can overcome the influence of nonlinear radiation differences between optical and SAR remote sensing images and the SAR speckle noise， improving the repeatability rate of feature points effectively. The registration results on the real optical and SAR images showed that， compared with the other three algorithms， the matching points increased by 23， 26， and 35 pairs and the root mean square error decreased by 12.6%， 37.2%， and 40.8%， respectively. The performance of registration algorithm was improved effectively. 可见光和SAR影像 非线性辐射差异 斑点噪声 相位一致性 特征点提取 optical and SAR images nonlinear radiation difference speckle noise phase congruency feature points detection 1 引 言 1 随着传感器技术的发展，传感器的种类逐渐增加。可见光传感器技术相对成熟，使得可见光图像具有成像效果好，分辨率高和易于解读等优点，但是，可见光传感器成像容易受到天气的影响。合成孔径雷达（Synthetic Aperture Radar，SAR）属于主动发出微波的成像系统，具有穿透力强、对人造目标尤其是金属目标十分敏感、能够有效地识别伪装和揭露掩盖等特点，但是，SAR图像的可读性不如可见光图像 ［ 1 1 - 2 2 1-2 ］ 。将两种图像配合使用有助于形成信息互补，可应用于诸多场景，例如，图像融合，地物目标分类与识别等 ［ 3 3 - 5 5 3-5 ］ 。为了融合异源遥感数据实现对地观测，需要对图像进行高精度配准，而特征提取作为图像配准的前提步骤，其提取精度和特征点的重复率严重影响配准结果，因此，研究适用于可见光和SAR图像的特征提取方法，提升特征提取精度和重复率，已成为异源遥感图像处理领域的研究热点之一。 近年来，在计算机视觉领域出现了许多特征点检测算子，典型的算子如SIFT ［ 6 6 ］ 、SURF ［ 7 7 ］ 和Harris ［ 8 8 ］ 等。由于可见光和SAR图像间的非线性辐射差异和SAR斑点噪声的影响，将常规的特征点提取算法直接用于可见光和SAR图像提取时效果并不理想，特征点重复率较低，进而降低了配准算法性能 ［ 9 9 ］ 。为了克服上述问题，诸多学者在特征提取上做出改进，Ye等学者结合了Harris-Laplace和DOG两种算子的优势，提出了Harris-DOG检测算子 ［ 10 10 ］ ，该算子能够同时检测影像间的角点和Blob点，为影像匹配提供更多、更稳健的特征点。Xiang等学者对可见光和SAR图像分别采用不同的梯度计算方式，提取出可见光和SAR图像的一致性梯度信息，在此基础上构造了两个不同的Harris尺度空间来提取高度重复的特征点 ［ 11 11 ］ 。Zhang等学者同样改进了梯度计算方式，利用多尺度Sobel和多尺度ROEWA算子分别计算可见光和SAR图像的梯度来获得一致性较好的梯度结果，并在此基础上进行特征点提取，有效地提高了特征点的重复率 ［ 12 12 ］ 。Fan等学者提出一种均匀分布的特征点检测方法UNDSS-Harris ［ 13 13 ］ ，UNDSS-Harris利用非线性扩散滤波来构建非线性扩散尺度空间，较好地保留了边缘梯度特征；在此基础上，采用多尺度Harris算子、比例系数和分块策略，在可见光和SAR图像中提取出均匀分布的特征点。保文星等学者提出了一种改进的基于信息熵约束和KAZE特征提取的预处理算法 ［ 14 14 ］ 。该方法利用滑动窗口对整幅遥感图像进行遍历，计算出每个窗口的信息熵并形成直方图，通过阈值设定保留高信息熵以便用于特征提取。总的来说，上述方法的本质仍是基于图像梯度信息来提取影像间的显著特征，但可见光和SAR影像间的灰度差异较大，加之SAR斑点噪声的影响，严重影响了梯度计算结果，导致特征提取的重复率较低，进而降低了算法的匹配性能。 近年来，频域特征检测技术在异源遥感领域中的优势逐渐突显出来。比较典型的如相位一致性（Phase Congruency， PC） ［ 15 15 ］ 。该方法不受光照、灰度差异的影响，可大范围的检测图像中的角点、边缘和纹理等特征 ［ 16 16 ］ 。鉴于PC在特征检测方面的优势，其已成为学者们的研究热点。孙明超等学者基于相位一致性最大矩和相位一致性最小矩的叠加图进行Harris特征点提取，获得了稳定的角点和边缘点，用于后续的特征匹配 ［ 17 17 ］ 。Paul和Pati两位学者结合Gabor奇滤波器和多尺度Harris函数进行特征点检测 ［ 18 18 ］ 。Li等学者在相位一致性最大矩和最小矩上分别进行FAST（或Harris）特征点检测 ［ 19 19 ］ 。 诸多学者基于相位一致性提取出影像间的同名特征，但上述方法未充分考虑影像间的非线性辐射差异和SAR斑点噪声对特征点提取结果的影响。相位一致性在特征提取和特征描述方面尚有待发掘，算法稳定性有待提升。针对此问题，本文提出了一种基于相位一致性矩特征的Harris特征点提取算法——MMPC-Harris。首先，对待处理图像进行包含重叠区域的分块操作，增加特征点的均匀分布特性；然后，结合相位一致性最大矩和最小矩构建相位一致性多矩图，通过Harris算法实现特征的精细提取；最后，通过对多矩图上特征点进行投票和统计，选取重复出现超过半数的特征点作为最终的特征点。最终的特征点集合中包含了Harris特征点和在边缘上提取的Harris点，该方法保证了特征点的数量和稳定性。 2 相位一致性特征检测理论 1 相位一致性是基于频域信息进行特征提取的理论，最早由Morrone和Owens两位学者提出 ［ 20 20 ］ 。早期的相位一致性理论仅适用于一维信号，而图像为二维信号，为此，Kovesi在该理论基础上做出了改进，将其扩展到二维，提出了二维相位一致性特征检测理论。Kovesi采用多尺度、多方向的Log Gabor Filter来计算图像的局部相位特征。二维Log Gabor Filter的表达方式为： <math display="block" id="M1"><mi>F</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><mi>i</mi><mo>×</mo><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">d</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995965&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995964&type= 58.75866318 5.16466665 ， 其中， <math id="M2"><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995972&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995969&type= 16.25600052 5.16466665 和 <math id="M3"><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">d</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995977&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995975&type= 15.40933418 5.16466665 分别表示在尺度 <math id="M4"><mi>s</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995983&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995980&type= 1.26999998 2.79399991 和方向 <math id="M5"><mi>o</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995988&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995986&type= 1.86266661 2.79399991 上的偶对称和奇对称滤波器。利用二维Log Gabor Filter对输入图像 <math id="M6"><mi>I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995994&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32995991&type= 1.26999998 2.79399991 滤波，其实质为利用二维Log Gabor Filter的奇对称和偶对称滤波器分别与输入图像进行卷积运算，则图像上任意像素点 <math id="M7"><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996020&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996017&type= 11.09133339 3.97933316 处的响应分量为： <math display="block" id="M8"><mtable><mtr><mtd><mfenced open="[" close="]" separators="|"><mrow><msubsup><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">d</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>,</mo><msubsup><mrow><mi>O</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced><mo>=</mo></mtd></mtr><mtr><mtd><mfenced open="[" close="]" separators="|"><mrow><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo>,</mo><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><msubsup><mrow><mi>F</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">d</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup></mrow></mfenced></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996000&type= 43.77266693 11.51466751 ， 其中： <math id="M9"><msubsup><mrow><mi>O</mi></mrow><mrow><mi mathvariant="normal">e</mi><mi mathvariant="normal">v</mi><mi mathvariant="normal">e</mi><mi mathvariant="normal">n</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996007&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996005&type= 16.59466553 5.16466665 为偶分量， <math id="M10"><msubsup><mrow><mi>E</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">d</mi><mi mathvariant="normal">d</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msubsup><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996014&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996011&type= 15.49400043 5.16466665 为奇分量。进一步可以得出图像上任意像素点 <math id="M11"><mi>I</mi><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996020&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996017&type= 11.09133339 3.97933316 对应的相位一致性值 <math id="M12"><msub><mrow><mi>V</mi></mrow><mrow><mi mathvariant="normal">P</mi><mi mathvariant="normal">C</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996026&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996023&type= 5.75733376 3.80999994 为： <math display="block" id="M13"><mtable><mtr><mtd><msub><mrow><mi>V</mi></mrow><mrow><mi mathvariant="normal">P</mi><mi mathvariant="normal">C</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo></mtd></mtr><mtr><mtd><mfrac><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>o</mi></mrow></munder><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>s</mi></mrow></munder><mrow><msub><mrow><mi>W</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mfenced open="⌊" close="⌋" separators="|"><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mi mathvariant="normal">Δ</mi><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>-</mo><mi>N</mi></mrow></mfenced></mrow></mrow></mrow><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>o</mi></mrow></munder><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>s</mi></mrow></munder><mrow><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>+</mo><mi>ς</mi></mrow></mrow></mrow></mfrac></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996032&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996029&type= 69.84999847 20.74333191 ， 其中： <math id="M14"><msub><mrow><mi>W</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996039&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996035&type= 14.81666565 4.91066647 为随着滤波器响应变化的加权函数； <math id="M15"><mfenced open="⌊" close="⌋" separators="|"><mrow><mo>⋅</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996045&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996042&type= 3.64066648 3.97933316 运算表示符号内代数式的值为正时，符号内的计算结果与符号内代数式的计算结果相等，否则结果为0； <math id="M16"><mi>N</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996053&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996050&type= 2.87866688 2.87866688 表示估计的噪声阈值； <math id="M17"><mi>ς</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996058&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996055&type= 1.60866666 3.64066648 是一个小常数，避免分母为零的情况； <math id="M18"><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996064&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996061&type= 15.49400043 4.91066647 和 <math id="M19"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996069&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996066&type= 14.73200035 4.91066647 分别表示PC的幅值和相角，计算方式如下： <math display="block" id="M20"><msub><mrow><mi>A</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mroot><mrow><msub><mrow><mi>E</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msub><mrow><mi>O</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><msup><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996075&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996072&type= 61.21400070 5.24933338 ， <math display="block" id="M21"><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">a</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">c</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><mi mathvariant="normal">n</mi><mtext>  </mtext><mfenced separators="|"><mrow><msub><mrow><mi>O</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>,</mo><msub><mrow><mi>E</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996083&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996078&type= 65.95532990 4.82600021 ， 其中， <math id="M22"><mi mathvariant="normal">Δ</mi><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996090&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996088&type= 17.27199936 4.91066647 为相位差，计算方式如下： <math display="block" id="M23"><mtable><mtr><mtd><mi mathvariant="normal">Δ</mi><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mtext> </mtext><mo>=</mo><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">s</mi><mtext> </mtext><mo stretchy="false">(</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mover><mrow><mi>ϕ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><mo>-</mo><mfenced open="|" close="|" separators="|"><mrow><mi mathvariant="normal">s</mi><mi mathvariant="normal">i</mi><mi mathvariant="normal">n</mi><mtext>  </mtext><mo stretchy="false">(</mo><msub><mrow><mi>ϕ</mi></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mover><mrow><mi>ϕ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow></mfenced><mo>,</mo></mtd></mtr></mtable></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996097&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996094&type= 65.70133209 13.29266548 其中， <math id="M24"><msub><mrow><mover><mrow><mi>ϕ</mi></mrow><mo>¯</mo></mover></mrow><mrow><mi>s</mi><mo>,</mo><mi>o</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996103&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996100&type= 15.15533447 5.07999992 为平均相位角度。 进一步，Kovesi分析了相位一致性矩随方向变化的情况，提出利用基于相位一致性度量进行角点特征和边缘特征的检测方法 ［ 16 16 ］ 。Kovesi计算任意方向 o 和尺度 s 上的 PC 2 值来获得相位一致性的最大矩 M 和最小矩 m ，最终得到影像的角点特征和边缘特征， M 和 m 具体形式如下： <math display="block" id="M25"><mi>M</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo stretchy="false">(</mo><mi>a</mi><mo>+</mo><mi>c</mi><mo>+</mo><mroot><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><mi>a</mi><mo>-</mo><mi>c</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996110&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996106&type= 51.81600189 7.87400007 ， <math display="block" id="M26"><mi>m</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mo stretchy="false">(</mo><mi>a</mi><mo>+</mo><mi>c</mi><mo>-</mo><mroot><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mo stretchy="false">(</mo><mi>a</mi><mo>-</mo><mi>c</mi><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996117&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996113&type= 51.22333527 7.87400007 ， 其中： a，b，c 是中间变量，其具体形式如下： <math display="block" id="M27"><mi>a</mi><mo>=</mo><msup><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>o</mi></mrow></munder><mrow><mo stretchy="false">(</mo><msub><mrow><mi>V</mi></mrow><mrow><mi>P</mi><mi>C</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mi>c</mi><mi>o</mi><mi>s</mi><mtext> </mtext><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996125&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996121&type= 37.16867065 7.36600018 ， <math display="block" id="M28"><mi>b</mi><mo>=</mo><mn mathvariant="normal">2</mn><munder><mo movablelimits="false">∑</mo><mrow><mi>o</mi></mrow></munder><mrow><mo stretchy="false">(</mo><msub><mrow><mi>V</mi></mrow><mrow><mi>P</mi><mi>C</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mi>c</mi><mi>o</mi><mi>s</mi><mtext> </mtext><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>⋅</mo><mo stretchy="false">(</mo><msub><mrow><mi>V</mi></mrow><mrow><mi>P</mi><mi>C</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mi>s</mi><mi>i</mi><mi>n</mi><mtext> </mtext><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996131&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996128&type= 63.66933060 6.60400009 ， <math display="block" id="M29"><mi>c</mi><mo>=</mo><msup><mrow><munder><mo movablelimits="false">∑</mo><mrow><mi>o</mi></mrow></munder><mrow><mo stretchy="false">(</mo><msub><mrow><mi>V</mi></mrow><mrow><mi>P</mi><mi>C</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mi>s</mi><mi>i</mi><mi>n</mi><mtext> </mtext><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996134&type= 36.15266418 7.36600018 ， 其中， <math id="M30"><msub><mrow><mi>θ</mi></mrow><mrow><mi>o</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996141&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996139&type= 2.96333337 3.72533321 表示方向 o 对应的角度。 M 反映了边缘强度信息，可以用于边缘特征的检测； m 相当于角点检测子中的角点率，当像素点处的 m 值较大时，则该像素点是角点的可能性较大，通过设定 m 的特定阈值来筛选影像的角点特征。 3 基于改进相位一致性的可见光和SAR影像特征提取 1 基于相位一致性的最大矩和最小矩信息能够获取图像的边缘和角点信息。可见，图像的相位一致性矩信息可以有效表示图像的特征。本文提出图像的相位一致性矩特征表示如下： <math display="block" id="M31"><msub><mrow><mi>M</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mfenced separators="|"><mrow><mi>a</mi><mo>+</mo><mi>c</mi></mrow></mfenced><mo>+</mo><mfrac><mrow><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><mroot><mrow><msup><mrow><mi>b</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><msup><mrow><mfenced separators="|"><mrow><mi>a</mi><mo>-</mo><mi>c</mi></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996147&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996145&type= 56.38800049 7.87400007 ， 其中： a，b，c 是式（9）~ 式（11） 式（11） 的中间变量， <math id="M32"><msub><mrow><mi>M</mi></mrow><mrow><mi>k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996154&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996151&type= 4.40266657 3.72533321 表示第 <math id="M33"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996160&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996157&type= 1.69333339 2.79399991 个矩特征图像。进一步可以根据最大矩和最小矩特征将 公式（12） 公式（12） 描述为： <math display="block" id="M34"><msub><mrow><mi>M</mi></mrow><mrow><mi>k</mi></mrow></msub><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn><mo>+</mo><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">5</mn></mrow></msub><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn><mo>-</mo><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msub><mrow><mi>M</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996166&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996161&type= 44.78866577 7.95866632 ， 其中， <math id="M35"><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996172&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996170&type= 2.45533323 3.72533321 用于控制相位一致性矩的值，取值范围为： <math id="M36"><mfenced open="[" close="]" separators="|"><mrow><mo>-</mo><mn mathvariant="normal">1,1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996178&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996175&type= 11.76866722 3.97933316 。若相位一致性多矩图数量为 <math id="M37"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996229&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996226&type= 1.86266661 2.79399991 ，则参数 <math id="M38"><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996257&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996252&type= 2.45533323 3.72533321 以步长 <math id="M39"><mi>h</mi><mo>=</mo><mrow><mrow><mn mathvariant="normal">2</mn></mrow><mo>/</mo><mrow><mi>n</mi><mo>-</mo><mn mathvariant="normal">1</mn></mrow></mrow></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996196&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996193&type= 18.37266731 3.30200005 变化。 可见光和SAR影像间存在明显的非线性辐射差异，但二者的边缘相似度极高 ［ 21 21 ］ ，因此，考虑先提取出一致性较好的边缘特征，在边缘上进行特征点检测，来增加特征点数量。而相位一致性是一种频域特征检测理论，相比于直接利用梯度来提取边缘，能够更好地抵抗影像间的非线性辐射差异。鉴于此，本文提出了一种基于相位一致性矩进行特征提取的方法（MMPC-Harris），该方法旨在抵抗SAR斑点噪声的影响和影像间的非线性辐射差异，来获得均匀分布且重复率较高的特征点对，进一步提升配准算法的性能。 图1 图1 给出了算法流程图。 MMPC-Harris的具体步骤如下： （1）分块处理：将输入图像分成 <math id="M40"><msub><mrow><mi>S</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>×</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996211&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996208&type= 12.02266598 3.72533321 个图像块，并在相邻图像块之间增加 <math id="M41"><msub><mrow><mi>n</mi></mrow><mrow><mi mathvariant="normal">o</mi><mi mathvariant="normal">p</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996217&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996213&type= 4.06400013 3.72533321 个像素的重叠区域，以防止在分块边界处丢失特征信息； （2）对每个图像块构建相位一致性多矩特征图。通过多次实验发现，随着相位一致性多矩图数量 <math id="M42"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996229&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996226&type= 1.86266661 2.79399991 的增加，特征点重复率有所提高， <math id="M43"><mi>n</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996229&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996226&type= 1.86266661 2.79399991 超过5之后增加的速度明显变慢。综合考虑到算法的计算量，本文实验选取 <math id="M44"><mi>n</mi><mo>=</mo><mn mathvariant="normal">5</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996238&type= 8.63599968 2.79399991 。在每个相位一致性矩特征图像上基于Harris提取特征点，记录特征点位置和数量。 （3）基于统计原理，投票确定最终的特征点。若某一特征点在多矩图中出现的次数超过多矩图数量的一半，则将该特征点放入最终特征提取结果的集合中，并利用其在多矩图中坐标取平均值作为最终该特征点的位置信息。 （4）将每个图像块的特征提取结果合并，构成本文算法最终的特征点提取集合。 为了便于程序实现，参数设置为 <math id="M45"><mi>n</mi><mo>=</mo><mn mathvariant="normal">5</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996241&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996238&type= 8.63599968 2.79399991 ， <math id="M46"><mi>h</mi><mo>=</mo><mn mathvariant="normal">0.5</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996250&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996246&type= 11.59933376 2.79399991 ， <math id="M47"><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996257&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996252&type= 2.45533323 3.72533321 表示的序列为 <math id="M48"><mfenced open="{" close="}" separators="|"><mrow><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><mo>-</mo><mn mathvariant="normal">0.5,0</mn><mo>,</mo><mn mathvariant="normal">0.5,1</mn></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996266&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996261&type= 33.69733429 3.97933316 ，给出MMPC-Harris程序逻辑如图2所示。 MMPC-Harris算法 输入图像 1） 分块策略：将图像分成 <math id="M49"><msub><mrow><mi>S</mi></mrow><mrow><mi>n</mi></mrow></msub><mo>×</mo><msub><mrow><mi>S</mi></mrow><mrow><mi>m</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996273&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996270&type= 10.41399956 3.21733332 个子图像； 2） 分别对每个子图像计算PC值； for k = 1， <math id="M50"><mo>⋅</mo><mo>⋅</mo><mo>⋅</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996296&type= 3.30200005 2.37066650 ，5 do ｛ <math id="M51"><msub><mrow><mi>k</mi></mrow><mrow><mi mathvariant="normal">t</mi></mrow></msub><mo>=</mo><mo>-</mo><mn>1</mn><mo>+</mo><mn>0.5</mn><mo stretchy="false">(</mo><mi>k</mi><mo>-</mo><mn>1</mn><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996286&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996283&type= 29.12533188 3.80999994 ，构建 <math id="M52"><msub><mrow><mi>M</mi></mrow><mrow><mi>k</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996292&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996289&type= 3.80999994 3.21733332 ，Harris特征点提取；｝ end for 3） for i = 1， <math id="M53"><mo>⋅</mo><mo>⋅</mo><mo>⋅</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996296&type= 3.30200005 2.37066650 ， （角点数） do 统计特征点 i 在 <math id="M54"><msub><mrow><mi>M</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>~</mo><msub><mrow><mi>M</mi></mrow><mrow><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996306&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996302&type= 10.92199993 3.21733332 上出现的次数； if （出现次数）≥3 do ｛计算该特征点在 <math id="M55"><msub><mrow><mi>M</mi></mrow><mrow><mn>1</mn></mrow></msub><mo>~</mo><msub><mrow><mi>M</mi></mrow><mrow><mi>n</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996313&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996307&type= 10.92199993 3.21733332 上的坐标的平均值，作为该特征点最终的坐标；｝ end if end for 4） 合并特征点 输出： MMPC-Harris检测结果 下面通过模拟可见光和SAR图像进行实验分析，如 图3 图3 所示。在可见光图像上分别加入高斯噪声和斑点噪声来生成仿真数据，利用仿真图像作为待处理的可见光和SAR图像。 由仿真结果可知，在最小矩上提取的特征点重复率较高，但容易受噪声的影响，出现了一定的错误点对，如 图3 图3 （b）。在最大矩上提取的边缘点的数量较多，能够保证特征点数量，如 图3 图3 （f）。将两者直接相加作为最终的特征点集，尽管结果中含有稳定的角点和数量较多的边缘点，但仍存在一定问题，如直接叠加方式难以抵抗噪声干扰，导致结果中仍存在一些不稳定的特征点，如 图3 图3 （g）。而本文算法基于相位一致性多矩特征和统计原理来获得最终的特征点，以此来抵抗噪声干扰，剔除了虚假点，进而得到稳定的特征点对，保证了特征点的重复率；再加上分块处理，获取均匀分布的特征点，如 图3 图3 （h）。 4 实验与分析 1 本小节对MMPC-Harris特征点检测算子的性能进行实验分析，实验内容包括三个部分：首先，验证MMPC-Harris特征点检测算子对SAR斑点噪声的鲁棒性；其次，验证MMPC-Harris特征点检测算子对非线性辐射差异的鲁棒性；最后，基于一组实测可见光和SAR图像验证MMPC-Harris特征点检测算子对整体配准算法性能提升的影响。 4.1　MMPC-Harris对噪声鲁棒性评估 2 4 . 1 . 1　实验数据 3 目前尚且没有公开的可见光和SAR数据集用于测试特征点检测算法的优劣。为了便于计算特征点重复率，验证本文算法的有效性，本小节实验部分仍采用仿真数据，仿真数据的生成方式为对高清可见光图像分别加高斯噪声和不同级别的斑点噪声来生成仿真可见光和SAR图像，原始高清光学图像如 图4 图4 所示。 高清图像为长光卫星科技公司拍摄的智利圣地亚哥的机场和港口，拍摄时间为2017年3月，影像分辨率为1 m/pixel，图像大小均为1 000 pixel <math id="M56"><mo>×</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996348&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996344&type= 2.53999996 2.79399991 1 000 pixel。对上述两组高清光学图像分别加入高斯噪声和乘性噪声所生成仿真光学和SAR图像，结果如 图5 图5 所示。 4 . 1 . 2　评价方法和评价准则 3 为验证本文算法的先进性，选取了典型的特征点检测算法Harris ［ 8 8 ］ ，SAR-Harris ［ 21 21 ］ 和m+M-Harris ［ 17 17 ］ 作为对比算法。通过定性和定量两种方式来评价算法的性能。定性分析通过观察特征点的分布情况，均匀分布的特征点有利于几何变换模型的求取。定量分析则以不同SAR噪声级别下，特征点重复率 R 作为评价指标。参考图像 A 和待配准图像 B 上的一对特征点属于同名点前提为： <math display="block" id="M57"><msub><mrow><mfenced open="‖" close="‖" separators="|"><mrow><msub><mrow><mi>p</mi></mrow><mrow><mi>A</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>p</mi></mrow><mrow><mi>B</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>≤</mo><mi>d</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996364&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996360&type= 43.34933472 5.41866684 ， 其中： <math id="M58"><msub><mrow><mi>p</mi></mrow><mrow><mi>A</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996370&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996367&type= 13.29266548 4.91066647 和 <math id="M59"><msub><mrow><mi>p</mi></mrow><mrow><mi>B</mi></mrow></msub><mo stretchy="false">(</mo><mi>x</mi><mo>,</mo><mi>y</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996376&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996374&type= 13.12333298 4.91066647 分别表示特征点在参考图像和待配准图像中的位置坐标，符号 <math id="M60"><msub><mrow><mfenced open="‖" close="‖" separators="|"><mrow><mo>⋅</mo></mrow></mfenced></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996383&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996381&type= 5.84200001 4.57200003 为欧式距离测度， d 为欧式距离阈值。本文 d 设置为2，表征若两点之间的欧式距离 <math id="M61"><mo>≤</mo><mn mathvariant="normal">2</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996390&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996387&type= 6.18066692 3.30200005 个像素，则这两个点为一对同名点。 <math id="M62"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996427&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996423&type= 2.70933342 2.79399991 定义为： <math display="block" id="M63"><mi>R</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">2</mn><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></mrow><mrow><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub><mo>+</mo><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></mrow></mfrac><mo>×</mo><mn mathvariant="normal">100</mn><mi mathvariant="normal">%</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996402&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996398&type= 35.05199814 8.38199997 ， 其中： <math id="M64"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">c</mi><mi mathvariant="normal">o</mi><mi mathvariant="normal">r</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996408&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996405&type= 5.50333309 3.80999994 为正确匹配点对数， <math id="M65"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">o</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996413&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996411&type= 4.06400013 3.80999994 和 <math id="M66"><msub><mrow><mi>N</mi></mrow><mrow><mi mathvariant="normal">s</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996420&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996417&type= 3.80999994 3.80999994 分别为在可见光和SAR图像中提取的特征点总数量， <math id="M67"><mi>R</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996427&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=32996423&type= 2.70933342 2.79399991 值越大，特征点检测算子的性能越好。 4 . 1 . 3　实验结果与分析 3 本实验通过调整阈值，使得所提取的特征点总数量保持一致，每种算法均提取约600对特征点，取10次计算的平均值作为最终实验结果，以减小误差的影响。 图6 图6 给出了四种算法对两组图像的特征点重复率变化曲线。 由结果可见，MMPC-Harris在任何SAR噪声级别下的特征点重复率均获得了最高值，说明MMPC-Harris能在一定程度上抵抗SAR斑点噪声的影响，获得较高的特征点重复率。SAR-Harris的效果弱于本文方法，但优于Harris算法，这是由于SAR-Harris是在ROEWA算子提取SAR信息的基础上进行的Harris特征提取，ROEWA可以有效地抵抗SAR斑点噪声，获得较为理想的边缘特征，因此，在此基础上进行Harris特征提取能够抵抗噪声干扰。但是，SAR-Harris仅考虑了在SAR图像上进行特征提取所遇到的问题，并未对可见光和SAR这两种图像的异源性进行处理，所以，效果略差于本文方法。Harris算子的性能是最差的，随噪声水平的增加，特征点重复率值下降较快，表明Harris算子容易受到SAR乘性噪声影响，不适合用于直接处理SAR图像。 4.2　MMPC-Harris对非线性辐射差异的鲁棒性评估 2 4 . 2 . 1　实验数据 3 对 图4 图4 （a）的列乘以可变系数来模拟影像间的非线性辐射差异，结果如 图7 图7 所示，由图可见，影像间的同名特征呈现出不同的灰度特性，这种差异增加了同名特征提取的难度。 4 . 2 . 2　实验结果与分析 3 对比算法和评价指标同4.1节实验一致。选取第一组实验数据的特征提取结果进行对比，结果如 图8 图8 所示。 1;2; Harris和SAR-Harris算子的特征提取结果分别如 图8 图8 （a）和 图8 图8 （b）所示，二者基于灰度梯度信息进行特征提取，易受到影像间非线性辐射差异的影响，在可见光图像中可以检测到的特征在模拟的非线性辐射差异的图像中未被检出，特征点分布不均匀，较容易集中在灰度差异较大的区域。 图8 图8 （c）和 图8 图8 （d）分别给出了m+M-Harris和MMPC-Harris的特征提取结果，二者均能克服影像间的非线性辐射差异，特征点分布较均匀，在光学图像中提取的特征，在对应图像上仍然能够被检测到，这是因为本文方法和m+M-Harris方法是基于PC进行提取的，而PC算子对灰度和光照变化不敏感。进一步给出定量分析的结果来验证本文算子的性能，如 表1 表1 所示。 Group Harris SAR- Harris m+M-Harris MMPC-Harris Group 1 72.85% 41.76% 64.51% 87.91% Group 2 60.50% 71.56% 65.13% 90.53% 四种算子中，MMPC-Harris的特征点重复率值最高，这是由于本文方法在相位一致性多矩图上进行投票，剔除了一些不稳定的特征点，故本文算子可以较好地抵抗影像间的非线性辐射差异的影响，获得更多的同名点对。 4.3　MMPC-Harris特征点检测算子的配准性能 2 MMPC-Harris可以提升特征点的重复率，进而提升配准算法性能。为了证实这一点，基于一组实测的可见光和SAR图像进行配准实验，选取了Harris、SAR-Harris和m+M-Harris特征点检测算子作为对比算法，所有描述符均采用HOSMI ［ 22 22 ］ ，以验证MMPC-Harris特征点检测算子对可见光和SAR图像配准性能的提升。 4 . 3 . 1　实验数据和评价指标 3 实验数据如 图9 图9 所示，左侧为在Google Earth上截取的同一地区的可见光影像，右侧为机载SAR影像，拍摄时间为2017年7月。可见光和SAR图像的尺寸均为1 741 pixel×1 075 pixel，分辨率为3 m/pixel。拍摄场景为郊区，包含河流、低矮建筑物、农田和道路等。 配准实验采用主观和客观两种评价方式。主观评价方式为直接给出配准结果，客观评价指标为均方根误差（Root Mean Square Error， RMSE）和正确匹配点对数（Number of Correct Matches， NCM）。通过调整阈值，使得本文方法和对比方法均提取约1 200对特征点。采用平方差之和（Sum of Square Differences， SSD）作为匹配准则，SSD阈值设置为3个像素，并利用FSC ［ 23 23 ］ 剔除误匹配。 4 . 3 . 2　实验结果与分析 3 利用四种检测算子的特征提取结果进行配准，结果如 图10 图10 所示， 表2 表2 显示了四种检测算子的NCM和RMSE。 Algorithm RMSE NCM Harris+HOSMI 2.82 23 SAR-Harris+HOSMI 2.66 32 m+M-Harris+HOSMI 1.91 35 MMPC-Harris+HOSMI 1.67 58 由配准结果可见，四种检测算子都能完成配准，但配准性能存在差异。Harris和SAR-Harris是直接利用灰度信息进行特征提取的，所获得的正确匹配点对分布不够均匀，这说明直接利用灰度信息进行特征提取难以抵抗影像间的非线性辐射差异。m+M-Harris和MMPC-Harris均取得了较好的效果，这说明基于PC进行特征提取有助于抵抗影像间的非线性辐射差异。 结果表明，与其他特征点检测算子相比，在描述符一致的情况下，MMPC-Harris获得了最多的NCM和最小的RMSE，这是因为MMPC-Harris是基于PC特征进行特征点检测的，在一定程度上抵抗了影像间的非线性辐射差异，同时，MMPC-Harris检测了角点和边缘上角点，保证了特征点的数量，通过建立相位一致性多矩图和设计的投票策略，在一定程度上剔除了噪声的影响，从而获得大量稳定特征点用于图像匹配。此外，MMPC-Harris使用了分块策略，能够获得均匀分布的同名点对，可以更好地适应影像间的几何形变。总之，MMPC-Harris可以显著提高在可见光和SAR影像间所提取的特征点重复率，进而提升配准精度。 5 结 论 1 本文针对可见光和SAR影像间的非线性辐射差异和SAR斑点噪声导致在影像间提取的特征点重复率较低的问题，提出了一种基于相位一致性的特征点提取算法MMPC-Harris。MMPC-Harris通过计算图像的相位一致性矩信息，构建了相位一致性多矩图，在所构建的相位一致性多矩图上提取Harris角点，并设计了投票策略，选取了重复出现超过半数的点作为最终的特征点，获得了稳定的角点和边缘点。仿真数据的测试结果表明，MMPC-Harris可以较好地抵抗影像间的非线性辐射差异和SAR斑点噪声的影响，有效地提高了特征点的重复率。一组实测可见光和SAR图像上的测试结果表明，与其他三种方法对比，MMPC-Harris的正确匹配点对数分别提高了23、26和35对，均方根误分别降低了12.6%、37.2%和40.8%，能够有效地提升配准算法性能。但是，MMPC-Harris不具备旋转和尺度不变性，下一步工作将针对此问题进行更加深入地研究。 参考文献 SHITOLE S ， SHARMA M ， DE S ， et al . Local contrast based adaptive SAR speckle filter ［J］. Journal of the Indian Society of Remote Sensing ， 2017 ， 45 （ 3 ）： 451 - 462 . doi: 10.1007/s12524-016-0613-2 http://dx.doi.org/10.1007/s12524-016-0613-2 Local contrast based adaptive SAR speckle filter 吴文达 ， 张葆 ， 洪永丰 ， 等 . 机载红外与合成孔径雷达共孔径天线设计 ［J］. 中国光学 ， 2020 ， 13 （ 3 ）： 595 - 604 . doi: 10.3788/CO.2019-0160 http://dx.doi.org/10.3788/CO.2019-0160 机载红外与合成孔径雷达共孔径天线设计 WU W D ， ZHANG B ， HONG Y F ， et al . Design of co-aperture antenna for airborne infrared and synthetic aperture radar ［J］. Chinese Optics ， 2020 ， 13 （ 3 ）： 595 - 604 . （in Chinese） . doi: 10.3788/CO.2019-0160 http://dx.doi.org/10.3788/CO.2019-0160 Design of co-aperture antenna for airborne infrared and synthetic aperture radar 眭海刚 ， 徐川 ， 刘俊怡 . 基于特征的光学与SAR遥感图像配准 ［M］. 北京 ： 科学出版社 ， 2017 . 基于特征的光学与SAR遥感图像配准 SUI H G ， XU C ， LIU J Y . Feature Based Optical and SAR Remote Sensing Image Registration ［M］. Beijing ： Science Press ， 2017 . （in Chinese） Feature Based Optical and SAR Remote Sensing Image Registration KULKARNI S C ， REGE P P . Pixel level fusion techniques for SAR and optical images： a review ［J］. Information Fusion ， 2020 ， 59 ： 13 - 29 . doi: 10.1016/j.inffus.2020.01.003 http://dx.doi.org/10.1016/j.inffus.2020.01.003 Pixel level fusion techniques for SAR and optical images： a review 吕宝林 ， 佟首峰 ， 徐伟 ， 等 . 基于配准的机载红外非均匀性校正技术应用 ［J］. 中国光学 ， 2020 ， 13 （ 5 ）： 1124 - 1137 . doi: 10.37188/CO.2020-0109 http://dx.doi.org/10.37188/CO.2020-0109 基于配准的机载红外非均匀性校正技术应用 LV B L ， TONG S F ， XU W ， et al . Non-uniformity correction of airborne infrared detection system based on inter-frame registration ［J］. Chinese Optics ， 2020 ， 13 （ 5 ）： 1124 - 1137 . （in Chinese） . doi: 10.37188/CO.2020-0109 http://dx.doi.org/10.37188/CO.2020-0109 Non-uniformity correction of airborne infrared detection system based on inter-frame registration LOWE D G . Distinctive image features from scale-invariant keypoints ［J］. International Journal of Computer Vision ， 2004 ， 60 （ 2 ）： 91 - 110 . doi: 10.1023/b:visi.0000029664.99615.94 http://dx.doi.org/10.1023/b:visi.0000029664.99615.94 Distinctive image features from scale-invariant keypoints BAY H ， TUYTELAARS T ， GOOL L V . SURF： Speeded Up Robust Features ［C］. European Conference on Computer Vision Springer， Berlin， Heidelberg ， 2006 ： 404 - 417 . doi: 10.1007/11744023_32 http://dx.doi.org/10.1007/11744023_32 SURF： Speeded Up Robust Features HARRIS C ， STEPHENS M . A combined corner and edge detector ［C］. Proceedings of the Alvey Vision Conference 1988. Manchester. Alvey Vision Club ， 1988 ： 147 – 152 . doi: 10.5244/c.2.23 http://dx.doi.org/10.5244/c.2.23 A combined corner and edge detector 樊建伟 . 基于特征点的SAR图像配准算法研究 ［D］. 西安 ： 西安电子科技大学 ， 2017 . doi: 10.1109/tgrs.2016.2626373 http://dx.doi.org/10.1109/tgrs.2016.2626373 基于特征点的SAR图像配准算法研究 FAN J W . Research on SAR Image Registration Based on Feature Point ［D］. Xi'an ： Xidian University ， 2017 . （in Chinese） . doi: 10.1109/tgrs.2016.2626373 http://dx.doi.org/10.1109/tgrs.2016.2626373 Research on SAR Image Registration Based on Feature Point YE Y X ， WANG M M ， HAO S Y ， et al . A novel keypoint detector combining corners and blobs for remote sensing image registration ［J］. IEEE Geoscience and Remote Sensing Letters ， 2021 ， 18 （ 3 ）： 451 - 455 . doi: 10.1109/lgrs.2020.2980620 http://dx.doi.org/10.1109/lgrs.2020.2980620 A novel keypoint detector combining corners and blobs for remote sensing image registration XIANG Y M ， WANG F ， YOU H J . OS-SIFT： a robust SIFT-like algorithm for high-resolution optical-to-SAR image registration in suburban areas ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2018 ， 56 （ 6 ）： 3078 - 3090 . doi: 10.1109/tgrs.2018.2790483 http://dx.doi.org/10.1109/tgrs.2018.2790483 OS-SIFT： a robust SIFT-like algorithm for high-resolution optical-to-SAR image registration in suburban areas ZHANG X T ， WANG Y H ， LIU H W . Robust optical and SAR image registration based on OS-SIFT and cascaded sample consensus ［J］. IEEE Geoscience and Remote Sensing Letters ， 2022 ， 19 ： 1 - 5 . doi: 10.1109/lgrs.2021.3069761 http://dx.doi.org/10.1109/lgrs.2021.3069761 Robust optical and SAR image registration based on OS-SIFT and cascaded sample consensus FAN J W ， WU Y ， LI M ， et al . SAR and optical image registration using nonlinear diffusion and phase congruency structural descriptor ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2018 ， 56 （ 9 ）： 5368 - 5379 . doi: 10.1109/tgrs.2018.2815523 http://dx.doi.org/10.1109/tgrs.2018.2815523 SAR and optical image registration using nonlinear diffusion and phase congruency structural descriptor 保文星 ， 桑斯尔 ， 沈象飞 . 基于信息熵约束和KAZE特征提取的遥感图像配准算法研究 ［J］. 光学 精密工程 ， 2020 ， 28 （ 8 ）： 1810 - 1819 . 基于信息熵约束和KAZE特征提取的遥感图像配准算法研究 BAO W X ， SANG S E ， SHEN X F . Remote sensing image registration algorithm based on entropy constrained and KAZE feature extraction ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 8 ）： 1810 - 1819 . （in Chinese） Remote sensing image registration algorithm based on entropy constrained and KAZE feature extraction KOVESI P . Image Features from Phase Congruency ［J］. Viderei： Journal of Computer Vision Research ， 1999 ， 1 （ 3 ）： 1 - 26 . Image Features from Phase Congruency KOVESI P . Phase Congruency Detects Corners and Edges ［C］. The Australian Pattern recognition society conference. Sydney， Australia ： DICTA 2003 ： 309 - 318 . Phase Congruency Detects Corners and Edges 孙明超 ， 马天翔 ， 宋悦铭 ， 等 . 基于相位特征的可见光和SAR遥感图像自动配准 ［J］. 光学 精密工程 ， 2021 ， 29 （ 3 ）： 616 - 627 . doi: 10.37188/OPE.20212903.0616 http://dx.doi.org/10.37188/OPE.20212903.0616 基于相位特征的可见光和SAR遥感图像自动配准 SUN M C ， MA T X ， SONG Y M ， et al . Automatic registration of optical and SAR remote sensing image based on phase feature ［J］. Opt. Precision Eng. ， 2021 ， 29 （ 3 ）： 616 - 627 . （in Chinese） . doi: 10.37188/OPE.20212903.0616 http://dx.doi.org/10.37188/OPE.20212903.0616 Automatic registration of optical and SAR remote sensing image based on phase feature PAUL S ， PATI U C . Automatic optical-to-SAR image registration using a structural descriptor ［J］. IET Image Processing ， 2020 ， 14 （ 1 ）： 62 - 73 . doi: 10.1049/iet-ipr.2019.0389 http://dx.doi.org/10.1049/iet-ipr.2019.0389 Automatic optical-to-SAR image registration using a structural descriptor LI J ， HU Q ， AI M . RIFT： multi-modal image matching based on radiation-variation insensitive feature transform ［J］. IEEE Transactions on Image Processing： a Publication of the IEEE Signal Processing Society ， 2019 ： 2019 Dec17. doi: 10.1109/tip.2019.2959244 http://dx.doi.org/10.1109/tip.2019.2959244 RIFT： multi-modal image matching based on radiation-variation insensitive feature transform MORRONE M C ， OWENS R A . Feature detection from local energy ［J］. Pattern Recognition Letters ， 1987 ， 6 （ 5 ）： 303 - 313 . doi: 10.1016/0167-8655(87)90013-4 http://dx.doi.org/10.1016/0167-8655(87)90013-4 Feature detection from local energy DELLINGER F ， DELON J ， GOUSSEAU Y ， et al . SAR-SIFT： a SIFT-like algorithm for SAR images ［J］. IEEE Transactions on Geoscience and Remote Sensing ， 2015 ， 53 （ 1 ）： 453 - 466 . doi: 10.1109/tgrs.2014.2323552 http://dx.doi.org/10.1109/tgrs.2014.2323552 SAR-SIFT： a SIFT-like algorithm for SAR images WANG L N ， SUN M C ， LIU J H ， et al . A robust algorithm based on phase congruency for optical and SAR image registration in suburban areas ［J］. Remote Sensing ， 2020 ， 12 （ 20 ）： 3339 . doi: 10.3390/rs12203339 http://dx.doi.org/10.3390/rs12203339 A robust algorithm based on phase congruency for optical and SAR image registration in suburban areas WU Y ， MA W P ， GONG M G ， et al . A novel point-matching algorithm based on fast sample consensus for image registration ［J］. IEEE Geoscience and Remote Sensing Letters ， 2015 ， 12 （ 1 ）： 43 - 47 . doi: 10.1109/lgrs.2014.2325970 http://dx.doi.org/10.1109/lgrs.2014.2325970 A novel point-matching algorithm based on fast sample consensus for image registration 10.37188/OPE.20223014.1738 TP751.1 A 1004-924X（2022）14-1738-11 王丽娜,梁怀丹,王中石等.面向可见光和SAR影像配准的特征点检测[J].光学精密工程,2022,30(14):1738-1748. 王丽娜,梁怀丹,王中石等.面向可见光和SAR影像配准的特征点检测[J].光学精密工程,2022,30(14):1738-1748. WANG Lina,LIANG Huaidan,WANG Zhongshi,et al.Feature point detection for optical and SAR remote sensing images registration[J].Optics and Precision Engineering,2022,30(14):1738-1748. WANG Lina,LIANG Huaidan,WANG Zhongshi,et al.Feature point detection for optical and SAR remote sensing images registration[J].Optics and Precision Engineering,2022,30(14):1738-1748. 国家自然科学基金项目（No. 62103396）；吉林省自然科学基金项目（No.YDZJ202101ZYTS048）；吉林省教育厅科学研究项目（No.JJKH20210811KJ） 2021-11-24 2022-01-17 2022-07-25 2022-08-08 光学精密工程 14 30