2022年 11月 13日

Python(arcpy) 重采样

任务描述:有高精度的数据,想要重采样降低精度以便跟其它数据匹配

import arcpy
from arcpy import env     #导入 env 类(env 类包含所有地理处理环境)
#定义工作环境
arcpy.env.workspace = "F:\\"   

#重采样时的捕捉栅格(定义环境),这一句可要可不要,取决于你有没有需要对齐的目标栅格                    
arcpy.env.snapRaster = "F:\\sample_std.tif"   

#输入栅格
in_raster = "F:\\sample_in.tif"

#输出栅格
out_raster = "F:\\sample_out.tif"

#获取重采样时对标栅格的像元大小
pixel_size = arcpy.GetRasterProperties_management(arcpy.env.snapRaster,"CELLSIZEX")           
print(pixel_size)

#重采样
arcpy.Resample_management (in_raster, out_raster,pixel_size)
  • 1
  • 2
  • 3
  • 4
  • 5
  • 6
  • 7
  • 8
  • 9
  • 10
  • 11
  • 12
  • 13
  • 14
  • 15
  • 16
  • 17
  • 18
  • 19
  • 20

总结:重采样涉及snapRaster 工具、GetRasterProperties_management工具和Resample_management工具,详细介绍请阅读
捕捉栅格
获取栅格属性
重采样

I hope this brings you some insight!