如何在python中将雷达图像与shapefile重叠(mach)和裁剪

import geopandas as gpd
import rasterio
from rasterio.mask import mask
from shapely.geometry import mapping

shapefile = 'path/to/shapefile.shp'
gdf = gpd.read_file(shapefile)

image_file = 'path/to/radar_image.tif'
raster = rasterio.open(image_file)

geoms = gdf.geometry.values
geometry = mapping(geoms[0])

out_image, out_transform = mask(raster, [geometry], crop=True)

out_meta = raster.meta.copy()
out_meta.update({"driver": "GTiff",
                 "height": out_image.shape[1],
                 "width": out_image.shape[2],
                 "transform": out_transform})
output_file = 'path/to/output_image.tif'
with rasterio.open(output_file, "w", **out_meta) as dest:
    dest.write(out_image)