arcgis使用arcpy计算土地利用LUCC变化转移矩阵、马尔科夫模型变化表,最好用的方法TabulateArea实现两期数据变化统计:
coding = 'utf-8' #utf8编码
import os #加载os库
import arcpy #加载arcpy库
from arcpy import env #导入环境变量
from arcpy.sa import * #导入Spatial Analyst模块
env.workspace = r"C:\Users\HP\Desktop\LUCC\test.gdb" #设置工作空间为数据文件夹!注意是文件夹!(r字母代表转义字符)
#设置变量参数
inZoneData = "T2008" #土地利用数据1期
zoneField = "LAND_USE_G" #区域字段
inClassData = "T2013" #土地利用数据2期
classField = "LAND_USE_G" #类字段
outTable = "ChangeTable_" + inZoneData + "_" + inClassData #输出表路径
arcpy.CheckOutExtension("Spatial") # 检查Spatial Analyst服务开启
print("--------开始Start!--------") #打印开始日志
# 开始计算土地利用转移矩阵
TabulateArea(inZoneData, zoneField, inClassData, classField, outTable)
print("--------完成Success!--------") #打印完成日志