rasterio 切割矢量

要使用 rasterio 库来切割矢量数据，您可以按照以下步骤进行操作：

1. 导入所需的库：

import rasterio
from rasterio.mask import mask
import geopandas as gpd

2. 加载矢量数据和栅格数据：

# 加载矢量数据
vector_data = gpd.read_file('path/to/vector_file.shp')

# 加载栅格数据
with rasterio.open('path/to/raster_file.tif') as src:
    raster_data = src.read(1)  # 读取第一个波段（如果有多个波段，可以选择其他波段）
    raster_meta = src.meta.copy()  # 复制栅格元数据

3. 对矢量数据进行裁剪：

# 将矢量数据转换为栅格坐标系
vector_data = vector_data.to_crs(raster_meta['crs'])

# 将矢量数据转换为GeoJSON格式
vector_geojson = vector_data.geometry.__geo_interface__

# 利用 rasterio 库的 mask 方法进行裁剪
masked_data, masked_transform = mask(dataset=raster_data, shapes=vector_geojson, crop=True)

在上述代码中，我们首先将矢量数据转换为与栅格数据相同的坐标系，然后将其转换为 GeoJSON 格式。然后，我们使用 rasterio 库的 `mask` 方法将栅格数据裁剪为与矢量数据重叠的部分。`masked_data` 变量将包含裁剪后的栅格数据，`masked_transform` 变量将包含裁剪后的栅格的变换信息。

请确保安装了 rasterio 和 geopandas 库，并将 "path/to/vector_file.shp" 和 "path/to/raster_file.tif" 替换为您实际的矢量和栅格数据文件路径。

希望能对您有所帮助！如果还有其他问题，请随时提问。