地形特征在石漠化程度遙感影像分类中的应用.docVIP

地形特征在石漠化程度遥感影像分类中的应用 蒙琳，胡宝清，段炼 （广西师范学院 资源与环境科学学院，广西 南宁 530001） 摘要：地形、地貌要素既是影响石漠化过程的重要因素，又是石漠化程度的重要表征。本文由DEM数据提取出地形特征，并将地形特征数据与SPOT5遥感影像进行叠加生成复合影像。接着，应用最大似然法对复合影像进行分类并提取出石漠化程度信息。结果表明，遥感影像和地形特征应用于石漠化程度信息的提取有着很强的实用性。 关键词：地形特征；石漠化；遥感 引言 喀斯特石漠化(karst rocky desertification)是指在亚热带脆弱的喀斯特环境背景下，以人类活动为主要驱动力，以降水、径流等运移作用造成的现代水体流失为主导，引起土壤严重侵蚀，基岩大面积出露，土地生产力严重下降，地表出现类似荒漠景观的土地退化过程 REF _Ref207010041 \r \h \* MERGEFORMAT [1]。地形、地貌要素是最基本的自然地理单元,地形的起伏与变化特征直接影响地表物质迁移与能量转换的方式、规模与速率,制约着地表过程的进程、发展方向及地面景观的形成 REF _Ref207890065 \r \h \* MERGEFORMAT [2]。因此，地形、地貌要素既是影响石漠化过程的重要因素，又是石漠化程度的重要表征。 现有的石漠化遥感调查研究中，或是在建立石漠化解译标志基础上进行目视解译和判读 REF _Ref207890065 \r \h \* MERGEFORMAT [3],或是在寻求不同石漠化程度下TM数据各波段像元DN值变化范围特征的基础上进行分类 REF _Ref207890065 \r \h \* MERGEFORMAT [4]，或是设计出单一的表征石漠化程度的指数进行分类 REF _Ref207890065 \r \h \* MERGEFORMAT [5] REF _Ref207890065 \r \h \* MERGEFORMAT [6]。这些方法都存在着不足：有的易受人的主观因素的干扰，有的工作效率不高，有的有一定的片面性；而且它们都忽略了地形特征这一重要因素。本文在1：5万地形图的基础上生成DEM并提取出地形特征，将地形特征信息参与到10米分辨率的SPOT遥感影像的石漠化程度信息的提取。 研究区 研究区处于都安县北部，位于广西壮族自治区中部偏西。地处云贵高原向广西盆地过渡的斜坡地带，都阳山脉东段，东经107°51′～108°30′，北纬23°47′～24°35′之间，南北长100公里，东西宽90公里，幅员面积4095.2平方公里。境内除澄江流域、刁江流域谷地及红水河流域的丘陵地带地势较为平坦外，其余均属大石山区，峰峦丛错，叠嶂连绵，洼地密布。 基于DEM数据提取地形特征 Wood, J.D把地形特征要素定义为凸点(peak)、凹点(pit)、山脊点(ridge)、谷底点(chan-nel)、平地点(plane)和交线点(pass)等六类，然后采用3×3的固定窗口计算DEM的每个像素点的一阶导数和二阶导数,从而提取出以上各个特征要素 REF _Ref207890065 \r \h \* MERGEFORMAT [7]。如表1所示： 表1 wood提取地形特征要素的原理 地形特征名 地形特征的导数表示 描述 Peak凸点 δ2z δx20,δ2zδy20, 位于凸面上的点(各方向上均高于相临各点) Ridge山脊点 δ2z δx20,δ2zδy2=0, 位于凸线上的点 pass交线点 δ2z δx20,δ2zδy20, 位于成正交的凹、凸线的交点处 Plane平地点 δ2z δx2=0,δ2zδy2=0, 位于平面上的点(各方向上均没有凹凸性) Channel谷底点 δ2z δx20,δ2zδy2=0, 位于凹线上的点 Pit凹点 δ2z δx20,δ2zδy20, 位于凹面上的点(各方向上均低于相临各点) 其算法的实施是根据固定窗口中的DEM数据衍生出的slope、 crosc、 maxic 、minic数据组合起来进行地形特征的判别。Envi软件的地形工具的子菜单TOPOGRAPHIC FEATURES提供了实现该提取地形特征的功能。其中可以通过设置topographic kernel size、slope tolerance、curvature tolerance三个参数来调节实现的效果，设置第一个参数可以调节内核窗口的大小，通过设置坡度、曲率容限值避免地形特征都判断为山脊点和谷底点。一般而言，由DEM提取地形特征要素的结果得到山脊点、谷底点、平地点三类的占了绝大部分。跟据这六类地形特征的描述以及提取石漠化程度的最终目的，本文将凹点合并到谷底点，交线点合并到脊线点，平地点与凸点不变。