基于特征密度聚类与正态分布变换的侧扫声呐配准方法技术

本发明专利技术涉及图像配准技术领域，特别涉及基于特征密度聚类与正态分布变换的侧扫声呐配准方法，主要解决目前的方法难以快速有效地配准侧扫声呐图像的问题。主要包括以下步骤：首先使用均值滤波器对待配准的侧扫声呐图像进行滤波，在此基础上计算图像的梯度，阈值化梯度图像得到全部特征点后，使用基于密度的空间聚类方法DBSCAN对梯度图像中的像素点进行聚类得到图像特征点，最后使用正态分布变换NDT算法优化特征点，得到两幅图像的配准关系。本发明专利技术克服了基于特征点方法中由于图像特征点少而配准失败的困难，也不存在两幅图像旋转角度大而配准失败的问题,可广泛应用于海底地形地貌探测、AUV自主导航和图像处理等领域。

【技术实现步骤摘要】

本专利技术涉及图像配准
，特别涉及基于特征密度聚类与正态分布变换的侧扫声呐配准方法。
技术介绍
图像配准是将传感器在不同位置获取的图像转换到同一坐标系下的技术。通过图像配准获取的信息可用于水下航行器感知周围环境，能有效提高水下航行器的自主感知能力。水下环境由于其能见度低、无光照、会对水下照明系统发生散射等特点，传统的光学传感器用于水下图像获取存在许多挑战。由于侧扫声呐不受这些光学因素影响，使用它进行水下图像获取越来越受到重视，侧扫声呐配准技术也相应发展起来。目前，用于声呐图像配准的方法主要有两种：基于特征的方法和基于区域的方法。基于特征的方法用于图像中特征点较多且图像是连续帧的情况，例如Harris角点检测算法(NegahdaripourS，Onprocessingandregistrationofforward-scanacousticvideoimagery，ComputerandRobotVision，2005)和SIFT(ScaleInvariantFeatureTransform)结合RANSAC(RANdomSAmpleConsensus)的算法(NegahdaripourS，Dynamicsceneanalysisandmosaicingofbenthichabitatsbyfssonarimaging—issuesandcomplexities，InProceedingsofOCEANS，2011)。但是这些方法在特征点数较少或者匹配成功的特征数较少时配准结果较差。为了解决基于特征的方法在配准声呐图像时存在的问题，基于区域的配准方法被用于声呐图像配准。基于区域的配准方法把整幅图像作为处理对象，如相位相关法。但是这种方法只适用于旋转角小的情况。AUV(AutonomousUnderwaterVehicle)搭载侧扫声呐提取周围环境信息，通过配准声呐图像获取水下全局地形地貌，提高AUV对环境的感知能力，而如何快速有效地配准侧扫声呐图像是目前亟待解决的主题。
技术实现思路
本专利技术针对现有技术中的不足，提出一种基于特征密度聚类与正态分布变换的侧扫声呐配准方法。首先使用均值滤波器对待配准的侧扫声呐图像进行滤波，在此基础上计算图像的梯度，并使用基于密度的空间聚类方法DBSCAN(Density-BasedSpatialClusteringofApplicationswithNoise)对梯度图像中的像素点进行聚类得到图像特征点，最后使用正态分布变换NDT算法(NormalDistributedTransform)优化特征点，得到两幅图像的配准关系。该技术方法可应用于海底地形地貌探测、AUV自主导航和图像处理等领域。基于特征密度聚类与正态分布变换的侧扫声呐配准方法具体包括下列步骤：步骤一，中值滤波器处理原始侧扫声呐图像。原始侧扫声呐图像中有噪声，首先要使用中值滤波器对声呐图像去噪处理。步骤二，计算梯度幅值图像。使用Sobel算子(Sobeloperator，索贝尔算子)计算得到梯度幅值图像。步骤三，阈值化梯度图像得到全部特征点。步骤四，使用DBSCAN聚类得到图像特征点。DBSCAN是一种基于密度的空间聚类算法，该算法将具有足够密度的区域划分为簇，并在具有噪声的空间数据库中发现任意形状的簇，它将簇定义为密度相连的点的最大集合。DBSCAN算法的目的在于过滤低密度区域，发现稠密度样本点。步骤五，使用NDT算法(NormalDistributedTransform，正态分布变换)得到两幅图像的配准结果。NDT算法将二维特征点转换为分段连续可微的概率分布。对一幅有聚类特征点的图像进行NDT处理时，首先将图像分成相同大小的栅格图像，对包含特征点的栅格单元，以概率密度形式对其进行分段连续可微，并用正态分布表示每个栅格单元中特征点的概率分布：p(x)∝exp(-(x-q)TΣ-1(x-q)2)]]>其中，x是一个栅格单元中i个点的二维坐标，q是栅格单元中包含的特征点位置的平均值，即q=1nΣi=1nxi;]]>Σ是这些点的协方差矩阵，Σ=1nΣi=1n(xi-q)(xi-q)T.]]>两幅图像的空间变换关系T可以表示为：T:x′y′=cosθ-sinθsinθcosθxy+txty]]>其中，txty]]>是两幅图像的平移关系，θ是两幅图像的旋转关系。其中，步骤四具体聚类过程如下：a.扫描整个梯度图像，找到任意一个核心点，对该核心点进行扩充。扩充的方法是寻找从该核心点出发的所有密度相连的数据点。b.遍历该核心点邻域内的所有核心点，寻找与这些像素点密度相连的点，直到没有可以扩充的像素点为止。最后聚类成的簇的边界节点都是非核心数据点。c.之后就是重新扫描梯度图像(不包括之前寻找到的簇中的任何像素点)，寻找没有被聚类的核心点，再重复上面的步骤，对该核心点进行扩充直到数据集中没有新的核心点为止。数据集中没有包含在任何簇中的数据点就构成异常点。最后的数据集就是聚类后得到的图像特征点。步骤五中，配准的目的就是为了找到两幅声呐图像之间的这些变换参数。具体包括如下步骤：a.计算一幅声呐图像的NDT；b.初始化平移参数和旋转参数(使用0或者里程计数据进行初始化)；c.对于另一幅声呐图像，根据初始化参数将其变换到第一幅图像的坐标系中；d.计算坐标变换后声呐图像中每个点的正态分布；e.通过评价每个点的正态分布之和计算平移参数和旋转参数的分数；f.使用Hessian矩阵法对这些分数进行优化，计算新的参数估计值；g.返回步骤c，直到满足收敛要求。参数txty]]>和θ的分数可以表示为：score(p)=Σiexp(-(xi′-qi)TΣi-1(xi′-qi)2)]]>其中，p＝(tx,ty,θ)T，xi'是根据坐标变换参数将第二幅图像中的点xi变换到第一幅图像坐标系中坐标。qi和分别是x′i的均值和协方差。上述第e步中对分数进行优化，就是使分数score(p)最大，即-score(p)最小。牛顿迭代算法的详细操作如下：为了使函数f最小，在每次的迭代中要解方程：HΔp＝-g计算出的Δp与当前参数相加，得到新的参数p←p+Δpp用0值初始化，所以此时的p＝Δp。g是f的转置梯度：gi=∂f∂pi]]>H是f的Hessian矩阵：本文档来自技高网...

【技术保护点】
一种基于特征密度聚类与正态分布变换的侧扫声呐配准方法，包括下列步骤：步骤一，中值滤波器处理原始侧扫声呐图像：原始侧扫声呐图像中有噪声，首先要使用中值滤波器对声呐图像去噪处理；步骤二，计算梯度幅值图像：使用Sobel算子计算得到梯度幅值图像；步骤三，阈值化梯度图像得到全部特征点；步骤四，使用DBSCAN聚类得到图像特征点：将具有足够密度的区域划分为簇，并在具有噪声的空间数据库中发现任意形状的簇，它将簇定义为密度相连的点的最大集合；DBSCAN算法的目的在于过滤低密度区域，发现稠密度样本点；步骤五，使用NDT算法得到两幅图像的配准结果：NDT算法将二维特征点转换为分段连续可微的概率分布；对一幅有聚类特征点的图像进行NDT处理时，首先将图像分成相同大小的栅格图像，对包含特征点的栅格单元，以概率密度形式对其进行分段连续可微，并用正态分布表示每个栅格单元中特征点的概率分布：p(x)∝exp(-(x-q)TΣ-1(x-q)2)]]>其中，x是一个栅格单元中i个点的二维坐标，q是栅格单元中包含的特征点位置的平均值，即q=1nΣi=1nxi;]]>Σ是这些点的协方差矩阵，Σ=1nΣi=1n(xi-q)(xi-q)T;]]>两幅图像的空间变换关系T可以表示为：T:x′y′=cosθ-sinθsinθcosθxy+txty]]>其中，txty]]>是两幅图像的平移关系，θ是两幅图像的旋转关系。...

【技术特征摘要】
1.一种基于特征密度聚类与正态分布变换的侧扫声呐配准方法，包括下列步骤：
步骤一，中值滤波器处理原始侧扫声呐图像：原始侧扫声呐图像中有噪声，首先要使用
中值滤波器对声呐图像去噪处理；
步骤二，计算梯度幅值图像：使用Sobel算子计算得到梯度幅值图像；
步骤三，阈值化梯度图像得到全部特征点；
步骤四，使用DBSCAN聚类得到图像特征点：将具有足够密度的区域划分为簇，并在具
有噪声的空间数据库中发现任意形状的簇，它将簇定义为密度相连的点的最大集合；DBSCAN
算法的目的在于过滤低密度区域，发现稠密度样本点；
步骤五，使用NDT算法得到两幅图像的配准结果：
NDT算法将二维特征点转换为分段连续可微的概率分布；对一幅有聚类特征点的图像进
行NDT处理时，首先将图像分成相同大小的栅格图像，对包含特征点的栅格单元，以概率密
度形式对其进行分段连续可微，并用正态分布表示每个栅格单元中特征点的概率分布：
p(x)∝exp(-(x-q)TΣ-1(x-q)2)]]>其中，x是一个栅格单元中i个点的二维坐标，q是栅格单元中包含的特征点位置的平均
值，即q=1nΣi=1nxi;]]>Σ是这些点的协方差矩阵，Σ=1nΣi=1n(xi-q)(xi-q)T;]]>两幅图像的空间变换关系T可以表示为：
T:x′y′=cosθ-sinθsinθcosθxy+txty]]>其中，txty]]>是两幅图像的平移关系，θ是两幅图像的旋转关系。
2.根据权利要求1所述的基于特征密度聚类与正态分布变换的侧扫声呐配准方法，其特
征在于，步骤四中算法的具体聚类过程如下：
a.扫描整个梯度图像，找到任意一个核心点，对该核心点进行扩充；扩充的方法是寻找
从该核心点出发的所有密度相连的数据点；
b.遍历该核心点邻域内的所有核心点，寻找与这些像素点密度相连的点，直到没有可以
扩充的像素点为止；最后聚类成的簇的边界节点都是非核心数据点；
c.重新扫描梯度图像，扫描内容不包括之前寻找到的簇中的任何像素点，寻找没有被聚
类的核心点，再重复上面的步骤，对该核心点进行扩充直到数据集中没有新的核心点为止；

\t数据集中没有包含在任何簇中的数据点就构成异常点；最后的数据集就是聚类后得到的图像
特征点。
3.根据权利要求1所述的基于特征密度聚类与正态分布变换的侧扫声呐配准...

【专利技术属性】
技术研发人员：何波，宋艳，张论娟，年睿，沈钺，沙启鑫，高强，冯晨，
申请(专利权)人：中国海洋大学，
类型：发明
国别省市：山东;37