如何在python打开shp文件，并用来掩膜

要在Python中打开shp文件并用来掩膜，可以使用`geopandas`库和`rasterio`库。以下是一个示例代码：

import geopandas as gpd
import rasterio.mask

# 打开shp文件
mask = gpd.read_file('path/to/shapefile.shp')

# 打开要掩膜的栅格图像
with rasterio.open('path/to/raster.tif') as src:
    # 将shp文件的几何图形转换为栅格图像的坐标系
    mask_geojson = mask.to_crs(src.crs.data)
    # 掩膜栅格图像
    masked_image, masked_transform = rasterio.mask.mask(src, mask_geojson.geometry, crop=True)
    # 将掩膜后的图像保存到文件
    with rasterio.open('path/to/masked_raster.tif', 'w', driver='GTiff', width=masked_image.shape[2], height=masked_image.shape[1], count=masked_image.shape[0], dtype=masked_image.dtype, crs=src.crs, transform=masked_transform) as dst:
        dst.write(masked_image)

在这个示例中，`geopandas`库用于打开shp文件，并将其几何图形转换为栅格图像的坐标系。`rasterio`库用于打开栅格图像、掩膜和保存掩膜后的图像。请注意，这个示例假设shp文件只包含一个几何图形，如果shp文件包含多个几何图形，则需要进行适当的更改。