用Python写一个等分切割矢量多边形为100米*100米的小多边形的脚本

要将矢量多边形等分切割为 100 米 * 100 米的小多边形，需要考虑以下步骤：

1. 将矢量多边形转换为 Shapely 的 Polygon 对象。你可以使用 Shapely 库来处理矢量几何对象。

2. 计算矢量多边形的边界框（bounding box），得到矢量多边形的最小外接矩形。

3. 根据最小外接矩形的范围确定切割的网格大小。根据需求，这里将网格大小设为 100 米 * 100 米。

4. 创建一个网格，将最小外接矩形划分为等大小的单元格。可以使用 numpy 来创建网格。

5. 对于每个单元格，检查其是否与多边形相交。可以使用 Shapely 的 intersects() 方法来判断是否相交。

6. 如果单元格与多边形相交，则将相交区域作为一个小多边形保存起来。

下面是一个示例脚本，用于等分切割矢量多边形为 100 米 * 100 米的小多边形：

import numpy as np
from shapely.geometry import Polygon

def split_polygon(polygon, grid_size):
    min_x, min_y, max_x, max_y = polygon.bounds
    
    x_steps = int((max_x - min_x) / grid_size)
    y_steps = int((max_y - min_y) / grid_size)
    
    polygons = []
    
    for i in range(x_steps):
        for j in range(y_steps):
            x_min = min_x + i * grid_size
            y_min = min_y + j * grid_size
            x_max = min_x + (i + 1) * grid_size
            y_max = min_y + (j + 1) * grid_size
            
            cell_polygon = Polygon([(x_min, y_min), (x_min, y_max), (x_max, y_max), (x_max, y_min)])
            
            if polygon.intersects(cell_polygon):
                intersection = polygon.intersection(cell_polygon)
                
                if intersection.area > 0:
                    polygons.append(intersection)
    
    return polygons

# 示例用法
polygon = Polygon([(0, 0), (0, 5), (5, 5), (5, 0)])  # 定义一个矩形多边形
grid_size = 100  # 设置网格大小为 100 米

result = split_polygon(polygon, grid_size)

for sub_polygon in result:
    print(sub_polygon)

这个示例中，我们使用了 Shapely 库来处理矢量几何对象。`split_polygon` 函数接受一个多边形和网格大小作为输入，并返回切割后的小多边形列表。示例用法中的多边形是一个矩形，网格大小设置为 100 米。程序将切割结果打印出来。

请注意，这个示例假设输入的多边形是平面几何对象，并且没有考虑地球曲率等因素。如果需要考虑地理空间坐标系，可以使用适当的地理空间库（如 GeoPandas）来处理坐标转换和距离计算。