土地利用情况分析（土地利用指数）

一、土地利用程度综合指数简介

  土地利用程度综合指数是针对一个具体年份内所有土地利用类型整体反映的土地利用集约化程度。土地利用程度综合指数通过对土地利用程度的分级,量化人类活动对于土地系统的影响程度,定量地描述该地区土地利用的综合水平和变化趋势，是衡量指定区域内土地利用情况的一个重要指标。

  那在实际中我们如何应用土地利用程度综合指数来进行研究呢？

二、数据简介

  本例子使用的数据为2020年山东省30米空间分辨率的土地利用数据、2020年山东省行政区划数据，数据来源于中国国家基础地理信息中心提供的全球30米地表覆盖数据（GlobeLand30）。

三、方法流程

（1）栅格转矢量
  由于GlobeLand30提供的是tif格式的栅格数据，后续需要进行矢量分割、属性计算、分级等操作，因此我们首先需要在ArcGIS里面进行栅格转矢量的操作，将栅格数据转化为面要素数据。这里我们使用ArcMap中的工具Raster to Polygon来实现。

（2）土地利用情况分级
  接下来我们需要将转化后的面要素图层进行数据分级，参考GlobeLand30产品中给出的分类系统信息，然后按照我们自己的研究需要进行级别划分，本例子中划分为四个级别：
1：冰川和永久积雪、裸地、无数据
2：林地、草地、灌木地、湿地、水体、苔原
3：耕地
4：人造地表
  注意：所用的数据中不仅仅有10至100等十种数值，还有个别255的数值，这种一般是无数据的情况。

  为新的面要素图层添加新的一个属性Class，并对应赋值1、2、3、4四个等级数。这里我们使用ArcGIS中的字段计算器来实现。

  在字段计算器中，我们可以编写预逻辑代码，以此来依据旧的属性情况，对新的属性进行赋值。解析程序可以使用Python或者VB脚本，本文将分别提供这两种方式的具体使用方法。
  ①VB脚本：如果解析程序选择VB脚本的话，则不需要预逻辑代码。首先我们依据上述的四个等级，分别筛选（右击图层，打开属性表，在属性表内即可进行排序筛选等操作），例如先筛选等级1的要素出来，然后直接在字段计算器中对其Class属性进行赋值，Class=1即可，以此类推，完成赋值。这样的方式相对写一个VB预逻辑函数进行筛选来说要较为简单。
  ②Python：如果解析程序选择Python，那就更简单了，我们可以在预逻辑代码里面写一个函数来实现我们的需要。

def level(GRIDCODE):
    if GRIDCODE >=90:
        return 1
    elif GRIDCODE >= 20 and GRIDCODE <= 70:
        return 2
    elif GRIDCODE == 10:
        return 3
    else:
        return 4

  然后在字段计算器输出结果那里直接使用Class=level(!GRIDCODE!)。注意，GRIDCODE这个变量前后必须有感叹号，这是ArcGIS的特殊要求，可以通过双击字段计算器中的字段框内的字段，来实现使用，它会自动按照正确格式输入的。
  同时注意，在使用字段计算器时候，上述图中红框内都是需要注意的选项，根据自己操作时候的实际情况来进行选择，如果属性类别选错了，比如创建的是字符串类型，计算器输入时候选的是数字类型，那么之后计算时候就会报错，错误内容如下：
There was a failure during processing,check the Geoprocessing Results window for detail

  下面是字段计算器计算后的结果：

（3）面要素分割
  由于所用的土地利用数据中会存在个别类型占地较大，从而形成较大的图斑，出现横跨多个区县的情况，对于之后进行各区县的情况分析不利，因此我们需要对转化后的面要素图层进行分割，按照山东曲线的矢量边界数据来进行分割。这里使用ArcMap中的Split工具，按照山东省各区划来分割，确保数据中跨区划的图斑分割至各区。

（4）合并
  为了之后统计、分析数据的方便，我们在分割之后，还需要将所有文件合并为一个新的面要素图层。这里使用ArcMap工具箱中的Merge工具。

（5）空间链接
  这一步是将合并后的新图层与区县边界数据进行空间链接，将分割并合并后的土地利用数据与区县数据进行链接，使其依照位置关系赋值到相同地点。
（6）属性表导出Excel
  接下来需要将赋值好的各个区县的新属性数据进行导出。
  ArcMap中提供了Table To Excel工具，可以实现该目的，但是当数据量过大的时候，就会出现报错情况，如下所示。Excel表格中允许的最大行数为65535，当属性表中数据量超过时，就无法使用该工具直接导出Excel文件了。

  那面对这种情况时，我们可以换个方法，不使用ArcMap中的工具，绕路而行。分三步走，首先，把属性表转存至txt文件，这一步通过属性表上方菜单中的Export即可实现。

  第二步是将txt文件转存为csv文件。使用记事本打开txt文件，然后另存为csv文件，值得注意的是，在输出csv文件时，需要将编码方式设置为ANSI，而不能像往常一样使用UTF-8的编码方式，否则文件中的中文会出现乱码。

  最后将csv文件另存为Excel文件即可。
（7）Excel数据可视化
  Excel作为一个数据存储与处理的高级软件，妙用无穷。本文就将在Excel中对导出数据进行可视化的举例。数据透视图是Excel中常用的一种可视化方式，通过插入图、筛选条件等步骤，将我们需要的数据从上万行数据中提取出来，统计绘图。

（8）Excel中进行土地利用程度综合指数计算
  这一步按照本文中提到的土地利用程度综合指数公式，在Excel中进行计算。首先，我们需要将划分后的各个类分别制作一个sheet存放，这里仍使用数据透视表筛选，例如统计出各区县的1类土地面积，然后将筛选后数据存放好。同时，也需要统计各个区县的土地总面积，单独存放于一个sheet中。然后在新的sheet中进行所有数据的汇总。
  说到这里，大家可能就会疑惑，为什么不一开始就干脆存放在一个sheet中呢？
  其实是这样的，个别区县可能并不存在某一类土地，例如博兴县就没有1类土地，也就是它的1类土地面积为0。但是Excel的数据透视表在进行数据统计筛选时，并不会将面积为0的项列出，也就是说在不同类型土地时，可能筛选结果项数就不一样了，比如1类土地筛选得到100个有面积的区县，2类土地筛选得到130个有面积的区县。因此，如果一开始就用一个sheet存放所有数据，那么就会导致数据项匹配错误。因此，我们采用多个sheet进行数据存放。
  那么在用于汇总的那个sheet当中，我们又如何匹配复制列数据，使得每一个区县都有四种类型土地的数据并且正确赋值呢？可以使用如下函数进行筛选计算。

IFERROR(VLOOKUP(A2,Class1!$A$1:$B$80,2,),0)

  VLOOKUP函数用于遍历查找指定区域内指定项的指定值；
  $这个字符主要表示绝对引用，使用之后，在下拉单元格进行批量赋值时，该字符后的单元格位置不会改变。使用它能够保证我们的函数在批量赋值时候参数不发生变化。
  IFERROR函数：当函数内的公式返回错误值（#N/A）时，将单元格赋值为0，而不是继续显示错误值。

  完成上述步骤后，用各个类型的统计面积除以总面积，最后再计算土地利用程度综合指数即可。

100*SUM(G2*1,H2*2,I2*3,J2*4)

（9）将计算得到的指数作为新的属性字段连接到ArcGIS中区县数据图层
  右击该图层，选择Join即可。

（10）分级显示
  最后只需要在该图层中调整属性，选择分级色彩显示即可。