修改稿基于GIS的大连旅顺口区地形起伏度分析-w.doc

1基于 GIS 技术的大连旅顺口区地形起伏度测算研究（修改稿） 王慧鹏 1,王 利 1，2（1.辽宁师范大学城市与环境学院，中国 大连 116029；2.辽宁师范大学海洋经济与可持续发展研究中心，中国 大连 116029）摘要:依据地形起伏度研究的理论和方法，充分运用地理信息系统软件的分析功能和制图技术，以 1:5 万DEM 数据为基础，利用 0.01km2 的动态格网以街道为评价单元进行空间信息的处理和分析，计算大连旅顺口区的地形起伏度，制成地形起伏度专题地图,划分相应起伏区分析得出旅顺口区局部地区地形起伏度较大，起伏度指数以中低值为主，并测算出各街道实际用地情况，为本区域经济发展和城市建设提供参考，也为区域制定科学的发展规划提供依据关键词:GIS；地形起伏度；旅顺口区Measurement and Research of Relief Amplitude of Lyvshunkou District based on GISWanghui Peng1, Wang Li 1, 2,(1. College of Urban and Environment, Liaoning Normal University, Dalian, China, 116029; 2.Research Centerfor Marine Economy and Sustainable Development, Liaoning Normal University, Dalian, China, 116029)Abstract:This article fully uses analytical function and cartographic technique of geographic information system software according to theory and method for relief amplitude research. On the basis of Lyvshunkou District 1:50000 data, proceed processing and analysis of spatial information making use of 0.01km 2dynamic grid and township as statistical unit. The relief amplitude of Lvshunkou district is extracted and thematic map is drawn and divide the relief amplitude region . A conclusion is reached through analysis of relief amplitude index of local fluctuation is bigger which mainly focus on middle and low value. Kinds of topographic relief area and corresponding land using situation in various villages and towns are measured. This measured data can be reference in aspects of economic development, urban construction and establishment of regional scientific development planning.Key words: GIS; Relief amplitude; Lyvshunkou district1 引言地形分析是地理信息系统空间分析中的重要组成部分，也是各类空间规划的基础。

近年来随着 GIS 技术的不断发展及《全国主体功能区规划》的公布实施，关于地形起伏度的研究逐渐得到了许多专家学者的关注与重视，对地形起伏度的概念也存在多种表述，如朱红春等认为地形起伏度是指在所指定的分析区域内所有栅格中最大高程与最小高程的差，是反映地形起伏的宏观地形因子 [1]；刘新华等认为地形起伏度是指地面一定距离范围内最大的高程差，其实质仍是坡度概念的一个延伸，可以作为区域水土流失评价的地形指标 [2]；封志明等把地形起伏度又称为地表起伏度，是指区域海拔高度和地表切割程度的综合表征[3]由此可见，地形起伏度仍然需要继续深入地研究旅顺口区位于辽东半岛最南端,属大连市辖区东临黄海，西濒渤海，南与山东半岛隔海相望，北依大连，地理坐标为东经 120°8′至 121°35′之间，北纬 38°7′至 39°之间，是北抵俄罗斯、南达我国三亚的重要欧亚大陆桥，境内港口众多，区域位置优越旅顺口区陆地南北纵距 26.1 公里，东西横距 31.2 公里，总面积 506.8 平方公里，海岸线长 169.7 公里基金项目：国家科技支撑计划(2008BAH31B06) [Foundation item: national science and technogy support program,No.2008BAH31B06]第一作者简介:王慧鹏(1988—） ，男，山西省阳曲人，硕士研究生，研究方向为城市与区域规划。

通讯作者：王利(1965—)，男，辽宁省建平人，博士，教授，博士生导师主要研究方向为区域发展规划与地理信息系统应用2全境属长白山余脉构成的沿海丘陵地带，多山地丘陵，少平原低地，地势东高西低(参见图1)当前对旅顺口区的研究主要集中在区域产业生态链构建研究 [4]、自然旅游资源开发研究 [5]、生态农业建设规划研究 [6]及区域生态环境可持续发展研究 [7]等方面，还没有专门关于旅顺口区地形起伏度方面的研究所以本文基于 GIS 技术的支持，应用 1:5 万 DEM 数据，对大连旅顺口区地形进行详细研究，分析旅顺口区的地形起伏度特征，为在较小区域内的地形起伏度研究提供一定的方法，同时也为区域规划的制定提供参考图1 旅顺口区行政区划图 图2 旅顺口区60米等高线地形图Fig.1 The administrative division map of LyvShunKou Fig.2 The 60 meters terrain contour map of LyvShunKou 2 研究思路与方法2.1 研究思路地形起伏度的研究最早起源于 1948 年苏联科学院地理研究所提出的割切深度 [8]。

近年来，许多学者从各自的研究领域出发，在较大范围内提取大尺度的 DEM 数据分析其地形起伏度的空间分布特征及其对区域人口分布的影响 [9-11]；把地形起伏度同滑坡、泥石流等地质灾害相联系，为区域防灾减灾提供方法 [12-14]；将地形起伏度与城市土地利用类型结合起来，探讨地形对城市空间结构发展的影响 [15]；通过对区域地形起伏度的计算和分析，把区域划分为不同的起伏区，指出区域经济发展的方向 [16]此外，把地形起伏度作为地图绘制和地貌特征分析的重要组成部分正日益成为地形研究的基本发展方向 [17]从国内外的研究进展来看，关于地形起伏度研究主要集中在对较大区域宏观地形状况的研究上，而在小范围内提取高分辨率的地形起伏度并进行分析研究的还为数不多 根据大连城市规划、土地规划和发展规划等的用地要求，即“沿海地区 60 米等高线，坡度 25 度” ，作为建设用地限制区，但把沿海和内地两部分综合起来看，这项指标并不太适宜，需要用一个更加客观的指标——地形起伏度来衡量( 参见图 2)目前为止所获得的DEM 数据及地形分析基础不能满足大连旅顺口区区域规划编制的要求本文基于 GIS 技术的支持，应用 1:5 万 DEM 数据，利用高程、高差、坡度与坡向以及地形起伏度等参数对本区域地形做详细分析，为发展区域经济和编制区域规划提供参考。

2.2 研究模型与方法地形起伏度的计算方法有多种：早期计算地形起伏度的公式是从地形起伏度的定义出发，进行宏观上的地形起伏度计算，其公式如下：(1.1)maxinR=H-3(1.1)式中R—地形起伏度， —单位面积内最大高程值， —单位面积内最小高程值maxHminH牛文元等考虑了地形起伏度对生态环境脆弱的影响因素，把平地比重 [9]引入地形起伏度的公式中，但利用此公式计算的结果较为粗糙，各分析单元差异也不显著，其公式如下：(1.2)=ax(h)-in/ax()-in1-(A)/RDLS P(1.2)式中RDLS—地形起伏度指数， max(h)—地区的最高海拔高度，min(h)—地区的最低海拔高度，max(H)—全国陆地最高海拔高度， min(H)—全国陆地最低海拔高度，P(A) —地区平地所占面积，A—地区总面积 [18]封志明等将500m作为中国基准山体高度，在原有基础上更新了计算模型，分析了地形起伏度与人口分布的关系，但在小范围的区域里没有得到运用其公式如下：(1.3)=max(H)-in1-(A)/50RDLSP(1.3)式中RDLS—地形起伏度指数 ,max(H)—区域内最高海拔，min(H)—区域内最低海拔，P(A)—区域内的平地面积,A—区域总面积 [9]。

杜鹏等根据国家人口与计划生育委员会计算人居环境适宜性评价中地形起伏度指数的计算方法，结合了营口市的具体情况进一步改进了计算地形起伏度的公式，使其适合于较小区域内的地形起伏度分析其公式如下：(1.4)=/(H)+()-Minh1-(A)/(H)-VG(h)RDLSTVGaxP(1.4)式中：ALT—以某一栅格单元为中心一定区域内的平均海拔，Max(H)—该区域最高海拔，Min(h)—该区域最低海拔 ,P(A)—区域内平地面积，A—区域总面积，AVG(H)—区域平均高程，AVG(H)-AVG(h)—区域内平均高程差 [16]通过比较和分析关于计算地形起伏度的各种方法的优缺点以及区域适宜性，并且充分考虑到大连旅顺口区区域范围较小的实际情况，本文将采用(1.4)式来计算旅顺口区各评价单元的地形起伏度指数3 技术运用3.1 研究过程对地形起伏度的研究主要包括地形数据的收集和整理、数据空间运算以及对研究结果的分析与讨论等过程，具体研究流程如下：基础数据准备 插值(IDW) 等高线地形图 DEM 转出高程点数据空间运算 制作格网 提取等高线节点 起伏度专题地图 制作街道边界 行政区划图 分析地形起伏度 划分地形起伏区 4图3 地形起伏度分析流程图Fig. 3 The flow diagram based on the relief amplitude3.2 研究步骤地形起伏度的研究步骤主要是对MapInfo、Vertical Mapper等地形分析软件的配合使用；是对插值运算、赋值运算等数据运算的综合应用以及专题地图的绘制和分析。

具体步骤如下：第一步：数据准备采用MapInfo软件进行数据格式转换及等高线矢量化，匹配坐标系为大连旅顺口区坐标系第二步：点集提取Vertical Mapper(VM)是一套具备完整系统架构的3D(三维) 分析系统工具集，主要采用MapInfo Professional的宏语言MapBasic 开发而成的MBX扩展工具，密切配合MapInfo Professional使用VM主要用于分析空间上连续数据的趋势，通过多重插值路径(TIN 、IDW、Kriging 等定制点估计) ，可以生成渲染效果的图像(Grid格式)，也可以生成3D效果图，进而实现高程分析以及坡度、坡向等地形参数提取 [19]利用VM3.5制作大连旅顺口区数字高程模型，设定精度为10m ，并利用其POLY—TO—POINT 命令，提取等高线节点第三步：插值运算插值的方。