ArcGIS 空间分析基本操作

空间分析基本操作

这里笔者使用的软件版本是 ArcGIS 10.5

一、前提条件

当进行空间数据分析会报错，进行空间分析，需要先打开扩展模块，就是Customize→Extensions，检查第二个图扩展模块的扩展是否打开。


----------------------------分隔线-------------------------------

二、空间分析的基本操作

1. 了解栅格数据

加载栅格数据，我们可以查看栅格数据的属性。

2. 用任意多边形裁剪栅格数据（空间数据的叠加，或者说两个不同图层的叠加怎么操作？）

两个不同数据的叠加怎么操作？ 【两种方法】
首先创建一个面状要素并绘制一个面

（1） 第一种方法：Extract by Msak（掩膜提取）

Extract by Msak 这个方法不需要统一图层属性 我觉得这一点很好。
我们可以看到图层【landuse】是栅格数据。

而图层【ClipPoly】是矢量数据

好啦，！我们接下来看一下怎么操作咯~

在这里插入代码片

Spatial Analyat Toola 
     →→→ Extraction 
         →→→ Extract by Msak

输出结果：系统会自动为其命名【Extract_land1】

（2）第二种方法：Raster Calculator （栅格计算器）

这个方法需要我们先统一数据类型，我们需要先将矢量数据转换成栅格数据。

Step 1：数据类型转换 → 从矢量到栅格

Conversion Tools
     →→→ To Raster
         →→→ Polygon to Raster

第3个框，输入像元大小是必须要做的，因为我们要做一个两个图层的叠加，如果像元大小不统一，那么就算转换了数据类型，做出的两个图层的叠加也是没什么用的。（语言表述可能不太准确）

这里请注意：像元大小要根据“被剪栅格数据”的大小为准。

像元大小可以通过属性查出来，然后输入在3框中。


当然也可以通过 3 框后方的文件夹，输入“被裁减的栅格图层”的路径，程序会自动提取该栅格图层的像元大小。

Step 2 ：为属性表添加一个 “ IDD ” 字段，并赋值为1（该字段可以自己命名）

Step 3：给进行栅格转换后的图层的属性表添加一个新的 “ IDD ” 字段，并赋值为 “ 1 ”

Step 4：使用 Raster Calculator（栅格计算器）

Spatial Analyst Tools
   →→→ Map Algebra
       →→→ Raster Calculator

将要叠加的两个图层相乘 在我的这个操作中，我的公式就是【“ClipPoly_PolygonToRaster1” * “landuse”】
同理，这个输出的结果，程序也会自动命名，也可以根据自己的想法修改。

得到结果：

----------------------------分隔线-------------------------------

3. 栅格重分类（Raster Reclassify)

Step 1：载入坡度栅格数据

Step 2：选择要重分类的数据，使用Raster Reclassify工具，并根据实际需要对其进行重新分类

Spatial Analyst
    →→→Reclass
        →→→Raster Reclassify

我们把坡度数据重新分为5类：0-8，8-15，15-25，35度以上，如上图所示。
分类结果如下，我们可以看到左方TOC中的【Reclass_slop2】已被分为五类。不同级别的颜色，程序会自动添加。

Step 3：查看结果

----------------------------分隔线-------------------------------

4. 查询符合条件的栅格，（提取符合条件的栅格）

（1）第一种方法：栅格计算——（Raster Calculator）

step 1：载入数据

step 2：使用Raster Calculator（栅格计算器）工具

Spatial Analyst Tools
   →→→ Map Algebra
       →→→ Raster Calculator

例如： 我想查询 坡度≥25° 的坡度数据， 【“slope1” >= 25】，也可以将输出的结果添加入自己的路径

（2）第二种方法：Less than Equal 工具