CN104992431B - 多光谱图像配准的方法及装置 - Google Patents
多光谱图像配准的方法及装置 Download PDFInfo
- Publication number
- CN104992431B CN104992431B CN201510345294.6A CN201510345294A CN104992431B CN 104992431 B CN104992431 B CN 104992431B CN 201510345294 A CN201510345294 A CN 201510345294A CN 104992431 B CN104992431 B CN 104992431B
- Authority
- CN
- China
- Prior art keywords
- intersection point
- straight line
- intersection
- mtd
- multispectral image
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Active
Links
- 238000000034 method Methods 0.000 title claims abstract description 27
- 230000009466 transformation Effects 0.000 claims description 55
- 239000011159 matrix material Substances 0.000 claims description 20
- 238000006243 chemical reaction Methods 0.000 abstract 2
- 238000010586 diagram Methods 0.000 description 2
- 230000008569 process Effects 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 239000003814 drug Substances 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000012986 modification Methods 0.000 description 1
- 230000004048 modification Effects 0.000 description 1
- 230000003287 optical effect Effects 0.000 description 1
- 230000004044 response Effects 0.000 description 1
- 238000011524 similarity measure Methods 0.000 description 1
- 238000006467 substitution reaction Methods 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/10—Image acquisition modality
- G06T2207/10032—Satellite or aerial image; Remote sensing
- G06T2207/10036—Multispectral image; Hyperspectral image
Landscapes
- Image Analysis (AREA)
Abstract
本发明提供一种多光谱图像配准的方法及装置。本发明多光谱图像配准的方法，包括：获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；确定所述参考多光谱图像上与所述第一交点匹配的第二交点；根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准。本发明提高了多光谱图像配准的准确度。
Description
技术领域
本发明涉及图像处理技术领域，具体涉及一种多光谱图像配准的方法及装置。
背景技术
图像配准(Image registration)是指依据一些相似性度量决定图像间的变换参数，使从不同探测器、不同视角、不同时间获取的同一场景的两幅或多幅图像，变换到同一坐标系下，在像素层上得到最佳匹配的过程，它是图像融合、目标识别、目标变化检测、计算机视觉等过程的一个重要前期步骤，在军事、遥感、医学等领域有着广泛的应用。
现有的一种图像配准方法是采用尺度不变特征转换(Scale-invariant featuretransform，SIFT)方法，具体步骤如下：首先生成图像尺度空间，然后检测尺度空间中的局部极值点，再通过剔除局部极值点中的低对比度点和边缘响应点，得到特征点；对特征点进行描述时，对以特征点为中心的区域进行梯度方向直方图统计，得到特征点的方向分配，并生成特征描述子；最后通过特征描述子寻找图像之间匹配的特征点对，然后通过匹配的特征点对得到图像空间坐标变换参数，最后根据坐标变换参数进行图像配准。
但是上述方法用于多光谱图像配准时，由于多光谱图像的内容比较复杂，而且像素值之间没有对应关系，纹理结构信息有很大差别。因而，根据上述方法确定出的特征点不稳定，不具有代表性，通过特征点进行配准会导致配准的准确度较低。
发明内容
本发明提供一种多光谱图像配准的方法及装置，能够解决现有技术中配准的准确度较低的问题。
第一方面，本发明提供一种多光谱图像配准的方法，包括：
获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准。
可选地，所述确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体包括：
分别对所述第一交点和所述第二交点计算边缘方向直方图EOH描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点。
可选地，所述对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配之前，还包括：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔所述预设距离选取第二特征点；所述第二特征点包括所述第一交点。
可选地，所述对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体包括：
计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；所述第一特征点与所述第二特征点的位置对应；
若有超过预设的数量阈值的所述欧氏距离小于所述预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；所述第三距离阈值小于所述第二距离阈值。
可选地，所述根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数，包括：
根据如下的公式(1)确定坐标变换矩阵H，将所述坐标变换矩阵H作为坐标变换参数；
其中，为所述第二交点的坐标；为所述第一交点的坐标；坐标变换矩阵
可选地，所述根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，包括：
将所述待配准多光谱图像中的像素点的坐标与所述坐标变换矩阵相乘，使得所述待配准多光谱图像变换到与所述参考多光谱图像相同的坐标系下。
第二方面，本发明实施例通过一种多光谱图像配准的装置，包括：
获取模块，用于获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
确定模块，用于确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
所述确定模块，还用于根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
处理模块，用于根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准。
可选地，所述确定模块，具体用于：
分别对所述第一交点和所述第二交点计算边缘方向直方图EOH描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点。
可选地，所述处理模块，还用于：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔所述预设距离选取第二特征点；所述第二特征点包括所述第一交点。
可选地，所述确定模块，具体用于：
计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；所述第一特征点与所述第二特征点的位置对应；
若有超过预设的数量阈值的所述欧氏距离小于所述预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；所述第三距离阈值小于所述第二距离阈值。
本发明提供的多光谱图像配准的方法及装置，通过获取待配准多光谱图像和参考多光谱图像中多组相交的直线对的交点，将直线的交点作为特征点进行匹配，确定匹配的交点对，并根据确定出的匹配的交点对确定坐标变换参数，最后根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，与现有技术相比，由于直线特征在多光谱图像中具有更好的稳定性，因此选择直线的交点作为特征点进行匹配提高了多光谱图像配准的准确度。
附图说明
为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例或现有技术描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其它的附图。
图1为本发明多光谱图像配准的方法一实施例的流程图；
图2为本发明多光谱图像配准的装置一实施例的结构示意图。
具体实施方式
下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有作出创造性劳动前提下所获得的所有其它实施例，都属于本发明保护的范围。
图1为本发明多光谱图像配准的方法一实施例的流程图。如图1所示，本实施例的方法，包括：
步骤101、获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
步骤102、确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
步骤103、根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
步骤104、根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准。
具体来说，进行多光谱图像配准时，选择至少两幅多光谱图像，一幅作为参考多光谱图像，另外的至少一幅作为待配准多光谱图像；获取至少两幅多光谱图像中的直线，具体可以通过Hough变换检测直线；然后根据直线的方程得到多组相交的直线对，以及交点；待配准多光谱图像中的相交的直线对的交点作为第一交点，参考多光谱图像中的相交的直线对的交点作为第二交点。
确定所述参考多光谱图像上与所述第一交点匹配的第二交点，例如本发明实施例中可以通过计算第一交点的EOH描述符与第二交点的EOH描述符之间的欧氏距离实现，上述确定与第一交点匹配的第二交点，需对参考多光谱图像上的所有第二交点进行遍历。
根据所述第一交点以及与所述第一交点匹配的第二交点，即确定出多组交点匹配对，根据该多组交点匹配对之间的坐标关系确定待配准多光谱图像的坐标变换参数；根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，即根据所述坐标变换参数将所述待配准多光谱图像进行坐标变换，实现所述待配准多光谱图像和所述参考多光谱图像配准。
本实施例提供的多光谱图像配准的方法，通过获取待配准多光谱图像和参考多光谱图像中多组相交的直线对的交点，将直线的交点作为特征点进行匹配，确定匹配的交点对，并根据确定出的匹配的交点对确定坐标变换参数，最后根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，与现有技术相比，由于直线特征在多光谱图像中具有更好的稳定性，因此选择直线的交点作为特征点进行匹配提高了多光谱图像配准的准确度。
在图1所示的实施方式的基础上，进一步的，在实际应用中，确定匹配的交点的方式可以有多种，可选的，作为一种可实施的方式，确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体可以采用如下方式实现：
对所述第一交点和所述第二交点计算边缘方向直方图(Edge OrientedHistogram，简称EOH)描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点。
具体来说，由于需要确定所述参考多光谱图像上与所述第一交点匹配的第二交点，本发明实施例中可以根据交点的EOH描述符之间的欧氏距离来确定匹配的第一交点和第二交点，因此需在匹配之前计算第一交点和第二交点的EOH描述符。
计算EOH描述符可以通过如下步骤实现：
以下采用第一交点举例说明，首先将第一交点的N×N邻域的像素区域等分为16块(N为正整数，例如N取80或100)，计算每一块的边缘方向直方图。每个边缘方向直方图记录着对应块内梯度方向分布信息。每个边缘方向直方图都分为5个子区域(bins)，其中有四个bins各自代表不同的方向(这四个方向的梯度值最大的点的数目大于其余方向的梯度值最大的点的数目)，剩下的一个bin代表没有方向。每个bin的数值就是对应块内该方向的梯度值最大的点的数目，没有方向的那个bin的数值就是对应块内梯度方向不明显的点的数目(即除了上述四个方向其余的方向的点)。这样，每个第一交点就可用一个80维的向量描述，该向量就是第一交点的EOH描述符。
计算出EOH描述符之后，计算第一交点的EOH描述符与第二交点的EOH描述符之间的欧氏距离；若判断出所述欧氏距离小于预设的第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点，若小于第一距离阈值的第二交点有多个，则可以选择任一个，或选择欧氏距离最小的第二交点。
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点。
具体的，若所述欧氏距离大于所述第一距离阈值且小于所述第二距离阈值，则确定出多个所述第二交点；其中，第二距离阈值可以是比第一距离阈值大的阈值；通过预设的第二距离阈值，先确定出多个匹配的第二交点，进一步采用直线匹配最终确定与第一交点匹配的第二交点；采用直线匹配即对多个第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配。
可选地，所述对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配之前，还包括：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第二特征点；所述第二特征点包括所述第一交点。
在上述实施例的基础上，进一步的，在实际应用中，进行直线匹配的方式可以有多种，可选的，作为一种可实施的方式，对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体可以采用如下方式：
计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；
若有超过预设的数量阈值的所述欧氏距离小于预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点。
具体来说，计算特征点的EOH描述符，具体可以通过如下步骤实现：
(1)首先把两条相交的直线分别设为一号直线和二号直线，规则如下：如果一条直线逆时针旋转一个锐角角度后与另一条直线重合，那么该直线就被标为一号直线，另一条直线就是二号直线；
(2)在一号直线上距离第一交点位置每隔一段预设距离d1取一个点，直到取的点超出待配准多光谱图像的范围；
(3)在一号直线的两侧隔一段预设距离d2，然后沿着平行于一号直线的方向继续取点，此时取的点的间隔可以是d1；
(4)在通过第一交点，方向与一号直线垂直的直线上继续取点，包括该直线的两侧，取点规则同上；
(5)对于二号直线，也进行(2)、(3)、(4)的操作；
(6)将取到的所有点作为特征点分别计算EOH描述符。
对于每一个第二交点所在的直线对的第一特征点，都分别计算第一特征点与对应位置的第二特征点的EOH描述符之间的欧氏距离；
若对于某一第二交点所在的直线对，有超过预设的数量阈值的所述欧氏距离小于预设的第三距离阈值，则确定该直线对的所述第二交点为与所述第一交点匹配的第二交点。若有多个第二交点所在的直线对，有超过预设的数量阈值的所述欧氏距离小于预设的第三距离阈值，则可以选择欧氏距离小于预设的第三距离阈值的数量最多的直线对，对应的第二交点作为匹配的第二交点。
上述具体实施方式中，由于是根据多个第二交点所在的相交直线对与所述第一交点所在的相交直线对进行匹配，计算时利用了较多点的描述符，与只根据两个交点进行匹配相比更准确。
在上述实施例的基础上，进一步的，在实际应用中，确定坐标变换参数的方式可以有多种，可选的，作为一种可实施的方式，根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数，具体可以采用如下方式：
根据如下的公式(1)确定坐标变换矩阵H，将坐标变换矩阵H作为坐标变换参数；
其中，为所述第二交点的坐标；为所述第一交点的坐标；坐标变换矩阵
具体来说，最终可以至少确定出三组匹配的第一交点和第二交点，将坐标代入上述公式(1)中计算坐标变换矩阵H的参数，确定出该坐标变换矩阵。
在上述实施例的基础上，进一步的，在实际应用中，作为一种可实施的方式，根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，具体可以采用如下方式：
将所述待配准多光谱图像中的像素点的坐标与所述坐标变换矩阵相乘，使得所述待配准多光谱图像变换到与所述参考多光谱图像相同的坐标系下。
具体来说，将所述待配准多光谱图像中的像素点的坐标与所述坐标变换矩阵相乘，得到变换后的待配准多光谱图像，变换后的待配准多光谱图像与参考多光谱图像在相同的坐标系下。
上述具体实施方式中，根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换矩阵，最终实现了对所述待配准多光谱图像和所述参考多光谱图像进行配准。
图2为本发明多光谱图像配准的装置一实施例的结构示意图。如图2所示，本实施例的装置，可以包括：获取模块201、确定模块202和处理模块203；
其中，获取模块201，用于获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
确定模块202，用于确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
所述确定模块202，还用于根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
处理模块203，用于根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准。
可选地，所述确定模块202，具体用于：
分别对所述第一交点和所述第二交点计算边缘方向直方图EOH描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点。
可选地，所述处理模块203，还用于：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔所述预设距离选取第二特征点；所述第二特征点包括所述第一交点。
可选地，所述确定模块202，具体用于：
计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；所述第一特征点与所述第二特征点的位置对应；
若有超过预设的数量阈值的所述欧氏距离小于所述预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；所述第三距离阈值小于所述第二距离阈值。
可选地，所述确定模块202，具体用于：
根据如下的公式(1)确定坐标变换矩阵H，将所述坐标变换矩阵H作为坐标变换参数；
其中，为所述第二交点的坐标；为所述第一交点的坐标；坐标变换矩阵
可选地，所述处理模块203，具体用于：
将所述待配准多光谱图像中的像素点的坐标与所述坐标变换矩阵相乘，使得所述待配准多光谱图像变换到与所述参考多光谱图像相同的坐标系下。
本实施例的装置，可以用于执行如图1所示方法实施例的技术方案，其实现原理和技术效果类似，此处不再赘述。
本领域普通技术人员可以理解：实现上述方法实施例的全部或部分步骤可以通过程序指令相关的硬件来完成，前述的程序可以存储于一计算机可读取存储介质中，该程序在执行时，执行包括上述方法实施例的步骤；而前述的存储介质包括：ROM、RAM、磁碟或者光盘等各种可以存储程序代码的介质。
最后应说明的是：以上实施例仅用以说明本申请的技术方案，而非对其限制；尽管参照前述实施例对本申请进行了详细的说明，本领域的普通技术人员应当理解：其依然可以对前述各实施例所记载的技术方案进行修改，或者对其中部分技术特征进行等同替换；而这些修改或者替换，并不使相应技术方案的本质脱离本申请各实施例技术方案的精神和范围。
Claims (4)
1.一种多光谱图像配准的方法，其特征在于，包括：
获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准；
其中，所述确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体包括：
分别对所述第一交点和所述第二交点计算边缘方向直方图EOH描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
其中，所述对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配之前，还包括：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔所述预设距离选取第二特征点；所述第二特征点包括所述第一交点；
其中，所述对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点，具体包括：
计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；所述第一特征点与所述第二特征点的位置对应；
若有超过预设的数量阈值的所述欧氏距离小于预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；所述第三距离阈值小于所述第二距离阈值；
若有多个所述第二交点所在的直线对，有超过预设的数量阈值的所述欧氏距离小于所述预设的第三距离阈值，则选择所述欧氏距离小于所述预设的第三距离阈值的数量最多的所述第二交点所在的直线对，确定所述第二交点为与所述第一交点匹配的第二交点。
2.根据权利要求1所述的方法，其特征在于，所述根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数，包括：
根据如下的公式(1)确定坐标变换矩阵H，将所述坐标变换矩阵H作为坐标变换参数；
<mrow>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msup>
<mi>x</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<msup>
<mi>y</mi>
<mo>&prime;</mo>
</msup>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>=</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mn>11</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>12</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>13</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<msub>
<mi>a</mi>
<mn>21</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>22</mn>
</msub>
</mtd>
<mtd>
<msub>
<mi>a</mi>
<mn>23</mn>
</msub>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>0</mn>
</mtd>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>&times;</mo>
<mfenced open = "[" close = "]">
<mtable>
<mtr>
<mtd>
<mi>x</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mi>y</mi>
</mtd>
</mtr>
<mtr>
<mtd>
<mn>1</mn>
</mtd>
</mtr>
</mtable>
</mfenced>
<mo>-</mo>
<mo>-</mo>
<mo>-</mo>
<mrow>
<mo>(</mo>
<mn>1</mn>
<mo>)</mo>
</mrow>
</mrow>
其中，为所述第二交点的坐标；为所述第一交点的坐标；坐标变换矩阵
3.根据权利要求2所述的方法，其特征在于，所述根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准，包括：
将所述待配准多光谱图像中的像素点的坐标与所述坐标变换矩阵相乘，使得所述待配准多光谱图像变换到与所述参考多光谱图像相同的坐标系下。
4.一种多光谱图像配准的装置，其特征在于，包括：
获取模块，用于获取待配准多光谱图像中多组相交的直线对的第一交点，以及获取参考多光谱图像中多组相交的直线对的第二交点；
确定模块，用于确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
所述确定模块，还用于根据所述第一交点以及与所述第一交点匹配的第二交点，确定坐标变换参数；
处理模块，用于根据所述坐标变换参数对所述待配准多光谱图像和所述参考多光谱图像进行配准；
其中，所述确定模块，具体用于：
分别对所述第一交点和所述第二交点计算边缘方向直方图EOH描述符；
计算所述第一交点的EOH描述符与所述第二交点的EOH描述符之间的欧氏距离；
判断所述欧氏距离是否小于预设的第一距离阈值，若小于所述第一距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；
若大于所述第一距离阈值，则判断所述欧氏距离是否小于预设的第二距离阈值，若小于所述第二距离阈值，则对所述第二交点所在的相交的直线对与所述第一交点所在的相交的直线对进行直线匹配，确定所述参考多光谱图像上与所述第一交点匹配的第二交点；
其中，所述处理模块，还用于：
在所述第二交点所在的相交的直线对上，对任一条直线，依次在所述第二交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第二交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔预设距离选取第一特征点；所述第一特征点包括所述第二交点；
在所述第一交点所在的相交的直线对上，对任一条直线，依次在所述第一交点所在的直线、与所述直线平行并在所述直线两侧的直线、与所述直线垂直并过所述第一交点的垂线、与所述垂线平行并在所述垂线两侧的直线上间隔所述预设距离选取第二特征点；所述第二特征点包括所述第一交点；
其中，所述确定模块，还用于：计算所述第一特征点的EOH描述符与所述第二特征点的EOH描述符之间的欧氏距离；所述第一特征点与所述第二特征点的位置对应；
若有超过预设的数量阈值的所述欧氏距离小于预设的第三距离阈值，则确定所述第二交点为与所述第一交点匹配的第二交点；所述第三距离阈值小于所述第二距离阈值；
若有多个所述第二交点所在的直线对，有超过预设的数量阈值的所述欧氏距离小于所述预设的第三距离阈值，则选择所述欧氏距离小于所述预设的第三距离阈值的数量最多的所述第二交点所在的直线对，确定所述第二交点为与所述第一交点匹配的第二交点。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510345294.6A CN104992431B (zh) | 2015-06-19 | 2015-06-19 | 多光谱图像配准的方法及装置 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201510345294.6A CN104992431B (zh) | 2015-06-19 | 2015-06-19 | 多光谱图像配准的方法及装置 |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| CN104992431A CN104992431A (zh) | 2015-10-21 |
| CN104992431B true CN104992431B (zh) | 2018-02-27 |
Family
ID=54304239
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201510345294.6A Active CN104992431B (zh) | 2015-06-19 | 2015-06-19 | 多光谱图像配准的方法及装置 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN104992431B (zh) |
Families Citing this family (5)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN106157306B (zh) * | 2016-06-30 | 2019-05-07 | 北京理工大学 | 一种基于多个光谱图像的配准方法以及装置 |
| CN108305277B (zh) * | 2017-12-26 | 2020-12-04 | 中国航天电子技术研究院 | 一种基于直线段的异源图像匹配方法 |
| CN109389628B (zh) * | 2018-09-07 | 2021-03-23 | 北京邮电大学 | 图像配准方法、设备和存储介质 |
| CN111681271B (zh) * | 2020-08-11 | 2020-10-30 | 湖南大学 | 一种多通道多光谱相机配准方法、系统及介质 |
| CN113096168B (zh) * | 2021-03-17 | 2024-04-02 | 西安交通大学 | 一种结合sift点和控制线对的光学遥感图像配准方法及系统 |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101315698A (zh) * | 2008-06-25 | 2008-12-03 | 中国人民解放军国防科学技术大学 | 基于直线特征图像配准中的特征匹配方法 |
| CN103679741A (zh) * | 2013-12-30 | 2014-03-26 | 北京建筑大学 | 一种基于三维线特征的激光点云数据自动配准方法 |
-
2015
- 2015-06-19 CN CN201510345294.6A patent/CN104992431B/zh active Active
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN101315698A (zh) * | 2008-06-25 | 2008-12-03 | 中国人民解放军国防科学技术大学 | 基于直线特征图像配准中的特征匹配方法 |
| CN103679741A (zh) * | 2013-12-30 | 2014-03-26 | 北京建筑大学 | 一种基于三维线特征的激光点云数据自动配准方法 |
Non-Patent Citations (3)
| Title |
|---|
| A Novel SIFT-Like-Based Approach for FIR-VS Images Registration;Cristhian Aguilera 等;《11th International Conference on Quantitative InfraRed Thermography》;20120614;摘要，第3节，图3 * |
| 基于干线对的红外与可见光最优图像配准算法;高峰 等;《计算机学报》;20070615;第30卷(第6期);第2.1、3.2.2、4节 * |
| 多光谱图像配准算法研究;臧丽;《中国优秀硕士学位论文全文数据库 信息科技辑》;20080115;I138-329 * |
Also Published As
| Publication number | Publication date |
|---|---|
| CN104992431A (zh) | 2015-10-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| CN104992431B (zh) | 多光谱图像配准的方法及装置 | |
| CN109829398B (zh) | 一种基于三维卷积网络的视频中的目标检测方法 | |
| CN103456022B (zh) | 一种高分辨率遥感图像特征匹配方法 | |
| JP5385105B2 (ja) | 画像検索方法およびシステム | |
| CN103700099B (zh) | 一种旋转和尺度不变的宽基线立体匹配方法 | |
| CN105427298B (zh) | 基于各向异性梯度尺度空间的遥感图像配准方法 | |
| CN105184801B (zh) | 一种基于多层次策略的光学与sar影像高精度配准方法 | |
| CN104200495B (zh) | 一种视频监控中的多目标跟踪方法 | |
| CN104200461B (zh) | 基于互信息图像选块和sift特征的遥感图像配准方法 | |
| CN107248159A (zh) | 一种基于双目视觉的金属工件缺陷检测方法 | |
| CN105184830B (zh) | 一种对称图像对称轴检测定位方法 | |
| CN102722887A (zh) | 一种图像配准方法及装置 | |
| CN105631872B (zh) | 基于多特征点的遥感图像配准方法 | |
| CN102722731A (zh) | 一种基于改进sift算法的高效图像匹配方法 | |
| CN103679702A (zh) | 一种基于图像边缘矢量的匹配方法 | |
| CN103489191B (zh) | 一种遥感图像显著目标变化检测方法 | |
| CN104978582B (zh) | 基于轮廓弦角特征的遮挡目标识别方法 | |
| Benseddik et al. | SIFT and SURF Performance evaluation for mobile robot-monocular visual odometry | |
| CN104992433B (zh) | 基于线段匹配的多光谱图像配准的方法及装置 | |
| CN101650784A (zh) | 一种利用结构上下文特征进行图像匹配的方法 | |
| CN110675442A (zh) | 一种结合目标识别技术的局部立体匹配方法及系统 | |
| CN105678720A (zh) | 一种全景拼接判断图像匹配方法及装置 | |
| CN103336964B (zh) | 一种基于模值差镜像不变性的sift图像匹配方法 | |
| CN107578029B (zh) | 计算机辅助图片认证的方法、装置、电子设备及存储介质 | |
| Wang et al. | LBP-based edge detection method for depth images with low resolutions |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| GR01 | Patent grant | ||
| GR01 | Patent grant |