土地覆被变化制图与土地利用转移矩阵

土地利用数据来源

1. 自己解译

哨兵数据易康面向对象SVM分类

eCognition（易康）面向对象分类与分类结果处理

1. 网上各种土地覆被数据

• 清华宫鹏课题组：http://data.ess.tsinghua.edu.cn/
• 中科院资源环境科学与数据中心：http://www.resdc.cn/data.aspx?DATAID=99
• GlobalLand30: http://www.webmap.cn/commres.do?method=globeIndex
• MCD12Q1：https://e4ftl01.cr.usgs.gov/MOTA/MCD12Q1.006/
• 国家地球系统科学数据中心土地覆被专题库：http://www.geodata.cn/data/index.html?categoryId=18

数据源与制图目的

已有数据

已有2010，2018年两个年份的矢量土地覆被图

图1. 黄河口2010，2018主要植物群落和土地覆被分布

分析结果

根据上面的土地覆被图，绘制下面的侵蚀变化图，植物群落与地类变化图。制作土地利用转移矩阵。

图2. 三角洲侵蚀变化图

图3. 黄河三角洲植物群落与地类变化图

表1. 黄河口群落与土地利用转移矩阵

三角洲侵蚀变化图制作

总体思路

• 提取水陆边界
• 提取2010，2018水陆边界交集
• 2010，2018土地覆被求并集
• 土地覆被并集-水陆边界交集，求侵蚀变化情况

具体做法

1. 根据图1绘制水陆分界

图4. 水陆分界图

1. 将图4两个矢量求交集，求得2010，2018均为陆地的区域（输出要素类：OLI水陆分界_Intersect）
2. 将图1两个矢量求交集（输出要素类：TM2010_Intersect1）
3. 在TM2010_Intersect1减去OLI水陆分界_Intersect（输出要素类：TM2010_Intersect1_Erase）
4. 新建一个字段，计算侵蚀淤积变化情况

使用python进行计算

需要注意的是，我前面用的是汉字，下面代码里面判断的是地类编码，由于ArcGIS中Python2对于汉字支持不太好，所以还是建议判断数字。

def Cal(a,b):
  if a==b:
    return "无变化"
  elif a ==1 and b !=a:
    return "淤积"
  else:
    return "侵蚀"
 
 
Cal(  !CODE! , !CODE_1! )

1. 最后对Change字段进行分类符号化即可。

变化制图

黄河三角洲植物群落与地类变化图

以水域为例，介绍地类变化图制作

• 按属性选择，筛选2010，2018年均为水域的所有要素，并导出要素类
• 新建一个Change字段（文本型），计算变化情况
• 计算代码和上面侵蚀淤积代码类似
• 转入转出分类符号化

根据转入转出分类符号化

土地利用转移矩阵

• 将前面TM2010_Intersect1要素类属性表导出

属性表导出为文本型

• NOTEPAD++改一下ANSI编码，防止EXCEL打开乱码
• 数据透视表制作转移矩阵，行2010年地类，列2018年地类，求和项为面积

制作数据透视表

• 由于数据使用的是投影坐标系，单位为米，对应的面积单位为平方米，因此土地利用转移矩阵需要进一步进行单位换算，换算为所需要的单位。