一种遥感图像分类方法技术

一种遥感图像分类方法，其步骤为：①在参考遥感图像上，分别在不同类型的区域中，截取小图作为训练样本；②设计加博（Ｇａｂｏｒ）滤波器，将频率域分解至不同的尺度与角度，对训练样本进行滤波，得到滤波图像集合；③计算每个训练样本在不同尺度下的排序直方谱特征向量；④通过对训练样本进行聚类，得到最佳分类尺度；⑤将待分类遥感图像分成子块，使用④中得到的尺度，计算子块的排序直方谱特征向量；⑥使用分类器依次对子块进行分类，最终得到待分类遥感图像的分类结果。本发明专利技术能对遥感图像有效分类，操作简便，分类效率高。

【技术实现步骤摘要】

本专利技术属于遥感图像分类
，具体涉及一种利用排序直方谱特征对遥感图像进行分类的方法。
技术介绍
遥感图像是按一定的比例尺，客观真实地记录和反映地表物体的电磁辐射的强弱信息，是遥感探测所获得的遥感信息资料的一种表现形式。所以遥感技术的核心问题是根据地物辐射电磁辐射强弱在遥感图像上表现的特征，判读识别地面物体的类属及其分布特征。在遥感图像技术的研究中，无论是专业信息提取，运动变化预测，还是专题地图制作和遥感数据库的建立等都离不开分类，对照地面地物类型，便可以从图像上地物的识别。随着近年来计算机技术的飞速发展，计算机识别分类成了遥感技术应用的一个重要组成部分。遥感图像数据类别多，合混度大，如何解决多类别分类识别并满足一定的分类精度，是当前遥感图像研究中的一个关键问题，也是人们关注的焦点。在目前的遥感图像分类应用中，用的较多的方法有共生矩阵，行程长度，加博(Gabor)的一阶二阶矩。其中，共生矩阵和行程长度由于只对图像在空间域的特性进行分析，而Gabor的一阶二阶矩的统计量过分简单，所以分类能力都非常有限。
技术实现思路
本专利技术的目的在于，提供，该方法能达到分类的高准确率。本专利技术按照以下步骤完成(1)在参考遥感图像上截取若干小块作为训练样本；(2)设计加博(Gabor)滤波器组，将频率域分解至S个尺度和T个角度，使用滤波器组对所有训练样本进行滤波，得到训练样本的滤波图像集合，记为P；(3)从P中提取基于不同尺度的滤波图像集合，通过分别计算该集合中滤波图像最小值和最大值的平均值，得到该集合对应的直方图统计范围(l＝1，2，...，S)，将该范围平均分成小区间，得到相应的直方图区间划分方式(z＝1，2，...，c)其中，range=Qmax(l)-Qmin(l)c,]]>c为区间个数；(4)按照以下步骤提取每个训练样本基于不同尺度的排序直方谱特征向量(4.1)从P中提取一个训练样本基于某尺度的滤波图像集合，根据基于该尺度的滤波图像集合的直方图区间划分方式，得到其对应的直方图集合和特征值集合；(4.2)对直方图集合里的元素按照特征值集合里对应元素的值的大小进行排序，得到训练样本基于该尺度的排序直方谱特征向量；(4.3)按照步骤(4.1)-(4.2)的方式，遍历所有训练样本和尺度，得到每个训练样本基于不同尺度的排序直方谱特征向量；(5)使用χ2统计距离，根据步骤(4)得到的排序直方谱特征向量对训练样本进行聚类，得到最佳分类尺度；(6)将待分类图像分成小块，按照步骤(4.1)-(4.2)的方式计算每个小块基于最佳分类尺度的排序直方谱特征向量；根据排序直方谱特征向量，使用分类器依次对小块分类并标记，最终得到整幅图像的分类结果。本专利技术同时对空间域和频率域进行分析，能够获取更多有助于分类的图像纹理信息；计算排序直方谱特征向量，并将其作为本专利技术所需统计量，使得分类结果更为精确。附图说明图1为本专利技术的总体流程图；图2为提取每个训练样本基于不同尺度的排序直方谱特征向量流程图；图3为提取训练样本基于不同尺度的排序直方谱特征向量的示例图 其中，图3a为农田训练样本；图3b为该训练样本在第4个尺度下的6个直方图；图3c为该训练样本在第4个尺度下的排序直方谱；图4为本专利技术对待分类图像分类的流程图；图5为本专利技术关于遥感大图子块划分方式的示意图；图6为本专利技术与一般方法的分类准确率比较图。具体实施例方式下面结合附图和实施例对本专利技术做详细说明。(1)截取训练样本。假定待分类遥感图像有M类区域需要划分。在参考遥感图像上，分别在每个类型所属区域中截取若干小图作为训练样本，令第m类事物的第n个训练样本为Imn。要求每一类的训练样本能够体现该类区域的所有基本特征。从理论上说，各训练样本大小可以不相同，但是为了编制程序的方便和达到更稳定的分类效果，本专利技术建议使用大小相同的小图，本专利技术取小图大小为50像素*50像素。(2)使用加博(Gabor)滤波器组得到滤波图像集合。(2.1)设计Gabor滤波器组。滤波器组的设计必须满足设计的滤波器组能够覆盖频率域。设计将频率域分解至S个尺度和T个角度的Gabor滤波器组，则该滤波器组一共有S*T个滤波器，令第l个尺度，第r个角度的Gabor滤波器为Glr。一般采用4个尺度和6个角度。例如，采用4个尺度，其对应的中心频率分别为0.05，0.1，0.2和0.4，6个角度分别为0，30，60，90，120，150度，一共24个滤波器。(2.2)得到训练样本的滤波图像集合。对应于每一个训练样本和每一个滤波器，可以得到一幅滤波图像。例如，对于训练样本Imn和滤波器Glr，通过卷积运算，得到滤波后的图像Fm,nl,r=con(Inm*Glr).]]>利用Gabor滤波器组对所有训练样本进行滤波，得到包括A*S*T(A为训练样本总数)幅滤波图像的滤波图像集合，记为P。(3)获取基于不同尺度的滤波图像集合，计算其对应的直方图区间划分方式。通过不同尺度的Gabor滤波器，可以得到差异很大的滤波图像。一幅滤波图像的区间分布就叫做这幅滤波图像的直方图，任何一幅图像通过一组滤波器，都可以得到一组滤波图像，以及相应的一组直方图。从P中提取基于不同尺度的滤波图像集合，对每个集合中的滤波图像采用相同的统计区间计算直方图，步骤如下(3.1)确定基于不同尺度的滤波图像集合的直方图统计范围。令Qmin(l)=1A*QΣm=1Σn=1Σr=1min(Fm,nl,r),]]>Qmax(l)=1A*QΣm=1Σn=1Σr=1max(Fm,nl,r),]]>(l＝1，2，...，S)为基于第l个尺度的滤波图像集合的直方图统计范围。(3.2)确定基于不同尺度的滤波图像集合的直方图统计区间。将步骤(3.1)得到的直方图统计范围均匀分成c个区间(一般取c为10)，则第l个尺度对应的直方谱区间长度range=Qmax(l)-Qmin(l)c,]]>其第z个直方谱统计区间是(z＝1，2，...，c)。(4)提取每个训练样本基于不同尺度的排序直方谱特征向量。训练样本Imn在第l个尺度下提取排序直方谱特征向量，按照以下步骤完成(4.1)计算训练样本Imn基于第l个尺度的滤波图像直方图集合。定义滤波图像Fm，nl，r的直方图为Hm,nl,r=(Hm,nl,r(1),Hm,nl,r(2),...,Hm,nl,r(c)),]]>其中Hm,nl,r(z)=sum(Qmin(l)+range*(z-1)<Fm,nl,r(x,y)<Qmin(l)+range*z)(z=1,2,...,c).]]>同时，定义Um,nl,r=max{Hm,nl,r(z),z=1,2...,c}]]>为Hm，nl，r的特征值。训练样本Imn通过第l个尺度下不同角度的滤波器滤波后，可以得到T幅滤波图像的集合为{Fm，nl，1，Fm，nl，2，...，Fm，nl，T}，这些滤波图像相应的直方图的集合为{Hm，nl，1，Hm，nl，2，...，Hm，nl，T}，对应的特征值的集合为{Um，nl，1，Um本文档来自技高网...

【技术保护点】
一种遥感图像分类方法，其特征是，该方法包括以下步骤：（１）在参考遥感图像上截取若干小块作为训练样本；（２）设计加博（Ｇａｂｏｒ）滤波器组，将频率域分解至Ｓ个尺度和Ｔ个角度，使用滤波器组对所有训练样本进行滤波，得到训练样本的滤波图像集合，记为Ｐ；（３）从Ｐ中提取基于不同尺度的滤波图像集合，通过分别计算该集合中滤波图像最小值和最大值的平均值，得到该集合对应的直方图统计范围［Ｑ↓［ｍｉｎ］（ｌ），Ｑ↓［ｍａｘ］（ｌ）］（ｌ＝１，２，…，Ｓ），将该范围平均分成小区间，得到相应的直方图区间划分方式：［Ｑ↓［ｍｉｎ］（ｌ）＋ｒａｎｇｅ＊（ｚ－１），Ｑ↓［ｍｉｎ］（ｌ）＋ｒａｎｇｅ＊ｚ］（ｚ＝１，２，…，ｃ）其中，＊＊＊，ｃ为区间个数；（４）按照以下步骤提取每个训练样本基于不同尺度的排序直方谱特征向量：（４．１）从Ｐ中提取一个训练样本基于某尺度的滤波图像集合，根据基于该尺度的滤波图像集合的直方图区间划分方式，得到其对应的直方图集合和特征值集合；（４．２）对直方图集合里的元素按照特征值集合里对应元素的值的大小进行排序，得到训练样本基于该尺度的排序直方谱特征向量；（４．３）按照步骤（４．１）－（４．２）的方式，遍历所有训练样本和尺度，得到每个训练样本基于不同尺度的排序直方谱特征向量；（５）使用χ↑［２］统计距离，根据步骤（４）得到的排序直方谱特征向量对训练样本进行聚类，得到最佳分类尺度；（６）将待分类图像分成小块，按照步骤（４．１）－（４．２）的方式计算每个小块基于最佳分类尺度的排序直方谱特征向量；根据排序直方谱特征向量，使用分类器依次对小块分类并标记，最终得到整幅图像的分类结果。...

【技术特征摘要】

【专利技术属性】
技术研发人员：曹治国，范绎，卓问，肖阳，
申请(专利权)人：华中科技大学，
类型：发明
国别省市：83[中国|武汉]