网上有很多教程,但是,一个简单的裁剪。。。搞得那么复杂,这儿我说下一般两种方法(下面的代码都可以直接copy使用):1)使用矢量文件生成掩膜。用掩膜裁剪TIFF(使用掩膜的方法谁愿意补充我可以加上);2)直接调用接口。下面且听我分说。
直接调用,你可以直接copy使用
from osgeo import gdal
input_shape = r"C:/Users/dad/Desktop/GDAL/area/area_WGS84.shp"
output_raster=r'C:UsersdadDesktopGDALtest2.tif'
# tif输入路径,打开文件
input_raster = r"C:/Users/dad/Desktop/GDAL/Landsat8.tif"
# 矢量文件路径,打开矢量文件
input_raster=gdal.Open(input_raster)
# 开始裁剪,一行代码,爽的飞起
ds = gdal.Warp(output_raster,
input_raster,
format = 'GTiff',
cutlineDSName = input_shape,
cutlineWhere="FIELD = 'whatever'",
dstNodata = 0)
# 关闭文件
ds=None
使用掩膜这种方法:这种方法我仅仅是裁剪了最小外接矩形,这还算不上掩膜&#