一种基于序列光谱曲面特征的滑坡识别方法
【技术领域】
[0001] 本发明涉及一种识别滑坡地质灾害的方法，具体的说，是涉及一种基于序列光谱 曲面特征的滑坡识别方法。
【背景技术】
[0002] 当前遥感技术监测与识别滑坡地质灾害的方法较多，但多数是从监测与识别其它 目标对象的方法中总借鉴而来的。现有技术中，识别滑坡地质灾害的方法主要有光谱距离 方法、光谱角度方法、决策树方法等;现有方法实现原理要么是基于单个像元的光谱信息与 选取方法，要么是必须有较多其它数据参与。单像元方法极易导致椒盐现象，产生无效破碎 斑块，有其它数据参与的方法依赖性较强，很难切实应用与推广。
【发明内容】
[0003] 本发明的目的在于克服上述缺陷，提供一种方法简单、实现方便、科学的基于序列 光谱曲面特征的滑坡识别方法。
[0004] 为了实现上述目的，本发明采用的技术方案如下：
[0005] -种基于序列光谱曲面特征的滑坡识别方法，包括以下步骤：
[0006] (1)将高光谱遥感图像计算成反射率图像；
[0007] (2)将反射率图像分割成若干单元，单元即为目标单元，选定部分目标单元为样本 单元；
[0008] (3)基于单元构建与单元大小一样的光谱曲面，在每个波段图像中每个单元构建 一个光谱曲面，将所有波段中，基于同一个单元构建的光谱曲面组合成曲面集合；目标单元 对应的为目标曲面集合，样本单元对应的为样本曲面集合；
[0009] (4)通过对比目标曲面集合与样本曲面集合的相似程度来确定滑坡灾害的可能 性；
[0010] 步骤(3)中，一个波段包含m个光谱曲面，m为分割成的单元数目；一个单元则对应η 个光谱曲面，其中，η为高光谱遥感波段数目。
[0011] 进一步的，所述步骤(2)中采用普适方法将反射率图像分割成若干单元，且分割的 单元尺寸与滑坡尺度匹配。
[0012] 进一步的，所述步骤(3)中的光谱曲面为以单元空间轮廓为边界，以对应的波段反 射率值为顶点的不规则曲面。
[0013] 进一步的，所述步骤(4)的具体方法如下：
[0014] (41)通过计算光谱曲面的特征参数，得到每个光谱曲面的特征值；
[0015] (42)计算样本单元的特征值与目标单元的特征值之间差异，得到差异指数；
[0016] (43)结合每个波段形成的差异指数，通过加权方法计算出样本单元与目标单元的 综合差异值；
[0017] (44)将最终的综合差异值与预设的滑坡识别阀值进行比较，如果差异值小于或者 等于阀值，就认为本单元为滑坡分布，否则，本单元为非滑坡分布。
[0018] 进一步的，所述步骤(43)中的波段为预设置的滑坡灾害的敏感光谱波段。
[0019] 进一步的，所述步骤(41)中光谱曲面的特征参数包括：曲面面积和投影面积的比 例、曲面凹凸点的波动区间、曲面凹凸点平均值和方差。
[0020] 与现有技术相比，本发明具有以下有益效果：
[0021] (1)本发明将滑坡灾害看成一个整体目标单元，而不是一个(系列)独立像素，综合 考虑了目标单元包含的像素之间的逻辑相关性，客观上符合对滑坡灾害现象的认识。
[0022] (2)本发明充分利用了光谱序列优势，通过计算每个光谱曲面细微差异特征，捕捉 滑坡灾害在遥感影像中的细微响应，通过序列特征计算综合特征，更具科学性。
[0023] (3)本发明在特征综合时，考虑到了滑坡体物质组成对遥感数据的影响贡献，并通 过加权的方法将其贡献进行量化，有效提升了方法可靠性。
【附图说明】
[0024] 图1为区域反射率图像。
[0025]图2为滑坡识别的效果图。
【具体实施方式】
[0026]下面结合实施例对本发明作进一步说明，本发明的实施方式包括但不限于下列实 施例。
[0027]实施例
[0028] 如图1、2所示，本实施例提供了一种基于序列光谱曲面特征的滑坡识别方法，该识 别方法结合高光谱序列光谱曲面特征，充分利用了光谱序列优势，通过计算每个光谱曲面 细微差异特征，捕捉滑坡灾害在遥感影像中的细微响应，通过序列特征计算综合特征，更具 科学性。本识别方法设计原理是:将图像分割成合适的目标单元，然后利用每个目标单元区 域内的一个波段的反射率值构建一个与目标单元大小一样的光谱曲面，一个波段包含m个 光谱曲面，m为分割成的单元数目；一个单元则对应η个光谱曲面，其中，η为高光谱遥感波段 数目；所有波段中，基于同一个单元构建的光谱曲面组合成曲面集合；目标单元对应的为目 标曲面集合，样本单元对应的为样本曲面集合。通过对比该目标单元的曲面集合与样本单 元的曲面集合的相似程度来确定滑坡灾害的可能性。具体的方法如下：
[0029] 首先，将高光谱遥感图像计算成反射率图像;对图像进行面向滑坡灾害目标的多 尺度分割(方法较多，本实施例采用普适方法），目的是利用光谱信息和空间信息划分出合 适单元，选择与滑坡尺度相适应的单元作为待识别的基本单元。
[0030] 然后利用每个目标单元区域内的一个波段的反射率值构建一个与目标单元大小 一样的光谱曲面，一个波段对应一个光谱曲面，一个单元则对应η(高光谱遥感波段数目）个 光谱曲面，形成曲面集合，为了便于区分，将其命名为目标单元的曲面集合。选取部分目标 单元(已确认为滑坡的单元）的曲面集合为样本单元的曲面集合。其中，光谱曲面为以单元 空间轮廓为边界，以对应的波段反射率值为顶点的不规则曲面。
[0031] 最后，通过对比该目标单元的曲面集合与样本单元的曲面集合的相似程度来确定 滑坡灾害的可能性。在本实施例中曲面的相似程度可以由曲面特征进行描述，比如曲面数 据集合值的方差、均值、起伏度等。最后，通过预设置的滑坡灾害的敏感光谱波段(敏感光谱 曲面），将每个光谱曲面的结果进行加权计算，计算结果为综合差异值，对比该差异值与预 设阀值进行滑坡识别。
[0032] 本实施例中，目标曲面集合与样本曲面集合相似程度的确定方法具体如下:首先， 通过计算曲面的特征参数，得到每个曲面的特征值;特征参数如下：
[0033]
[0035]为了实现方便，优先考虑以下参数:平均值：；
平
均凸点个数：平均凹点个数：
；综合加权： 还有特征的加权 值。
[0036]根据得到的特征值，计算样本单元的特征值与目标单元的特征值之间差异，得到 差异指数;然后，结合每个波段形成的差异指数，通过加权方法计算出样本单元与目标单元 的综合差异值;其中，加权值是基于滑坡体的物质组成确定，主要依据是滑坡体物质在高光 谱遥感波段范围内的敏感性，比如地质体敏感波段贡献大、水汽空气影响波段贡献小、近红 外波段贡献大(岩石诊断光谱分布带)等。
[0037]最后，将每个光谱曲面的结果进行加权计算，计算结果与预设阀值进行比较来识 别滑坡。
[0038]按照上述实施例，便可很好地实现本发明。值得说明的是，基于上述设计原理的前 提下，为解决同样的技术问题，即使在本发明所公开的结构基础上做出的一些无实质性的 改动或润色，所采用的技术方案的实质仍然与本发明一样，故其也应当在本发明的保护范 围内。
【主权项】
1. 一种基于序列光谱曲面特征的滑坡识别方法，其特征在于，包括以下步骤： (1) 将高光谱遥感图像计算成反射率图像； (2) 将反射率图像分割成若干单元，单元即为目标单元，选定部分目标单元为样本单 元； (3) 基于单元构建与单元大小一样的光谱曲面，在每个波段图像中每个单元构建一个 光谱曲面，将所有波段中，基于同一个单元构建的光谱曲面组合成曲面集合；目标单元对应 的为目标曲面集合，样本单元对应的为样本曲面集合； (4) 通过对比目标曲面集合与样本曲面集合的相似程度来确定滑坡灾害的可能性； 步骤(3)中，一个波段包含m个光谱曲面，m为分割成的单元数目；一个单元则对应η个光 谱曲面，其中，η为高光谱遥感波段数目。2. 根据权利要求1所述的一种基于序列光谱曲面特征的滑坡识别方法，其特征在于， 所述步骤(2)中采用普适方法将反射率图像分割成若干单元，且分割的单元尺寸与滑坡尺 度匹配。3. 根据权利要求1所述的一种基于序列光谱曲面特征的滑坡识别方法，其特征在于， 所述步骤(3)中的光谱曲面为以单元空间轮廓为边界，以对应的波段反射率值为顶点的不 规则曲面。4. 根据权利要求1所述的一种基于序列光谱曲面特征的滑坡识别方法，其特征在于， 所述步骤(4)的具体方法如下： (41) 通过计算光谱曲面的特征参数，得到每个光谱曲面的特征值； (42) 计算样本单元的特征值与目标单元的特征值之间差异，得到差异指数； (43) 结合每个波段形成的差异指数，通过加权方法计算出样本单元与目标单元的综合 差异值； (44) 将最终的综合差异值与预设的滑坡识别阀值进行比较，如果差异值小于或者等于 阀值，就认为本单元为滑坡分布，否则，本单元为非滑坡分布。5. 根据权利要求4所述的一种基于序列光谱曲面特征的滑坡识别方法，其特征在于， 所述步骤(43)中的波段为预设置的滑坡灾害的敏感光谱波段。6. 根据权利要求4所述的一种基于序列光谱曲面特征的滑坡识别方法，其特征在于， 所述步骤(41)中光谱曲面的特征参数包括：曲面面积和投影面积的比例、曲面凹凸点的波 动区间、曲面凹凸点平均值和方差。
【专利摘要】本发明公开了一种基于序列光谱曲面特征的滑坡识别方法，解决了现有滑坡识别方法实现复杂、识别精度欠佳等缺陷。该识别方法包括以下步骤：（1）将高光谱遥感图像计算成反射率图像；（2）将反射率图像分割成若干目标单元，选定部分目标单元为样本单元；（3）基于单元构建与单元大小一样的光谱曲面，在每个波段图像中每个单元构建一个光谱曲面，将所有波段中，基于同一个单元构建的光谱曲面组合成曲面集合；（4）通过对比目标曲面集合与样本曲面集合的相似程度来确定滑坡灾害的可能性。本发明将滑坡灾害看成一个整体目标单元，而不是一个独立像素，综合考虑了目标单元包含的像素之间的逻辑相关性，客观上符合对滑坡灾害现象的认识。
【IPC分类】G06K9/34, G06K9/00
【公开号】CN105469058
【申请号】CN201510855658
【发明人】叶成名, 赵磊, 李尧, 苗放, 邓地娟, 梁莉, 胡源, 徐松浦
【申请人】成都理工大学
【公开日】2016年4月6日
【申请日】2015年11月30日