第一次用python来处理栅格数据,虽然ArcGIS也可以裁剪,但是如果数据比较大的话,ArcGIS就可能会崩掉,所以尝试用python来处理数据。从网上看了一下别人的代码,虽然代码很简单,但是很有用。
from osgeo import gdal # 如果python中没有安装GDAL库,会报错,需要安装GDAL库
input_shape = r"E:/temp/Export_Output.shp"
output_raster=r'E:/code/temporary_clip/test1.tif'
# tif输入路径,打开文件
input_raster = r"E:/code/Hmp.tif"
# 矢量文件路径,打开矢量文件
input_raster=gdal.Open(input_raster)
# 开始裁剪
ds = gdal.Warp(output_raster,
input_raster,
format = 'GTiff',#设置格式
cutlineDSName = input_shape,
cropToCutline=True, # 保证裁剪后影像大小跟矢量文件的范围大小一致(当设置为False,结果大小会跟待裁剪影像大小一样,会出现大量的空值区域)
dstNodata = 0)