shapely 对一个多边形生成外包矩形，并对矩形进行投影，把矩形划分成100米*100米的多边形

要使用Shapely对一个多边形生成外包矩形，并对矩形进行投影，并将矩形划分成100米*100米的多边形，你可以使用GeoPandas、Shapely和pyproj库。下面是一个示例代码，展示了如何实现这个过程：

import geopandas as gpd
from shapely.geometry import Polygon
from shapely.ops import transform
from pyproj import Transformer

# 定义原始多边形
polygon = Polygon([(0, 0), (0, 10), (10, 10), (10, 0)])

# 计算多边形的外包矩形
bounding_box = polygon.bounds

# 定义原始坐标系（EPSG:4326）和目标坐标系（EPSG:3857）
src_crs = "EPSG:4326"
dst_crs = "EPSG:3857"

# 创建一个转换器，用于进行投影转换
transformer = Transformer.from_crs(src_crs, dst_crs)

# 对外包矩形进行投影转换
projected_bounding_box = transform(transformer.transform, bounding_box)

# 计算投影后的矩形的边界值
x_min, y_min, x_max, y_max = projected_bounding_box

# 计算每个小格子的大小（100米*100米）
grid_size = 100

# 创建一个空的GeoDataFrame来存储划分后的多边形
grid = gpd.GeoDataFrame()

# 划分矩形成小格子
for x in range(int(x_min), int(x_max), grid_size):
    for y in range(int(y_min), int(y_max), grid_size):
        poly = Polygon([(x, y), (x+grid_size, y), (x+grid_size, y+grid_size), (x, y+grid_size)])
        # 如果小格子与原始多边形相交，则将其添加到GeoDataFrame中
        if poly.intersects(polygon):
            grid = grid.append({'geometry': poly}, ignore_index=True)

# 设置GeoDataFrame的坐标系为目标坐标系
grid.crs = dst_crs

# 保存划分后的多边形为Shapefile文件
grid.to_file("grid.shp")

这段代码会对给定的多边形生成外包矩形，并将矩形投影到目标坐标系中。然后，它会将矩形划分成100米*100米的小格子，并将与原始多边形相交的小格子保存为一个Shapefile文件。请确保你已经安装了pyproj、geopandas和shapely库。你可以根据需要调整多边形的顶点坐标、目标坐标系和小格子的大小。