arcpy将1km精度栅格数据重采样到30m精度栅格数据

使用Arcpy中的Resample工具可以将1km精度的降水栅格数据降尺度到30m精度的栅格数据。以下是一个简单代码:


import arcpy

# 设置输入数据路径
input_raster = r"C:\data\precipitation_1km.tif"  # 1km精度的降水栅格数据

# 设置输出数据路径和空间参考
output_raster = r"C:\data\precipitation_30m.tif"  # 30m精度的栅格数据
output_cell_size = 30  # 输出栅格单元的大小(单位:米)

# 进行栅格数据降尺度
arcpy.Resample_management(input_raster, output_raster, str(output_cell_size), "NEAREST")

# 更新输出栅格数据的空间参考信息
arcpy.DefineProjection_management(output_raster, arcpy.Describe(input_raster).spatialReference)
 

在上述代码中,首先设置了输入数据的路径,即1km精度的降水栅格数据。然后,设置了输出数据的路径和栅格单元的大小(30m)。接下来,使用`Resample_management`函数将输入栅格数据进行降尺度处理,并保存为输出栅格数据。最后,通过使用`DefineProjection_management`函数更新输出栅格数据的空间参考信息,保持与输入栅格数据一致。

请注意,代码中的路径和文件名应根据你自己的数据进行相应修改。此外,`"NEAREST"`是插值方法参数,表示使用最近邻插值方法来进行降尺度处理。如果你需要使用其他的插值方法(如双线性插值),可以相应修改该参数。

猜你喜欢

转载自blog.csdn.net/weixin_58851039/article/details/131240166