土地覆被变化制图与土地利用转移矩阵
土地利用数据来源
- 自己解译
哨兵数据易康面向对象SVM分类
eCognition（易康）面向对象分类与分类结果处理
- 网上各种土地覆被数据
- 清华宫鹏课题组：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绘制水陆分界
【土地覆被变化制图与土地利用转移矩阵】图4. 水陆分界图
- 将图4两个矢量求交集，求得2010，2018均为陆地的区域（输出要素类：OLI水陆分界_Intersect）
- 将图1两个矢量求交集（输出要素类：TM2010_Intersect1）
- 在TM2010_Intersect1减去OLI水陆分界_Intersect（输出要素类：TM2010_Intersect1_Erase）
- 新建一个字段，计算侵蚀淤积变化情况
使用python进行计算
需要注意的是，我前面用的是汉字，下面代码里面判断的是地类编码 ， 由于ArcGIS中Python2对于汉字支持不太好，所以还是建议判断数字 。
def Cal(a,b):if a==b:return "无变化"elif a ==1 and b !=a:return "淤积"else:return "侵蚀"Cal(!CODE! , !CODE_1! )- 最后对Change字段进行分类符号化即可 。
变化制图
黄河三角洲植物群落与地类变化图
以水域为例，介绍地类变化图制作- 按属性选择，筛选2010，2018年均为水域的所有要素 ， 并导出要素类
- 新建一个Change字段（文本型），计算变化情况
- 计算代码和上面侵蚀淤积代码类似
- 转入转出分类符号化
根据转入转出分类符号化
土地利用转移矩阵
- 将前面TM2010_Intersect1要素类属性表导出
属性表导出为文本型
- NOTEPAD改一下ANSI编码，防止EXCEL打开乱码
- 数据透视表制作转移矩阵 ， 行2010年地类，列2018年地类，求和项为面积
制作数据透视表
- 由于数据使用的是投影坐标系，单位为米 ， 对应的面积单位为平方米，因此土地利用转移矩阵需要进一步进行单位换算，换算为所需要的单位 。
