一种融合DEM、光学遥感和形变信息的潜在滑坡识别方法与流程

本发明属于潜在滑坡识别领域，具体涉及一种融合dem、光学遥感和形变信息的潜在滑坡识别方法。

背景技术：

滑坡是一种常见的工程地质问题导致的自然灾害，现在进行滑坡的早期识别主要是基于常规野外图幅调查或者常规工程勘察，这些方法主要依赖于专业人员的经验，不仅耗时费力，覆盖面有限，还需要投入大量的人力资源。而光学遥感解译技术识别滑坡，尽管其覆盖范围大，且可以通过人机交互进行识别滑坡；但是，常规的光学遥感调查往往只能提供半定量滑坡识别成果，存在漏判、错判现象，自动化程度低，且多数是针对滑坡之后的定位，难以对潜在的滑坡进行早期识别，导致很多灾害难以得到有效的防治。

sar影像不仅含有强度信息，还包含相位信息，insar技术可以获取厘米级甚至毫米级的地表形变，从而提高滑坡识别与监测的可靠性与正确率；但是它只能获取一维形变而且存在阴影、顶底倒置和透视缩短等现象，因此也会存在对滑坡的漏判、错判现象，特别是对于空间尺度较小的滑坡，基于insar形变直接进行识别往往不成功。高精度dem数据可以客观定量的分析精细尺度下的地表特征进行滑坡识别，但也存在漏判、错判现象。insar技术和光学遥感融合的方法中，需要构建复杂的分类规则集，且其是在大范围形变的基础上探测疑似滑坡，然后以单个疑似滑坡为单位进行提取，不利于对大范围的潜在滑坡进行快速识别，仍旧错判、漏判的现象。

技术实现要素：

针对现有技术中存在错判、漏判的缺点，本申请的目的在于，提供了一种融合dem、光学遥感和形变信息的潜在滑坡识别方法。

为了实现上述目的，本发明采取以下技术方案予以实现：

一种融合dem、光学遥感和形变信息的潜在滑坡识别方法，具体包括以下步骤：

步骤1、对研究区域内的sar影像进行处理，得到研究区域的形变速率图；获取研究区域的光学遥感影像数据和dem数据，所述光学遥感影像数据包括全色波段数据和4波段多光谱数据；通过dem数据对光学遥感影像数据进行预处理，得到融合的光学遥感影像图；所述dem指数字高程模型；

步骤2、对研究区域的形变速率图定义投影和重采样，得到重采样后的形变速率图；对所述dem数据定义投影，得到具有投影信息的dem数据；采用多尺度分割方法对具有投影信息的dem数据和融合的光学遥感影像图中的4波段多光谱数据整体进行分割，得到多个对象；

还包括以下步骤：

步骤3、将重采样后的形变速率图与得到的多个对象进行叠加，得到融合的光学遥感影像图中各像元的形变速率；通过融合的光学遥感影像图中各像元的形变速率计算得到多个对象的形变速率；采用阈值分类的方法对多个对象的形变速率进行处理，得到疑似滑坡区；

步骤4、在研究区域内分别选取已知滑坡区和已知非滑坡区，分别沿已知滑坡区和已知非滑坡区的边界，对具有投影信息的dem数据进行剪裁，得到已知滑坡区的dem数据和已知非滑坡区的dem数据；给定多个尺度因子，计算多个尺度因子下已知滑坡区的dem数据的小波系数的方差和多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；通过多个尺度因子下已知非滑坡区的dem数据的小波系数的方差，对多个尺度因子下已知滑坡区dem数据的小波系数的方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差；判断得到的多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差的大小，得到特征尺度因子；采用特征尺度因子对dem数据进行小波变换，得到dem数据中每个节点在多个整数尺度因子下的小波系数的平方和；

步骤5、设定多个给定阈值，判断多个整数尺度因子下dem数据中每个节点的小波系数的平方和与多个给定阈值的大小，得到多个给定阈值对应的疑似滑坡区域；采用误差分析方法对得到的多个给定阈值对应的疑似滑坡区域进行处理，得到优选阈值；从多个给定阈值对应的疑似滑坡区域中选取优选阈值对应的疑似滑坡区域，将优选阈值对应的疑似滑坡区域与疑似滑坡区进行交运算，得到潜在滑坡区域。

进一步的，所述步骤3的具体操作如下：

将重采样后的形变速率图与得到的多个对象进行叠加，得到融合的光学遥感影像图中各像元的形变速率；分别计算多个对象中每个对象的所有像元形变速率的平均值，将得到的平均值分别作为多个对象的形变速率；采用阈值分类的方法对多个对象的形变速率进行分类，得到阈值范围内的形变速率，从多个对象中提取阈值范围内的形变速率所对应的对象，将提取的阈值范围内的形变速率所对应的对象作为疑似滑坡区。

进一步的，所述步骤4的具体操作如下：

步骤4.1、在研究区域内分别选取已知滑坡区和已知非滑坡区，分别沿已知滑坡区和已知非滑坡区的边界，对具有投影信息的dem数据进行剪裁，分别得到已知滑坡区的dem数据和已知非滑坡区的dem数据；给定多个尺度因子，计算得到多个尺度因子下已知滑坡区的dem数据的小波系数和多个尺度因子下已知非滑坡区的dem数据的小波系数；

步骤4.2、计算多个尺度因子下已知滑坡区的dem数据的小波系数的方差，和多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；

步骤4.3、通过多个尺度因子下已知非滑坡区的dem数据的小波系数的方差，对多个尺度因子下已知滑坡区的dem数据的小波系数的方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差；将多个尺度因子下已知滑坡区域dem数据归一化后小波系数方差中的最大值的1/2作为截断值；提取多个尺度因子下已知滑坡区域dem数据归一化小波系数方差中大于截断值的已知滑坡区域dem数据归一化小波系数方差，将大于截断值的已知滑坡区域dem数据归一化小波系数方差所对应的尺度因子作为特征尺度因子，所述特征尺度因子中包括多个整数尺度因子；

步骤4.4、采用特征尺度因子中的多个整数尺度因子对dem数据进行小波变换，得到多个整数尺度因子下的dem数据的小波系数；通过多个整数尺度因子下的dem数据的小波系数，计算dem数据中每个节点在多个整数尺度因子下的小波系数的平方和。

进一步的，所述步骤5的具体操作如下：

步骤5.1、设定多个给定阈值，对dem数据中的每个节点和多个给定阈值分别进行编号，得到dem数据节点序列和给定阈值序列；将dem数据节点序列中的第一个节点作为当前节点，将给定阈值序列中的第一个给定阈值作为当前阈值；

步骤5.2、从dem数据中每个节点在多个整数尺度因子下的小波系数的平方和中选取当前节点在多个整数尺度因子下的小波系数的平方和；

当当前节点在多个整数尺度因子下得到的小波系数的平方和大于等于当前阈值时，将当前节点作为当前阈值对应的疑似滑坡区域节点；将当前节点的下一个节点作为当前节点，重复执行与步骤5.2相同的操作，直至判断完最后一个节点结束，得到当前阈值对应的疑似滑坡区域的多个节点；根据得到的当前阈值的疑似滑坡区域的多个节点构成得到当前阈值对应的疑似滑坡区域；

当当前节点在多个整数尺度因子下得到的小波系数的平方和小于当前阈值时，表示当前节点属于非滑坡区域；将当前节点的下一个节点作为当前节点，重复执行与步骤5.2相同的操作，直至判断完dem数据中所有节点操作完成时操作结束；

步骤5.3、将当前阈值的下一个给定阈值作为当前阈值，重复执行与步骤5.2相同的操作，直至判断完最后一个给定阈值结束，得到多个给定阈值对应的疑似滑坡区域；

步骤5.4、采用误差分析方法对得到的多个给定阈值对应的疑似滑坡区域进行处理，得到优选阈值；从多个给定阈值对应的疑似滑坡区域中选取优选阈值对应的疑似滑坡区域，将优选阈值对应的疑似滑坡区域与疑似滑坡区进行交运算，得到潜在滑坡区域。

进一步的，所述步骤4.1中，采用公式(1)计算得到多个尺度因子下已知滑坡区的dem数据的小波系数，采用公式(2)计算得到多个尺度因子下已知非滑坡区的dem数据的小波系数；

其中，(x,y)表示已知滑坡区域的dem数据的节点，x和y均大于等于0；zfld(x,y)为已知滑坡区域的dem数据中节点为(x,y)的高程，zfld(x,y)为实数；cfld(si,a,b)为zfld(x,y)经过小波变换得到的小波系数，cfld(si,a,b)∈r；si表示第i个尺度因子，si>0；为已知滑坡区的dem数据中节点为(x,y)在第i个尺度因子下的小波函数值，a、b均表示小波变换的位移，a和b均为实数；

其中，(m,n)表示已知非滑坡区的dem数据中的节点，m和n均大于等于0；zunfld(m,n)为已知非滑坡区的dem数据中节点为(m,n)的高程，zunfld(m,n)为实数；cunfld(si,a,b)为zunfld(m,n)经过小波变换得到的小波系数，cunfld(si,a,b)∈r；si表示第i个尺度因子，si>0；为已知非滑坡区的dem数据中节点为(m,n)在第i个尺度因子下的小波函数值，a、b均表示小波变换的位移，a和b均为实数。

进一步的，所述步骤4.2的具体操作如下：

通过多个尺度因子下已知滑坡区的dem数据的小波系数，采用公式(3)计算得到多个尺度因子下已知滑坡区的dem数据的小波系数的方差；

其中，vcwt_fld(si)表示第i个尺度因子时已知滑坡区域的dem数据的小波系数的方差，vcwt_fld(si)≥0；na×nb表示dem数据中的节点个数，na,nb分别表示dem数据中每一行和每一列的节点个数，na,nb均大于0；

通过多个尺度因子下已知非滑坡区的dem数据的小波系数，采用公式(4)计算得到多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；

其中，vcwt_unfld(si)表示第i个尺度因子时已知非滑坡区域的dem数据的小波系数的方差，vcwt_unfld(si)≥0；na,nb表示dem数据中的节点个数，na,nb分别表示dem数据中每一行和每一列的节点个数，na,nb均大于0。

进一步的，所述步骤4.3中，通过多个尺度因子下已知非滑坡区的dem数据的小波系数方差，采用公式(5)对多个尺度因子下已知滑坡区dem数据的小波系数方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差；

vcwt_norm(si)＝vcwt_fld(si)/vcwt_unfld(si)(5)

其中，vcwt_norm(si)表示第i个尺度因子的已知滑坡区的dem数据的归一化小波系数方差，vcwt_fld(si)表示第i个尺度因子的已知滑坡区的dem数据的小波系数的方差，vcwt_unfld(si)表示第i个尺度因子的已知非滑坡区的dem数据的小波系数的方差，vcwt_norm(si),vcwt_fld(si)，vcwt_unfld(si)均大于等于0。

进一步的，所述步骤4.4中，采用特征尺度因子中的多个整数尺度因子对dem数据进行小波变换，具体采用公式(6)计算得到多个整数尺度因子下的dem数据的小波系数；

其中，(xj,yj)表示dem数据中第j个节点，且j>0，xj与yj均为实数；z(xj,yj)为dem数据中第j个节点的高程，z(xj,yj)为实数；c(s'i,a,b)为z(xj,yj)经过小波变换得到的小波系数，c(s'i,a,b)∈r；s'i表示特征尺度因子中第i个整数尺度因子，s'i>0，i>0；为dem数据中第j个节点在第i个整数尺度因子下的小波函数值，且a、b均表示小波变换的位移，且a,b均为实数。

进一步的，所述步骤1的具体操作如下：

步骤1.1、对研究区域内的sar卫星采集得到的sar影像，采用小基线集技术对sar影像进行处理，得到研究区域的的形变速率图；

步骤1.2、通过遥感卫星获取覆盖研究区域的光学遥感影像数据，所述光学遥感影像数据包括分辨率为0.5m的全色波段数据和分辨率为2m的4波段多光谱数据，所述4波段包括红、绿、蓝和近红外波段；

步骤1.3、利用无人机获取研究区域的分辨率为3m的dem数据，利用光学遥感影像数据自带的rpc文件和rpc模型得到完整的rpc模型；利用dem数据和完整的rpc模型分别对光学遥感影像数据进行正射校正，得到正射光学遥感影像，所述rpc指有理多项式系数；

步骤1.4、采用nndiffusepansharpening算法对正射光学遥感影像进行融合，得到融合的光学遥感影像图。

进一步的，所述步骤2的具体操作如下：

步骤2.1、采用地图投影方法对研究区域的形变速率图定义投影，得到具有投影信息的形变速率图，使具有投影信息的形变速率图与融合的光学遥感影像图的投影坐标系统一致；采用三次卷积的方法对具有投影信息的形变速率图进行重采样，得到重采样后的形变速率图，使重采样后的形变速率图与融合的光学遥感影像图的分辨率保持一致；

步骤2.2、采用地图投影方法对所述dem数据定义投影，得到具有投影信息的dem数据，使具有投影信息的dem数据与融合的光学遥感影像图的投影坐标系统保持一致；

步骤2.3、将具有投影信息的dem数据作为一个波段数据，采用多尺度分割方法，对具有投影信息的dem数据和融合的光学遥感影像图中的4波段多光谱数据整体进行分割，得到多个对象，所述对象是指具有同质性的像元集合。

与现有技术相比，本发明的有益效果如下：

1、本发明融合了高精度dem、光学遥感和insar形变信息自动化识别潜在滑坡区域，以形变信息作为一特征因子参与光学遥感对潜在滑坡的识别提取，得到利用面向对象分类技术获取的疑似滑坡区，并利用二维连续小波变化基于高精度dem自动绘制研究区域的疑似滑坡区域，将利用面向对象分类技术获取的疑似滑坡区与优选阈值对应的疑似滑坡区域进行融合，得到潜在滑坡区域，能够有效避免沉降区域，减少错判现象；

2、本发明无需野外调查，对于人员的地质专业知识要求较低，对于人员难以到达的区域也能进行研究；相较于传统的方法，本发明不需要建立规则集，更加简便、客观，自动化程度更高，得到的结果空间上更加连续、完整，有利于进一步的滑坡分析，为及时预防滑坡灾害的发生提供技术支撑。

附图说明

图1为本发明流程图；

图2为基于光学遥感影像及insar形变信息并利用obia技术获取的疑似滑坡的结果图；

图3为利用投影后的dem数据通过小波变换自动绘制的疑似滑坡区域的结果图；

图4为本发明提取的潜在滑坡区域的结果图。

具体实施方式

如图1所示，本发明提供了一种融合dem、光学遥感和形变信息的潜在滑坡识别方法，包括以下步骤：

步骤1、对研究区域内的sar影像进行处理，得到研究区域的形变速率图；获取研究区域的光学遥感影像数据和dem数据，通过dem数据对光学遥感影像数据进行预处理，得到融合的光学遥感影像；所述dem指数字高程模型；所述光学遥感影像数据包括全色波段数据和4波段多光谱数据；

步骤2、对研究区域的形变速率图定义投影和重采样，得到重采样后的形变速率图；对所述dem数据定义投影，得到具有投影信息的dem数据；采用多尺度分割方法对具有投影信息的dem数据和融合的光学遥感影像中的4波段多光谱数据整体进行分割，得到多个对象；

步骤3、将重采样后的形变速率图与得到的多个对象进行叠加，得到融合的光学遥感影像图中各像元的形变速率；通过融合的光学遥感影像图中各像元的形变速率计算得到多个对象的形变速率根据；采用阈值分类的方法对多个对象的形变速率进行处理，得到疑似滑坡区(即利用面向对象分类技术获取的疑似滑坡区)；

步骤4、在研究区域内分别选取已知滑坡区和已知非滑坡区，分别沿已知滑坡区和已知非滑坡区的边界，对具有投影信息的dem数据进行剪裁，得到已知滑坡区的dem数据和已知非滑坡区的dem数据；给定多个尺度因子，计算多个尺度因子下已知滑坡区的dem数据的小波系数的方差和多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；通过多个尺度因子下已知非滑坡区的dem数据的小波系数方差，对多个尺度因子下已知滑坡区dem数据的小波系数方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差；判断得到的多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差的大小，得到特征尺度因子；采用特征尺度因子对dem数据进行小波变换，得到dem数据中每个节点在多个整数尺度因子下的小波系数的平方和；

步骤5、设定多个给定阈值，判断多个整数尺度因子下dem数据中每个节点的小波系数的平方和与多个给定阈值的大小，得到多个给定阈值对应的疑似滑坡区域；采用误差分析方法对得到的多个给定阈值对应的疑似滑坡区域进行处理，得到优选阈值；从多个给定阈值对应的疑似滑坡区域中选取优选阈值对应的疑似滑坡区域，将优选阈值对应的疑似滑坡区域与疑似滑坡区进行交运算，得到潜在滑坡区域。

本发明融合了高精度dem、光学遥感和insar形变信息自动化识别潜在滑坡区域，得到利用面向对象分类技术获取的疑似滑坡区；并利用二维连续小波变化基于高精度dem自动绘制研究区域的疑似滑坡区域，得到优选阈值对应的疑似滑坡区域，将利用面向对象分类技术获取的疑似滑坡区与优选阈值对应的疑似滑坡区域进行融合，得到潜在滑坡区域，能够有效避免沉降区域，减少错判现象；本发明无需野外调查，对于人员的地质专业知识要求较低，对于人员难以到达的区域也能进行研究；相较于传统的方法，不需要建立规则集，更加简便、客观，自动化程度更高，得到的结果空间上更加连续、完整，有利于进一步的滑坡分析，为及时预防滑坡灾害的发生提供技术支撑。

具体的，所述步骤1的具体操作如下：

步骤1.1、对研究区域内的sar卫星采集得到的sar影像，采用sbas(小基线集，smallbaselinesubsets，简称sbas)技术对sar影像进行处理，得到研究区域的的形变速率图。

步骤1.2、通过遥感卫星获取覆盖研究区域的光学遥感影像数据，所述光学遥感影像数据包括分辨率为0.5m的全色波段数据和分辨率为2m的4波段多光谱数据，所述4波段包括红、绿、蓝和近红外波段；

步骤1.3、利用无人机获取研究区域的分辨率为3m的dem数据，利用光学遥感影像数据自带的rpc文件和rpc模型得到完整的rpc模型；利用dem数据和完整的rpc模型分别对光学遥感影像数据中的全色波段数据与4波段多光谱数据进行正射校正，得到包括全色波段数据和4波段多光谱数据的正射光学遥感影像，所述rpc指有理多项式系数；

步骤1.4、采用nndiffusepansharpening算法对包括全色波段数据和4波段多光谱数据的正射光学遥感影像进行融合，得到融合的光学遥感影像图。

该方式通过采集sar影像、研究区域的光学遥感影像数据和dem数据，处理得到研究区域的形变速率图与融合的光学遥感影像图。

具体的，所述步骤2的具体操作为：

步骤2.1、采用地图投影方法对研究区域的形变速率图定义投影，得到具有投影信息的形变速率图，使具有投影信息的形变速率图与融合的光学遥感影像图的投影坐标系统一致；采用三次卷积的方法对具有投影信息的形变速率图进行重采样，得到重采样后的形变速率图，使重采样后的形变速率图与融合的光学遥感影像图的分辨率保持一致；

步骤2.2、采用地图投影方法对所述dem数据定义投影，得到具有投影信息的dem数据，使具有投影信息的dem数据与融合的光学遥感影像图的投影坐标系统保持一致；

步骤2.3、将具有投影信息的dem数据作为一个波段数据，采用多尺度分割方法，对具有投影信息的dem数据和融合的光学遥感影像图中的4波段多光谱数据整体进行分割，得到多个对象，所述对象是指具有同质性的像元集合。

该方式中通过对dem数据及研究区域的形变速率图进行处理，得到多个对象及重采样后的形变速率图；能够采用形变信息这一要素进行分类，提取疑似滑坡区域。

具体的，所述步骤3的具体操作为：

将重采样后的形变速率图与得到的多个对象进行叠加，得到融合的光学遥感影像图中各像元的形变速率；分别计算多个对象中每个对象的所有像元的形变速率的平均值，将得到的平均值作为多个对象的形变速率；设置形变速率的两个阈值m、n，其中m<0,n>0，通过阈值分类的方法提取得到多个对象中对象的形变速率k在k≤m或k≥n的范围内的对象，将提取的对象作为利用obia技术获取的疑似滑坡区；

所述obia表示object-basedimageanalysis，面向对象分类。

该方式中，只利用形变信息这一要素进行分类，设置形变阈值直接提取疑似滑坡区域，不需要计算一系列的特征要素(地形特征、纹理特征、光谱特征等)，既不需要专家经验来决定分类特征要素，也不需要建立十分复杂的规则集对非滑坡地类进行逐一剔除，对工作人员的解译经验及专业知识没有苛刻的要求，操作步骤简洁明了，自动化程度高，十分适合于在大范围内进行疑似滑坡的识别和提取。

具体的，所述步骤4的具体操作为：

步骤4.1、在研究区域内分别选取已知滑坡区和已知非滑坡区，分别沿已知滑坡区和已知非滑坡区的边界，对具有投影信息的dem数据进行剪裁，分别得到已知滑坡区的dem数据和已知非滑坡区的dem数据；给定多个尺度因子，采用公式(1)计算得到多个尺度因子下已知滑坡区的dem数据的小波系数，采用公式(2)计算得到多个尺度因子下已知非滑坡区的dem数据的小波系数；

其中，(x,y)表示已知滑坡区域的dem数据的节点，x和y均大于等于0；zfld(x,y)为已知滑坡区域的dem数据中节点为(x,y)的高程，zfld(x,y)为实数；cfld(si,a,b)为zfld(x,y)经过小波变换得到的小波系数，cfld(si,a,b)∈r；si表示第i个尺度因子，si>0；为已知滑坡区的dem数据中节点为(x，y)在第i个尺度因子下的小波函数值，a、b均表示小波变换的位移，a和b均为实数；

其中，(m,n)表示已知非滑坡区的dem数据中的节点，m和n均大于等于0；zunfld(m,n)为已知非滑坡区的dem数据中节点为(m,n)的高程，zunfld(m,n)为实数；cunfld(si,a,b)为zunfld(m,n)经过小波变换得到的小波系数，cunfld(si,a,b)∈r；si表示第i个尺度因子，si>0；为已知非滑坡区的dem数据中节点为(m,n)在第i个尺度因子下的小波函数值，a、b均表示小波变换的位移，a和b均为实数；

步骤4.2、通过多个尺度因子下已知滑坡区的dem数据的小波系数，采用公式(3)计算得到多个尺度因子下已知滑坡区的dem数据的小波系数的方差；通过多个尺度因子下已知非滑坡区的dem数据的小波系数，采用公式(4)计算得到多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；

其中，vcwt_fld(si)表示第i个尺度因子时已知滑坡区域的dem数据的小波系数的方差，vcwt_fld(si)≥0；na,nb表示dem数据中的节点个数，na,nb分别表示dem数据中每一行和每一列的节点个数，na,nb均大于0；

其中，vcwt_unfld(si)表示第i个尺度因子时已知非滑坡区域的dem数据的小波系数的方差，vcwt_unfld(si)≥0；na,nb表示dem数据中的节点个数，na,nb分别表示dem数据中每一行和每一列的节点个数，na,nb均大于0；

步骤4.3、通过多个尺度因子下已知非滑坡区的dem数据的小波系数方差，采用公式(5)将多个尺度因子下已知滑坡区dem数据的小波系数方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差，及多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差与尺度因子之间的对应关系；将多个尺度因子下已知滑坡区域dem数据归一化后小波系数方差中的最大值的1/2作为截断值，提取多个尺度因子下已知滑坡区域dem数据归一化小波系数方差中大于截断值的已知滑坡区域dem数据归一化小波系数方差，将大于截断值的已知滑坡区域dem数据归一化小波系数方差所对应的尺度因子作为特征尺度因子，所述特征尺度因子中包括多个整数尺度因子；

vcwt_norm(si)＝vcwt_fld(si)/vcwt_unfld(si)(5)

其中，vcwt_norm(si)表示第i个尺度因子的已知滑坡区的dem数据的归一化小波系数方差，vcwt_fld(si)表示第i个尺度因子的已知滑坡区的dem数据的小波系数的方差，vcwt_unfld(si)表示第i个尺度因子的已知非滑坡区的dem数据的小波系数的方差，vcwt_norm(si),vcwt_fld(si)，vcwt_unfld(si)均大于等于0。

步骤4.4、根据特征尺度因子中的多个整数尺度因子，采用公式(6)对dem数据进行小波变换，得到多个整数尺度因子下的dem数据的小波系数；通过多个整数尺度因子下的dem数据的小波系数，计算得到dem数据中每个节点在多个整数尺度因子下的小波系数的平方和；

其中，(xj,yj)表示dem数据中第j个节点，且j>0，xj与yj均为实数；z(xj,yj)为dem数据中第j个节点的高程，z(xj,yj)为实数；c(s'i,a,b)为z(xj,yj)经过小波变换得到的小波系数，c(s'i,a,b)∈r；s'i表示特征尺度因子中第i个整数尺度因子，s'i>0，i>0；为dem数据中第j个节点在第i个整数尺度因子下的小波函数值，且a、b均表示小波变换的位移，且a,b均为实数。

具体的，所述步骤5具体操作如下：

步骤5.1、设定多个给定阈值，对dem数据中的每个节点和多个给定阈值分别进行编号，得到dem数据节点序列和给定阈值序列；将dem数据节点序列中的第一个节点作为当前节点，将给定阈值序列中的第一个给定阈值作为当前阈值；

步骤5.2、从dem数据中每个节点在多个整数尺度因子下的小波系数的平方和中选取当前节点在多个整数尺度因子下的小波系数的平方和；

当当前节点在多个整数尺度因子下得到的小波系数的平方和大于等于当前阈值时，将当前节点作为当前阈值对应的疑似滑坡区域节点；将当前节点的下一个节点作为当前节点，重复执行与步骤5.2相同的操作，直至判断完最后一个节点结束，得到当前阈值对应的疑似滑坡区域的多个节点；根据得到的当前阈值的疑似滑坡区域的多个节点构成得到当前阈值对应的疑似滑坡区域；

当当前节点在多个整数尺度因子下得到的小波系数的平方和小于当前阈值时，表示当前节点属于非滑坡区域；将当前节点的下一个节点作为当前节点，重复执行与步骤5.2相同的操作，直至判断完dem数据中所有节点操作完成时操作结束；

步骤5.3、将当前阈值的下一个给定阈值作为当前阈值，重复执行与步骤5.2相同的操作，直至判断完最后一个给定阈值结束，得到多个给定阈值对应的疑似滑坡区域；执行步骤5.4；

步骤5.4、根据研究区域已知的滑坡编目图，采用误差分析方法对得到的多个给定阈值对应的疑似滑坡区域进行处理，得到优选阈值；从多个给定阈值对应的疑似滑坡区域中选取优选阈值所对应的疑似滑坡区域；将优选阈值所对应的疑似滑坡区域与疑似滑坡区进行交运算，得到潜在滑坡区域。

步骤4与步骤5中，基于高精度的dem数据，利用二维连续小波变换进行疑似滑坡区域的绘制，得到优选阈值对应的疑似滑坡区域，即研究区域的疑似滑坡区域；高精度的dem数据能够对滑坡的精细地表特征进行客观定量的分析；将利用面向对象分类技术获取的疑似滑坡区与优选阈值对应的疑似滑坡区域进行融合，得到潜在滑坡区域，能够有效避免沉降区域，减少错判现象。

实施例

本发明实验数据采用了真实的terrasar-x降轨数据和worldview-02光学遥感影像数据。其中，terrasar-x降轨数据为19景，worldview-02光学遥感影像数据为一景，包括分辨率为0.5m的全色数据与分辨率为2m的多光谱数据，采集日期为2015年11月29日；

实验过程：

步骤1、对研究区域内的sar影像(terrasar-x降轨数据)进行处理，得到研究区域的形变速率图；获取研究区域的光学遥感影像数据和dem数据，所述光学遥感影像数据包括全色波段数据和4波段多光谱数据；根据dem数据对光学遥感影像数据进行预处理，得到融合的光学遥感影像图；所述dem指数字高程模型；

步骤2、对研究区域的形变速率图定义投影和重采样，得到重采样后的形变速率图；对所述dem数据定义投影，得到具有投影信息的dem数据；采用多尺度分割方法对具有投影信息的dem数据和融合的光学遥感影像图中的4波段多光谱数据整体进行分割，得到多个对象；

步骤3、将重采样后的形变速率图与得到的多个对象进行叠加，得到融合的光学遥感影像图中各像元的形变速率；通过融合的光学遥感影像图中各像元的形变速率计算多个对象的形变速率；采用阈值分类的方法对多个对象的形变速率进行处理，得到疑似滑坡区；

步骤4、在研究区域内分别选取已知滑坡区和已知非滑坡区，根据已知滑坡区和已知非滑坡区的边界，对具有投影信息的dem数据进行剪裁，得到已知滑坡区的dem数据和已知非滑坡区的dem数据；给定多个尺度因子，计算多个尺度因子下已知滑坡区的dem数据的小波系数的方差和多个尺度因子下已知非滑坡区的dem数据的小波系数的方差；通过多个尺度因子下已知非滑坡区的dem数据的小波系数的方差，对多个尺度因子下已知滑坡区dem数据的小波系数的方差进行归一化处理，得到多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差；判断得到的多个尺度因子下已知滑坡区的dem数据归一化后的小波系数方差的大小，得到特征尺度因子；采用特征尺度因子对dem数据进行小波变换，得到dem数据中每个节点在多个整数尺度因子下的小波系数的平方和；

步骤5、设定多个给定阈值，判断多个整数尺度因子下dem数据中每个节点的小波系数的平方和与多个给定阈值的大小，得到多个给定阈值对应的疑似滑坡区域；采用误差分析方法对得到的多个给定阈值对应的疑似滑坡区域进行处理，得到优选阈值；从多个给定阈值对应的疑似滑坡区域中选取优选阈值对应的疑似滑坡区域，将优选阈值对应的疑似滑坡区域与疑似滑坡区进行交运算，最终的得到如图4中用黑线绘制的潜在滑坡区域。

利用光学遥感影像及insar形变信息设置形变阈值的方法得到如图2中用黑线绘制的疑似滑坡区域；利用投影后的dem数据通过小波变换自动绘制得到如图3中用黑线绘制的疑似滑坡区域。

通过图2、图4的对比可以看出，insar技术监测得到的形变不仅包含有潜在滑坡的形变，也包含有地面沉降的形变。仅利用光学遥感影像进行分割得到对象后利用形变信息提取得到的疑似滑坡区域会将沉降区包含在内，这是明显的错判现象；通过图3、图4的对比可以看出，基于高精度dem数据利用小波变换进行滑坡区域的自动绘制所获取的疑似滑坡区域，可以有效避免沉降区域，但其覆盖范围广泛，存在大量的错判现象。

从图4中可以看出，本发明融合了高精度dem数据、光学遥感影像和insar形变信息自动化识别潜在滑坡区域，以形变信息作为一特征因子参与光学遥感影像对潜在滑坡的识别提取，并利用二维连续小波变化基于高精度dem数据自动绘制研究区域的疑似滑坡区域，以二者的交集作为最终的潜在滑坡区域，与滑坡编目图中识别出的滑坡具有很高的一致性，并可以在大范围内自动化识别潜在滑坡区域，不需要建立规则集，更加简便、客观，自动化程度更高，得到的结果空间上更加连续、完整，有利于进一步的滑坡分析，为及时预防滑坡灾害的发生提供技术支撑。

以上公开的仅为本发明的具体实施例，但是，本发明实施例并非局限于此，任何本领域的技术人员能思之的变化都应落入本发明的保护范围。