发明名称 一种遥感图像分类方法
摘要 一种遥感图像分类方法,其步骤为:①在参考遥感图像上,分别在不同类型的区域中,截取小图作为训练样本;②设计加博(Gabor)滤波器,将频率域分解至不同的尺度与角度,对训练样本进行滤波,得到滤波图像集合;③计算每个训练样本在不同尺度下的排序直方谱特征向量;④通过对训练样本进行聚类,得到最佳分类尺度;⑤将待分类遥感图像分成子块,使用④中得到的尺度,计算子块的排序直方谱特征向量;⑥使用分类器依次对子块进行分类,最终得到待分类遥感图像的分类结果。本发明能对遥感图像有效分类,操作简便,分类效率高。
申请公布号 CN101067659B 申请公布日期 2010.08.04
申请号 CN200710052417.2 申请日期 2007.06.08
申请人 华中科技大学 发明人 曹治国;范绎;卓问;肖阳
分类号 G01S17/89(2006.01)I 主分类号 G01S17/89(2006.01)I
代理机构 华中科技大学专利中心 42201 代理人 曹葆青
主权项 一种遥感图像分类方法,其特征是,该方法包括以下步骤:(1)在参考遥感图像上截取若干小块作为训练样本;(2)设计加博(Gabor)滤波器组,将频率域分解至S个尺度和T个角度,使用滤波器组对所有训练样本进行滤波,得到训练样本的滤波图像集合,记为P;(3)从P中提取基于不同尺度的滤波图像集合,通过分别计算该集合中滤波图像最小值和最大值的平均值,得到该集合对应的直方图统计范围[Qmin(l),Qmax(l)],l=1,2,…,S,将该范围平均分成小区间,得到相应的直方图区间划分方式:[Qmin(l)+range*(z-1),Qmin(l)+range*z],z=1,2,…,c其中, <mrow> <mi>range</mi> <mo>=</mo> <mfrac> <mrow> <msub> <mi>Q</mi> <mi>max</mi> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>-</mo> <msub> <mi>Q</mi> <mi>min</mi> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> </mrow> <mi>c</mi> </mfrac> <mo>,</mo> </mrow>c为区间个数, <mrow> <msub> <mi>Q</mi> <mi>min</mi> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>A</mi> <mo>*</mo> <mi>T</mi> </mrow> </mfrac> <munder> <mi>&Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <munder> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <munder> <mi>&Sigma;</mi> <mrow> <mi>r</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <mi>min</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> <mrow> <mi>l</mi> <mo>,</mo> <mi>r</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow> <mrow> <msub> <mi>Q</mi> <mi>max</mi> </msub> <mrow> <mo>(</mo> <mi>l</mi> <mo>)</mo> </mrow> <mo>=</mo> <mfrac> <mn>1</mn> <mrow> <mi>A</mi> <mo>*</mo> <mi>T</mi> </mrow> </mfrac> <munder> <mi>&Sigma;</mi> <mrow> <mi>m</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <munder> <mi>&Sigma;</mi> <mrow> <mi>n</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <munder> <mi>&Sigma;</mi> <mrow> <mi>r</mi> <mo>=</mo> <mn>1</mn> </mrow> </munder> <mi>max</mi> <mrow> <mo>(</mo> <msubsup> <mi>F</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> <mrow> <mi>l</mi> <mo>,</mo> <mi>r</mi> </mrow> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow> <mrow> <msubsup> <mi>F</mi> <mrow> <mi>m</mi> <mo>,</mo> <mi>n</mi> </mrow> <mrow> <mi>l</mi> <mo>,</mo> <mi>r</mi> </mrow> </msubsup> <mo>=</mo> <mi>con</mi> <mrow> <mo>(</mo> <msubsup> <mi>I</mi> <mi>m</mi> <mi>n</mi> </msubsup> <mo>*</mo> <msubsup> <mi>G</mi> <mi>l</mi> <mi>r</mi> </msubsup> <mo>)</mo> </mrow> <mo>,</mo> </mrow>Imn表示第m类事物的第n个训练样本,Glr表示第l个尺度的第r个角度的Gabor滤波器,con()为卷积运算,A为训练样本总数;(4)按照以下步骤提取每个训练样本基于不同尺度的排序直方谱特征向量:(4.1)从P中提取一个训练样本基于某尺度的滤波图像集合,根据基于该尺度的滤波图像集合的直方图区间划分方式,得到其对应的直方图集合和特征值集合;(4.2)对直方图集合里的元素按照特征值集合里对应元素的值的大小进行排序,得到训练样本基于该尺度的排序直方谱特征向量;(4.3)按照步骤(4.1)-(4.2)的方式,遍历所有训练样本和尺度,得到每个训练样本基于不同尺度的排序直方谱特征向量;(5)使用χ2统计距离,根据步骤(4)得到的排序直方谱特征向量对训练样本进行聚类,得到最佳分类尺度;(6)将待分类图像分成小块,按照步骤(4.1)-(4.2)的方式计算每个小块基于最佳分类尺度的排序直方谱特征向量;根据排序直方谱特征向量,使用分类器依次对小块分类并标记,最终得到整幅图像的分类结果。
地址 430074 湖北省武汉市洪山区珞瑜路1037号