基于DEM与DLG的福建省地貌形态自动分类

第１５卷第１期２０１３年２月

地球信息科学学报

ＪＯＵＲＮＡＬＯＦＧＥＯＩＮＦＯＲＭＡＴＩＯＮＳＣＩＥＮＣＥ　　－　

Ｖｏｌ．１５，Ｎｏ．１

，２０１３Ｆｅｂ．

基于ＤＥＭ与ＤＬＧ的福建省地貌形态自动分类

陈加兵，李　慧，陈文惠，郑达贤

（）福建师范大学地理研究所，福州３５０００７

，，摘要：以福建省１数字线划图）和Ｄ数字高程模型）∶１０万的ＤＬＧ（ＤｉｉｔａｌＬｉｎｅＧｒａｈｉｃＥＭ（ＤｉｉｔａｌＥｌｅｖａｔｉｏｎＭｏｄｅｌ　　　　ｇｐｇ为基础，在Ａ依据坡度和面积组合指标，将福建省地貌形态单元分为平缓地和山丘地；提取地形ｒｃＧＩＳ软件支持下，构建了福建省地方相对基准面，并以流域水文分析方法提取了山丘地地貌形态实图中水系与等高线的高程点交集，

，体单元；把平缓地和山丘地貌形态实体单元分别与地方侵蚀基准面叠加，参照“中国１将福∶１００万地貌制图规范”台地、低丘、高丘、低山、中山。研究结果表明：该方法符合福建传统地貌分类体系，能建省地貌形态类型分为平地、

够较好快速实现平缓地与山丘地、平地（平原）与台（阶）地的自动划分，有效地提取山丘地貌形态实体单元的界线。关键词：地貌形态类型；福建省ＤＬＧ；ＤＥＭ；：／ＤＯＩ１０．３７２４ＳＰ．Ｊ．１０４７．２０１３．０００７５

１　引言

地貌是陆域表面各自然地理要素相互作用相互影响最为活跃的界面之一。在很大程度上控制着其他自然地理要素的发生、发展与变化，在一定区域尺度内，地貌与其他自然地理要素，如气候、水

１］

。文、植被、土壤等相互关系中，其起着重要作用［

形统计单元进行计算以减少在计算小起伏地形时的误差。此后Ｄｒａｕｔ在上述基础上进一步修正与ｇ发展，从Ｄ坡度分级、剖面曲率、平ＥＭ中提取海拔、面曲率等地形因子，采用图像分割方法来完成分

９］

。此后，类［国内外许多专家学者在规则统计单元］１０２０－

，的地貌自动分类方面进行了深入研究［这些研

地貌的长期、系统和深入研究，积累了大量的珍贵资料，在农业生产、生态环境建设与保护和国防建设等方面都发挥着巨大的作用量的人力、物力、财力和时间。

当前，随着计算机、遥感和ＧＩＳ等技术的飞速

５］

，发展，地貌制图工作可在室内完成［加之野外调能在较短的时间内推陈出新。因此，人们查辅助，

开始尝试研究各种地貌类型数字制图的新思路和方法。

１９６４年Ｈａｍｍｏｎｄ最早提出地貌类型的自动

划分方法，指出可通过统计区域单元内的坡度、相

５］

，对起伏度进行判定［这种思路由Ｄｉｋａｕ等于１９９１６］

，年实现［之后这种方法得到Ｂｒａｂｎ、Ｍｏｒａｎ等的ｙｇ］７８－

。Ｄ不断尝试、修正和发展［ｉｋａｕ等利用正方形格

［］１４－

究的共同点在于利用规则单元统计区域内各类地形因子（如地形起伏度、地表粗糙度、高程变异系坡度变化率、平均高程、剖面曲率、平面曲率数、

，等）实现了地貌类型的快速分类。其缺点是采用提取的地貌类型界线往往与传统规则的统计单元，

人工勾绘的地貌界线相差较大，在地形起伏较大的山区尤为明显。传统地貌制图大多强调地貌实体

２１］

为此，肖飞等［提出以Ｄ单元的完整性，ＥＭ进行

。但该方法

工作量大，进行大范围地貌信息提取，需要花费大

完整地貌实体的尝试划分研究：首先，提取平原、山然后，针对山地地貌类型按照约定的规则地２类，

进一步划分。但上述分类标准与方法主要是针对整个中国大陆区域提出的，故其并不完全适用于福建传统的地貌分类体系。为此，本研究以福建省为例，提出一种以Ｄ数字线划图）为基ＥＭ与ＤＬＧ（结合约定的地貌形态划分标准，进行以地理实础，

体为单元的地貌形态类型自动划分。

网作为统计单元进行划分，而Ｂｒａｂｎ则利用了圆ｙ

；修回日期：收稿日期：２０１２０９２０２０１２１２１０．－－－－

；；基金项目：国家自然科学基金项目（福建省教育厅Ａ类项目（福建省自然科学基金项目４１００１０６４）ＪＡ１００８４）

（）。２０１２Ｊ０５０５１

，：作者简介：陈加兵（男，博士，副研究员，研究方向为资源与环境信息系统应用。Ｅ－ｍ１９７７ａｉｌｃｂｆｎｕ６３．ｃｏｍ－）＠１ｊｊ

７６

地球信息科学学报　　　　　　　　　　　　　　　　　　　２０１３年

［２］

均较大。经１中获得的∶１０万ＤＬＧ生成的ＤＥＭ２

２　研究区地理背景

福建省位于中国东南部，东海之滨。陆域介于东西最大宽１１５°５０′２０°４３′Ｅ，２３°３２′８°１９′Ｎ，～１～２

南北最大长度约５度约５４０ｋｍ，５０ｋｍ。西南与广东省接壤、西邻江西省，北靠浙江省，东南隔台湾海峡

２

。全省陆域总面积１与台湾省相望，２．１４万ｋｍ

高程和坡度数据统计分析：福建省大部分县（市、区）海拔最高点与最低点的相对高差在１０００ｍ以陆域坡度≤２的平地和微坡面仅占１坡上；°３．１％，

度≤１区域面积占４而＞１的陡坡、峻坡５°１．１％，５°其中，的峻坡和峭坡和峭坡的面积占５８．９％，５°≥２面积占２４．３％。

福建地势西北高、东南低，与海岸大体平行的两列大山带呈ＮＮ西部一列ＥＳＳＷ走向斜贯全境，－为武夷山脉，东部一列为鹫峰山－戴云山－博平岭山。在两列山脉之间，脉（闽中大山带）是一条由纵横交错的丘陵山地和大大小小的盆谷地镶嵌组成的地势逐渐山地盆谷带。闽中大山带向东直至海岸，下降，构成福建东部由低山、丘陵、台地、平原和半主要也是丘陵台地）与海湾交错组成的岛岛屿（

闽中、闽西地区ＮＥＳＷ走向的沿海地带。闽东北、－

以中山、低山和丘陵为主，分布众多串珠状河谷盆东南沿海地区地势低缓，以平原、台地和丘陵为地；

主；海岸线漫长曲折，拥有众多港湾和半岛，岛屿星）

。图１罗棋布（

３　福建省地貌形态的分析及划分

３．１　地貌形态分类标准及定义３．１．１　基本地貌面单元

在进行福建省地貌形态分类前，需对地貌相关称地概念进行定义。属于以下三类地貌面之一者，貌形态分析的基本地貌面单元：

水域：包括河（溪）流、湖泊、水库、水塘等。平缓面：大于一定面积且坡度≤２的地面及其°相邻的坡度≤５的区域。°

锥形面：具有一个山峰点（高地点）且为一条封闭的汇水线包围的一个区域。

单斜面：不属于上述二类，而位于上述二类划或与水域之间的区域，包括峰点位出的区域之间、于边界线上的区域。３．１．２　地方相对基准面

区内所有溪河水系线（水系与等高线的交点的）高程所构成的一个虚拟地形面，用于计算地貌单元的相对高度。福建现代地貌发育的外营力主目前，已达到一个相对稳定阶要以流水侵蚀为主，

段，因而本地方相对基准面近似于侵蚀基准面。３．１．３　地貌形态分类及其标准

地貌形态分类可以是一个基本地貌面单元，或是由几个相邻的基本地貌面单元合并组成的分类评价单元，也可由一个基本地貌面单元根据一定原则分割而成的新的分类评价单元。本文参照“中国

［２３］

，结合福建省地貌类型１∶１００万地貌制图规范”

特点，按其形态分为：平地、台地、低丘、高丘、低山、

图１　福建省地势图

Ｆｉ．１　ＴｈｅｈｓｏｍｅｔｒｉｃｍａｏｆＦｕｉａｎＰｒｏｖｉｎｃｅ　　　　ｇｙｐｐｊ　

中山６个类型。具体划分标准及定义如下：（）平缓地：坡度≤２的区域和其相邻的坡度≤１°

周边的台阶地前坡面所构成的区域。５°

平地（河谷平地、山间平地，包括堆积①平原、

：，平地和剥蚀平地）平缓地中坡度≤２且与地方侵°蚀基准面的高差≤５ｍ的区域。

阶）地：平缓地中地面与地方侵蚀基准面②台（

福建山高坡陡，最高峰为武夷山脉的黄岗山，海拔２地形破碎，地貌形态类型多样，并以山１５８ｍ；地丘陵为主，占全省土地总面积的８５％以上。由于中低山构成全省的主体，因而各地的地面相对高差

１期陈加兵等：基于ＤＥＭ与ＤＬＧ的福建省地貌形态自动分类

７７　　　

高差５～５以及与其相邻、海拔低于该区０ｍ的地面，域的单斜面合并组成的区域。

（）山丘地：平缓地和水域以外的其他区域。２

相对高度＜２００ｍ的山丘地。①丘陵：

高丘：相对高度１００～２００ｍ。

低丘：相对高度＜１００ｍ。

相对高度≥２００ｍ的山丘地。②山地：低山：相对高度≥２最高峰海拔＜１００ｍ，０００ｍ的山丘地。

中山：相对高度≥２最高峰海拔≥１００ｍ，０００ｍ的山丘地。

３．２　基于ＤＥＭ与ＤＬＧ的地貌形态分类实验３．２．１　基础数据及软件基本数据：地图投影和１∶１０万的福建省ＤＬＧ（北京１坐标系统采用Ｇａｕｓｓｒｕｅｒ投影，９５４坐标－Ｋｇ系，单位为ｍ）和由此生成的栅格ＤＥＭ。

所采用的软件：美国ＥＳＲＩ公司开发的ＡｒｃＧＩＳ）桌面版（及工作版（ＡｒｃＧＩＳＡｒｃＩｎｆｏＡｒｃＩｎｆｏＷｏｒｋ　　－）。ｓｔａｔｉｏｎ

福建省栅格Ｄ采用三角化不规则网ＥＭ建立：（方法，将矢量线划图转换成ＴＴＩＮ）ＩＮ数据结构，进行空间插值；对ＴＩＮ数据结构转换成格网结构；，并格网数据进行采样形成栅格ＤＥＭ（ＧＲＩＤ格式）用福建省边界进行裁切，生成福建３０ｍ×３０ｍ的

［２］

。ＤＥＭ２

图２　福建省地方相对基准面

Ｆｉ．２　ＴｈｅｒｅｌａｔｉｖｅｒｅｆｅｒｅｎｃｅｓｕｒｆａｃｅｏｆＦｕｉａｎＰｒｏｖｉｎｃｅ　　　　　　ｇｊ

缓面基本地貌单元和面状水域后的其余区域的地形倒置，形成负地形。可将其Ｄ就可ＥＭ乘以－１，把山坡把所有山峰变成一个个锥形漏斗的最低点，间的汇水线变成锥形漏斗的边交界线。②从每一在Ａ个锥底点出发，ｒｃＧＩＳ流域水文分析模块支持

下，按划分流域的方法划出流域区，即提取山丘地

［２］的沟谷线。本文采用Ｄ提取山脊线的方８算法２法进行流域划分，即通过填洼、消峰、流向分析、汇

３．２．２　实验与结果分析

（））坡度图生成：利用Ｓ函数命令从１ＬＯＰＥ（ＤＥＭ提取坡度图。（）地方相对基准面构建：从１利２∶１０万ＤＬＧ，用Ａ共获得ｒｃＧＩＳ软件寻找水系与等高线的交点，交点３采用三角化不规１８０８２个。利用这些交点，　则网（方法进行空间插值，形成地方相对基准ＴＩＮ））。面（图２

（）提取面状水域：从１３∶１０万ＤＬＧ中提取面状水域，并栅格（化。３０ｍ×３０ｍ）

（）提取平缓地基本地貌单元：４ＣＧＩＳ①在ＡＲ支持下，从坡度图中去除水面水域，后提取坡度≤

２

，（面积≥０．最小制图单元）的区域，这种平２°１ｋｍ

缓地只能出现在河谷、山间盆地和平原上，不会在

流分析和流域识别等流程，形成流域单元区（面积

２

，在０．以上）也是每一个漏斗的边缘交界线、１ｋｍ

边缘线所包围的区域，即为山丘地锥形面基本地貌单元。

（）提取单斜面基本地貌单元：提取过平缓面、６

面状水系和锥形面之后的其余区域为单斜面基本地貌单元。

（）平原、平地、台阶地的划分和归并：７①在把平缓面与地方侵蚀基准面作差ＡｒｃＧＩＳ支持下，

值运算，确定平缓面的相对高度，将相对高度≤５ｍ的区域提取出来，即为平原或平地；②把平缓面的其余区域及其相邻的海拔比它低的单斜面归并一起，即为台阶地。

（）山地和丘陵地的判别和划分：利用Ｚ８ＯＮ－

））函数和Ｚ函数分别计ＡＬＲＡＮＧＥ（ＯＮＡＬＭＡＸ（

算各锥形面的相对高度和山峰点海拔高度。判别方法为：相对高度≤１００ｍ为低丘；１００～２００ｍ为高

山丘地上。②把上述所提取的区域与其周边坡度的区域合并，形成平缓地基本形态地貌单元。°≤５

）（提取锥形面基本地貌单元：５①将提取过

平

７８

地球信息科学学报　　　　　　　　　　　　　　　　　　　２０１３年

表１　福建省地貌各形态类型面积和所占比例ｒｏｏｒｔｉｏｎｅｏｍｏｒｈｏｌｏｉｃａｌＴａｂ．１　Ｔｈｅａｒｅａａｎｄｏｆ　　　　　ｐｐｇｐｇ

ｔｅｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅ　　　ｙｐｊ

地貌类型

水域

平原、平地

台地

丘陵低丘

高丘

山地低山

中山

丘；相对高度＞２其中，海拔≤１００ｍ为山地，０００ｍ为低山，０００ｍ为中山。＞１

（））生成数字地貌图：通过ＭＥ函数或９ＲＧＥ（

把具有相同投影和坐标系统ＧＲＩＤＩＮＳＥＲＴ命令，

并转成矢量图，且消除细的各类型地貌单元合并，

小碎多边形，可形成福建省以完整地貌实体为单元。经统计分析，的地貌形态类型图（图３）福建省地貌各形态类型面积和所占比例如表１。从表１可以看出，福建省地貌形态类型中丘陵面积占３４．１％，山地面积占５两都合计占８平原、平１．８％，５．９％，地和台地合计只占１３．１％。本文研究方法与传统

２４］

（的手工制图方法制作的福建省地貌类型图［１∶

２）面积（万ｋ０．１１２３３７９９１５９５３３ｍ　１．　０．　０．　３．　４．　１．

百比例（％）０．９０．１０２５．９０．８１．０　１　３．　８．　２　４　１

合计

０．９０．１０　１　３．　

３４．１　

５１．８

４　结论与讨论

经实验认为：

）（综合利用坡度和面积两个指标，计算坡度１

小于一定阈值的连通栅格的面积，能够较好快速实现平缓地与山丘地的自动划分；通过构建福建地方并与平缓地地貌单元叠加，可实现平地（平基准面，

原）与台（阶）地的快速划分；以Ｄ８算法为基础的流域水文分析方法能够有效地通过山丘地负地形提取福建省山丘地地貌形态实体单元的界线。该方法从中国地貌基本形态分类的相关定义出发，所提取统计单元能够保证地貌实体的完整性，并较好符合于福建传统地貌分类体系，分类结果从中观或宏观上反映出福建省各级地貌体的空间分异。（）通过专家知识修正，一定程度上可提高地２

为更大尺度范围地貌信貌基本形态类型分类精度，

息的快速、准确提取与集成研究奠定基础。但以文ＤＥＭ的地貌实体单元提取尚有许多不足之处：中构建的地方侵蚀基准面并不是真正地貌意义上的侵蚀基准面；因Ｄ特ＥＭ的精度和划分指标问题，别是地貌类型急剧变化区域，山地与丘陵间的划分等方面还存在不足，需在不同地貌区域进行人工调整。鉴此，在以后的研究中拟增加其他指标（如山，丘地锥形面判别指标等）以推进中国传统地貌分类体系下的地貌形态类型的快速分类与制图研究。随着空间对地观测技术的不断发展，ＤＥＭ的精度将不断提高，本文研究方法不仅适用于宏观和中观省域尺度，甚至于更大尺度的区域，应用前景广阔。

参考文献：

［］斯皮里顿诺夫Ａ　北京：地质出版１Ｎ．地貌制图学［Ｍ］．

社，１９５６．

［］沈玉昌．］中国地貌的类型与区划问题的商榷［中国第２Ｊ．

（）：四纪研究，１９５８，１１３３４１．－

相比，结合１∶１研究发现：本１６０万）０万的ＤＬＧ，

不文研究方法获得的地貌实体单元边界更加精细，受制图人员的认识水平及技术的影响，实现了以地貌实体为单元的省域地貌形态类型快速自动划分。从分类结果看：平原、平地和台地分类结果在两种图上的范围基本一致，而且本方法与基于规则单元

１４］

所获得的地貌图［相比，能更好地识别出分出河谷

（山间）平地和台（阶）地，而这恰恰是福建的特色所在，即大山带和山区间夹杂着许多大小不一的串珠状河谷、盆地。但丘陵、山地之间的分类与传统分类相比还存在一定的误判，使丘陵范围和面积有所增加，还需增加判别指标进一步加以区分

。

图３　福建省地貌形态类型图

Ｆｉ．３　ＧｅｏｍｏｒｈｏｌｏｉｃａｌｔｅｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅ　　　　ｇｐｇｙｐｊ

１期陈加兵等：基于ＤＥＭ与ＤＬＧ的福建省地貌形态自动分类

７９　　　

［］中国自然区划工作委员会．［中国地貌区划（初稿）３Ｍ］．

北京：科学出版社，１９５９．

［］４ＣｈｅｎＷ　Ｍ，ＺｈｏｕＣ　Ｈ，ＬｉＢＹ，ｅｔａｌ．Ｓｔｒｕｃｔｕｒｅａｎｄ　　　　　ｇ　

ｃｏｎｔｅｎｔｓｅｏｏｆｌａｅｒｅｄｃｌａｓｓｉｆｉｃａｔｉｏｎｓｓｔｅｍｏｆｄｉｉｔａｌ　　　　　　　－ｇｙｙｇ［］ｍｏｒｈｏｌｏｆｏｒＣｈｉｎａＪ．ＪｏｕｒｎａｌｏｆＧｅｏｒａｈｉｃａｌＳｃｉ　　　　－ｐｇｙｇｐ　，（）：ｅｎｃｅｓ２０１１，２１５７７１７９０．－

［］Ｈｒｏｅｒｔｉｅｓｅ５ａｍｍｏｎｄＥ　Ｈ．Ａｎａｌｓｉｓｏｆｉｎｌａｎｄｆｏｒｍ　　　　　　－ｐｐｇｙ

：Ａｏｒａｈｎａｌｉｃａｔｉｏｎｔｏｂｒｏａｄｓｃａｌｅｌａｎｄｆｏｒｍ　ｍａ　　　－　－ｇｐｙｐｐｐ］ｉｎＪ．ＡｎｎａｌｓｏｆＡｓｓｏｃｉａｔｉｏｎｏｆＡｍｅｒｉｃａｎＧｅｏｒａ　　　　　－ｐｇ［ｇ，（）：ｈｅｒｓ１９６４５４１１１９．－ｐ

［］６ＤｉｋａｕＲ，ＢｒａｂｂＥＥ，ＭａｒｋＲ　Ｍ．Ｌａｎｄｆｏｒｍｃｌａｓｓｉｆｉｃａ　　　　　－

［ｔｉｏｎｏｆＮｅｗ　ＭｅｘｉｃｏｂｃｏｍｕｔｅｒＭ］．Ｕ．Ｓ．Ｇｅｏｌｏｉｃａｌ　　　ｙｐｇ　，Ｓｕｒｖｅ１９９１．ｙ

［／７］ＢｒａｂｎＬ．ＧＩＳａｎａｌｓｉｓｏｆｍａｃｒｏｌａｎｄｆｏｒｍ［Ｍ］．／　　　　　ｙｙ

ＳＩＲＣ，Ｔｈｅ１０ｔｈＣｏｌｌｏｕｉｕｍｏｆｔｈｅＳａｔｉａｌＩｎｆｏｒｍａｔｉｏｎ　　　　　　ｑｐ：，ＲｅｓｅａｒｃｈＣｅｎｔｒｅ．ＮｅｗＺｅａｌａｎｄＵｎｉｖｅｒｓｉｔｏｆＯｔａｏ　　　ｙｇ　１９９８，３５４８．－

［］Ｍ８ｏｒａｎＪＭ，ＬｅｓｈＡ　Ｍ．Ｄｅｖｅｌｏｉｎｌａｎｄｆｏｒｍ　ｍａｓｕ　　　　－ｇｐｇｐ　

ｓｉｎＥＳＲＩ ｓＭｏｄｅｌＢｕｉｌｄｅｒ［Ｃ］．ＥＳＲＩＵｓｅｒＣｏｎｆｅｒｅｎｃｅ　　　ｇ　：／／／／２００５Ｐｒｏｃｅｅｄｉｎｓ．ｈｔｔｉｓ．ｅｓｒｉ．ｃｏｍｌｉｂｒａｒｕｓｅｒ－ｇｐｇｙ／／／，ｃｏｎｆｆｒｏｃ０５２００５．ａｅｒｓａ２２０６．ｄｆｐｐｐｐｐ

［９］ＤｒａｕｔＬ，ＢｌａｓｃｈｋｅＴ．Ａｕｔｏｍａｔｅｄｃｌａｓｓｉｆｉｃａｔｉｏｎｏｆ　　　　ｇ

ｌａｎｄｆｏｒｍｅｌｅｍｅｎｔｓｕｓｉｎｏｂｅｃｔｂａｓｅｄｉｍａｅａｎａｌｓｉｓ　　－　　ｇｊｇｙ　［］（）：Ｊ．Ｇｅｏｍｏｒｈｏｌｏ．２００６８１３３０３４４．－ｐｇｙ

［］１０ＰｒｉｍａＯＤＡ，ＥｃｈｉｏＡ，ＹｏｋｏａｍａＲ，ｅｔａｌ．Ｓｕｅｒ　　　　　　－ｇｙｐ

ｖｉｓｅｄｌａｎｄｆｏｒｍｃｌａｓｓｉｆｉｃａｔｉｏｎｏｆＮｏｒｔｈｅａｓｔＨｏｎｓｈｕｆｒｏｍ　　　　　　，ＤＥＭ－ｄｅｒｉｖｅｄｔｈｅｍａｔｉｃｍａｓ［Ｊ］．Ｇｅｏｍｏｒｈｏｌｏ　　ｐｐｇｙ（）：２００６７８３７３３８６．－

［］，１１ＩｗａｈａｓｈｉＪＰｉｋｅＲＪ．Ａｕｔｏｍａｔｅｄｃｌａｓｓｉｆｉｃａｔｉｏｎｓｏｆｔｏ　　　　　　－

ｏｒａｈｆｒｏｍ　ＤＥＭｓｂａｎｕｎｓｕｅｒｖｉｓｅｄｎｅｓｔｅｄｅａｎｓ　　　－ｍｐｇｐｙｙｐ　　［］ａｌｏｒｉｔｈｍａｎｄａｔｈｒｅｅａｒｔｅｏｍｅｔｒｉｃｓｉｎａｔｕｒｅＪ．Ｇｅ　　　－ｐ　　－ｇｇｇ

，（）：ｏｍｏｒｈｏｌｏ２００７８６４０９４４０．－ｐｇｙ

［］朱红春，陈楠民，刘海英，等．自１１２∶１００００比例尺ＤＥＭ

——以陕北黄土高原的实验为例［］提取地形起伏度—Ｊ．（）：测绘科学，２００５，３０４８６８８．－

［］刘爱利，汤国安．中国地貌基本形态Ｄ１３ＥＭ的自动划分

］（）：研究［地球信息科学，Ｊ．２００６，８４８１４．－

［］郎玲玲，程维明，朱启疆，等．多尺度Ｄ１４ＥＭ提取地势起

——以福建低山丘陵区为例［］地球伏度的对比分析—Ｊ．（）：信息科学，２００７，９６１８．－

［］程维明，周成虎，柴慧霞，等．中国陆地地貌基本形态类１５

］型定量提取与分析［地球信息科学学报，Ｊ．２００９，１１（）：６７２５７３６．－

［］曹伟超，陶和平，孔博，等．基于Ｄ１６ＥＭ数据分割的西南

］地区地貌形态自动识别研究［中国水土保持，Ｊ．２０１１（）：３３８４１．－

［］曹伟超，陶和平，孔博，等．利用最佳地形特征空间进行１７

］，地貌形态自动识别［武汉大学学报（信息科学版）Ｊ．（）：２０１１，３６１１１３７６１３８０．－

［］韦金丽，王国波，凌子燕．基于高分辨率Ｄ１８ＥＭ的地形

］特征提取与分析［测绘与空间地理信息，Ｊ．２０１２，３５（）：１３３３６．－

［］汪禹芹，李艳，刘爱利．基于面向对象思想的中国地貌１９

］（）：形态类型划分［遥感信息，Ｊ．２０１２１１３１８．－［］马士彬，安裕伦．基于Ａ２０ＳＴＥＲＧＤＥＭ数据喀斯特区域　

］：地貌类型划分与分析［地理科学，Ｊ．２０１２，３２（３）３６８－３７３．

［］肖飞，张百平，凌峰，等．基于Ｄ２１ＥＭ的地貌实体单元自

］（）：地理研究，动提取方法［Ｊ．２００８，２７２４５９４６６．－［］陈加兵，励惠国，郑达贤，等．基于Ｄ２２ＥＭ的福建省小流

］（）：地球信息科学，域划分研究［Ｊ．２００７，９２７４７７，９５．－［］中国科学院地理科学与资源研究所．中国１∶１２３００００００　　

］地貌制图规范［北京：科学出版社，Ｓ．１９８７．

［］张寒龙．福建省自然地图集［福州：福建省地图出２４Ｍ］．

版社，２０００．

ＡｕｔｏＣｌａｓｓｉｆｉｃａｔｉｏｎｏｆＧｅｏｍｏｒｈｏｌｏｉｃａｌＴｅｓＢａｓｅｄｏｎＤＬＧａｎｄ－　　　　　　　ｐｇｙｐ

ＤＥＭｆｏｒＦｕｉａｎＰｒｏｖｉｎｃｅ　　　ｊ

，，ＣＨＥＮＪｉａｂｉｎＬＩＨｕｉＣＨＥＮ　ＷｅｎＨｕｉａｎｄＺＨＥＮＧＤａｘｉａｎ　　　　　ｇ

（ＩｎｓｔｉｔｕｔｅｏＧｅｏｒａｈＦｕｉａｎ　ＮｏｒｍａｌＵｎｉｖｅｒｓｉｔＦｕｚｈｏｕ３５０００７，Ｃｈｉｎａ）　　ｆ　ｇｐｙ，ｊｙ，

：ＡｂｓｔｒａｃｔＢａｓｅｏｎＤｉｉｔａｌＬｉｎｅＧｒａｈｉｃ（ＤＬＧ）ａｔ１∶１００，０００ｓｃａｌｅａｎｄＤｉｉｔａｌＥｌｅｖａｔｉｏｎＭｏｄｅｌ（ＤＥＭ：　　　　　　　　　ｇｐｇ

），，３０×３０ｍｒｅｓｏｌｕｔｉｏｎｆｉｒｓｔｌｔｈｅｅｏｍｏｒｈｏｌｏｉｃａｌｔｅｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅｃａｎｂｅｄｉｖｉｄｅｄｉｎｔｏｔｗｏｂａｓｉｃ　　　　　　　　　　　ｙｇｐｇｙｐｊ：ａｃｃｏｒｄｉｎｔｏｔｈｅｔｗｏｍａｉｎｃｏｍｂｉｎａｔｉｏｎａｌｉｎｄｉｃａｔｏｒｓｓｌｏｅａｎｄａｒｅａｓｈａｌｌｏｗｓｌｏｅｒｅｉｏｎａｎｄｈｉｌｌｔｅｓ　　　　　　　　　－　　　－ｇｐｐｇｙｐ　

），ｍｏｕｎｔａｉｎ（ｓｔｅｅｓｌｏｅｒｅｉｏｎ．ＳｅｃｏｎｄｌｔｈｅｅｌｅｖａｔｉｏｎａｌｏｉｎｔｓｔｈａｔｃｏｎｔｏｕｒｓｔｈｅｒｉｖｅｒｓｓｔｅｍｉｎＤＬＧａｒｅ－　　　　　　　　　　ｐｐｇｙｐｙ，ｂＡｒｃＧＩＳｓｏｆｔｗａｒｅ．ＯｎｔｈｉｓｂａｓｉｓｔｈｅｒｅｌａｔｉｖｅｒｅｆｅｒｅｎｃｅｓｕｒｆａｃｅｏｆＦｕｉａｎＰｒｏｖｉｎｃｅｉｓｃｏｎｓｔｒｕｃｅｘｔｒａｃｔｅｄ　　　　　　　　　　　　－ｙｊ　

８０

地球信息科学学报　　　　　　　　　　　　　　　　　　　２０１３年

，ｔｅｄｔｈｒｏｕｈｓａｃｉａｌｉｎｔｅｒｏｌａｔｉｏｎｍｅｔｈｏｄ．Ｔｈｉｒｄｌｅｏｍｏｒｈｏｌｏｉｃａｌｅｎｔｉｔｉｅｓｉｎｓｔｅｅｓｌｏｅｒｅｉｏｎａｒｅｄｉ　　　　　　　－　　　－ｇｐｐｙｇｐｇｐｐｇ

，ｖｉｄｅｄｉｎｔｏｔａｅｒｕｎｉｔｓｗｉｔｈｔｈｅｍｅｔｈｏｄｏｆｗａｔｅｒｓｈｅｄｈｄｒｏｌｏｉｃａｌｍｏｄｅｌ．Ｆｉｎａｌｌａｃｃｏｒｄｉｎｔｏｔｈｅ１：１，　　　　　　　　　　　　ｐｙｇｙｇ　，０００，０００ＧｅｏｍｏｒｈｏｌｏｉｃａｌＭａｉｎＳｅｃｉｆｉｃａｔｉｏｎｆｏｒＣｈｉｎａｃｏｍｂｉｎｉｎｔｈｅｃｈａｒａｃｔｅｒｉｓｅｏｍｏｒｈｏｌｏｉｃａｌ　　　　　－ｐｇｐｐｇｐｇｇｐｇ　　

，ｔｉｃｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅａｎｄｏｖｅｒｌａｉｎｔｈｅａｂｏｖｅｅｎｔｉｏｎｅｄｅｎｔｉｔｉｅｓａｎｄｒｅｌａｔｉｖｅｒｅｆｅｒｅｏｍｏｒｈｏｌｏｉｃａｌ　　　　　－ｍ　　　　　－ｊｙｇｇｐｇ　，ｅｎｃｅｓｕｒｆａｃｅｏｆＦｕｉａｎＰｒｏｖｉｎｃｅｒｅｓｅｃｔｉｖｅｌｔｈｅｅｏｍｏｒｈｏｌｏｉｃａｌｅｎｔｉｔｉｅｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅｃａｎｂｅｄｉｖｉｄ　　　　　　　　　　　　　－ｊｐｙｊｇｐｇ，，，，，ｅｄｉｎｔｏｓｉｘｔｅｓｉ．ｅ．ｌｏｗｒｅｌｉｅｆｈｉｌｌｈｉｈｒｅｌｉｅｆｈｉｌｌｌｏｗｒｅｌｉｅｆｍｏｕｎｔａｉｎａｎｄｉｎｔｅｒｍｅｄｉｌａｉｎｌａｔｆｏｒｍ，　　　　　　　　　　　－ｙｐｇｐｐ：，ａｔｅｒｅｌｉｅｆｍｏｕｎｔａｉｎ．ＴｈｅｒｅｓｕｌｔｓｈｏｗｓａｍｏｎｅｏｍｏｒｈｏｌｏｉｃａｌｔｅｓｏｆＦｕｉａｎＰｒｏｖｉｎｃｅｈｉｌｌａｃｃｏｕｎｔｓ　　　　　　　　　ｇｇｐｇｙｐｊ　

ｆｏｒ３４．１％，ｍｏｕｎｔａｉｎａｃｃｏｕｎｔｓｆｏｒ５１．８％，ｂｏｔｈｔｏｅｔｈｅｒａｃｃｏｕｎｔｆｏｒ８５．９％，ａｎｄｌａｉｎ，ｆｌａｔａｎｄｌａｔｆｏｒｍ　　　　　　　　　　　ｇｐｐｔｏｅｔｈｅｒａｃｃｏｕｎｔｆｏｒｏｎｌ１３．１％．Ｔｈｅｒｅｓｕｌｔｉｓｂａｓｉｃａｌｌｉｄｅｎｔｉｃａｌｔｏｔｈａｔｂｔｒａｄｉｔｉｏｎａｌｍａｎｕａｌｍａｉｎ．　　　　　　　　　　　ｇｙｙｙｐｐｇ　　　Ｔｈｅｍｅｔｈｏｄｓｄｅｖｅｌｏｅｄｉｎｔｈｉｓａｒｔｉｃｌｅａｒｅｃｏｒｒｅｓｏｎｄｉｎｔｏｔｒａｄｉｔｉｏｎａｌｃｌａｓｓｉｆｉｃａｔｉｏｎｓｓｔｅｍｏｆｅｏｍｏｒｈｏｌ　　　　　　　　　　　　－ｐｐｇｙｇｐ　

，ｏｆｏｒＦｕｉａｎ，ａｎｄｃａｎａｕｔｏｍａｔｉｃａｌｌｃｌａｓｓｉｆｔｈｅｅｏｍｏｒｈｏｌｏｉｃａｌｅｎｔｉｔｉｅｓｆａｓｔａｎｄｗｅｌｌｓｕｃｈａｓｔｈｅ　　　　　　　　　　ｇｙｊｙｙｇｐｇ　　　，ｓｈａｌｌｏｗｓｌｏｅｒｅｉｏｎａｎｄｈｉｌｌｏｕｎｔａｉｎｒｅｉｏｎｌａｉｎａｎｄｌａｔｆｏｒｍ，ｗｈｉｃｈａｒｅｆｅａｓｉｂｌｅａｃｃｏｒｄｉｎｔｏｒｅｃｉ－　　　－ｍ　　　　　　　－ｐｇｇｐｐｇｐ　

，，ｍ，ｓｉｏｎａｎｄｅｆｆｉｃｉｅｎｃｓａｖｉｎｍａｎｏｗｅｒａｔｅｒｉａｌｒｅｓｏｕｒｃｅｓａｎｄｆｉｎａｎｃｉａｌｒｅｓｏｕｒｃｅｓｒｅａｔｌａｎｄｏｆｆｅｒｉｎａ　　　　　　　　ｙｇｐｇｙｇ　　ｎｅｗｔｈｉｎｋｉｎａｎｄｍｅｔｈｏｄｏｌｏｆｏｒｄｏｍｅｓｔｉｃａｎｄｉｎｔｅｒｎａｔｉｏｎａｌｓｉｍｉｌａｒｒｅｓｅａｒｃｈｅｓ．　　　　　　　ｇｇｙ　　：；ＫｅｗｏｒｄｓｅｏｍｏｒｈｏｌｏｉｃａｌｔｅｓＤＬＧ；ＤＥＭ；ＦｕｉａｎＰｒｏｖｉｎｃｅ　　ｇｐｇｙｐｊｙ