专利名称：一种采用图像分类技术从遥感图像中进行滑坡检测的方法
技术领域：
本发明涉及一种采用图像分类技术从遥感图像中进行滑坡检测的方法，可以应用于各种遥感图像的滑坡检测。
背景技术：
滑坡是一类破坏力极强的地质灾害，它常常给工农业生产以及人民生命财产造成巨大损失。因此，滑坡研究也越来越受到人们的重视。早期的滑坡检测主要采用传统的地面调查方法，该方法野外工作强度很大，特别在危险及气候恶劣地区，工作效率十分低下。 到20世纪90年代末，立体镜航空照片解译配合一定的地面验证成为滑坡检测和制图最常用的方法，然后将解译结果通过手工转绘到相应比例尺的地形图上，制作滑坡分布图，该方法使部分野外工作转移到室内，在一定程度上提高了工作效率，减轻了野外工作强度。近年来，开始利用ArcGIS、CorelDRAW等软件平台结合数字高程模型进行滑坡检测和分析，但该方法需要获取三维立体地形学信息，计算量很大。
遥感作为一种滑坡调查和监测手段，随着传感器的不断发展和遥感图像分辨率的不断提高，越来越受到地质灾害研究人员的关注。从遥感图像中进行滑坡检测是滑坡分析、 评价、预测和监测的基础，它利用滑坡发生处的影像光谱信息、地形地貌和形态特征等，对遥感图像进行分析和判别，获取滑坡灾害发生范围，从而达到灾害调查和制图的目的。近年来，随着高分辨率遥感图像在资源环境监测、灾害管理等领域的广泛应用，高分辨率遥感图像信息提取方法以及遥感图像分类方法已成为近年来研究的热点。现有的采用遥感图像进行滑坡检测的方法主要有基于对象的高分辨率遥感图像滑坡检测方法和基于监督分类的遥感图像滑坡检测方法。但是，基于对象的高分辨率遥感图像滑坡检测方法需要对多期遥感图像进行多尺度、多层次分割，检测效果直接依赖于图像分割的结果，总体检测率不高；基于监督分类的遥感图像滑坡检测方法只能实现将图像划分为滑坡解译所需要的简单几种类型，且检测过程中仍需借助于遥感专业处理软件ENVI，检测效率和检测精度都有待提闻。发明内容
要解决的技术问题
为了避免现有技术的不足之处，本发明提出了一种采用图像分类技术从遥感图像中进行滑坡检测的方法，能够从遥感图像中进行滑坡检测，且具有很高的检测正确率。
本发明的思想在于(I)采用基于重叠的面积均分法对预处理后的遥感图像进行分块，得到面积相等的正方形图像块，然后将得到的图像块分成两个集合——训练集和测试集；⑵提取训练集和测试集中所有图像块的SIFT特征，然后对训练集中的SIFT特征采用k-means聚类方法得到单词和词典；(3)用BoVW模型表示训练集和测试集中的每个图像块；(4)利用pLSA模型提取每个图像块的主题，最后用KNN分类器将测试集中所有的图像块分成滑坡和非滑坡两类，从而实现对遥感图像的滑坡检测。
技术方案一种采用图像分类技术从遥感图像中进行滑坡检测的方法，其特征在于步骤1 预处理采用加权平均法对遥感图像的RGB三个分量进行加权平均得到灰度图像，然后利用灰度线性变换函数将灰度图像的灰度范围映射到[/mm,/max]的灰度区间，得到预处理后的图像；其中加权平均法计算公式为f(i, j) = 0. 3R(i, j)+0. 59G(i, j)+0. llB(i, j)， f(i，j)为加权平均法得到的(i，j)像素点的灰度值，R(i，j)、G(i，j)和B(i，j)分别为(i，
j)像素点的RGB三个分量值；灰度线性变换函数为/α/>=(L· -LjxffJ)jfmm +/mm，
J max J min
为(υ)像素点经过灰度映射后的灰度值，fmin、f_分别为灰度映射前的灰度图像的最小灰度值和最大灰度值；步骤2 图像分块采用anxaii的滑动窗口对预处理后的灰度图像进行重叠分块，相邻图像块的重叠区域为m个像素，得到图像块的集合D = U1, L，IJ ；从集合D中选取Nlandslide个具有滑坡特征的滑坡类图像块和Nnrarfandslide个不具有滑坡特征的非滑坡类图像块组成训练集M = L ’ _J ；将集合D中剩余的图像块组成测试集万=·^,L ,^teJ;其中1表示集合D中的图像块；下标N为集合D中图像块的个数；m的取值范围为 150 300
T T，Y为遥感图像的空间分辨率，单位为米；a表示训练集中的图像块；下标Nteaining
为训练集中图像块的个数，且Ntraining = Nlandslide+Nnon_landslide ‘Niandsiide/Nnon_landslide ^ 0. 5 ;b ^ 示测试集中的图像块；下标Ntest为测试集中图像块的个数；步骤3 单词和词典的形成共分两个步骤，具体过程如下步骤al 提取集合D中每个图像块的SIFT特征，具体过程如下首先计算每个图像块的图像卷积‘ =/*^Τ^(χ2+/)/2σ"2唭中ο = 0. 5，*表示卷积运算；然后计算每个图像块的图像金字塔GX^7 =/σ ^-J^e-^2^2，图像金字塔共0
“2πσ
组，每组S层，下一组的第一幅图像由上一组的最后一幅图像降采样得到，采样因子为2 ；其中σ = σ 0 . 2。+s/s，σ 0 = 1. 6 · 21/s, ο = 0，K 0-1, s = 0，K S-I, S = 3，0 = Iog2 2m ；再对每个图像块中同一组相邻的GSS。求差分得到DOG。，将DOG。每个像素点的值分别和它同尺度的8个相邻点和上下相邻尺度对应的9X2个点共沈个点比较，如果该像素点的值为极小值或者极大值，则该像素点为图像显著点，显著点周围以ο为半径的区域为显著区域，由此可以得到一系列的图像显著点；最后计算每个图像块的梯度图像卷积1 %eW)/2.0,^，其中 7&为允的梯度图像；然后在t上计算以图像显著点为中心、。为半径的圆形区域的梯度
°nuηuη
方向直方图，梯度方向直方图共分36个区间，每个区间10度，梯度方向直方图峰值的方向区域为该显著点的主方向；最后在t上将以每个图像显著点为中心、σ为半径的圆形区
u η域，按主方向及其垂直方向等分成16个区域，在每个小区域中分别统计梯度方向直方图， 梯度方向直方图共分8个区间，每个区间45度，并将每个方向直方图的幅值量化到
区间，得到16X8 = 128维的SIFT特征；由此，每一个图像块I可以得到一系列的SIFT特征；
步骤bl :设定聚类中心数目M，利用K-means聚类算法对训练集中所有图像块的 SIFT特征进行聚类，得到含有M个单词的词典W= {Wl，L，wM}；
其中M表示词典中单词的个数，取值范围为500 2500 表示词典中的单词；
步骤4 :用BoVW模型表示图像块将训练集和测试集中所有图像块的每个SIFT 特征映射为词典中的一个单词；然后统计每一个图像块中所有单词的频率，得到一个 NtrainingXM 的矩阵Ntrammg =(Kli^)y.)Wt mmgXM和一个 NtestXM 的矩阵Ntest =(/7( ,'%.) xM ；
其中mh，Wj)为图像块Bi的BoVW模型表示，代表单词Wj在图像块Bi中出现的频率；下标ij表示元素Wapwj)在Ntraining中处于第i行、第j列，i = I, 2, L, Ntraining, j = 1，2，L，M ；n(bk, Wj)为图像块bk的BoVW模型表示，代表单词Wj在图像块bk中出现的频率； 下标kj表示元素n(bk, W」)在Ntejst中处于第k行、第j列，k = 1,2, L, Ntest ；
步骤5 :用pLSA模型提取图像块的主题共分两个步骤，具体过程如下
步骤a2 :计算训练集图像块的P (W11 zt)和P (zt | a,)，具体过程如下
权利要求
1. 一种采用图像分类技术从遥感图像中进行滑坡检测的方法，其特征在于 步骤1 预处理采用加权平均法对遥感图像的RGB三个分量进行加权平均得到灰度图像，然后利用灰度线性变换函数将灰度图像的灰度范围映射到[/mm,/max]的灰度区间，得到预处理后的图像；其中加权平均法计算公式为 f(i，j) = 0. 3R(i, j)+0. 59G(i, j)+0. llB(i，j)，f (i，j) 为加权平均法得到的(i，j)像素点的灰度值，R(i，j)、G(i，j)和B(i，j)分别为(i，j)像素点的RGB三个分量值；灰度线性变换函数为/α _；+)=cL -Ljxf fJ)jfmm +L·，J max J minJ \i J)为(i，j)像素点经过灰度映射后的灰度值，fmin、fmax分别为灰度映射前的灰度图像的最小灰度值和最大灰度值；步骤2 图像分块采用anXail的滑动窗口对预处理后的灰度图像进行重叠分块，相邻图像块的重叠区域为m个像素，得到图像块的集合D= {Ii;L,In}；从集合D中选取Nlandslide 个具有滑坡特征的滑坡类图像块和Nn。n_landslide个不具有滑坡特征的非滑坡类图像块组成训练集X = P1,L , _」；将集合D中剩余的图像块组成测试集万=·^,L ,^teJ;其中1表示集合D中的图像块；下标N为集合D中图像块的个数；m的取值范围为 150 300~y 7，γ为遥感图像的空间分辨率，单位为米；a表示训练集中的图像块；下标Nteaining为训练集中图像块的个数，且Ntraining = Nlandslide+Nnon_landslide ‘Niandsiide/Nono_landslide ^ 0. 5 ;b ^ 示测试集中的图像块；下标Ntest为测试集中图像块的个数； 步骤3 单词和词典的形成共分两个步骤，具体过程如下 步骤al 提取集合D中每个图像块的SIFT特征，具体过程如下首先计算每个图像块的图像卷积‘ =/*^Τ^(χ2+/)/2σ"2 ；其中σ n = 0. 5，*表示卷积运算；然后计算每个图像块的图像金字塔^-J^e-^2^2，图像金字塔共0组，每“2πσ组S层，下一组的第一幅图像由上一组的最后一幅图像降采样得到，采样因子为2;其中σ = σ0. 2°+s/s，σ 0 = 1. 6 · 21/s, ο = 0，K 0-1, s = 0，K S-I, S = 3，0 = Iog2 2m ；再对每个图像块中同一组相邻的GSS。求差分得到DOG。，将DOG。每个像素点的值分别和它同尺度的8个相邻点和上下相邻尺度对应的9X2个点共沈个点比较，如果该像素点的值为极小值或者极大值，则该像素点为图像显著点，显著点周围以ο为半径的区域为显著区域，由此可以得到一系列的图像显著点；最后计算每个图像块的梯度图像卷积1 =ι0σ %一2”2)/2.(1.⑷2，其中:1“Ιπ-^.^σ)1Gt7n为/σ的梯度图像；然后在t上计算以图像显著点为中心、。为半径的圆形区域的梯度方向u ηu η直方图，梯度方向直方图共分36个区间，每个区间10度，梯度方向直方图峰值的方向区域为该显著点的主方向；最后在t上将以每个图像显著点为中心、。为半径的圆形区域，按u η主方向及其垂直方向等分成16个区域，在每个小区域中分别统计梯度方向直方图，梯度方向直方图共分8个区间，每个区间45度，并将每个方向直方图的幅值量化到
区间， 得到16X8 = 128维的SIFT特征；由此，每一个图像块I可以得到一系列的SIFT特征； 步骤bl :设定聚类中心数目M，利用K-means聚类算法对训练集中所有图像块的SIFT 特征进行聚类，得到含有M个单词的词典W= (W1, L, wM}；其中M表示词典中单词的个数，取值范围为500 2500 ;w表示词典中的单词；步骤4 :用BoVW模型表示图像块将训练集和测试集中所有图像块的每个SIFT特征映射为词典中的一个单词；然后统计每一个图像块中所有单词的频率，得到一个NtrainingXM的矩阵Ntrammg =( (%，％HmmgxM和一个 NtestXM 的矩阵Ntest xM ；其中为图像块Bi的BoVW模型表示，代表单词Wj在图像块Bi中出现的频率； 下标 ij 表示元素 Wai, Wj)在 Ntraining 中处于第 i 行、第 j 列，i = 1,2, L, Ntraining, j = 1,2, L, M ；n(bk, Wj)为图像块bk的BoVW模型表示，代表单词Wj在图像块bk中出现的频率；下标 kj表示元素n(bk, W」)在Ntejst中处于第k行、第j列，k = 1,2, L, Ntest ；步骤5 :用pLSA模型提取图像块的主题共分两个步骤，具体过程如下步骤a2 :计算训练集图像块的P (W11 zt)和P (zt I a,)，具体过程如下初始化P(Zi)=丄、/1(w, I zt)=丄和I Zt) T“计算 P(AlH):MP(Zt)Piwi I z^Piai I zt)training将 P (Zt,■)代入公式尸O. I zt)P(ar \Zt)'-iHiai7WW(ZtIai7Wj)Yjs^nias,W])P{zt\as,W])miH,)计算P (Wj I Zt)、P (at I zt)和P (zt),然后再次计算P (zt I ai; Wj)并重复上述过程，直到Izf)]的对数似然函数期望值的增加量小于指定的阈值S时停止迭代，得到P(w」Zt);然后，得到P(Zt I ai) =,---;其中Z表示主题；P (Wj I Zt)表示单词Wj在主题Zt中出现的概率；P (Zt I Sii)表示主题Zt 在图像块31中出现的概率；P(zt)表示主题Zt在所有主题中出现的概率；P(ai|zt)表示图像块ai在主题Zt中出现的概率；P (Zt I BijWj)表示主题Zt在图像块Bi和单词Wj条件下的后验概率；下标t和m表示主题的序号，t = I, 2, L, T, m = I, 2, L, T, T表示主题个数，T的取值范围为10 25 ;下标r表示单词的序号，r = 1,2, L, M ;下标s表示训练集图像块的序号，,2，L,'training;log表示对数运算符；S为设定的阈值，取值范围为0. I 0. 5步骤b2 :计算测试集图像块的P (zt I bk)，具体过程如下
全文摘要
本发明涉及一种采用图像分类技术从遥感图像中进行滑坡检测的方法，首先采用基于重叠的面积均分法对预处理后的遥感图像进行分块，得到面积相等的正方形图像块，并将得到的图像块分成两个集合——训练集和测试集；其次，提取训练集和测试集中所有图像块的SIFT特征，对训练集中的SIFT特征采用k-means聚类方法得到单词和词典；然后，用BoVW模型表示训练集和测试集中的每个图像块；最后，利用pLSA模型提取每个图像块的主题，并用KNN分类器将测试集中的图像块分成滑坡和非滑坡两类，从而实现对遥感图像的滑坡检测。该方法运算量小、检测效率高，同时又具有很高的检测正确率。
文档编号G06K9/62GK102542295SQ20121000475
公开日2012年7月4日 申请日期2012年1月8日 优先权日2012年1月8日
发明者任波, 余博, 程塨, 赵天云, 郭雷, 钱晓亮, 韩军伟 申请人:西北工业大学