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)