一种面向海洋数据可视化的图像序列匹配方法
【专利摘要】本发明公开一种面向海洋数据可视化的图像序列匹配方法，其包括步骤：融合点和区域的关键区域检测与描述；以及关键区域匹配。该关键区域匹配的步骤进一步包括基于欧式距离的匹配；以及排除误检匹配。该排除误检匹配的步骤包括基于平均运动距离的粗匹配和基于随机抽样一致性的精确匹配。本发明能适应于多样化的图像数据，具有明显视觉特征，精确度高，且匹配效率高。
【专利说明】一种面向海洋数据可视化的图像序列匹配方法
【技术领域】
[0001]本发明涉及图像处理领域，特别涉及一种面向海洋数据可视化的图像序列匹配方法。
【背景技术】
[0002]图像序列数据在海洋数据中占有较大比例，直接通过图像序列数据展示海洋信息不够直观。因此，研究者通常使用基于全景图技术或基于三维重建技术对海洋信息进行展示。在由图像序列数据获得全景图和三维信息的过程中，图像序列的匹配问题是其关键。
[0003]早期的技术人员直接使用图像的颜色信息来匹配图像，但是由于光照的变化，视角的多样性等影响，使得其匹配准确率不高。近年来技术人员提出了基于局部特征的图像匹配方法，该方法的核心技术主要包括两部分:关键区域的检测和关键区域的描述。在检测部分，通常仅使用点的检测或区域的检测，不能满足多样性的海洋数据的需求；在基于点的检测方法中，现在应用最广泛的是尺度不变特征变换(SIFT, Scale Invariant FeatureTransform)检测子，该类方法使用的特征点不具有明显的视觉特性，在海洋数据可视化的应用中受到一定的限制；在误匹配的排除过程中，通常直接使用随机抽样一致性(RANSAC，RANdom Sample Consensus)方法，时间复杂度较高。
[0004]现有技术有如下缺点:
[0005](I)关键区域仅使用点检测或区域检测得到，较单一，影响匹配准确率且不能适应多样化的图像数据。
[0006](2)所使用的特征点不具有视觉特性。
[0007](3)精确匹配过程耗时较长。
【发明内容】
[0008]本发明主要解决的技术问题是提供一种面向海洋数据可视化的图像序列匹配方法，解决现有基于特征的图像匹配方法中所使用关键区域较单一，不能适应于多样化的图像数据，不具有明显视觉特征，精确度不高，耗时较长的问题。
[0009]为解决上述技术问题，本发明采用的技术方案是:提供一种面向海洋数据可视化的图像序列匹配方法，其包括步骤:融合点和区域的关键区域检测与描述；以及关键区域匹配。该关键区域匹配的步骤进一步包括基于欧式距离的匹配；以及排除误检匹配。该排除误检匹配的步骤包括基于平均运动距离的粗匹配和基于随机抽样一致性的精确匹配。
[0010]在本发明一个较佳实施例中，所述融合点和区域的关键区域检测包括如下步骤:输入图像；对输入图像进行点检测，所使用的点检测方法为Harris角点，得到角点集合；对角点集合进行筛选，取特征最明显的前1/4角点构成的角点集作为构成关键区域的部分之一；对输入图像进行区域检测，所使用的区域检测方法为最稳定极值区域检测，得到区域集合；对区域集合进行筛选，将周围区域纹理不明显的区域排除，得到区域集作为构成关键区域的部分之一；以及由角点集和区域集构成最终的关键区域集合。[0011]在本发明一个较佳实施例中，所述融合点和区域的关键区域描述包括Harris角点方向描述和最稳定极值区域描述。
[0012]在本发明一个较佳实施例中,所述Harris角点方向描述方法包括如下步骤:选择Harris角点周围3*3邻域的点；计算邻域内点的梯度大小和方向；通过直方图统计法对邻截至内的点进行投影；以及累加值最大的三个方向描述角点方向。
[0013]在本发明一个较佳实施例中，选择Harris角点周围3*3邻域的9个点；分别对该9 个点利用公式 m(x, y) = {(L(x+1), y)-L(x_l, y))2+(L(x, y+l)_L(x, y_l))2}1/2 和 θ (x, y)=tan-1 (((L(x, y+l)-L(x, y-1))/(L(x+1, y)-L(χ-l, y)))计算梯度大小和方向，其中，m(x,y)为梯度大小，θ (χ, y)为梯度方向，(χ, y)为角点对应的坐标，L(x, y)为对应于点(x, y)的灰度值；通过直方图统计法，将角点3*3邻域内的9个点投影到0-360度平均分成的具有36个柱的直方图中。
[0014]在本发明一个较佳实施例中，所述最稳定极值区域通过椭圆拟合方法对其进行拟合，然后进行椭圆归一化将其归一化为圆；且以梯度累加值最大的第一个方向作为最稳定极值区域的方向。
[0015]在本发明一个较佳实施例中，所述关键区域匹配过程包括:基于欧式距离的匹配；以及排除误检匹配；其中，排除误检匹配步骤进一步包括基于平均运动距离的粗匹配和基于随机抽样一致性的精确匹配。
[0016]在本发明一个较佳实施例中，所述基于欧式距离的匹配步骤中，通过欧式距离对关键区域进行匹配得到最原始的匹配结果，欧式距离的计算方法如下式所示:
[0017]
【权利要求】
1.一种面向海洋数据可视化的图像序列匹配方法，其包括步骤: 融合点和区域的关键区域检测与描述；以及 关键区域匹配，该步骤进一步包括: 基于欧式距离的匹配；以及 排除误检匹配，该步骤包括基于平均运动距离的粗匹配和基于随机抽样一致性的精确匹配。
2.如权利要求1所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述融合点和区域的关键区域检测包括如下步骤: 输入图像； 对输入图像进行点检测，所使用的点检测方法为Harris角点，得到角点集合； 对角点集合进行筛选，取特征最明显的前1/4角点构成的角点集作为构成关键区域的部分之一； 对输入图像进行区域检测，所使用的区域检测方法为最稳定极值区域检测，得到区域集合; 对区域集合进行筛选，将周围区域纹理不明显的区域排除，得到区域集作为构成关键区域的部分之一；以及 由角点集和区域集构 成最终的关键区域集合。
3.如权利要求1所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述融合点和区域的关键区域描述包括Harris角点方向描述和最稳定极值区域描述。
4.如权利要求3所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述Harris角点方向描述方法包括如下步骤: 选择Harris角点周围3*3邻域的点； 计算邻域内点的梯度大小和方向； 通过直方图统计法对邻截至内的点进行投影；以及 累加值最大的三个方向描述角点方向。
5.如权利要求4所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述步骤中，选择Harris角点周围3*3邻域的9个点；分别对该9个点利用公式 m(x，y) = {(L(x+1), y)-L(x_l, y))2+(L(x, y+l)_L(x, y_l))2}1/2 和 θ (χ，y)=tarT1 (((L(x, y+l)-L(x, y-1))/(L(x+1, y)-L(χ-1, y)))计算梯度大小和方向，其中，m(x, y)为梯度大小，Θ (17)为梯度方向，(1，7)为角点对应的坐标，1(1，50为对应于点(x，y)的灰度值；通过直方图统计法，将角点3*3邻域内的9个点投影到0-360度平均分成的具有36个柱的直方图中。
6.如权利要求4所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述最稳定极值区域通过椭圆拟合方法对其进行拟合，然后进行椭圆归一化将其归一化为圆；且以梯度累加值最大的第一个方向作为最稳定极值区域的方向。
7.如权利要求5所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述最稳定极值区域通过椭圆拟合方法对其进行拟合，然后进行椭圆归一化将其归一化为圆；且以梯度累加值最大的第一个方向作为最稳定极值区域的方向。
8.如权利要求1所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述关键区域匹配过程包括如下步骤: 基于欧式距离的匹配；以及 排除误检匹配；其中，排除误检匹配步骤进一步包括基于平均运动距离的粗匹配和基于随机抽样一致性的精确匹配。
9.如权利要求8所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述基于欧式距离的匹配步骤中，通过欧式距离对关键区域进行匹配得到最原始的匹配结果，欧式距离的计算方法如下式所示:
10.如权利要求8所述的面向海洋数据可视化的图像序列匹配方法，其特征在于，所述基于平均运动距离的粗匹配包括如下子步骤: 使用欧氏距离公式计算图像中关键区域的运动距离；
【文档编号】G06T7/00GK103544706SQ201310517989
【公开日】2014年1月29日 申请日期:2013年10月28日 优先权日:2013年10月28日
【发明者】黄磊, 魏志强, 纪筱鹏, 殷波, 丛艳平, 盛艳秀, 洪博为 申请人:中国海洋大学