arcpy实现LUCC土地利用变化转移矩阵,2期影像或者矢量数据土地利用变化统计,生成统计表。使用“面积制表TabulateArea”轻松实现变化统计。
coding = 'utf-8' #utf8编码
import os #加载os库
import arcpy #加载arcpy库
from arcpy import env #导入环境变量
from arcpy.sa import * #导入Spatial Analyst模块
env.workspace = r"D:\Ex1" #设置工作空间为数据文件夹!注意是文件夹!(r字母代表转义字符)
#设置变量参数
inZoneData = "zonedata.shp" #土地利用数据1期
zoneField = "IDStr" #区域字段
inClassData = "valueraster.shp" #土地利用数据2期
classField = "VALUE" #类字段
outTable = r"D:\Ex1\areatable01.dbf" #输出表路径
arcpy.CheckOutExtension("Spatial") # 检查Spatial Analyst服务开启
print("--------开始Start!--------") #打印开始日志
# 开始计算土地利用转移矩阵
TabulateArea(inZoneData, zoneField, inClassData, classField, outTable)
print("--------完成Success!--------") #打印完成日志