中国地貌基本形态DEM的自动划分研究_刘爱利

第８卷第４期

２００６年１２月地球信息科学ＧＥＯ－ＩＮＦＯＲＭＡＴＩＯＮＳＣＩＥＮＣＥＶｏｌ．８，Ｎｏ．４Ｄｅｃ．，２００６

收稿日期：２００５－１２－１５；修回日期：２００６－０１－１４．

资助项目：国家自然科学基金项目（４０２７１０８９），（４０５７１１２０）．作者简介：刘爱利（１９７９－），女，陕西华阴人，南京师范大学在读博士研究生，主要从事地图学与ＧＩＳ空间分析研究。

Ｅ－ｍａｉｌ：ａｉｌｉｉ６６＠１２６．ｃｏｍ

中国地貌基本形态ＤＥＭ的自动划分研究

刘爱利，汤国安

（南京师范大学虚拟地理环境实验室，南京２１００９７）

摘要：我国１∶１００万的数字高程模型，是在１∶５万及１∶１０万基本地形图上，高精度采集方里网交点高程所构建的１ｋｍ分辨率地面高程数字矩阵。本文利用该ＤＥＭ数据及其所派生的多种地貌信息进行地貌形态类型自动划分的技术方法。实验提取地形起伏度、地表切割度、地表粗糙度、高程变异系数、平均坡度、平均高程６个地形因子，并将各因子置于不同的信息层面中，通过主成分分析，ＩＳＯＤＡＴＡ非监督分类法与Ｂａｙｅｓｉａｎ最大似然监督分类法相结合，对中国地貌的基本形态进行了多维信息综合分类。研究结果表明：①我国１∶

１００万比例尺ＤＥＭ在宏观地貌分类方面具有重要的价值和应用潜力；②所提取的地形因子能宏观地反映我国地形的起伏特征，

为地貌形态分类提供重要的依据；③采用ＩＳＯＤＡＴＡ非监督分类法与Ｂａｙｅｓｉａｎ最大似然监督分类法，

能有效地实现我国地貌基本形态类型的定量化、自动化划分；④依据数据的统计特征进行分类，

较合理地解决了类型模糊的形态实体的归类问题。实验结果不仅揭示了此项技术在地貌形态分类中的巨大潜力，同时对于完善ＤＥＭ数字地形分析的理论与方法也具有重要的意义。

关键词：ＤＥＭ；

中国；地貌；分类１引言

传统地貌分类与制图主要采用地形图量测，野

外实地勘测的方法，制图质量的好坏很大程度上依

赖于专家的个人经验［１，２］。该方法工作量大，

耗时，已越来越难满足现代社会发展的需要［３］。

数字高程模型（ＤｉｇｉｔａｌＥｌｅｖａｔｉｏｎＭｏｄｅｌ，简称

ＤＥＭ

）作为地形表面的一种数字表达形式，与传统地形图相比，其易实现地形因子的自动提取；容易

实现地形分析的定量化、自动化；容易与其他数字

地图或影像进行叠加分析等。然而，目前对ＤＥＭ的

应用还普遍局限在地形指标的提取方面，较少进行

地貌分类的研究［３］。综合分析前人研究成果，以高精

度ＤＥＭ为主要信息源进行地貌分类研究国内较

少；基于专家经验方法较多，而基于统计等数学方

法的较少；严格分类法较多，而模糊分类法较少；微

观地貌研究多，而宏观地貌研究少。

中国１∶１００万ＤＥＭ是在我国１∶５万及１∶１０万

基本比例尺地形图上，高精度采集方里网交点高程

所构建的１ｋｍ分辨率地面高程数字矩阵，其采样精度远高于常规的１∶１００万地形图。本文拟采用基于统计分析的模糊地貌分类法（ＩＳＯＤＡＴＡ非监督分类法与Ｂａｙｅｓｉａｎ最大似然监督分类法相结合），探讨１∶１００万ＤＥＭ定量化、自动化的地貌分类应用。２信息源分析中国１∶１００万比例尺ＤＥＭ是利用全国１∶５万（８７４０幅）、１∶１０万（３８６１幅）地形图，按照２８．１２″（经差）×１８．７５０″（纬差）的格网间隔，采集格网交叉点的高程值，经过查错修改、编辑处理，以１∶１００万图幅范围为单位入库。原始数据的高程允许最大误差为１０￣２０ｍ。为保证数据质量，人工读点两遍，进行对照检查和录入计算机自动对检，最大限度地消除判读误差。利用该模型数据内插国内任一点高程值的中误差（表１），其内插精度符合全国１∶１００万地形图的精度要求，其ＤＥＭ图见图版Ⅲ彩图１。以此数据为基本信息源，再进行重采样

４期刘爱利等：中国地貌基本形态ＤＥＭ的自动划分研究（１０００ｍ×１０００ｍ）

、投影转换（转为等积圆锥投影）、边界修正、值域规范化等预处理，使其能够满足区

域地形分析的需要。该数据的特点在于既能反映我

国地形的宏观特征，又具有较高的采样精度。３地貌类型划分原理与方法３．１地貌类型识别原理地形因子是描述地貌形态特征的量化指标。由于相同地貌类型实体具有最大的相似性，最小的差异性；不同地貌类型实体具有最小的相似性和最大的差异性。因此从相同地貌类型所提取的同种地形因子，其值也更为接近，而从不同类型的地貌实体提取的地形因子则存在较大的差异。这种相似性和差异性即是地貌形态类型划分的依据。假设表示的地面特征有ｎ个地形因子，将每个

地形因子看作单波段，则地面上的点（ｉ，ｊ）

在每个波段上的归一化值可以构成一个矢量Ｘ＝（ｘ１，ｘ２∧，ｘｎ）Ｔ，

Ｘ称作该点的特征值，

包含Ｘ的ｎ维空间称为特征空间，这样地面上的点便可用ｎ维特征空间中的一

系列点来表示。通常情况下，同一类地貌实体其同

种地形因子的值较接近，在特征空间上的对应点就

较为集中于该类中心附近。而多种地貌类型在特征

空间就形成多个点集，地貌类型划分等价于在特征

空间中找到若干条曲线或若干个曲面将特征空间

中的点集分开。因此，地貌类型划分的核心就是确

定判别函数和相应的判别准则。

常用的判别函数很多，本文采用Ｂａｙｅｓｉａｎ最大

似然法进行监督分类。它是求出各点对于各类别的

似然度（ｌｉｋｅｌｉｈｏｏｄ），然后将该点分到似然度最大的

类别中。可表示为公式（１）

Ｐ（ｗｉ／Ｘ）＝Ｐ（ｗｉ）Ｐ（Ｘ／ｗｉ）

ｍｉ＝１"Ｐ（Ｘ／ｗｉ

）Ｐ（ｗｉ）

（１）Ｐ（ｗｉ／Ｘ）为Ｘ属于ｗｉ的概率，

也称后验概率。Ｐ（ｗｉ）为先验概率，也就是在被分类的图像中类别ｗｉ出现的

概率。Ｐ（Ｘ／ｗｉ）为似然度，指当观测到像元数据ｘ时，

它是从分类类别ｋ中得到的（后验）概率。从式（１）中可以看出分母与类别ｗｉ无关，因此作判别时可将Ｐ（Ｘ）去掉。应用最大可能性判别规则，再加上贝氏（Ｂａｙｅｓｉａｎ）的使平均损失最小的原则，都表明ｇｉ（ｘ）＝Ｐ（Ｘ／ｗｉ）Ｐ（ｗｉ）ｉ＝１，２，……，ｍ

（２）是一组理想的判别函数。判别规则为：若Ｐ（ｗｉ）Ｐ（ｘ／

ｗｉ）≥Ｐ（ｗｊ）Ｐ（ｘ／ｗｊ），

则ｘ∈ｗｉ。在地貌形态划分中，其实质是计算地面上的点与各地貌类型的相似程度，相

似程度最大的类，即为该点所属地貌类型。

３．２地形因子选取及地貌分类方法

在一定空间尺度下，参与分类的地形因子能否

准确反映研究区域的地形特征，是地貌形态分类的

关键，其流程见图１。

本研究的目的是对我国宏观地貌形态类型的

自动划分，因此选择宏观地形因子作为地貌形态类

型划分的指标。研究中选取地形起伏度、地表切割

度、高程变异系数、地表粗糙度、平均高程、平均坡

度６个宏观地形因子。各因子的提取均采用ＧＩＳ窗

口分析，窗口大小为１３×１３的矩形

（在ＤＥＭ栅格分辨率为５００ｍ的空间尺度下）。各因子提取结果如表

２、图版Ⅲ彩图２所示。图１地貌类型划分技术流程图Ｆｉｇ．１Ｃｌａｓｓｉｆｉｃａｔｉｏｎｐｒｏｃｅｄｕｒｅ３．３数据预处理海拔高度是地貌类型划分的主要依据。参考“１∶４００万中国及其毗邻地区地貌图”中地貌类型的海拔高度分级，将该层面０￣５００ｍ、５００￣８００ｍ、８００￣２０００ｍ、２０００￣３０００ｍ、３０００￣５５００ｍ、５５００￣８８４８ｍ各高程段分别乘以权重１￣６，使高差增大以便于图像

９・・