生成落入多边形内的经纬度坐标格网

# 导入必要的库
import numpy as np
from shapely.geometry import Point, Polygon

# 定义多边形区域的顶点坐标
polygon_coords = [(30, 10), (40, 40), (20, 40), (10, 20)]

# 定义格网的经纬度间隔
grid_interval = 1

# 定义格网范围
min_lon = min(polygon_coords, key=lambda x: x[0])[0]
max_lon = max(polygon_coords, key=lambda x: x[0])[0]
min_lat = min(polygon_coords, key=lambda x: x[1])[1]
max_lat = max(polygon_coords, key=lambda x: x[1])[1]

# 划分格网
lon_range = np.arange(min_lon, max_lon, grid_interval)
lat_range = np.arange(min_lat, max_lat, grid_interval)

# 判断格网是否在多边形内
polygon = Polygon(polygon_coords)
grid_points = []
for lon in lon_range:
    for lat in lat_range:
        point = Point(lon, lat)
        if polygon.contains(point):
            grid_points.append((lon, lat))

# 打印落入多边形内的经纬度坐标格网
for point in grid_points:
    print(point)