arcpy数据处理1 批量统计每年的区域均值【自用代码备份】

原创

已于 2024-11-27 20:22:59 修改 · 322 阅读

0 ·

CC 4.0 BY-SA版权

文章标签：

#arcgis

于 2024-11-27 16:42:30 首次发布

要画GPP NPP CUE年均值变化图要把数据投影变成等面积投影再用以表格显示分区统计最后把表转为excel Arcpy代码如下【暂时还不会批量读取表格数据汇总】

import arcpy
from arcpy import env
from arcpy.sa import *
import os

# 设置工作环境
env.overwriteOutput = True

# 定义文件夹路径
tif_folder = "路径/"  # tif 文件所在文件夹
border_layer = "路径/pro_xiboliya_border.shp"  # 分区边界文件
output_folder = "路径/sta_npp_results/"  # 输出文件夹
projection = arcpy.SpatialReference(102013)  # #后来用了102025 Asia_North_Albers_Equal_Area_Conic  投影坐标系，Europe_Albers_Equal_Area_Conic 的 EPSG 代码

# 定义函数：将 DBF 文件转换为 Excel
def dbf_to_excel(dbf_path):
    excel_path = "{}.xls".format(dbf_path.replace(".dbf", ""))
    arcpy.conversion.TableToExcel(dbf_path, excel_path)
    print("已将 {} 转换为 {}".format(dbf_path, excel_path))
    return excel_path

# 获取 tif 文件夹下的所有 .tif 文件
tif_files = [f for f in os.listdir(tif_folder) if f.endswith(".tif")]
print("找到的 tif 文件: {}".format(tif_files))

# 处理每个 tif 文件
for tif_file in tif_files:
    # 获取完整路径和文件名（不含扩展名）
    tif_path = os.path.j