基于MeanShift的高分辨率遥感影像分割距离度量优化方法
【专利摘要】基于Meanshift的高分辨率遥感图像分割中合并区域过程的距离度量优化方法，充分考虑遥感图像数据的特点，利用光谱匹配距离度量代替传统的欧式距离，具体有光谱角度匹配度量，光谱相似性度量和核映射光谱匹配度量，使分割结果更加精确。
【专利说明】基于MeanShift的高分辨率遥感影像分割距离度量优化方法
【技术领域】
[0001]本发明是在高分辨率遥感影像处理领域中，针对基于MeanShift算法的遥感图像分割技术，对其中区域合并过程所涉及到的距离度量计算进行优化，从而进一步得到精度更高的图像分割效果。
【背景技术】
[0002]随着遥感技术的发展，高分辨率遥感影像逐渐在许多领域中得到广泛的应用。由于高分辨率遥感影像拥有更丰富的光谱波段、更详细的地形结构等数据信息，所以遥感图像处理技术的要求也随之不断提高，影像分割技术是面向对象的高分辨率遥感图像信息提取和分析的重要步骤之一，其分割质量直接决定了后续影像处理的精度，因此遥感图像分割具有重大的研究意义。
[0003]MeanShift均值漂移算法是一种无参的统计迭代方法，利用核密度估计求取一定范围内的梯度最大点，由于其简单的形式、较快的收敛速度以及较高的稳定性、抗噪性等，在计算机模式识别、聚类分析、图像跟踪和图像分割等方面有广泛的应用。
[0004]在基于MeanShift的遥感影像分割过程中，区域合并是其重要的步骤之一，它是在图像进行滤波聚类得到众多细分的同质区域小块之后，再进一步根据一定准则得到图像的最终分割结果，因此合并准则的精确对分割结果有重要的影响，而距离度量作为合并准则的重要环节，其具体的策略需要根据遥感图像数据的特点来进行设置。本发明据此提出了一种针对区域合并时的距离度量优化策略，使高分辨率遥感影像的分割在理论上得到更高精度的效果。
【发明内容】
[0005]本发明要克服现有技术分割图像时精度低的缺点，提供基于MeanShift的高分辨率遥感图像分割距离度量优化方法。
[0006]本发明为了在基于MeanShift的遥感图像分割中得到精度更高的分割效果，在区域合并过程中对距离度量的计算采用光谱匹配技术代替传统的欧式距离的计算，充分考虑遥感图像的数据特征，使区域之间的相似度量的计算更加精准，从而得到更精确的分割图像。
[0007]基于MeanShift的高分辨率遥感图像分割距离度量优化方法，包括以下步骤:
[0008]I)、高分辨率遥感图像输入，转换成栅格数据进行处理；
[0009]2)、利用MeanShift算法对遥感图像进行滤波处理，得到以模点为中心的大量同质区域；
[0010]3)、对滤波之后的大量同质区域进行区域合并，计算区域之间的相似性，将传统的欧式距离度量计算方法由光谱匹配度量计算方法或核光谱映射度量计算方法代替；
[0011]4)、设置恰当阈值来评判两区域的相似性度量，初步形成分割结果，在进一步的尺度区域合并过程中，通过区域之间相似度的比较，融入新的区域来达到规定的尺度。最终得到精度较高的分割结果。
[0012]进一步，步骤2)中的MeanShift算法是无参的核密度估计，完全依靠训练数据样本，其具体的原理主要依赖以下公式:
[0013]设样本集{xj，i = l,..., 11，为Rd维空间中的任意点，求其中χ点的密度函数估计，以κ(χ)为核函数，h为窗宽，利用单位矩阵，得:
[。。14]
【权利要求】
1.基于MeanShift的高分辨率遥感图像分割距离度量优化方法，包括以下步骤: 1)、高分辨率遥感图像输入，转换成栅格数据进行处理； 2)、利用MeanShift算法对遥感图像进行滤波处理，得到以模点为中心的大量同质区域； 3)、对滤波之后的大量同质区域进行区域合并，计算区域之间的相似性，将传统的欧式距离度量计算方法由光谱匹配度量计算方法或核光谱映射度量计算方法代替； 4)、设置阈值来评判两区域的相似性度量，初步形成分割结果，在进一步的尺度区域合并过程中，通过区域之间相似度的比较，融入新的区域来达到规定的尺度；最终得到精度较高的分割结果。
2.如权利要求1所述的基于MeanShift的高分辨率遥感图像分割距离度量优化方法，其特征在于:步骤2)中的MeanShift算法是无参的核密度估计，完全依靠训练数据样本，其具体的算法如下: 设样本集{xj，i = I,...,11，为Rd维空间中的任意点，求其中X点的密度函数估计，以K(X)为核函数，h为窗宽，利用单位矩阵，得:
求概率密度最大点，对上式求导并引用核函数G(X)，以及标准化系数ck，d，整理化简可得:
将等式右边最后面的部分作为MeanShift向量，具体如下:
从核概率密度式子推导得出MeanShift向量mh,e(X)，可知其总是指向概率密度增加最大的方向，利用此性质，可以通过对图像像素点的迭代处理，经过滤波过程，得到以模点值为中心的大量同质区域。
3.如权利要求1所述的基于MeanShift的高分辨率遥感图像分割距离度量优化方法，其特征在于:步骤3)中所涉及到的区域相似度的具体计算方式，传统的欧式距离的计算具体为以下式子:
其中，X，y为遥感图像中任意像素点，i = 1，...，n，其中n为波段值，因此两变量可看作是η维向量，此距离度量计算简单，几乎适用于任意情况下的距离度量计算，但在η值较大的高分辨率影像中，如此计算的距离度量没有充分考虑光谱值信息，因此合并效果并不能达到令人满意的效果。 光谱匹配是高光谱影像中像元光谱在每一个波段的变化量与方向的分析，能够突出特征谱段，提取光谱维信息，充分考虑到遥感图像数据的特点，从而能够较精确的定量表示像元间的相似性或者差异性。利用常见的几种光谱匹配技术代替欧式距离计算方法，具体公式如下: 光谱角度匹配:
x，y为η维光谱矢量，Θ为两矢量余弦夹角，计算两者的余弦角度与两矢量的长度无关，因此不受增益因素影响。 光谱相似度:
上式中，μ、O分别为均值和标准差，rxy为皮尔逊系数。光谱相似度综合考虑了光谱间的形状和距离关系，通常值越小，两矢量间的光谱相似性越大。 核映射光谱匹配度量:
d, H(x, y) = arccos exp (_ λ d2 (x, y)) 上式中λ为大于O的自定义参数，d(x，y)为光谱相似度量，此度量方式是光谱相似度结合核函数的作用，在样本计算过程中增加权重因素来计算区域相似度。 具体执行中将遥感图像不同方式下的距离度量值求出，根据具体数值可分析不同距离度量的精度效果。
【文档编号】G06T7/00GK104200458SQ201410367705
【公开日】2014年12月10日 申请日期:2014年7月30日 优先权日:2014年7月30日
【发明者】王卫红, 徐文涛, 夏列钢 申请人:浙江工业大学