arcgis10：中学选址-栅格重分类、叠合.docx

1 中学选址的依据本练习要求对现有中学、人口分布、规划土地使用进行分析，为新建中学的选址提供依据，可利用的数据有三种：（1）现有中学点状空间数据，新建中学不应离现有中学太近，为此建立离开现有中学的距离图2）人口分布点状空间数据，新建中学应在人口比较密集的地方，为此利用点状人口调查数据产生密度图3）规划土地使用矢量多边形空间数据，某些用地不应该建中学，如工业用地，某些用地不太适合建中学，如商业用地，而居住用地适合建中学根据上述数据和基本依据，将离开现有中学的距离分为 4 类，分别赋予评定指标：0～500 米：0，500～1000 米：1，1000～1500 米：2，>1500 米：3 栅格图层① 人口密度分为 4 类（人/公顷），分别赋予评定指标：0～50：0， 50～100：1， 100～200：2， >200：3 栅格图层② 土地使用再分为 3 类，分别赋予评定指标：工业或绿地：0， 商业：1， 居住：2 栅格图层③ 初步的是综合评定指标：图层① + 图层② + 图层③ （栅格图层④）得到栅格图层④，单元取值越大的位置，越适合选为中学但是离现有中学太近（<500m），人口密度太低（<50人/公顷），规划土地为工业、绿地，肯定不和合适，一旦符合这3个条件之一，就要排除，无论其他条件如何，为此需要特别处理。

单项指标相乘：图层① × 图层② × 图层③ （栅格图层⑤），只要某单项指标出现零，在栅格⑤中该单元肯定为零将栅格⑤分为两类，0：=0；1：>0（栅格⑥），栅格图层⑤中取零的保留，大于0的都取1综合评价：栅格图层④×栅格⑥（栅格⑦）通过乘法计算，栅格数据集⑦中取值为零的单元为明显不合适、要排除的位置，其他单元的取值依然是相加的结果，综合取值较大的位置，重新建立中学2 产生离开现有中学的距离图、再分类启动地图文档\ex12\ex12.mxd，激活 data frame1，有 4 个图层（见图6-1）：点状图层“现有中学”、点状图层“人口调查”、面状图层“规划土地使用”、 线状图层“道路”（只起限定栅格边界和背景图的作用）鼠标双击 data frame1， 在弹出的 Data Frame Properties 对话框的 General 标签中将 Map Unites 和 Display Units 从 Unknown改为 Meters（米），完成后按“确定”键关闭主菜单中选用“地理处理-环境…”，进一步设置：工作空间-当前工作空间：\ex12\temp 工作空间-临时工作空间：\ex12/temp输出坐标系：与输入相同处理范围：与图层道路相同栅格分析/像元大小/如下面的指定：50确定返回，完成初始设置。

选用菜单 Customize / Extension…，加载 Spatial Analyst，在菜单 Customize / Toolbars 下勾选 Spatial Analyst，调出空间分析工具条图 6-1 Data frame1 的显示ArcToolbox中选用菜单 Spatial Analyst Tools / Distance / 欧氏距离Euclidean distance…，出现 Euclidean distance参数设置对话框：Distance to：现有中学 下拉选择图层名，消防站为距离的参照点Output raster：school 键盘输入栅格数据名称，存放路径为初始化设置Maximum distance： 不限定最大的计算范围，保持空白Output cell size：50 默认，使用初始化设置按 OK 键，产生离开现有中学的距离分布图，自动加载（图6-2）图6-2 离现有中学的距离栅格选点击图层名 school，ArcToolbox中选用菜单 Spatial Analysis Tools / Reclass / 重分类Reclassify…，系统出现 Reclassify 对话框，点击 Classify…按钮， 出现 Classification 对话框：Method: Equal Interval 按等间距方式分类Classes：4 共分 4 类按 OK 键，返回 Reclassify 对话框，栅格数据被分为 4 类，还要按分析之前 确定的要求修改分类的间距：Old Values（原来值） New Value（修改值）0 – 500 0500 – 1000 11000 – 1500 21500 – 2200 3在 Output raster 中规定文件名及输出路径，这里为 R_school，按 OK 键返回， 系统产生新的分类图层 R_school（见图 6-3）。

图 6-3 离开现有中学的距离再分类图3 产生人口分布密度图、再分类ArcToolbox- Spatial Analysis Tools- 密度分析Density – 核密度分析Kernel Density，继续设置：Input data：人口调查 选择图层名，得到人口调查样本点Population Field：POPU 选择，该字段在“人口调查”要素属性表中Output raster：\GIS\ex12\temp\popu_den 键盘输入栅格数据名称，路径按初始化设置Output cell size：50 已经在初始化时设定，默认Search Radius：800 计算密度时的搜寻半径Area Units: SQUARE_MAP_UNITS 面积单位用地图单位的平方按 OK 键，产生估计的人口密度分布图 popu_denArcToolbox中选用菜单 Spatial Analysis Tools / Reclass /Reclassify…，出现 Reclassify 对话框，在 Input raster 下拉菜单中选择 popu_den，点击 Classify…按钮，出现 Classification 对话框：Method: Equal IntervalClasses：4按 OK 键后，返回 Reclassify 对话框，栅格数据被自动分为 4 类，再按要求 修改：Old Values（原来值） New Value（修改值）0 – 0.005 00.005 – 0.01 10.01 – 0.02 20.02 – 0.03 3在 Output raster 中指定产生栅格数据的名称及路径，这里名称为 R_popu，路径按初始设定，按 OK 键返回，系统产生新的分类图层“R_popu”（见图 6-4）。

本章开始，人口密度的面积按公顷（万平方米）考虑，练习中地图单位是米， 计算密度时按地图单位的平方即平方米计算，因此栅格单元的取值和当初确定的指标之间是 1：10000 的关系，再分类后和原来的定义一致图 6-4 重新分类后的人口密度图4 规划土地使用矢量图层转换为栅格图、再分类鼠标右键打开 “Attribute of 土地使用” ，可以看到该要素属性表有字段 LANDUSE，为每个多边形地块的规划土地使用性质，按城市规划专业习惯，土地使用分类编码的意义为：C：商业，G：绿地，M：工业，R1：一类居住，R2： 二类居住选用菜单 ArcToolbox/转换工具Conversion Tools/ 转为栅格To Raster / Polygon to Raster…（矢量要素转换为栅格），系统出现 Polygon to Raster 参数设置对话框：Input features：土地使用 选择图层名Value Field：LANDUSE 选择字段名，该字段在“土地使用”要素属性表中， 决定栅格单元取值Output raster：/temp/ld_use 输入栅格数据名称，路径按初始设置Cell assigment type：CELL_CENTER Cellsize：50 按 OK 键后，产生栅格状土地使用图层 ld_use（见图 6-5），打开图层属性表 “Attribute of ld_use”，该表的内容如下：Value Count landuse栅格单元取值 取该值的单元数累计 所对应的转换前的属性1 1568 M （工业用地）2 854 C （商业用地）3 4120 R2 （一类居住用地）4 1079 R1 （二类居住用地）5 658 G （绿地）图 6-5 土地使用矢量图转化成栅格图 按评价的需要对 LANDUSE字段再分类。

ArcToolbox中选用菜单Spatial Analysis Tools / Reclass /Reclassify…，系统出现 Reclassify 对话框，在 Input raster 下拉菜单中选择 ld_use，在 Reclass field中下拉选择 landuse，表示对 landuse 字段进行重新赋值注意，这里应直接修改New Values，具体操作如下：Old Values（原来值） New Value（修改值）C （商业） 1G （绿地） 0M （工业） 0R1（一类居住） 2R2（二类居住） 2在 Output 。