def clip_raster(in_raster, out_raster, mask_shp):
"""
:param in_raster: 輸入柵格
:param out_raster: 輸出柵格
:param mask_shp: 裁剪矢量
:param wkid: wkid
:return:
"""
gdal.Warp(out_raster,
in_raster,
format='GTiff',
dstSRS='EPSG:4326',
cutlineDSName=mask_shp,
cropToCutline=True, # 按掩膜圖層范圍裁剪
dstNodata=-9999,
outputType=gdal.GDT_Float64)