基于小波包变换和支持向量机的滑坡敏感性地图绘制方法与流程

本发明涉及遥感技术领域，具体来说，涉及一种基于小波包变换和支持向量机的滑坡敏感性地图绘制方法。

背景技术：

随着激光雷达技术的发展，激光雷达技术可以实现空间三维坐标的同步、快速、精确地获取，再现客观事物的实时的、真实的形态特性，为快速获取空间信息提供了简单有效手段，可用于生成精确的dem数据。卫星成像技术的发展，高分遥感影像的获取成为了可能。这些重要的数据来源，为之后的进一步数据提取、应用与分析打下基础。

虽然雷达数据和遥感数据为应用提供了更加丰富的空间数据信息，但滑坡的致灾因素复杂，致灾因子的选取较为困难。

技术实现要素：

针对相关技术中的上述技术问题，本发明提出一种基于小波包变换和支持向量机的滑坡敏感性地图绘制方法，能够克服现有技术的上述不足。

为实现上述技术目的，本发明的技术方案是这样实现的：

一种基于小波包变换和支持向量机的滑坡敏感性地图绘制方法，该方法包括以下步骤：

s1：从30m空间分辨率的dem数据和aster传感器的可见光近红外波段数据中计算并提取滑坡致灾敏感因子，归一化处理滑坡致灾敏感因子，其中，dem数据表示数字高程数据，滑坡致灾敏感因子包括高程、坡度、坡向、坡长、水流力指数、地形湿度指数、归一化植被指数；

s2：采用小波包变换，分解每个滑坡敏感因子，其中，小波包变换的函数为

x(t)为原始信号，是时间或位置的函数，此处是位置的函数，s为比例参数，μ为平移参数，*代表复共轭，母小波ψ(t)为变换函数。

s3：采用多项式核函数的支持向量机，有监督地学习滑坡敏感因子特征，累次向分类器输入敏感因子，得到分类结果，输出滑坡敏感因子集；

s4：依据得到的滑坡敏感因子集，生成滑坡敏感性地图。

进一步的，所述步骤s1包括以下步骤：

s11：采用地理信息系统软件arcgis和遥感图像处理平台envi，从30m空间分辨率的dem数据和aster传感器的可见光近红外波段数据中计算并提取7个滑坡致灾敏感因子；

s12：归一化处理7个滑坡致灾敏感因子，其中，公式为

v＇是归一化的数据矩阵，v是原始数据矩阵，而u和l是上下归一化边界。

进一步的，所述步骤s1包括以下步骤：

s31：建立多项式核函数分类器模型；

s32：输入模型中i+1个敏感因子的组合，其中，最初模型中敏感因子的组合数量为0；

s32：计算n-i个svm模型分类精度值，依次向模型中加入新的敏感因子，得到滑坡敏感因子集，其中，svm为支持向量机。

进一步的，所述步骤s2中，精细地处理滑坡敏感因子局部特征。

进一步的，离散小波包变换的比例尺和位置基于两个幂，函数表达式为

其中，两个幂为二进位比例尺、位置，j和k为整数，so表示固定的扩张步骤且大于零。

进一步的，所述步骤s2中，二维小波分解描述为沿水平轴和垂直轴连续应用的一维离散小波变换，其中，水平轴表示地图或图像的行，垂直轴表示地图或图像的列。

进一步的，所述步骤s32中，依次向模型中加入新的敏感因子，直至得到的分类精度不再增加为止。

本发明的有益效果：通过计算滑坡敏感因子，得到能够表征滑坡敏感性的特征；通过小波包变换，精细地处理敏感因子的局部特征描述；通过支持向量机对滑坡敏感因子特征进行有监督学习，能够能快速进行样本训练和测试，在数据量小、样本少的条件下获得较高的分类精度。

附图说明

为了更清楚地说明本发明实施例或现有技术中的技术方案，下面将对实施例中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本发明的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1是根据本发明实施例所述的基于小波包变换和支持向量机的滑坡敏感性地图绘制方法的流程框图；

图2是根据本发明实施例所述的svm分类器模型处理多因子高维特征分类的流程示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员所获得的所有其他实施例，都属于本发明保护的范围。

如图1-2所示，根据本发明实施例所述的基于小波包变换和支持向量机的滑坡敏感性地图绘制方法，包括以下步骤：

s1：从30m空间分辨率的dem数据和aster传感器的可见光近红外波段数据中计算并提取滑坡致灾敏感因子，归一化处理滑坡致灾敏感因子，其中，dem数据表示数字高程数据，滑坡致灾敏感因子包括高程、坡度、坡向、坡长、水流力指数、地形湿度指数、归一化植被指数；

s2：采用小波包变换，分解每个滑坡敏感因子，其中，小波包变换的函数为

x(t)为原始信号，是时间或位置的函数，此处是位置的函数，s为比例参数，μ为平移参数，*代表复共轭，母小波ψ(t)为变换函数。

s3：采用多项式核函数的支持向量机，有监督地学习滑坡敏感因子特征，累次向分类器输入敏感因子，得到分类结果，输出滑坡敏感因子集；

s4：依据得到的滑坡敏感因子集，生成滑坡敏感性地图。

步骤s1包括以下步骤：

s11：采用地理信息系统软件arcgis和遥感图像处理平台envi，从30m空间分辨率的dem数据和aster传感器的可见光近红外波段数据中计算并提取7个滑坡致灾敏感因子；

s12：归一化处理7个滑坡致灾敏感因子，其中，公式为

v＇是归一化的数据矩阵，v是原始数据矩阵，而u和l是上下归一化边界。

步骤s1包括以下步骤：

s31：建立多项式核函数分类器模型；

s32：输入模型中i+1个敏感因子的组合，其中，最初模型中敏感因子的组合数量为0；

s32：计算n-i个svm模型分类精度值，依次向模型中加入新的敏感因子，得到滑坡敏感因子集，其中，svm为支持向量机。

在本发明的一个具体实施例中，所述步骤s2中，精细地处理滑坡敏感因子局部特征。

在本发明的一个具体实施例中，离散小波包变换的比例尺和位置基于两个幂，函数表达式为

其中，两个幂为二进位比例尺、位置，j和k为整数，so表示固定的扩张步骤且大于零。

在本发明的一个具体实施例中，离散小波包变换执行两组功能，其中，功能为高通和低通滤波器。

在本发明的一个具体实施例中，所述步骤s2中，二维小波分解描述为沿水平轴和垂直轴连续应用的一维离散小波变换，其中，水平轴表示地图或图像的行，垂直轴表示地图或图像的列。

在本发明的一个具体实施例中，在小样本条件下，敏感因子高维特征问题通过支持向量机能有效解决，得到较好的分类效果，筛选出敏感因子集。

在本发明的一个具体实施例中，所述步骤s32中，依次向模型中加入新的敏感因子，直至得到的分类精度不再增加为止。

综上所述，借助于本发明的上述技术方案，通过计算滑坡敏感因子，得到能够表征滑坡敏感性的特征；通过小波包变换，精细地处理敏感因子的局部特征描述；通过支持向量机对滑坡敏感因子特征进行有监督学习，能够能快速进行样本训练和测试，在数据量小、样本少的条件下获得较高的分类精度。

以上所述仅为本发明的较佳实施例而已，并不用以限制本发明，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。