1、 南阳师范学院 20XX 届毕业生 毕业论文（设计） 题 目： 基于 SRTM-DEM的河南省基本地貌类型 提取方法研究 完 成 人： 班 级： 学 制： 专 业： 测绘工程 指导教师： 完成日期： 目 录 摘要 (1) 0 引言 (1) 1 研究区地理位置及地形地貌概述 (1) 2 数据准备 (2) 2.1 数据介绍 (2) 2.2 数据下载 (2) 2.2.1DEM 数据的下载 (2) 2.2.2 矢量数据的下载 (3) 2.2.2 河南地形数据的获取 (3) 2.3 我国三大常用坐标系区别（北京 54、西安 80 和 WGS-84） (4) 2.4 我国常用高程坐标系 (5) 3 地貌专
2、题图的制作 (6) 3.1 高程专题图 (6) 3.2 坡度专题图及分析 (7) 3.3 地形起伏度专题图及分析 (8) 4 地貌类型划分 (9) 4.1 河南省地貌类型划分 (9) 4.2 河南省县 级 行政地貌类型图 (11) 5 结语 (11) 5.1 河南地形 (12) 5.2 研究近况 (12) 参考文献 (12) Abstract (13) 第 1 页 (共 13 页 ) 基于 SRTM-DEM 的河南省基本地貌类型提取方法研究 摘要： 本文是基于美国 SRTM DEM 数据，在 GIS 平台下，通过对 SRTM DEM 数据的处理，制作 了河南省数字高程模型。用高程模型及起伏度模
3、型相结合，可初步实现基本地貌类型的自动提取，其结果可很好的反映区域地貌宏观规律。综合考虑这些制图信息，然后基于最新的中国基本地貌类型划分标准，可将河南省按地貌划分为低海拔平原、低海拔台地、低海拔丘陵、小起伏低山、中起伏低山、小起伏中山、中起伏中山、大起伏中山、中起伏中高山、大起伏中高山 10 类。研究表明，这种对地貌类型的自动提取方法在一定程度上改进了传统手工和野外调查为主的提取方法，实现了地貌信息的数字化。 关键词： 河南省；地貌特征；地貌类型；地理信息系统；数字高程模型 0.引言 地貌是自然地域综合体的主导因素，直接影响甚至决定着其他要素特征，与经济建设、社会发展、防灾减灾等方面关系密切
4、1，对其准确、快速地提取工作，不仅关系到更高级别地貌信息的准确提取，同时还对进一步的地貌格局分析、地貌区划等有着重要意义，并且对地球表面宏观地域分异、水分、能量的再分配等有着重要影响 ,并间接影响着土壤、 植被以及物质迁移和生态系统的总体演替和发展 2。地貌学的发展经历了一个定性描述到定量分析研究的阶段，传统地貌特征分析方法常以定性为主，主要是通过纸质地形图或野外调查获取数据资料进行，这种方法数据庞大、计算繁琐、成果精度低、主观性强、使得地貌学的定量研究发展缓慢。 20 世纪 90 年代以来，随着计算机、空间技术、3S(GIS、 RS、 GPS)、技术的发展和成熟，人们开始基于数字高程模型等数
5、据源尝试地貌特征提取方法，如张会平等利用 DEM 数据自动提取坡度等地形要素来研究岷江上游地貌形态 3；汤国安等利用 DEM来研究地 形定量因子的挖掘及地形起伏度的研究等 4,5；李军峰等利用 DEM 自动提取沟谷网络及节点技术等来研究水流累计量均值分级统计规律 6等。对基于 DEM 的河南省基本地貌类型类型提取方法的研究，以期实现地貌信息数字化的提取及基于领域分析的地势起伏度地貌形态的提取。 1.研究区地理位置及地形地貌概述 第 2 页 (共 13 页 ) 虽然国内外已将数字高程模型运用于地貌信息提取中，但由于地貌本身的区域性和复杂性，其方法还并不成熟，有由于河南省位于我国中部偏东，在黄河中
6、下游，全境处于北纬 31 23 -36 22和东经 110 21 -116 39之 间，南北跨越约540km，东西横跨约 580km，全省总面积约 16.51 万 km。河南分属四大流域，北部为海河流域，中部为黄河流域，西南为长江流域，东南为淮河流域，河南地形复杂，地势基本西高东低，北部为太行山余脉，西部为秦岭余脉，南部为大别山，中南部自西北向东南为横亘 400km 的伏牛山，东部为广阔的黄淮海平原，根据中国大陆地貌自西向东呈现出的 3个巨大地貌台阶逐渐急剧降低的特点，河南在全国地貌中的位置是跨第二和第三两级地貌台阶，河南地貌条件十分复杂，类型多种多样，形态结构和区域差异性显著 7。 2.数据
7、准备 2.1 数据介绍 我们常用的免费的 DEM数据就是 SRTM的 DEM 数据。 SRTM 数据主要是由美国太空总署（ NASA） 和国防测绘局（ NIMA） 联合测量的， SRTM 的全称是 Shuttle Radar Topography Mission，即航天飞机雷达地形测绘使命， 2000 年 2月 11 日上午 11 时 44 分，美国 “ 奋进 ” 号航天飞机在佛罗里达州卡那维拉尔角的航天发射中心发射升空， “ 奋进 ” 号上搭载的 SRTM 系统共计进行了 222 小时 23 分钟的数据采集工作，获取北纬 60 度至南纬 56 度之间，面积超过 1.19 亿平方公里的 9.8
8、万亿 字节的雷达影像数据，覆盖全球陆地表面的 80%以上，该计划共耗资 3.64 亿美元，获取的雷达影像数据经过两年多的处理，制成了数字地形高程模型，该测量数据覆盖中国全境。 SRTM 数据每经纬度方格提供一个文件，精度有 1 arc-second 和 3 arc-seconds 两种，称作 SRTM1和 SRTM3,或者称作 30M 和 90M 数据， SRTM1的文件里面包含 3601*3601个采样点的高度数据， SRTM3 的文件里面包含 1201*1201个采样点的高度数据。目前能够免费获取 中国境内的 RTM3 文件，是 90m的数据，每个 90m 的数据点是由 9个 30m 的数
9、据点算术平均得来的。 2.2 数据下载 2.2.1 DEM 数据的下载 SRTM 是目前几种重要的全球可获得性（ Global available） DEM 之一。 我第 3 页 (共 13 页 ) 们此次下载的数据就是从最原始的网址下载的 1。 它将全球分为 2472 的大小网格。根据河南的面积大小和经纬度可以知道，我们需要下载的有 SRTM_59_05、SRTM_59_06、 SRTM_60_05、 SRTM_59_06 四块数据。存储到一个文件夹里，解压即可。 2.2.2矢量数据的下载 要想获得河南省的高程地图仅有 DEM 数据是不够的，我们还需要得到河南省省界数据。这就需要先获得全国省
10、级行政区划，在从其中提取出河南省省界图，这个数据需要从下面这个网站获得 2。下载得到的是一个压缩数据，解压省界 bou2_4m.zip 即可得到全国省界地图。 2.2 3 河南地形数据的获取 地形主要分为平原、高原、山地、丘陵、盆地五中类型 8,9。而河南的地形中又包括其中的四种（除高原外）。 （ 1）平原：平原是陆地上最平坦的地域，它好像铺在大地上的绿色地毯，坦荡千里，辽阔无垠。平原地貌宽广低平，起伏很小，海拔多在 200 米以下。世界平原总面积约占全 球陆地面积的 1/4.平原可以分成堆积平原和侵蚀平原两类，堆积平原中的冲积平原，主要由河流冲积而成。它的特点是地面平坦，面积广大，多分布在大
11、江、大河的中、下游两岸地区。另一类侵蚀平原，主要由海水、风、冰川等外力的不断剥蚀、切割而成。这种平原地面起伏较大。 （ 2）高原：雄伟挺拔、险峻奇峭、蜿蜒起伏的高原，是在长期、连续、大面积的地壳抬升运动中形成的，海拔高度一般是在 1000 米以上。有的高原表面宽广平坦，地势起伏不大；有的高原则奇峰峻岭、山峦起伏，地势变化很大。 （ 3）山地：地球陆地的表面，有许多蜿蜒 起伏、巍峨奇特的群山。山由山顶、山坡和山麓三个部分组成，平均高度都在海拔 500 米以上。它们以较小的峰顶面积区别于高原，又以较大的高原区别于丘陵。这些群山层峦叠嶂，群居一起，形成一个山地大家族。 （ 4）丘陵：丘陵是陆地上起伏
12、和缓、连绵不断的高地。根据起伏高度，相对高度小于 100米者一般称为低丘陵，而 100-200米者一般称为高丘陵。 （ 5） 盆地：盆地四周高、中间低，整个地形象一个大盆。盆地的四周一般有高原或山地围绕，中间是平原或丘陵。盆地主要有两种类型。一种是地壳构1数据链接： http:/srtm.csi.cgiar.org/SELECTION/inputCoord.asp 2数据连接： http:/ 第 4 页 (共 13 页 ) 造运动形成的盆地，称为构造盆地， 如我国新疆的吐鲁番盆地、江汉平原盆地。另一种是由冰川、流水、风和岩溶侵蚀形成的盆地，称为侵蚀盆地，如我国云南西双版纳的景洪盆地，主要由澜沧
13、江及其支流侵蚀扩展而成 。 而我们所需要的最终数据需要经过一下几步操作。 （ 1） 需要拼接四块 DEM 数据，用 Data Management Tools / Raster / Mosaic To New Raster 进行拼接。 （ 2） 对已拼接的数据进行投影变化，在 Mosaic To New Raster 中选Coordinate system for the raster (optional) / Projected Coordinate Systems / UTM / WGS 1984 / 49N.prj 完成投影变换。 （ 3） 用 Editor工具从省界图中提取出河南省。
14、（ 4） 对数据进行掩膜处理，在工具箱中选 Spatial Analyst Tools / Extraction / Extract by Mask 可在 DEM数据中提取出河南的 DEM 地形图。 2.3 我国三大常用坐标系区别（北京 54、西安 80和 WGS-84） （ 1）北京 54 坐标系（ BJZ54） : 北京 54 坐标系为参心大地坐 标系，大地上的一点可用经度 L54、纬度 M54和大地高 H54定位，它是以克拉索夫斯基椭球为基础，经局部平差后产生的坐标系。 1954 年北京坐标系的历史：新中国成立以后，我国大地测量进入了全面发展时期，再全国范围内开展了正规的，全面的大地测量
15、和测图工作，迫切需要建立一个参心大地坐标系。由于当时的 “ 一边倒 ” 政治趋向，故我国采用了前苏联的克拉索夫斯基椭球参数，并与前苏联 1942年坐标系进行联测，通过计算建立了我国大地坐标系，定名为 1954年北京坐标系。因此， 1954 年北京坐标系可以认为是前苏联 1942年坐标系的延伸。 它的原点不在北京而是在前苏联的普尔科沃。北京 54 坐标系，属三心坐标系，长轴6378245m，短轴 6356863，扁率 1/298.3。 (2)西安 80 坐标系： 1978 年 4 月在西安召开全国天文大地网平差会议，确定重新定位，建立我国新的坐标系。为此有了 1980 年国家大地坐标系。 198
16、0 年国家大地坐标系采用地球椭球基本参数为 1975年国际大地测量与地球物理联合会第十六届大会推荐的数据，即 IAG 75 地球椭球体。该坐标系的大地原点设在我国中部的陕西省泾阳县永乐镇，位于西安市西北方向约 60 公里，故称 1980年西安坐标系 ，又简称西安大地原点。基准面采用青岛大港验潮站 1952 1979年确定的黄海平均海水面（即 1985 国家高程基准）。西安 80 坐标系，属三心坐标系，长轴 6378140m，短轴 6356755，扁率 1/298.25722101。 第 5 页 (共 13 页 ) （ 3） WGS-84 坐标系： WGS 84 坐标系（ World Geode
17、tic System）是一种国际上采用的地心坐标系。坐标原点为地球质心，其地心空间直角坐标系的 Z 轴指向国际时间局（ BIH） 1984.0 定义的协议地极（ CTP）方向， X 轴指向 BIH1984.0的协议子午面和 CTP 赤道的交点， Y轴与 Z轴、 X轴垂直构成右手坐标系，称为1984 年世界大地坐标系。这是一个国际协议地球参考系统（ ITRS），是目前国际上统一采用的大地坐标系。 GPS 广播星历是以 WGS-84 坐标系为根据的。 WGS84坐标系，长轴 6378137.000m，短轴 6356752.314，扁率 1/298.257223563。 由于采用的椭球基准不一样，并
18、且由于投影的局限性，使的全国各地并不存在一至的转换参数。对于这种转换由于量较大，有条件的话，一般都采用 GPS联测已知点，应用 GPS 软件自动完成坐标的转换。当然若条件不许可，且有足够的 重合点，也可以进行人工解算。 我们此次用的是 WGS-84 坐标系。 2.4 我国常用高程系 （ 1） 1956 年黄海高程系： 1956 年黄海高程系是在 1956 年确定的。 它是根据青岛验潮站 1950 年到 1956 年的黄海验潮资料，求出该站验潮井里横按铜丝的高度为 3.61 m，所以就确定这个钢丝以下 3.61 米处为黄海平均海水面。从这个平均海水面起，于 1956年推算出青岛水准原点的高程为
19、72.289m。 （ 2） 1985国家高程基准： 国家 85 高程基准其实也是黄海高程基准，只不过老的叫 “1956 年黄海高程系统 ” ，新 的叫 “1985 国家高程基准 ” ，新的比旧的低 0.029m。 我国于 1956年规定以黄海 (青岛 )的多年平均海平面作为统一基面，为中国第一个国家高程系统，从而结束了过去高程系统繁杂的局面。但由于计算这个基面所依据的青岛验潮站的资料系列（ 1950年 1956 年）较短等原因，中国测绘主管部门决定重新计算黄海平均海面，以青岛验潮站 1952 年 1979 年的潮汐观测资料为计算依据，并用精密水准测量接测位于青岛的中华人民共和国水准原点，得出
20、1985年国家高程基准高程和 1956年黄海高程的关系为： 1985年国家高程基准高程 =1956 年黄海高程 -0.029m。 1985年国家高程基准已于 1987年 5 月开始启用， 1956 年黄海高程系同时废止。我国目前通用的高程基准是：85 高程基准 。 我国各高程系统之间的关系如表 1所示。 第 6 页 (共 13 页 ) 表 1 各高程系统之间的关系 56 黄海高程基准 +72.289m 85 高程基准（最新的黄海高程） 56 高程基准 -0.029 吴淞高程系统 56 高程基准 +1.688 珠江高程系统 56 高程基准 -0.586 3.地貌专题图的制作 地貌即地球表面各种形
21、态的总称，也叫地形。地表形态是各种多样的，成因也不尽相同，是内、外力地质作用对地壳综合作用的结果。内力地质作用造成了地表的起伏，控制了海陆分布的轮廓及山地、高原、盆地和平原的地域配置，决定了地貌的构造格架。而外营力（流水、风力、太阳辐射能、大气和生物的生长和活动）地质作用，通过多种方式，对地壳表层物质不断进行风化、剥蚀、搬运和堆积，从而形成了现代地面的各种形态。 地表起伏的形态，如陆地上的山地、平原、河谷、沙丘，海底的大陆架、大陆坡、深海平原、海底山脉等。根据地表形态规模的大小，有全球地 貌，有巨地貌，有大地貌、中地貌、小地貌和微地貌之分。大陆与洋盆是地球表面最大的地貌单元，较小的地貌形态如有
22、在流水和风力作用下形成的沙垄和沙波等。地貌是自然地理环境的重要要素之一，对地理环境的其他要素及人类的生产和生活具有深刻的影响。 而我们此次主要介绍在栅格 DEM 的基础上，利用 GIS 软件得出研究区的高程专题图、坡度专题图、地形起伏度专题图等的方法与过程。 3.1 高程专题图及分析 DEM 所提供的原始数据就是高程数据，是 DEM 最基本的信息，根据研究区实际可知河南的最高的高程是 2392m(误差数据 ),最低的高程是 -65m（误差数据）。河南实际的最高峰是位于灵宝市小秦岭的老鸦贫 ,高程是 2413.8m,最低处为淮河出省出 23.2m。 我们在 Layer Properties 中选
23、 Classified 将高程分为 2000m。然后作图，步骤是：右击图层单击 Layer Properties/Classied/Classify 进行分类，完成河南省高程的自动提取，如下所示图 1 为河南省的高程分级结果图。由此我们可以看出，第 7 页 (共 13 页 ) 河南省高程空间分布特点 ：总体地形西部高，东部低，由西向向东呈阶梯状逐渐降低，平原主要分布于中部与东部，丘陵主要分布于东北部和南部，低山分布中部和东北部小部分区域，中高山分布于西部分地区。 图 1 河南省高程专题图 3.2 坡度专题图及分析 将研究区分为平原至微倾斜坡（ 45） 6 个等级，得到需要的专题图。 坡度的提取
24、过程为： （ 1）单击 Spatial Analyst 下拉箭头，单击 Surface Analysis，在弹出的下一菜单中单击 Slope，打开 Slope对话框。 （ 2）单击 Input surface 下拉箭头，选择输入栅格数据集。 （ 3）选择坡度表示方法，默认选项为度。 （ 4）在 Z factor文本框中输入高程变换系数。 （ 5）在 Output cell size 文本框中输入栅格大小。 （ 6）在 Output raster 文本框中输入输出坡度的存放路径与文件名。 第 8 页 (共 13 页 ) （ 7）单击 OK按钮，完成操作。结果如图 2所示。 图 2 河南省坡度专题
25、图 3.3 地形起伏度专题图及分析 地形起伏度又叫做地势或相对 高程，是指地面一定范围内最大的高程差，其实质是坡度在宏观区域的一个延伸，在宏观的区域内反映了地面的起伏度特征。陈志明 10对中国地貌形态按地形起伏度分为微缓起伏（ 0 20m） 、小起伏（ 20 75m） 、中起伏（ 75 200m）、山地起伏（ 200 600m） 和高山起伏（ 600m）等五级，我们依据此作河南省地形起伏度专题图。 地形起伏度的提取方法如下： （ 1）选中 DEM 数据，在 Spatial Analyst 下使用栅格邻域计算工具Neighborhood Statistic。设置 Statistic type 为最大值，邻域的类型为矩形（也可以是圆），邻域的大小为 7*7，则可得到一个邻域为 7*7 矩形的最大值层面，记为 A。 （ 2）重复第 1步，只是把 Statistic type 值设置为最小值，即可得到 DEM数据的最小值层面，记为 B。 （ 3）在 Spatial Analyst下使用栅格计算器 Raster Calculator，公式为 ：
