arcgis10中学选址栅格重分类叠合.docx
《arcgis10中学选址栅格重分类叠合.docx》由会员分享，可在线阅读，更多相关《arcgis10中学选址栅格重分类叠合.docx（9页珍藏版）》请在冰豆网上搜索。
arcgis10中学选址栅格重分类叠合
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\，激活dataframe1，有4个图层（见图6-1）：
点状图层“现有中学”、点状图层“人口调查”、面状图层“规划土地使用”、线状图层“道路”（只起限定栅格边界和背景图的作用）。
鼠标双击dataframe1，在弹出的DataFrameProperties对话框的General标签中将MapUnites和DisplayUnits从Unknown改为Meters（米），完成后按“确定”键关闭。
主菜单中选用“地理处理-环境…”，进一步设置：
工作空间-当前工作空间：
\ex12\temp
工作空间-临时工作空间：
\ex12/temp
输出坐标系：
与输入相同
处理范围：
与图层道路相同
栅格分析/像元大小/如下面的指定：
50
确定返回，完成初始设置。
选用菜单Customize/Extension…，加载SpatialAnalyst，在菜单Customize/Toolbars下勾选SpatialAnalyst，调出空间分析工具条。
图6-1Dataframe1的显示
ArcToolbox中选用菜单SpatialAnalystTools/Distance/欧氏距离Euclideandistance…，出现Euclideandistance参数设置对话框：
Distanceto：
现有中学下拉选择图层名，消防站为距离的参照点
Outputraster：
school键盘输入栅格数据名称，存放路径为初始化设置
Maximumdistance：
不限定最大的计算范围，保持空白
Outputcellsize：
50默认，使用初始化设置
按OK键，产生离开现有中学的距离分布图，自动加载（图6-2）。
图6-2离现有中学的距离栅格
选点击图层名school，ArcToolbox中选用菜单SpatialAnalysisTools/Reclass/重分类Reclassify…，系统出现Reclassify对话框，点击Classify…按钮，出现Classification对话框：
Method:
EqualInterval按等间距方式分类
Classes：
4共分4类
按OK键，返回Reclassify对话框，栅格数据被分为4类，还要按分析之前确定的要求修改分类的间距：
OldValues（原来值）NewValue（修改值）
0–5000
500–10001
1000–15002
1500–22003
在Outputraster中规定文件名及输出路径，这里为R_school，按OK键返回，系统产生新的分类图层R_school（见图6-3）。
图6-3离开现有中学的距离再分类图
3产生人口分布密度图、再分类
ArcToolbox-SpatialAnalysisTools-密度分析Density–核密度分析KernelDensity，继续设置：
Inputdata：
人口调查选择图层名，得到人口调查样本点
PopulationField：
POPU选择，该字段在“人口调查”要素属性表中
Outputraster：
\GIS\ex12\temp\popu_den键盘输入栅格数据名称，路径按初始化设置
Outputcellsize：
50已经在初始化时设定，默认
SearchRadius：
800计算密度时的搜寻半径
AreaUnits:
SQUARE_MAP_UNITS面积单位用地图单位的平方
按OK键，产生估计的人口密度分布图popu_den。
ArcToolbox中选用菜单SpatialAnalysisTools/Reclass/Reclassify…，出现Reclassify对话框，在Inputraster下拉菜单中选择popu_den，点击Classify…按钮，出现Classification对话框：
Method:
EqualInterval
Classes：
4
按OK键后，返回Reclassify对话框，栅格数据被自动分为4类，再按要求修改：
OldValues（原来值）NewValue（修改值）
0–0
–1
–2
–3
在Outputraster中指定产生栅格数据的名称及路径，这里名称为R_popu，路径按初始设定，按OK键返回，系统产生新的分类图层“R_popu”（见图6-4）。
本章开始，人口密度的面积按公顷（万平方米）考虑，练习中地图单位是米，计算密度时按地图单位的平方即平方米计算，因此栅格单元的取值和当初确定的指标之间是1：
10000的关系，再分类后和原来的定义一致。
图6-4重新分类后的人口密度图
4规划土地使用矢量图层转换为栅格图、再分类
鼠标右键打开“Attributeof土地使用”，可以看到该要素属性表有字段LANDUSE，为每个多边形地块的规划土地使用性质，按城市规划专业习惯，土地使用分类编码的意义为：
C：
商业，G：
绿地，M：
工业，R1：
一类居住，R2：
二类居住。
选用菜单ArcToolbox/转换工具ConversionTools/转为栅格ToRaster/PolygontoRaster…（矢量要素转换为栅格），系统出现PolygontoRaster参数设置对话框：
Inputfeatures：
土地使用选择图层名
ValueField：
LANDUSE选择字段名，该字段在“土地使用”要素属性表中，决定栅格单元取值
Outputraster：
/temp/ld_use输入栅格数据名称，路径按初始设置
Cellassigmenttype：
CELL_CENTER
Cellsize：
50
按OK键后，产生栅格状土地使用图层ld_use（见图6-5），打开图层属性表“Attributeofld_use”，该表的内容如下：
ValueCountlanduse
栅格单元取值取该值的单元数累计所对应的转换前的属性
11568M（工业用地）
2854C（商业用地）
34120R2（一类居住用地）
41079R1（二类居住用地）
5658G（绿地）
图6-5土地使用矢量图转化成栅格图
按评价的需要对LANDUSE字段再分类。
ArcToolbox中选用菜单SpatialAnalysisTools/Reclass/Reclassify…，系统出现Reclassify对话框，在Inputraster下拉菜单中选择ld_use，在Reclassfield中下拉选择landuse，表示对landuse字段进行重新赋值。
注意，这里应直接修改NewValues，具体操作如下：
OldValues（原来值）NewValue（修改值）
C（商业）1
G（绿地）0
M（工业）0
R1（一类居住）2
R2（二类居住）2
在Outputraster中指定栅格数据名为R_ld_use，按OK键，产生新的栅格图层R_ld_use（见图6-6）。
图6-6重新分类后的土地使用图层
5确定明显不合适的位置
本练习设定离开现有中学太近（<500m）、人口密度太低（<人/m）、土地使用为工业、绿地的范围内，不适合中学选址，这些位置要排除。
选用ArcToolbox/SpatialAnalysisTools/地图代数MapAlgebra/栅格计算器RasterCalculator，双击鼠标实现以下操作：
"R_popu"*"R_ld_use"*"R_school"
输出栅格：
/GIS/ex12/temp/calc01，确定，产生栅格数据集calc01，ArcToolbox中选用菜单SpatialAnalysisTools/Reclass/Reclassify…，继续设置：
输入栅格：
calc01
重分类字段：
Value
点击Classify…按钮，出现Classification对话框：
Method:
EqualInterval
Classes：
2
按OK键后，返回Reclassify对话框，栅格数据被自动分为2类，再按要求修改：
OldValues（原来值）NewValue（修改值）
00
1–121
在Outputraster中指定产生栅格数据的名称及路径，这里名称为r_calcu，路径按初始设定，按OK
键返回，系统产生新的分类图层“r_calcu”。
单元取零的位置，明显不合适选址（图6-7）。
图6-7选址范围被分成0和1两类
6计算综合评定指标
选用菜单ArcToolbox/SpatialAnalysisTools/地图代数MapAlgebra/RasterCalculator，出现RasterCalculator对话框，可供计算的图层名列在左侧图层选择框内，双击鼠标实现如下操作：
"R_popu"+"R_ld_use"+"R_school"
输出栅格：
/GIS/ex12/temp/calc02，确定，产生综合评价栅格数据集calc02，自动加载（图6-8）。
经加法计算，可以初步判定单元取值较大的位置，相对适合新建中学，反之不适合。
图6-8中学选址初步评价，单元取值大的位置相对合适
选用ArcToolbox/SpatialAnalysisTools/MapAlgebra/RasterCalculator，双击鼠标实现如下操作：
“calc02”*“r_ccalcu”
输出栅格：
/GIS/ex12/temp/calculation，确定，产生综合评价栅格calculation，自动加载，调出calculation图层属性，进入Symbology，调整符号如下：
显示：
Uniquevalues
Valuefield：
value
确定，除道路、现有中学、calculation外，关闭其他图层（图6-9）。
可以看出，颜色偏深的位置相对适合新建学校，颜色偏浅的不适合增设学校。
图6-9排除明显不合适的位置，颜色深的位置适合建中学
7分类计算面积
栅格单元为整数时，自动产生栅格取值属性表，打开AttributeofCalculation，可以看到栅格单元的汇总情况，tableoptions/AddField，
名称：
SUM_AREA
类型：
长整型
精度：
8
OK，右击SUM_AREA，选择FieldCalculator，用鼠标和键盘在文本框内输入：
[COUNT]*50*50，OK返回，得到评价值的汇总结果。
练习结束，选用菜单File/Exit，退出ArcMap。
升级会员 