使用ARCPY分析年度变更调查地类变化情况，比如变化年份，变化地类等

使用ARCPY ( arcpy ) 来分析年度变更调查地类的变化情况，通常涉及以下几个步骤：

1. **加载库和数据**:
- 导入arcpy模块并加载年度变更调查的数据集，可能包含多个年份的地类图层。

import arcpy

# 加载数据
base_layer_1 = r"C:\path\to\base_year_1.shp" # 第一年地类图层
base_layer_2 = r"C:\path\to\base_year_2.shp" # 第二年地类图层
change_layer = r"C:\path\to\difference_layer.shp" # 变化地类图层（如果有）

2. **创建临时图层**:
- 创建临时图层以方便处理和分析。

temp_layer = "in_memory/temp"
arcpy.CopyFeatures_management(base_layer_1, temp_layer)

3. **差异分析**:
- 使用`arcpy.sa.Difference`函数计算地类变化。

if change_layer:
    arcpy.sa_difference(temp_layer, base_layer_2, out_raster="in_memory/diff_results")
else:
    arcpy.sa_difference(base_layer_1, base_layer_2, out_raster="in_memory/diff_results")

4. **提取属性**:
- 通过查询工具或字段计算器，获取变化地类的属性，包括变化年份、变化类型等信息。

with arcpy.da.SearchCursor("in_memory/diff_results", ["SHAPE@", "YEAR"]) as cursor:
    for row in cursor:
        change_shape, change_year = row
        print(f"变化地类: {change_shape}, 变化年份: {change_year}")

5. **汇总结果**:
- 根据需要，你可以对这些数据进行汇总统计，比如按年份划分变化总数。

result_stats = arcpy.Statistics_analysis("in_memory/diff_results", "in_memory/statistics_table", [["YEAR", "SUM"]])