Terra Climate 数据预处理【Arcpy 代码分享】

twelve127

已于 2024-10-29 14:33:01 修改

阅读量316

点赞数 3

文章标签： arcgis python

于 2024-10-29 14:31:38 首次发布

本文链接：https://blog.youkuaiyun.com/twelve127/article/details/143328402

版权

具体步骤有按月份提取每个nc文件的数据，根据研究区按掩膜提取，对提取的数据进行重新投影（因为后续分析步骤需要等面积投影），重新投影后按照双线性插值法重采样。

import arcpy
import os

# 启用 Spatial Analyst 扩展
arcpy.CheckOutExtension("Spatial")

# 设置基础输入和输出路径
base_input_folder = "D:/Terra Climate data"  # 基础输入路径
base_output_folder = "D:/Terra Climate data/output_files"  # 基础输出文件夹路径
mask_file = "D:/Terra Climate data/mask/roi.shp"  # 掩膜文件路径（研究区矢量文件）
target_projection = arcpy.SpatialReference(102013)  # 目标投影为 Europe_Albers_Equal_Area_Conic 根据自己研究选择 这个102013代码可以在arcgis里用投影栅格选投影的时候看到 每个投影代码不一样

# 定义要处理的变量和相应文件夹名称，确保按处理顺序排列 变量名称一般就是下载数据的名称 具体信息也可以用panoply软件去看
variables = [
    ("ppt", "1_ppt"),
    ("soil", "2_soil"),
    ("srad", "3_srad"),
    ("vpd", "4_vpd")
]

# 遍历每个变量对应的文件夹
for variable, folder_name in variables:
    input_folder = os.path.join(base_