基于光谱图像空间的F-SIFT特征提取与匹配

基于光谱图像空间的F-SIFT特征提取与匹配 Feature extraction and matching of F-SIFT based on spectral image space 1 2 first-author 丁国绅（1992-），男，安徽怀宁人，博士研究生，2018年于安徽大学获得硕士学位，主要研究方向为计算机视觉与遥感图像处理算法。E-mail：guoshenahu@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848779&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848787&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848783&type= 22.01330566 32.00401306 丁国绅 （1992-），男，安徽怀宁人，博士研究生，2018年于安徽大学获得硕士学位，主要研究方向为计算机视觉与遥感图像处理算法。E-mail： guoshenahu@163.com guoshenahu@163.com 1 乔延利（1958-），男，河南洛阳人，研究员，博士生导师，主要研究方向为遥感辐射定标与偏振遥感图像探测机理。E-mail：ylqiao@aiofm.ac.cn http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848790&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848795&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848792&type= 22.01330566 32.00401306 乔延利 （1958-），男，河南洛阳人，研究员，博士生导师，主要研究方向为遥感辐射定标与偏振遥感图像探测机理。E-mail： ylqiao@aiofm.ac.cn ylqiao@aiofm.ac.cn 1 3 1 corresp E-mail：maria.lily@163.com E-mail：maria.lily@163.com 杜丽丽（1982-），女，副研究员，2005年于郑州大学获学士学位，2010年于中国科学院研究生院获博士学位，主要研究方向为遥感器辐射和光谱定标。E-mail： maria.lily@163.com http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848798&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848805&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848801&type= 22.01330948 32.00379944 杜丽丽 （1982-），女，副研究员，2005年于郑州大学获学士学位，2010年于中国科学院研究生院获博士学位，主要研究方向为遥感器辐射和光谱定标。E-mail： maria.lily@163.com maria.lily@163.com 中国科学院 安徽光学精密机械研究所 通用光学定标与表征技术重点实验室， 安徽 合肥 230031 中国科学院 安徽光学精密机械研究所 通用光学定标与表征技术重点实验室， 安徽 合肥 230031 Key Laboratory of Optical Calibration and Characterization， Anhui Institute of Optics and Fine Mechanics， Chinese Academy of Sciences， Hefei 230031， China Key Laboratory of Optical Calibration and Characterization， Anhui Institute of Optics and Fine Mechanics， Chinese Academy of Sciences， Hefei 230031， China 中国科学技术大学， 安徽 合肥 230026 中国科学技术大学， 安徽 合肥 230026 University of Science and Technology of China， Hefei 230026， China University of Science and Technology of China， Hefei 230026， China 西安科技大学 安全科学与工程学院， 陕西 西安 710054 西安科技大学 安全科学与工程学院， 陕西 西安 710054 Xi’an University of Science and Technology， Xi’an 710054， China Xi’an University of Science and Technology， Xi’an 710054， China 针对传统SIFT算法提取到的图像特征点数量稀少的问题，可利用高光谱图像的特性，在光谱维度上取差重构了图像的尺度空间，使得提取到的特征点数量得到了极大地提高。但特征点数量的大幅增加导致算法的时间开销随之增大，并且有效特征点的占比较低。为了解决特征点的冗余问题，提高匹配效率，提出了一种基于光谱图像空间的F-SIFT算法。利用FAST算法能够在像素层面上快速判断的特性，构建了以当前像素为中心的八邻域准则，在提取图像的特征点之前对差分金字塔的对应位置的像素点做预筛选，使得特征点数量降低到了原来的10%以下；另外传统的匹配方法只统计了目标象元邻域内的像素信息，而忽略了象元的几何位置信息，因此本文扩展了特征描述符向量，首先利用最近邻与次近邻之比对特征描述符做一次粗匹配，记录可靠性程度并将其纳入描述符向量，接着按照相似性程度的高低从前20组匹配中迭代选取4组匹配点用于构造三角平面，利用象元的位置信息进行精确匹配。实验结果表明本文方法能够有效降低冗余特征点的数量，剔除误匹配。 To solve the problem of the scarcity of the number of image feature points extracted by the traditional scale-invariant feature transform （SIFT） algorithms， in the previous work， we used the characteristics of a hyperspectral image and reconstructed the scale space of the image by considering the difference in the spectral dimension， so that the number of feature points extracted is greatly increased. However， the large increase in the number of feature points leads to the increase in the time cost of the algorithm， and the proportion of effective feature points is low. To solve the redundancy problem of feature points and improve the matching efficiency， we propose a novel flip-SIFT （F-SIFT） algorithm based on the spectral image space. Based on the ability of the proposed algorithm to perform matching quickly at the pixel level， an eight-neighborhood criterion centered on the current pixel is constructed. Before the feature points of the image are extracted， the pixels in the corresponding position of the difference pyramid are pre-filtered， so that the number of feature points is reduced to less than one tenth of the original. In addition， the traditional matching method only counts the pixel information in the neighborhood of the target pixel， but ignores the geometric location information of the pixel. Therefore， in this study， we extended the feature descriptor vector. First， the feature descriptor is roughly matched by the comparison between the nearest neighbor and the next nearest neighbor， and the similarity degree is recorded and included in the descriptor vector. Then， according to the reliability degree， four groups of matching points are selected iteratively from the first 20 sets of matching points to construct the triangle plane and the pixel position information is used for accurate matching. Experimental results showed that the proposed method can effectively reduce the number of redundant feature points and eliminate mismatches. 尺度不变特征变换 光谱图像空间 八邻域准则 双重位置迭代匹配 scale-invariant feature transform spectral image space eight neighborhood criterion double position iterative matching 1 引　言 1 随着机器视觉的兴起，图像处理技术已经成为机器视觉不可缺少的组成部分。同时，图像中含有大量的信息，因此需要采用图像处理算法来提取图像的显著特征提高处理速度。图像特征是计算机视觉领域最重要的概念之一，从几何上可以将图像特征分为三大类：点特征、线特征和块特征 ［ 1 1 - 2 2 1-2 ］ 。基于不同的准则衍生出了诸多图像特征处理方法，并被广泛应用于目标检测、视频压缩、影像配准、三维模拟等领域 ［ 3 3 - 5 5 3-5 ］ 。由于传统的图像特征提取算法不能同时满足实时、高效、高精度的要求，学者们从不同的角度提出了许多优秀的特征处理算法。HOG （Histogram of Orientation gradient） ［ 6 6 ］ ，PCA （Principle component analysis） ［ 7 7 ］ 和Krawtchouk矩阵法 ［ 8 8 ］ 是几种典型的区域块特征提取算法，由于块中包含更多的图像细节信息，因此在平行空间下能够取得良好的效果，但现实世界中的很多应用场景中的图像并不是通过单个相机拍摄得到，而是基于多角度、多维度汇总处理，这种对尺度和视角敏感的缺陷在很大程度上限制了此类方法的拓展空间。线特征能够规避区域特征存在的一些缺陷，但其对场景的要求也较为严格，所提取的目标线段应尽可能不被遮挡、不变形、不断裂，线段的空间完整性和尺度一致性是保证其最佳性能的前提 ［ 9 9 ］ 。点特征对于图像的旋转、拼接、尺度等维度上的变换适应性较前两种特征其鲁棒性更优，常见的点特征提取方法包括SURF （Speed Up Robust Feature） ［ 10 10 ］ ，SIFT （Scale Invariant Feature Transform） ［ 11 11 ］ ，ORB （Oriented Fast and Rotated Brief） ［ 12 12 ］ 等。由于SIFT算法是最早提出的基于尺度不变性的方法，其在图像旋转、视点和仿射变换条件下依旧能保持优良的稳定性，使得该算法被广泛应用于工业和实际生产中，也吸引了大批研究者的关注，从不同的方向对算法的性能进行了提升。 Yan等在SIFT方法中引入了2DPCA算子重构低维度的描述子，将改进后的PCA方法用于归一化特征描述子向量，大大压缩了算法的运行时间 ［ 13 13 ］ 。Rathgeb等提出了基于SIFT的虹膜识别方法，该方法增强了从近红外虹膜图像中获得的虹膜纹理归一化的效果，并将二值化技术纳入到整个系统中，对提取到的特征点做了二值化，形成了二值SIFT描述子，相对传统SIFT方法，系统最终的识别精度有了很大提高 ［ 14 14 ］ 。Xu等建立了基于图像局部信息的自适应分数阶微分的数学模型，建立最优阶数与图像局部信息的关系，并根据图像的特点自动计算出每个像素点的最优阶数，接着根据最优阶数和Riemann-Liouville分数阶定义构造了自适应分数阶微分动态掩模，提出了基于自适应分数阶微分的SIFT算子 ［ 15 15 ］ 。 针对SIFT算法提取的特征点数量较少的问题，在上一步的工作中我们开创性地提出了利用高光谱图像不同波段象元的DN值差异重构图像金字塔的方法，使得提取到的特征点数量呈现了指数式增长 ［ 16 16 ］ 。但特征点数量的增加必然导致运行开销的增大，因此本文从两个方面对算法进行了优化。首先利用FAST算法在像素级上快速筛选的特性，建立八邻域判断准则，将特征点数量降低到了初始值的1/10以下，极大地降低了冗余特征点的数量；其次提出了双重位置迭代匹配方法，将描述符向量的扩展维进行排序，从可靠性程度最高的20组匹配中迭代选择4组构建双重三角平面，利用像素点的初始位置关系进行精确判断，实验结果表明没有误匹配出现。 2 传统SIFT算法 1 SIFT算法是1999年由British Columbia大学的David G.Lowe教授提出并于2004年加以完善的基于图像点特征的提取算法。该算法引入了尺度空间的概念，对光照、旋转、仿射以及尺度的变化都具有良好的稳定性，在学术界和工业界都产生了深远的影响。概括来说，SIFT算法可以分为三个步骤。 2.1　关键点检测 2 给定输入目标图像 <math id="M1"><mi>I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848850&type= 1.26999998 2.79399991 ，对 <math id="M2"><mi>I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848850&type= 1.26999998 2.79399991 分别进行一次升采样和 <math id="M3"><mi>O</mi><mo>-</mo><mn mathvariant="normal">1</mn></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848821&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848818&type= 9.65200043 2.87866688 次降采样，每次采样都以当前图像为基准做高斯变换。其中，高斯核函数 <math id="M4"><mi>G</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848829&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848826&type= 2.79399991 2.79399991 的表达式为： <math display="block" id="M5"><mi>G</mi><mo>=</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn><mi mathvariant="normal">π</mi><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac><mi mathvariant="normal">e</mi><mi mathvariant="normal">x</mi><mi mathvariant="normal">p</mi><mfenced separators="|"><mrow><mo>-</mo><mfrac><mrow><mo stretchy="false">(</mo><mi>x</mi><mo>-</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>t</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><mi>y</mi><mo>-</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>t</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mn mathvariant="normal">2</mn><msubsup><mrow><mi>σ</mi></mrow><mrow><mi>t</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msubsup></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848835&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848832&type= 61.72200012 9.39799976 ， 其中： <math id="M6"><msub><mrow><mi>x</mi></mrow><mrow><mi>t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848840&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848838&type= 2.79399991 3.72533321 ， <math id="M7"><msub><mrow><mi>y</mi></mrow><mrow><mi>t</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848848&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848843&type= 2.62466669 4.57200003 分别是图像 <math id="M8"><mi>I</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848853&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848850&type= 1.26999998 2.79399991 的第 <math id="M9"><mi>t</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848859&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848857&type= 1.10066664 2.79399991 个像素点的像素值， <math id="M10"><mi>σ</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848865&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848862&type= 1.69333339 2.79399991 是采样尺度，每次采样之后的图像 <math id="M11"><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>t</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>t</mi></mrow></msub><mo>,</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mi>G</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>t</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>t</mi></mrow></msub><mo>,</mo><msub><mrow><mi>σ</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">)</mo><mi mathvariant="normal">*</mi><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>t</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>t</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848873&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848868&type= 53.93266296 4.74133301 。在构建完图像的尺度空间后，利用高斯差分的性质建立图像差分空间，即用图像金字塔中的每一组图像中的第 <math id="M12"><mi>k</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848880&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848875&type= 1.69333339 2.79399991 层减去上一层即可得到。最后比较差分空间每组图像中除去首尾两层剩下的每层图像的各像素点与其相邻的26个像素点的差异，若该点在此27个像素点中是极大值或极小值，则认为该点是候选关键点。 2.2　关键点定位 2 通过差分空间检测可以得到一系列离散的极值点，但实际上关键点可能存在于两个极值点之间，因此这就需要通过插值的方式得到真正的关键点。根据泰勒函数的性质，对尺度空间函数进行泰勒拟合，得到亚象元级别的DOG函数，其表达式为： <math display="block" id="M13"><mi mathvariant="bold-italic">D</mi><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><msub><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mn mathvariant="normal">0</mn></mrow></msub><mo>+</mo><mfrac><mrow><mo>∂</mo><msup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi>T</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow></mfrac><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>+</mo><mfrac><mrow><mn mathvariant="normal">1</mn></mrow><mrow><mn mathvariant="normal">2</mn></mrow></mfrac><msup><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mi mathvariant="bold-italic">D</mi></mrow><mrow><mo>∂</mo><msup><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848886&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848883&type= 63.41533279 8.97466660 . 对 式（2） 式（2） 进行求导并令导数为零可得极值点位置 <math id="M14"><msub><mrow><mover><mrow><mi mathvariant="bold-italic">X</mi></mrow><mo>¯</mo></mover></mrow><mrow><mi>o</mi></mrow></msub><mo>=</mo><mo>-</mo><mfrac><mrow><msup><mrow><mo>∂</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><msup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><msup><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac><mfrac><mrow><mo>∂</mo><mi mathvariant="bold-italic">D</mi><mo stretchy="false">(</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow><mrow><mo>∂</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848893&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848889&type= 44.02666473 8.97466660 ，代入 式（2） 式（2） 可得插值结果 <math id="M15"><mi mathvariant="bold-italic">D</mi><mo stretchy="false">(</mo><mover><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow><mo>¯</mo></mover><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="bold-italic">D</mi><mo>+</mo><mn mathvariant="normal">0.5</mn><mfrac><mrow><mo>∂</mo><msup><mrow><mi mathvariant="bold-italic">D</mi></mrow><mrow><mi mathvariant="normal">T</mi></mrow></msup></mrow><mrow><mo>∂</mo><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow></mfrac></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848899&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848896&type= 35.56000137 8.72066593 ，若 <math id="M16"><mi mathvariant="bold-italic">D</mi><mo stretchy="false">(</mo><mover><mrow><msub><mrow><mi mathvariant="bold-italic">X</mi></mrow><mrow><mi>o</mi></mrow></msub></mrow><mo>¯</mo></mover><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848906&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848901&type= 11.09133339 4.23333359 大于设定的阈值则被保留，否则被剔除。最后通过计算主曲率的值，删除曲率阈值范围之外的极值点以降低边缘效应。 2.3　描述符生成 2 由于拍摄的角度不同，同一目标在不同图像中会存在不同程度的旋转以及仿射变换，因此在确定关键点的位置之后还需能唯一表征此描述符的方向信息。SIFT算法依据当前关键点的4×4邻域内的像素信息统计其像素梯度分布，其梯度值和对应的方向如下： <math display="block" id="M17"><mi>M</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mroot><mrow><mtable columnalign="left"><mtr><mtd><mo stretchy="false">(</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo></mtd></mtr><mtr><mtd><mo stretchy="false">(</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo><mo>-</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><msup><mrow><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mtd></mtr></mtable></mrow><mrow/></mroot></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848913&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848910&type= 73.23667145 11.00666618 ， <math display="block" id="M18"><mi>θ</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mi mathvariant="normal">t</mi><mi mathvariant="normal">a</mi><msup><mrow><mi mathvariant="normal">n</mi></mrow><mrow><mo>-</mo><mn mathvariant="normal">1</mn></mrow></msup><mfenced separators="|"><mrow><mfrac><mrow><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo><mo>-</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn><mo stretchy="false">)</mo></mrow><mrow><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>+</mo><mn mathvariant="normal">1</mn><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><mi>I</mi><mo stretchy="false">(</mo><msub><mrow><mi>x</mi></mrow><mrow><mi>o</mi></mrow></msub><mo>-</mo><mn mathvariant="normal">1</mn><mo>,</mo><msub><mrow><mi>y</mi></mrow><mrow><mi>o</mi></mrow></msub><mo stretchy="false">)</mo></mrow></mfrac></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848925&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848917&type= 74.33734131 9.73666668 . 对于邻域内的每个子区域，将其梯度分布平均划分为8个方向，统计每个方向区间内的梯度值之和，最终形成128维的梯度向量。取梯度值最高的方向为描述符的主方向，最后对所有梯度向量进行归一化处理，形成SIFT特征描述符。 3 基于光谱图像空间的F-SIFT特征提取与匹配 1 在本节中，优化并进一步完善了基于光谱图像空间的SIFT算法 ［ 16 16 ］ ，在特征点的数量上取得突破的同时降低冗余特征点的占比，极大地提高了特征点的有效利用率，并且采用独特的双重迭代优化的方法来剔除误匹配，整个方法的流程概括如下。 3.1　构建光谱图像空间 2 传统SIFT算法的尺度空间是由一张输入图像经过升采样和降采样并且对于同一分辨率的图像采用不同模糊系数进行高斯模糊得到不同模糊程度的子图像集合。这种做法能够解决两张图像因为尺度不同而难以配准的问题，同时对图像的仿射变换具有一定的鲁棒性，因此在后续诸多改进的SIFT方法中都保留了这一做法 ［ 13 13 - 14 14 13-14 ］ 。虽然SIFT算法在不同尺度的图像上稳定性较好，但其尺度空间仅由一幅图像经系列变换得到，所带有的有效图像特征信息十分有限，因而提取到的图像特征点数量较少。高光谱图像在不同波段上的图像具有良好的相似性，目标在不同波段图像上的位置相同而DN（Digital Number）值不同，也即传统意义上的灰度值不同，但其所表征的同一目标的DN值趋势相同。利用这个特性，用可见光范围内的高光谱图像作为尺度空间的基准图像，并对其进行直接降采样得到图像的尺度空间，示意图如 图1 图1 所示。 3.2　FAST八邻域准则 2 在重构了图像的尺度空间后，可提取到的图像特征点的数量较传统SIFT算法提升了几十至几百倍，克服了原SIFT算法提取的特征点数量少的缺陷。但通过实验验证发现，虽然有效特征点的数量有了很大提升，但相较于特征点总数的指数级规模的增长，冗余特征点的数量明显具有相当的规模，因而如何降低冗余特征点的占比是本文的研究内容之一。 FAST算法是特征点检测中最直观地用原始图像像素点与其邻域内像素点的差异值来判断是否为特征点的方法之一。虽然该方法检测到的特征点不一定是最优特征点，但其快速选择的特性在本文中被融合到了差分金字塔的判定中。原FAST算法是在较大的邻域空间上对目标位置的像素值进行比较，虽然大部分的情况下像素值的变化趋势相似，但在像素值有明显突变的时候邻域范围较小更能较好地表征目标像素点的周围像素变化。因此本文采用了3+2+3的目标像素八邻域的方式，在构造差分金字塔中的每一层图像时对该图像上的像素做预筛选，若邻域内有 M 个点的像素值大于阈值 T 则保留该点的像素值，否则将该点的像素值置为 - 1。在特征点的判断阶段则首先判断该点的像素值是否为 - 1，若是则直接进入下一个特征点的判断阶段。经实验验证，加入FAST八邻域准则之后，F-SIFT算法所提取到的特征点数量降到了文献［ 16 16 ］中算法的1/10以下，且最终成功匹配的特征点数量不受大的影响。 3.3　双重位置迭代匹配 2 在完成特征点的提取之后要对两幅图像之间的特征点进行匹配，常见的方法是采用最近邻与次近邻的比值来选取最优匹配。这种方法是基于像素梯度值的分布来衡量特征点的相似程度的，图像中往往会存在一个或多个相似的特征点，这就容易造成误匹配现象，而且该方法对于比值阈值的确定异常敏感。传统SIFT算法在进行特征点匹配时常将比值设定为0.8，涉及到具体图像时还会有小范围的波动。这是因为当比值增大时匹配到的特征点数量会随之增加，有效匹配和误匹配数量也会增大；当比值减小时匹配到的特征点数量会随之降低，虽然误匹率降低但成功匹配的数量也明显衰减。 针对这种情况，本文提出了双重位置迭代匹配方法，将整个特征点的匹配过程分为两步。首先采用最近邻与次近邻之比对特征点进行一次粗匹配，不过在此阶段的比值设定较大，目的是为了尽可能多地包含有效匹配对，经实验验证此阶段的理想比值 H 设定为1.1。接着进行精细匹配操作，在粗匹配阶段依据最近邻与次近邻的相似性对特征点的描述符向量进行了一维扩展，虽然这一阶段的匹配中包含了大量的误匹配，但精细匹配的作用就是剔除这些误匹配。将匹配到的特征点对按照相似性程度进行排序，然后从相似性程度最高的20组匹配对中随机选择4组构建源-目标图像双三角平面。在之前的工作中选取的是相似性程度最高的3组匹配对来构建三角平面，但经过后来的实验发现这种情况下依旧可能出现误匹配的现象，因为只有在没有噪声的理想情况下才能通过三个点来构建一个完整的平面。为了排除噪声的干扰，本文采用了4组匹配对来构建双三角平面，即用排列组合的方式随机从4组中选取两个三维组合用来构建两个三角平面，最后根据所有匹配对的特征点在双三角平面中的位置关系进行是否为正确匹配对的判定。另一方面，没有采用相似性程度最高的4组匹配对是因为图像上会存在相似的特征点，将其作为基准点构建三角平面可能会包含更多的正确匹配信息，因而在一定的区间上取出4组匹配对会更有益。在兼顾算法效率的基础上将此区间设定为1~20。由于图像是用二维来映射三维空间位置关系的，因此用二维平面上的位置关系来阐述双三角平面的构建原理。 如 图2 图2 所示，特征点集合 <math id="M19"><mo stretchy="false">(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>D</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848947&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848943&type= 19.55800056 3.55599999 与 <math id="M20"><mo stretchy="false">(</mo><mi>A</mi><mi mathvariant="normal">'</mi><mo>,</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848958&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848953&type= 7.19666624 3.89466691 <math id="M21"><mi>B</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>C</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>D</mi><mi mathvariant="normal">'</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848966&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848962&type= 15.23999977 3.89466691 分别是源-目标图像中选取的4对基准特征点对，对应关系分别是： <math id="M22"><mi>A</mi><mo>↔</mo><mi>A</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848975&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848970&type= 11.09133339 3.30200005 ， <math id="M23"><mi>B</mi><mo>↔</mo><mi>B</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848985&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848979&type= 10.49866676 3.30200005 ， <math id="M24"><mi>C</mi><mo>↔</mo><mi>C</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848994&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848989&type= 10.75266743 3.30200005 ， <math id="M25"><mi>D</mi><mo>↔</mo><mi>D</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849003&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16848998&type= 11.17599964 3.30200005 ， <math id="M26"><mo>↔</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849012&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849008&type= 3.72533321 2.79399991 表示匹配关系。 <math id="M27"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849364&type= 2.87866688 3.21733332 点是源图像中的某个特征点，假设 <math id="M28"><msup><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849031&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849027&type= 3.38666677 3.47133350 是 <math id="M29"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849364&type= 2.87866688 3.21733332 点在目标图像中匹配到的特征点。因为 <math id="M30"><mo stretchy="false">(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>D</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849047&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849043&type= 19.55800056 3.55599999 与 <math id="M31"><mo stretchy="false">(</mo><mi>A</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>B</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>C</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>D</mi><mi mathvariant="normal">'</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849058&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849051&type= 22.35199928 3.89466691 是基准特征点对，故三角平面 <math id="M32"><mi>O</mi><mo stretchy="false">(</mo><mi>A</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>B</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>C</mi><mi mathvariant="normal">'</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849208&type= 20.48933220 3.89466691 可由 <math id="M33"><mi>O</mi><mo stretchy="false">(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849200&type= 18.45733452 3.55599999 经系列平面变换操作生成。直线 <math id="M34"><mi>L</mi><mo stretchy="false">(</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849081&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849077&type= 13.03866673 3.55599999 中存在如下等式： <math display="block" id="M35"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mo stretchy="false">(</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo><mi>x</mi><mtext> </mtext><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mi>y</mi><mtext> </mtext><mo>+</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mtext> </mtext><mo>-</mo></mtd></mtr><mtr><mtd><mtext>             </mtext><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr><mtr><mtd><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo><mi>x</mi><mtext> </mtext><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mi>y</mi><mtext> </mtext><mo>+</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mtext> </mtext><mo>-</mo></mtd></mtr><mtr><mtd><mtext>              </mtext><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo><mo>=</mo><mn mathvariant="normal">0</mn></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849087&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849084&type= 73.15200043 21.42066574 . 设 <math id="M36"><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849109&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849106&type= 21.67466736 4.14866638 ， <math id="M37"><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849120&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849113&type= 21.33600044 4.14866638 ， <math id="M38"><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>=</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849128&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849124&type= 53.84799957 4.48733330 ， <math id="M39"><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849137&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849132&type= 20.99733353 3.72533321 ， <math id="M40"><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849145&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849141&type= 19.47333336 3.72533321 ， <math id="M41"><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>=</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">)</mo><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849156&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849151&type= 49.86866760 4.48733330 ，可得点 <math id="M42"><mi>P</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849163&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849159&type= 2.45533323 2.79399991 坐标： <math display="block" id="M43"><mfenced open="{" close="" separators="|"><mrow><mtable><mtr><mtd><msub><mrow><mi>P</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>P</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>=</mo><mo stretchy="false">(</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo>-</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>T</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo stretchy="false">)</mo><mo>/</mo><mo stretchy="false">(</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><mo>-</mo><msub><mrow><mi>R</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub><msub><mrow><mi>S</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub><mo stretchy="false">)</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849171&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849167&type= 58.08132935 10.32933331 . 设 <math id="M44"><mi>M</mi><mo>=</mo><mo stretchy="false">|</mo><mi>B</mi><mi>P</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>B</mi><mi>C</mi><mo stretchy="false">|</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849181&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849176&type= 22.94466591 3.64066648 ， <math id="M45"><mi>N</mi><mo>=</mo><mo stretchy="false">|</mo><mi>A</mi><mi>P</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>A</mi><mi>Q</mi><mo stretchy="false">|</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849189&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849185&type= 23.45266724 3.64066648 ，有如下等式成立： <math display="block" id="M46"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mi>M</mi><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>x</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>y</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mtd></mtr><mtr><mtd><mi>N</mi><mo>=</mo><mfrac><mrow><mo stretchy="false">(</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>x</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>P</mi></mrow><mrow><mi>y</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow><mrow><mo stretchy="false">(</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup><mo>+</mo><mo stretchy="false">(</mo><msub><mrow><mi>A</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><msup><mrow><mo stretchy="false">)</mo></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msup></mrow></mfrac></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849197&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849193&type= 49.44533539 20.40466690 . 因为平面 <math id="M47"><mi>O</mi><mo stretchy="false">(</mo><mi>A</mi><mo>,</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849204&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849200&type= 18.45733452 3.55599999 与 <math id="M48"><mi>O</mi><mo stretchy="false">(</mo><mi>A</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>B</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>C</mi><mi mathvariant="normal">'</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849212&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849208&type= 20.48933220 3.89466691 中各点的相对位置是相同的，因而存在如下关系： <math display="block" id="M49"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><mo stretchy="false">|</mo><mi>B</mi><mi>P</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>B</mi><mi>C</mi><mo stretchy="false">|</mo><mo>=</mo><mo stretchy="false">|</mo><mi>B</mi><mi mathvariant="normal">'</mi><mi>P</mi><mi mathvariant="normal">'</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>B</mi><mi mathvariant="normal">'</mi><mi>C</mi><mi mathvariant="normal">'</mi><mo stretchy="false">|</mo></mtd></mtr><mtr><mtd><mo stretchy="false">|</mo><mi>A</mi><mi>P</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>A</mi><mi>Q</mi><mo stretchy="false">|</mo><mo>=</mo><mo stretchy="false">|</mo><mi>A</mi><mi mathvariant="normal">'</mi><mi>P</mi><mi mathvariant="normal">'</mi><mo stretchy="false">|</mo><mo>/</mo><mo stretchy="false">|</mo><mi>A</mi><mi mathvariant="normal">'</mi><mi>Q</mi><mi mathvariant="normal">'</mi><mo stretchy="false">|</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849222&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849216&type= 39.62400055 10.32933331 . 求得 <math id="M50"><mi>P</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849230&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849226&type= 3.13266683 3.30200005 坐标： <math display="block" id="M51"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>P</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>=</mo><mi>M</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>C</mi></mrow><mrow><mi>x</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>x</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></mtd></mtr><mtr><mtd><msub><mrow><mi>P</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>=</mo><mi>M</mi><mo stretchy="false">(</mo><msubsup><mrow><mi>C</mi></mrow><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo>-</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup><mo stretchy="false">)</mo><mo>+</mo><msubsup><mrow><mi>B</mi></mrow><mrow><mi>y</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849237&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849234&type= 37.33800125 10.15999985 . 进而可得到理想情况下的 <math id="M52"><mi>Q</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849381&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849376&type= 3.55599999 3.72533321 点的坐标，接着计算 <math id="M53"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849364&type= 2.87866688 3.21733332 点实际匹配到的点 <math id="M54"><msubsup><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849295&type= 7.70466709 4.40266657 与 <math id="M55"><mi>Q</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849381&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849376&type= 3.55599999 3.72533321 的位置误差 <math id="M56"><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849345&type= 3.13266683 3.72533321 ： <math display="block" id="M57"><mfenced open="{" close="" separators="|"><mrow><mtable columnalign="left"><mtr><mtd><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>,</mo><mi>x</mi></mrow></msub><mo>=</mo><mo stretchy="false">|</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>x</mi></mrow></msub><mo>-</mo><mi>Q</mi><msub><mrow><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow><mrow><mi>x</mi></mrow></msub><mo stretchy="false">|</mo></mtd></mtr><mtr><mtd><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">1</mn><mo>,</mo><mi>y</mi></mrow></msub><mo>=</mo><mo stretchy="false">|</mo><msub><mrow><mi>Q</mi></mrow><mrow><mi>y</mi></mrow></msub><mo>-</mo><mi>Q</mi><msub><mrow><msub><mrow><mi mathvariant="normal">'</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow></msub></mrow><mrow><mi>y</mi></mrow></msub><mo stretchy="false">|</mo></mtd></mtr></mtable></mrow></mfenced></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849290&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849285&type= 31.83466530 10.92199993 . 同理可得到 <math id="M58"><msubsup><mrow><mi>Q</mi></mrow><mrow><mi mathvariant="normal">T</mi><mi mathvariant="normal">r</mi><mi mathvariant="normal">u</mi><mi mathvariant="normal">t</mi><mi mathvariant="normal">h</mi></mrow><mrow><mi mathvariant="normal">'</mi></mrow></msubsup></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849300&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849295&type= 7.70466709 4.40266657 与 <math id="M59"><mi>Q</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849381&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849376&type= 3.55599999 3.72533321 在平面 <math id="M60"><mi>O</mi><mo stretchy="false">(</mo><mi>B</mi><mo>,</mo><mi>C</mi><mo>,</mo><mi>D</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849319&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849313&type= 18.28800011 3.55599999 与 <math id="M61"><mi>O</mi><mo stretchy="false">(</mo><mi>B</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>C</mi><mi mathvariant="normal">'</mi><mo>,</mo><mi>D</mi><mi mathvariant="normal">'</mi><mo stretchy="false">)</mo></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849329&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849324&type= 20.31999969 3.89466691 中的误差 <math id="M62"><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849359&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849354&type= 3.13266683 3.72533321 ，若 <math id="M63"><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">1</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849350&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849345&type= 3.13266683 3.72533321 与 <math id="M64"><msub><mrow><mi>δ</mi></mrow><mrow><mn mathvariant="normal">2</mn></mrow></msub></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849359&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849354&type= 3.13266683 3.72533321 均在允许误差范围内则认为 <math id="M65"><mi>Q</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849371&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849364&type= 2.87866688 3.21733332 与 <math id="M66"><mi>Q</mi><mi mathvariant="normal">'</mi></math> http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849381&type= http://html.publish.founderss.cn/rc-pub/api/common/picture?pictureId=16849376&type= 3.55599999 3.72533321 是一对正确匹配对，否则剔除这对匹配。 4 实验结果与分析 1 因为基于高光谱图像的SIFT算法是由我们首次提出的，并未出现此方面的公用数据集，结合实验需求，本文所用数据集都是基于IS210-0.4-1.0-L型号的高光谱成像相机拍摄的，数据集链接为： https：//pan.baidu.com/s/1qYEymTvj_wjVZEoj6wLfvQ https://pan.baidu.com/s/1qYEymTvj_wjVZEoj6wLfvQ 。实验的硬件环境为Windows10（64bit），Intel（R） Core（TM） i7-9700CPU @ 3.00 GHz处理器，编程环境为MatlabR2020a。本次实验分为两组，在可见光波段范围内选取的波段间隔为10~50，精细匹配阶段的迭代次数设置为1 000，验证F-SIFT算法对于简单和复杂目标情况下的基于不同波段间隔的算法性能。 图3 图3 和 图4 图4 分别是传统SIFT算法与F-SIFT算法在单一目标和复杂目标场景下的图像特征点匹配结果，图像波段间隔为25，其中图 3（a） 3（a） ~ 3（b） 3（b） 图是传统SIFT算法在最近邻与次近邻的比值为0.7和1.1时得到的，图 3（c） 3（c） ， 3（d） 3（d） 图时F-SIFT算法在最近邻与次近邻的比值为0.7和1.1时得到的。由图中可以看出，对于传统SIFT，在最近邻与次近邻的比值设定较小时算法提取到的特征点数量相对较少，误匹率较低；增大比值能够较大提升特征点提取能力，但误匹率也随之增加，这种结果往往是不能忍受的。而对于F-SIFT算法，虽然在最近邻与次近邻的比值较小时算法在某些场景下提取到的特征点数量有限，但能保证剔除所有误匹配，而且随着比值的增加这种零误匹率的优势更加明显，在比值为1.1时算法能同时兼顾特征点数量与正确匹配率两大优势，具有良好的应用潜力。 1,2;3,4; 1,2;3,4; 表1 表1 和 表2 表2 分别记录了在单一和复杂目标场景下传统SIFT算法与F-SIFT算法在不同最近邻与次近邻比值的情况下特征点匹配对数。其中SIFT-O表示传统SIFT算法只在最近邻与次近邻比值规则下得到的匹配对数目（包含不同程度的误匹配数目），SIFT-F表示在双重位置迭代匹配准则下传统SIFT算法得到的匹配对数目（不包含误匹配）。由表中数据可以看出F-SIFT算法较传统SIFT算法其有效特征点匹配对数提升了几十至上百倍，而且由于在匹配阶段采用的双重位置迭代匹配准则，这些匹配对中均没有误匹配。而且随着最近邻与次近邻的比值不断增大，F-SIFT算法的性能优势更为突出。另外，从波段间隔的变化上可以发现，在用高光谱图像重构尺度空间时，并不是波段间隔越大或越小算法的性能越好，这从图像层面上也很好理解。在波段间隔过小时，相邻波段的图像之间的差异很小，提取到的特征点信息就更接近，也即多样化信息就越少，因而能提取到的特征点数量就越少；在波段间隔过大时，不同波段之间的图像差异太大会对差分金字塔的特征点的确定造成负面的影响，因而也会导致提取到的特征点数量降低。 2 Dist-Ratio 2 SIFT-O 2 SIFT-F 9 Band-Interval（F-SIFT） 10 15 20 25 30 35 40 45 50 0.7 30 2 19 19 29 18 35 15 19 8 22 0.8 96 3 40 62 72 69 97 43 40 17 59 0.9 457 4 104 151 190 154 225 94 128 32 96 1.0 2 186 5 236 320 436 335 502 200 343 71 167 1.1 2 192 5 237 323 490 406 465 213 354 71 178 2 Dist-Ratio 2 SIFT-O 2 SIFT-F 9 Band-Interval（F-SIFT） 10 15 20 25 30 35 40 45 50 0.7 128 108 110 484 841 1 049 972 742 486 461 332 0.8 139 111 196 657 1 192 1 501 1 354 1 027 642 584 410 0.9 157 112 305 911 1 829 2 158 1 940 1 464 852 765 518 1.0 266 116 513 1 207 2 688 3 054 2 673 1 943 1 098 954 679 1.1 266 117 513 1 290 2 691 3 060 2 787 1 949 1 102 1 021 681 图5 图5 （a）表示F-SIFT算法与文献［ 16 16 ］中算法在不同波段下提取到的特征点数量曲线，纵轴的单位为万，可以看到相较于文献［ 16 16 ］，本文算法提取的特征点数量下降极为明显，这主要得益于FAST八邻域准则的加入，有效剔除了大量的冗余特征点。 1,2; 图5 图5 （b）表示F-SIFT算法与文献［ 16 16 ］中算法在不同波段下得到的最终正确匹配对数，由于文献［ 16 16 ］中可能含有误匹配，因而采用了相似性程度最高的4组匹配对用来构建双三角平面。由此也可以看出相似性程度越高并不一定就意味着是最佳的选择，经实验反复验证，从前20对相似性程度高的匹配对中随机选择基准特征点对能扩大算法的性能优势。 5 结　论 1 图像处理是机器视觉领域最基础的研究工作之一，特征点提取与匹配在图像处理中有着举足轻重的位置，良好的特征点获取与检校方法对于后期的工作具有重要的意义。因此本文从SIFT算法出发，对先前的工作进行了优化和总结，提出了一种基于光谱图像空间的F-SIFT特征提取与匹配方法。该方法从图像的尺度空间的构建、冗余特征点的删减、误匹配的剔除三个方面进行了改良，构造了一种新颖有效的图像特征点处理方法。接下来准备从高光谱图像本身出发，研究不同分辨率下的高光谱图像构造出的尺度空间对算法最终结果的影响，以及将本文方法引入图像的三维重构中，拓展算法的应用领域。 参考文献 BROWN M ， LOWE D G . Automatic panoramic image stitching using invariant features ［J］. International Journal of Computer Vision ， 2007 ， 74 （ 1 ）： 59 - 73 . Automatic panoramic image stitching using invariant features GUO C ， JIA F ， TANG W ， et al . A fast method for image matching and registration based on SIFT algorithm and image pyramid ［J］. Journal of Physics： Conference Series ， 2020 ， 1449 （ 1 ）： 012119 . A fast method for image matching and registration based on SIFT algorithm and image pyramid SHU Y ， TAN Z . 3D reconstruction based on spatial vanishing information ［J］. Chinese Optics Letters ， 2005 ， 3 （ 3 ）： 146 - 148 . 3D reconstruction based on spatial vanishing information ZHIXIN Z ， SHUHAO J . Design of incomplete 3D information image recognition system based on SIFT algorithm and wireless network ［J］. EURASIP Journal on Wireless Communications and Networking ， 2020 （ 1 ）. Design of incomplete 3D information image recognition system based on SIFT algorithm and wireless network WANG Y ， YUAN Y ， LEI Z . Fast SIFT feature matching algorithm based on geometric transformation （January 2020） ［J］. IEEE Access ， 2020 ， （ 99 ）： 1 - 1 . Fast SIFT feature matching algorithm based on geometric transformation （January 2020） 王广龙 ， 田杰 ， 朱文杰 ， 等 . 特征融合和自适应权重更新相结合的运动模糊目标跟踪 ［J］. 光学 精密工程 ， 2019 ， 27 （ 5 ）： 1158 - 1166 . 特征融合和自适应权重更新相结合的运动模糊目标跟踪 WANG G L ， TIAN J ， ZHU W J ， et al . Feature fusion and weight adaptive updating based motion blur object tracking ［J］. Opt. Precision Eng. ， 2019 ， 27 （ 5 ）： 1158 - 1166 . （in Chinese） Feature fusion and weight adaptive updating based motion blur object tracking KAHLON C S ， LI B ， BOARD J ， et al . Cluster and principle component analysis of soybean grown at various row spacings， planting dates and plant populations ［J］. Open Agriculture ， 2018 ， 3 （ 1 ）： 110 - 121 . Cluster and principle component analysis of soybean grown at various row spacings， planting dates and plant populations YANG T ， MA J ， MIAO Y ， et al . Privacy-preserving krawtchouk moment feature extraction over encrypted image data ［J］. Information ences ， 2020 ， 536 . Privacy-preserving krawtchouk moment feature extraction over encrypted image data 韩亚君 ， 杨德东 ， 李勇 ， 等 . 基于多特征融合相关滤波的红外目标跟踪 ［J］. 液晶与显示 ， 2019 ， 34 （ 2 ）： 177 - 187 . 基于多特征融合相关滤波的红外目标跟踪 HAN Y J ， YANG D D ， LI Y ， et al . Infrared target tracking based on correlation filter with multi-features fusion ［J］. Chinese Journal of Liquid Crystals and Displays ， 2019 ， 34 （ 2 ）： 177 - 187 . （in Chinese） Infrared target tracking based on correlation filter with multi-features fusion HE Q ， HE B ， ZHANG Y ， et al . Multimedia based fast face recognition algorithm of speed up robust features ［J］. Multimedia Tools and Applications ， 2019 ， 78 （ 17 ）： 24035 - 24045 . Multimedia based fast face recognition algorithm of speed up robust features LOWE D G . Distinctive Image Features from Scale-Invariant Keypoints ［J］. International Journal of Computer Vision ， 2004 ， 60 （ 2 ）： 91 - 110 . Distinctive Image Features from Scale-Invariant Keypoints RUBLEE E ， RABAUD V ， KONOLIGE K ， et al . ORB： An efficient alternative to SIFT or SURF ［C］. International Conference on Computer Vision . IEEE ， 2012 . ORB： An efficient alternative to SIFT or SURF 颜雪军 ， 赵春霞 ， 袁夏 . 2DPCA-SIFT：一种有效的局部特征描述方法 ［J］. 自动化学报 ， 2014 （ 4 ）： 675 - 682 . 2DPCA-SIFT：一种有效的局部特征描述方法 YAN X J ， ZHAO CH X ， YUAN X . 2DPCA-SIFT：An efficient local feature descriptor ［J］. Acta Automatica Sinica ， 2014 ， 40 （ 4 ）： 675 - 682 . （in Chinese） 2DPCA-SIFT：An efficient local feature descriptor RATHGEB C ， WAGNER J ， BUSCH C . SIFT-based iris recognition revisited： prerequisites， advantages and improvements ［J］. Pattern Analysis and Applications ， 2019 ， 22 （ 3 ）： 889 - 906 . SIFT-based iris recognition revisited： prerequisites， advantages and improvements XU K ， LIU J X ， MIAO J ， et al . An improved SIFT algorithm based on adaptive fractional differential ［J］. Journal of Ambient Intelligence and Humanized Computing ， 2019 . An improved SIFT algorithm based on adaptive fractional differential 丁国绅 ， 乔延利 ， 易维宁 ， 等 . 基于高光谱图像的改进SIFT特征提取与匹配 ［J］. 光学 精密工程 ， 2020 ， 28 （ 4 ）： 954 - 962 . 基于高光谱图像的改进SIFT特征提取与匹配 DING G SH ， YI W N ， DU Lili ， et al . Improved SIFT feature extraction and matching technology based on hyperspectral image ［J］. Opt. Precision Eng. ， 2020 ， 28 （ 4 ）： 954 - 962 . （in Chinese） Improved SIFT feature extraction and matching technology based on hyperspectral image 10.37188/OPE.20212905.1180 TP391.41 A 1004-924X（2021）05-1180-10 丁国绅,乔延利,易维宁等.基于光谱图像空间的F-SIFT特征提取与匹配[J].光学精密工程,2021,29(05):1180-1189. 丁国绅,乔延利,易维宁等.基于光谱图像空间的F-SIFT特征提取与匹配[J].光学精密工程,2021,29(05):1180-1189. DING Guo-shen,QIAO Yan-li,YI Wei-ning,et al.Feature extraction and matching of F-SIFT based on spectral image space[J].Optics and Precision Engineering,2021,29(05):1180-1189. DING Guo-shen,QIAO Yan-li,YI Wei-ning,et al.Feature extraction and matching of F-SIFT based on spectral image space[J].Optics and Precision Engineering,2021,29(05):1180-1189. 国家自然科学基金资助项目（No.41601379） 2020-11-04 2020-12-15 2021-05-15 2021-06-24 光学精密工程 5 29