[发明专利]一种采用图像分类技术从遥感图像中进行滑坡检测的方法

摘要

本发明涉及一种采用图像分类技术从遥感图像中进行滑坡检测的方法，首先采用基于重叠的面积均分法对预处理后的遥感图像进行分块，得到面积相等的正方形图像块，并将得到的图像块分成两个集合——训练集和测试集；其次，提取训练集和测试集中所有图像块的SIFT特征，对训练集中的SIFT特征采用k-means聚类方法得到单词和词典；然后，用BoVW模型表示训练集和测试集中的每个图像块；最后，利用pLSA模型提取每个图像块的主题，并用KNN分类器将测试集中的图像块分成滑坡和非滑坡两类，从而实现对遥感图像的滑坡检测。该方法运算量小、检测效率高，同时又具有很高的检测正确率。

主权项

1.一种采用图像分类技术从遥感图像中进行滑坡检测的方法，其特征在于：步骤1：预处理：采用加权平均法对遥感图像的RGB三个分量进行加权平均得到灰度图像，然后利用灰度线性变换函数将灰度图像的灰度范围映射到的灰度区间，得到预处理后的图像；其中：加权平均法计算公式为f(i，j)＝0.3R(i，j)+0.59G(i，j)+0.11B(i，j)，f(i，j)为加权平均法得到的(i，j)像素点的灰度值，R(i，j)、G(i，j)和B(i，j)分别为(i，j)像素点的RGB三个分量值；灰度线性变换函数为为(i，j)像素点经过灰度映射后的灰度值，f_min、f_max分别为灰度映射前的灰度图像的最小灰度值和最大灰度值；步骤2：图像分块：采用2m×2m的滑动窗口对预处理后的灰度图像进行重叠分块，相邻图像块的重叠区域为m个像素，得到图像块的集合D＝{I₁，L，I_N}；从集合D中选取N_landslide个具有滑坡特征的滑坡类图像块和N_{non-landslide}个不具有滑坡特征的非滑坡类图像块组成训练集将集合D中剩余的图像块组成测试集B={b1,L,bNtest};]]>其中：I表示集合D中的图像块；下标N为集合D中图像块的个数；m的取值范围为γ为遥感图像的空间分辨率，单位为米；a表示训练集中的图像块；下标N_training为训练集中图像块的个数，且N_training＝N_landslide+N_{non-landslide}，N_landslide/N_{ono-landslide}≈0.5；b表示测试集中的图像块；下标N_test为测试集中图像块的个数；步骤3：单词和词典的形成：共分两个步骤，具体过程如下：步骤a1：提取集合D中每个图像块的SIFT特征，具体过程如下：首先计算每个图像块的图像卷积其中：σ_n＝0.5，*表示卷积运算；然后计算每个图像块的图像金字塔图像金字塔共O组，每组S层，下一组的第一幅图像由上一组的最后一幅图像降采样得到，采样因子为2；其中：σ＝σ₀·2^o+s/S，σ₀＝1.6·2^1/S，o＝0，K O-1，s＝0，K S-1，S＝3，O＝log₂ 2m；再对每个图像块中同一组相邻的GSS_σ求差分得到DOG_σ，将DOG_σ每个像素点的值分别和它同尺度的8个相邻点和上下相邻尺度对应的9×2个点共26个点比较，如果该像素点的值为极小值或者极大值，则该像素点为图像显著点，显著点周围以σ为半径的区域为显著区域，由此可以得到一系列的图像显著点；最后计算每个图像块的梯度图像卷积其中：为的梯度图像；然后在上计算以图像显著点为中心、σ为半径的圆形区域的梯度方向直方图，梯度方向直方图共分36个区间，每个区间10度，梯度方向直方图峰值的方向区域为该显著点的主方向；最后在上将以每个图像显著点为中心、σ为半径的圆形区域，按主方向及其垂直方向等分成16个区域，在每个小区域中分别统计梯度方向直方图，梯度方向直方图共分8个区间，每个区间45度，并将每个方向直方图的幅值量化到[0，255]区间，得到16×8＝128维的SIFT特征；由此，每一个图像块I可以得到一系列的SIFT特征；步骤b1：设定聚类中心数目M，利用K-means聚类算法对训练集中所有图像块的SIFT特征进行聚类，得到含有M个单词的词典W＝{w₁，L，w_M}；其中：M表示词典中单词的个数，取值范围为500～2500；w表示词典中的单词；步骤4：用BoVW模型表示图像块：将训练集和测试集中所有图像块的每个SIFT特征映射为词典中的一个单词；然后统计每一个图像块中所有单词的频率，得到一个N_training×M的矩阵和一个N_test×M的矩阵Ntest=(n(bk,wj)kj)Ntest×M;]]>其中：n(a_i，w_j)为图像块a_i的BoVW模型表示，代表单词w_j在图像块a_i中出现的频率；下标ij表示元素n(a_i，w_j)在N_training中处于第i行、第j列，i＝1，2，L，N_training，j＝1，2，L，M；n(b_k，w_j)为图像块b_k的BoVW模型表示，代表单词w_j在图像块b_k中出现的频率；下标kj表示元素n(b_k，w_j)在N_test中处于第k行、第j列，k＝1，2，L，N_test；步骤5：用pLSA模型提取图像块的主题：共分两个步骤，具体过程如下：步骤a2：计算训练集图像块的P(w_j|z_t)和P(z_t|a_i)，具体过程如下：初始化P(zt)=1T,]]>P(wj|zt)=1M]]>和P(ai|zt)=1Ntraining,]]>计算P(zt|ai,wj)=P(zt)P(wj|zt)P(ai|zt)Σm=1TP(zm)P(wj|zm)P(ai|zm),]]>将P(z_t|a_i，w_j)代入公式P(wj|zt)=Σi=1Ntrainingn(ai,wj)P(zt|ai,wj)Σr=1MΣi=1Ntrainingn(ai,wr)P(zt|ai,wr)P(ai|zt)=Σj=1Mn(ai,wj)P(zt|ai,wj)Σj=1MΣs=1Ntrainingn(as,wj)P(zt|as,wj)P(zt)=Σi=1NtrainingP(ai|zt)Σm=1TΣi=1NtrainingP(ai|zm)]]>计算P(w_j|z_t)、P(a_i|z_t)和P(z_t)，然后再次计算P(z_t|a_i，w_j)并重复上述过程，直到L=Σi=1NtrainingΣj=1Mn(ai,wj)log[Σt=1TP(zt)P(wj|zt)P(ai|zt)]]]>的对数似然函数期望值的增加量小于指定的阈值δ时停止迭代，得到P(w_j|z_t)；然后，得到P(zt|ai)=Σj=1Mn(ai,wj)P(zt|ai,wj)Σj=1Mn(ai,wj);]]>其中：z表示主题；P(w_j|z_t)表示单词w_j在主题z_t中出现的概率；P(z_t|a_i)表示主题z_t在图像块a_i中出现的概率；P(z_t)表示主题z_t在所有主题中出现的概率；P(a_i|z_t)表示图像块a_i在主题z_t中出现的概率；P(z_t|a_i，w_j)表示主题z_t在图像块a_i和单词w_j条件下的后验概率；下标t和m表示主题的序号，t＝1，2，L，T，m＝1，2，L，T，T表示主题个数，T的取值范围为10～25；下标r表示单词的序号，r＝1，2，L，M；下标s表示训练集图像块的序号，s＝1，2，L，N_training；log表示对数运算符；δ为设定的阈值，取值范围为0.1～0.5；步骤b2：计算测试集图像块的P(z_t|b_k)，具体过程如下：初始化P(zt)=1T,]]>P(bk|zt)=1Ntest,]]>计算P(zt|bk,wj)=P(zt)P(wj|zt)P(bk|zt)Σm=1TP(zm)P(wj|zm)P(bk|zm),]]>将P(z_t|b_k，w_j)代入公式P(bk|zt)=Σj=1Mn(bk,wj)P(zt|bk,wj)Σj=1MΣd=1Ntestn(bd,wj)P(zt|bd,wj)P(zt)=Σk=1NtestP(bk|zt)Σm=1TΣk=1NtestP(bk|zm)]]>计算P(b_k|z_t)和P(z_t)，然后再次计算P(z_t|b_k，w_j)并重复上述过程，直到L=Σk=1NtestΣj=1Mn(bk,wj)log[Σt=1TP(zt)P(wj|zt)P(bk|zt)]]]>的对数似然函数期望值的增加量小于指定的阈值δ时停止迭代，计算其中：P(z_t|b_k)表示主题z_t在图像块b_k中出现的概率；P(b_k|z_t)表示图像块b_k在主题z_t中出现的概率；P(z_t|b_k，w_j)表示主题z_t在图像块b_k和单词w_j条件下的后验概率；下标d表示测试集图像块的序号，d＝1，2，L，N_test；步骤6：KNN分类：计算P(z_t|b_k)与P(z_t|a_i)的欧式距离，得到一个N_test×N_training的欧式距离矩阵然后根据设定的近邻数目K，利用KNN分类器将测试集中所有的图像块分成滑坡和非滑坡两类；其中：表示P(z_t|b_k)与P(z_t|a_i)的欧式距离，计算公式为下标ki表示元素在Θ中处于第k行、第i列；K为KNN分类器的近邻数目，取值范围为7～15。