一、实验目的
了解栅格成本距离的基本原理,学会综合运用表面分析、重分类、栅格叠加、
成本距离等方法,制作成本栅格,并能够根据成本栅格进行路径分析。
二、实验内容
某旅游景区拟修建一条由停车场至观景台的步行道路,要求考虑地形坡度、
地表起伏度、土地利用状况等因素的影响,确定一条最低成本路径。
三、实验原理与方法
实验原理:基于栅格数据的路径选择源于路径影响因素的分析,其核心在于
计算成本距离和方向。它首先需要确定路径的影响因素,然后根据这些因素选择
原始数据或计算派生数据,在对它们进行量化和标准化的基础上,生成单因素的
成本数据。不同要素的影响强度有所差异,需要确定因素的权重,在此基础上对
单因素成本数据进行加权求和得到综合成本数据。相邻像元间的资源移动,既要
考虑距离,也要考虑成本,即成本距离。以起点为源点,计算每个像元至源点的
累积最小成本距离和成本方向,在确定了终点(目标点)后,即可根据成本方向
确定最低成本路径。
实验方法:首先,根据 DEM 数据派生出坡度数据和起伏度数据。其次,利
用重分类工具将坡度数据、起伏度数据和土地利用数据进行量化和标准化处理,
得到单因素的成本数据。再次,根据各因素的权重,采用栅格计算器对单因素成
本数据进行加权叠加,输出综合成本栅格数据。再其次,利用成本距离工具计算
栅格数据中各像元到源点的成本距离和方向,输出累积最低成本距离栅格和成本
方向栅格。最后,利用成本路径工具生成最低成本路径。
图 实验流程图
四、实验设备与数据
(1)实验设备:计算机。
(2)主要软件:ArcGIS
(3)实验数据:高程数据(DEM)、土地利用数据(TDLY)、停车场数据(StartPoint)、
观景台数据(EndPoint)。
五、实验步骤
1.设置分析环境
(1)运行 ArcMap,激活 SpatialAnalyst 模块。 (2)加载“data02”文件夹中的所有实验数据。
(3)在“data02”文件夹中,新建一个 FileGeodatabase 可命名为“实验结果”,
用于存放实验过程与结果数据。
(4)设置 ArcToolbox 分析环境。右键点击【ArcToolbox】,选择【环境】,在打
开的设置窗口中,将临时工作空间设置为上一步生成的 Geodatabase;将处理范
围中的范围和栅格分析中的像元大小及掩膜均设置为与 DEM 数据一致。
2.计算成本数据
(1)生成坡度数据。打开坡度工具(【ArcToolbox】→【SpatialAnalyst 工具】
→【表面分析】→【坡度】),将 DEM 作为输入栅格,输出数据以 podu 命名,
其余采用默认值,点击【OK】,得到坡度数据,图 1。
(2)生成起伏度数据。打开邻域统计工具(【ArcToolbox】→【Spatial Analyst 工具】→【邻域分析】→【焦点统计】),将 DEM 作为输入栅格,输出数据以 qifudu 命名,将邻域设置为 11ⅹ11 像元的矩形,统计类型选择“RANGE”,点击 【确定】,得到起伏度数据,图 2。
(3)坡度数据重分类。打开重分类工具(【ArcToolbox】→【Spatial Analyst 工 具】→【重分类】→【重分类】)。先将坡度数据作为输入数据,然后点击【分类】 按钮,采用相等间隔法将坡度分为 10 级,坡度最小的一级赋值为 1,最大的一 级赋值为 10(即将坡度数据进行 10 分制的标准化处理,坡度越大,通行成本越 高);最后将输出数据命名为“Reclasspodu”。点击【确定】,得到坡度成本数据 (图 3)。
(4)起伏度数据重分类。再次打开重分类工具,以起伏度数据作为输入数据, 仍采用等间距法将其分为 10 级,起伏度最小的一级赋值为 1,最大的一级赋值 为 10,输出命名为“ReclassQFD”点击【确定】,得到起伏度成本数据(图 4)。
(5)土地利用数据重分类。打开重分类工具,将土地利用数据作为输入数据, 分类字段选择“DLMC”(地类名称),根据土地利用类型对道路修建成本的影响 大小,将河流分值设为 10,建设用地分值设为 8,林地分值设为 7,园地分值设 为 6,耕地分值设为 4,草地分值设为 2,裸地分值设为 1,输出命名为 “ReclassTDLY”,点击【确定】,得到土地利用成本数据(图 5)。
(6)加权计算综合成本栅格数据。本实验中,考虑山地坡度和起伏度较土地利 用对修建道路的成本影响较大,因此设坡度的权重为 0.5,起伏度为 0.3,土地利 用为 0.2。打开栅格计算器(【ArcToolbox】→【Spatial Analyst 工具】→【地图代 数】→【栅格计算器】),如图 12.8 所示,在中间的文本框中输入"Reclasspodu" * 0.5 + "ReclassQFD" * 0.3 + "ReclassTDLY" * 0.2,输出命名为"Cost",点击【确 定】,得到综合成本栅格数据(图 6)。
3.计算成本距离和方向。
打开成本距离(【ArcToolbox】→【Spatial Analyst 工具】→【距离】→【成 本距离】,)以起点 ( StartPoint )作为输入的源数据, 以综合成本数据 ( Cost ) 作 为 输 入 的 成 本 数 据 , 输 出 的 成 本 距 离 数 和 成 本 方 向 数 据 分 别 命 名 为 “CostDistance” 和“Direction”,点击【确定】, 生成累积最低成本距离加权数据和 后向 连接(成本方向)数据, 结果如图 7 和图 8 所示。
4.计算最低成本路径 打开成本路径工具 (【ArcToolbox】→【Spatial Analyst 工具】→【距离】→ 【成本路径】,)如选择 EndPoint 作为输入的目标点数据, 分别选择上一步生成 的 CostDistance 和 Direction 数据作为输入的成本距离数据和成本方向数据, 输出数据命名为“Path”,路径类型设置为默认值,点击【确定】,生成最低成本 路径,图 9。
六、实验心得
熟悉ArcGIS栅格数据距离制图、成本距离加权、数据重分类、多层面合并等空间分析功能。熟悉掌握利用ArcGIS上述空间分析功能和结果类似学校选址的实际应用问题的基本流程和操作过程。