一种高空间分辨率遥感图像的混合分类方法技术

本发明专利技术公开了一种高空间分辨率遥感图像的混合分类方法，涉及遥感图像处理技术领域，包括如下步骤：利用均值漂移滤波方法对遥感图像进行低通滤波；利用Canny方法获得每一个卫星波段的边缘梯度结果；借助矢量场模型得到了多光谱边缘强度；然后利用分水岭变换对遥感图像进行分割；利用支持向量机方法对遥感影像中的土地覆盖进行逐像素分类；利用投票原理完成基于对象的遥感分类。本发明专利技术将图像分割和逐像素分类相结合，完成了高分辨率遥感图像的分类，对高空间分辨率遥感数据分类具有良好的处理效果，有效提高分割区域边界精度。有效提高分割区域边界精度。有效提高分割区域边界精度。

【技术实现步骤摘要】
一种高空间分辨率遥感图像的混合分类方法


[0001]本专利技术属于遥感图像处理
，具体涉及一种高空间分辨率遥感图像的混合分类方法。

技术介绍

[0002]随着遥感技术的发展，遥感数据的分辨率越来越高，如SPOT、IKONOS、Quickbird、Worldview和Geoeyes等。有些空间分辨率甚至达到0.5米以下；因此，它提供了更多关于地球上土地覆盖的信息，并有可能在较小的尺度上监测土地覆盖。此外，高空间分辨率使得我们可以获得更多的地物结构细节。然而，传统的逐像素分类方法不能满足土地覆盖分类和地物识别的要求。因此，学者们提出了许多新的方法，利用地物的形状、边界和纹理来辅助分类。
[0003]在高分辨率遥感图像分析中，最突出的方法是基于对象的图像分析(OBIA)方法。首先，基于光谱、纹理等信息遥将感图像分割成有意义的区域。其次，采用图像分类方法，利用分割结果完成分类。在OBIA方法的应用过程中，关键的一步是图像分割，它决定了基于对象分类的应用效率。图像分割算法是将图像分割成有意义的区域，这些区域的内部对于遥感数据的光谱、纹理和空间特征都是均匀的。图像分割在遥感领域有许多应用，如分水岭变换、区域增长和Markov随机场模型、模糊图像区域等。
[0004]分水岭变换是一种基于数学形态学的图像分割方法。分水岭的概念来源于地形，将真实图像中每个像素的灰度值视为物理高程。集水区从地表局部极小点开始浸入水。当来自不同集水区的水汇合时，水坝被建成分水岭。将分水岭变换应用于实际图像时，由于噪声和量化误差的影响，会产生过分割。

技术实现思路

[0005]本专利技术的目的是为了解决现有技术中存在对高分辨率遥感图像产生过分割以及分割区域边界精度低的缺点，而提出的一种高空间分辨率遥感图像的混合分类方法。该高空间分辨率遥感图像的混合分类方法对高空间分辨率遥感数据分类具有良好的处理效果。
[0006]为了实现上述目的，本专利技术采用了如下技术方案：
[0007]设计一种高空间分辨率遥感图像的混合分类方法，具体包括如下步骤：
[0008]S1、基于均值漂移滤波方法对遥感图像进行低通滤波，获取低通滤波的卫星波段：波段1、波段2、
……
波段n，对各滤波段进行均值漂移滤波；
[0009]S2、采用Canny边缘检测方法对步骤1中的低通滤波的卫星波段进行边缘增强，获得遥感图像边缘梯度结果：梯度1、梯度2、
……
梯度n，提取遥感图像的边缘强度；
[0010]S3、基于S2的结果，采用矢量场模型进行复合梯度，得到了多光谱边缘强度，获得遥感图像总的边缘强度信息；
[0011]S4、利用分水岭变换方法，在遥感图像中采用自动提取内部和外部标记的算法进行图像分割；
[0012]S5、采用支持向量机方法对S4分割完成后的遥感影像进行面向像元分类，对遥感影像中的土地覆盖进行逐像素分类，提取建筑物；
[0013]S6、利用投票原理，将S4分割和S5分类的结果相结合，通过遍历每一个分割区域中类别占优来判定该区域的类别归属，最终完成面向对象的遥感分类。
[0014]进一步的，在S1中，均值漂移滤波方法的平滑过程中，引入空间核和多光谱核，具体为：
[0015]设D维空间Rd中n个数据点xi,i＝1,2,
…
,n,用点x中的G(x)计算mG(x)表示的多元均值移位向量，均值漂移的移位为：
[0016][0017]其中x是内核的中心，h是内核的带宽；
[0018]m
h,G
(x)是加权均值和以x为中心的核之间的差值，差值为：
[0019]G(x)＝cg(|x|2)
ꢀꢀ
(2)
[0020]其中g(|x|2)称为核G(x)的轮廓；c为归一化常数；
[0021]空间和多光谱核函数表达式如下所示：
[0022][0023]其中h
Spec
，h
s
分别为光谱核函数和空间核函数的带宽。
[0024]进一步的，采用矢量场模型来获得遥感图像总的边缘强度信息的步骤包括：
[0025]S31：给定f(x，y)为图像，Canny边缘检测器表达式为：
[0026][0027]S32：设定多光谱图像I,则给定点(x0，y0)处的I值是RN中的N维向量，多光谱图像为矢量场，两点P＝(x0，y0)和Q＝(x1，y1)处的图像值的差为：
[0028]ΔI＝I(P)
–
I(Q)
ꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀꢀ
(5)
[0029]S33：当P和Q之间的欧几里得距离d(P,Q)趋于零时，该差值为弧元素：
[0030][0031]其平方范数为：
[0032][0033]S34：允许测量多光谱图像中的变化，在2
×
2矩阵的特征向量的方向上获得二次形式的极值：
[0034][0035]多光谱边缘梯度为：
[0036][0037]其中：
[0038][0039]步骤34：给定样本标记对的标记训练集(X
i
，y
i
){(x1，y1)，...，(x
n
，y
n
)}，其中y
i
∈{
‑
1，+1}，
[0040]是Hilbert空间中的非线性映射，
[0041]进一步的，在S4中，自动提取内部和外部标记的算法是根据边缘强度，通过自适应寻找阈值来进行标记内部和外部区域。
[0042]进一步的，在S6中，支持向量机方法对遥感影像进行面向像元分类的分类工作如下：
[0043][0044][0045]其中w和b定义特征空间中的线性分类器，根据Cover定理推导了非线性映射函数表达式，正则化参数C控制分类器的泛化能力，同时它必须由用户选择，并且是能够处理允许的错误的正松弛变量。
[0046]本专利技术提出的一种高空间分辨率遥感图像的混合分类方法，有益效果在于：本专利技术将图像分割和逐像素分类相结合，完成了高分辨率遥感图像的分类，有效提高分割区域边界精度。具体地：
[0047](1)、使用标记控制分水岭变换进行图像分割，基于多光谱信息，利用Canny边缘检测方法和矢量场模型得到多光谱波段的边缘强度。可有效避免对遥感图像进行分割时容易产生过度分割。
[0048](2)、采用分水岭变换
‑
面向像元分类方法对原始遥感图像进行分类，支持向量机方法比传统的最大似然方法具有更高的分类精度，此外，支持向量机在只有少数训练样本可用的异构类中特别有利，最后，在面向像元分类后，利用投票原则结合高空间分辨率遥感影像的最终分类结果。
附图说明
[0049]附图用来提供对本专利技术的进一步理解，并且构成说明书的一部分，与本专利技术的实施例一起用于解释本专利技术，并不构成对本专利技术的限制。在附图本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种高空间分辨率遥感图像的混合分类方法，包括，其特征在于，具体包括如下步骤：S1、基于均值漂移滤波方法对遥感图像进行低通滤波，获取低通滤波的卫星波段：波段1、波段2、
……
波段n，对各滤波段进行均值漂移滤波；S2、采用Canny边缘检测方法对步骤1中的低通滤波的卫星波段进行边缘增强，获得遥感图像边缘梯度结果：梯度1、梯度2、
……
梯度n，提取遥感图像的边缘强度；S3、基于S2的结果，采用矢量场模型进行复合梯度，得到了多光谱边缘强度，获得遥感图像总的边缘强度信息；S4、利用分水岭变换方法，在遥感图像中采用自动提取内部和外部标记的算法进行图像分割；S5、采用支持向量机方法对S4分割完成后的遥感影像进行面向像元分类，对遥感影像中的土地覆盖进行逐像素分类，提取建筑物；S6、利用投票原理，将S4分割和S5分类的结果相结合，通过遍历每一个分割区域中类别占优来判定该区域的类别归属，最终完成面向对象的遥感分类。2.根据权利要求1所述的一种高空间分辨率遥感图像的混合分类方法，其特征在于，在S1中，均值漂移滤波方法的平滑过程中，引入空间核和多光谱核，具体为：设D维空间Rd中n个数据点xi,i＝1,2,
…
,n,用点x中的G(x)计算mG(x)表示的多元均值移位向量，均值漂移的移位为：其中x是内核的中心，h是内核的带宽；m
h,G
(x)是加权均值和以x为中心的核之间的差值，差值为：G(x)＝cg(|x|2)
ꢀꢀꢀꢀꢀꢀꢀ
(2)其中g(|x|2)称为核G(x)的轮廓；c为归一化常数；空间和多光谱核函数表达式如下所示：其中h
Spec
，h
s
分别为光谱核函数和空间核函数的带宽。3.根据权利要求1所述的一种高空间...

【专利技术属性】
技术研发人员：乔延春，吴钰林，
申请(专利权)人：江苏方寸图信息技术有限公司，
类型：发明
国别省市：