题目:

python对栅格加权叠加 栅格叠加分析_数据集

python对栅格加权叠加 栅格叠加分析_图层_02

 

python对栅格加权叠加 栅格叠加分析_数据集_03

 

 准备工作:

①新建数据库,将实验数据导入(利于后期矢量转栅格操作,未建数据库此步骤会一直出错);

②在工具条自定义窗口勾选spacial analyst模块,在菜单栏→扩展模块,勾选spacial analyst;确保实验后面的操作可以成功。

详细操作步骤:

1.对三个方案的规划路线(line.shp)建立50m缓冲区;打开ArcMap,在ArcToolBox→分析工具→邻域分析→缓冲区,打开缓冲区工具窗口,界面设置如图:线性单位填写50米,其余默认,单击“确定”。

python对栅格加权叠加 栅格叠加分析_数据分析_04

 缓冲区

python对栅格加权叠加 栅格叠加分析_栅格_05

 50m缓冲结果图

2.将生成的“缓冲区.shp”图层和土地利用类型图“landuse.shp”进行叠加分析;在ArcToolBox→分析工具→叠加分析→相交路径下,打开相交操作窗口,界面设置如下图所示;

python对栅格加权叠加 栅格叠加分析_图层_06

 相交

python对栅格加权叠加 栅格叠加分析_图层_07

 相交图层属性表

3.找出占耕地面积最少的方案;在菜单栏点击“选择”→“按属性选择”,窗口设置如下图所示;选择语句为:"方案" = '方案1' AND "LAND_USE" = '耕地';单击确定;打开“相交.shp”的属性表,选中“AREA”列,进行统计,统计图如下图所示;同理,对方案二方案三均进行此操作,结果如下。统计图显示,方案三所占用的耕地面积最小。

python对栅格加权叠加 栅格叠加分析_栅格_08

 按属性选择

python对栅格加权叠加 栅格叠加分析_图层_09

 方案一统计图

python对栅格加权叠加 栅格叠加分析_python对栅格加权叠加_10

 方案二统计图

python对栅格加权叠加 栅格叠加分析_数据分析_11

 方案三统计图

4.对方案3建立500m的缓冲区;首先按属性选择方案为“方案3”的路线;右击“line1.shp”图层导出方案3数据(“line方案三.shp”);附图如下;接着在ArcToolBox→分析工具→邻域分析→缓冲区,打开缓冲区工具窗口,界面设置如图:线性单位填写500米,其余默认,单击“确定”。

python对栅格加权叠加 栅格叠加分析_数据集_12

 按属性选择

python对栅格加权叠加 栅格叠加分析_python对栅格加权叠加_13

 导出数据

python对栅格加权叠加 栅格叠加分析_栅格_14

 缓冲区

python对栅格加权叠加 栅格叠加分析_数据集_15

 500m缓冲区图

5.将500m缓冲区与土地利用类型图进行叠加分析,得到500m缓冲区所覆盖的土地利用类型图;在ArcToolBox→分析工具→叠加分析→相交路径下,打开相交操作窗口,界面设置如下图所示;

python对栅格加权叠加 栅格叠加分析_栅格_16

 叠加分析“相交”

6.打开“landuse500.shp”的属性表,添加字段“权值”,设置短整型,根据实验指导书中的表1附权值。首先按属性对“landuse500.shp”进行选择,选中土地利用类型为“建设用地”的,再打开图层属性表,右击权值字段,选择“字段计算器”,赋权值为“1”;其余土地利用类型操作同上。相关截图如下。

python对栅格加权叠加 栅格叠加分析_图层_17

 添加字段

python对栅格加权叠加 栅格叠加分析_数据集_18

 按属性选择

python对栅格加权叠加 栅格叠加分析_数据分析_19

 字段计算器

python对栅格加权叠加 栅格叠加分析_图层_20

 Landuse500属性表

7.将“landuse500.shp”数据转化为栅格单元为50m*50m的栅格数据集;在ArcToolBox→转换工具→转为栅格→面转栅格,打开界面,相关设置如图;

python对栅格加权叠加 栅格叠加分析_数据分析_21

面转栅格

python对栅格加权叠加 栅格叠加分析_数据集_22

 Landuse500栅格数据图层

8.对DEM数据(ch)进行重采样操作;在ArcToolBox→数据管理工具→栅格→栅格处理→重采样,打开重采样窗口,设计如下,像元大小为50m,重采样技术默认即可,单击确定。

python对栅格加权叠加 栅格叠加分析_栅格_23

 重采样

9.生成坡度图(slope.shp),在ArcToolBox→spacial analyst工具→表面分析→坡度路径下,打开坡度窗口,进行如下设置:

python对栅格加权叠加 栅格叠加分析_栅格_24

 坡度

python对栅格加权叠加 栅格叠加分析_python对栅格加权叠加_25

 Slope栅格数据图层

10.按照实验指导书的表二对slope.shp重分类,在ArcToolBox→spacial analyst工具→重分类→重分类路径下,进行重分类,类别选择3,生成栅格数据集slopereclass;

python对栅格加权叠加 栅格叠加分析_图层_26

 重分类

11.将landuse500和slopereclass两个栅格数据集进行空间叠置分析,利用栅格计算器,将2个栅格数据集做乘运算,得到新的栅格数据集landuseslope。在ArcToolBox→spacial analyst工具→地图代数→栅格计算器路径下,打开栅格计算器窗口,设置如下图;

python对栅格加权叠加 栅格叠加分析_python对栅格加权叠加_27

 栅格计算器

python对栅格加权叠加 栅格叠加分析_数据分析_28

 Landuseslope栅格数据图层

12.按照实验指导书生态敏感性评价表对landuseslope.shp进行重分类,设置如下图; 

python对栅格加权叠加 栅格叠加分析_数据分析_29

 重分类landuseslope

python对栅格加权叠加 栅格叠加分析_图层_30

重分类栅格数据图层

13.制图输出

python对栅格加权叠加 栅格叠加分析_数据集_31