python gdal 基於shp檔案裁剪geotif圖並以最小外接矩形形式儲存

from osgeo import gdal
import os
import shapefile
#要裁剪的原圖
input_raster = r’raw.tif’
input_raster=gdal.Open(input_raster)

#shp檔案所在的資料夾
path=r’D:/shp_file/’

#裁剪結果儲存的資料夾
savepath=r’D:/result/’

#讀取shp檔案所在的資料夾
files= os.listdir(path)
for f in files: # 迴圈讀取路徑下的檔案並篩選輸出
if os.path.splitext(f)[1] == “.shp”:
name=os.path.splitext(f)[0]
input_shape=path+f
r = shapefile.Reader(input_shape)
output_raster=savepath+name+’.tif’
ds=gdal.Warp(output_raster,
input_raster,
format = ‘GTiff’,
outputBounds=r.bbox,
cutlineDSName = input_shape,
cutlineWhere=“FIELD = ‘whatever’”,
dstNodata = -1000)
ds=None

剛開始寫部落格，不大會用，發現程式碼間的縮排沒有了，那就直接再上傳張圖片吧，懶得改了。