在GeoPandas中将纬度/经度点转换为栅格面

import geopandas as gpd
from shapely.geometry import Point
from rasterio.transform import from_origin
from rasterio.enums import Resampling

# 创建一个包含纬度/经度点的GeoDataFrame
points = gpd.GeoDataFrame(geometry=[Point(lon, lat)])

# 将纬度/经度点转换为栅格坐标
points['x'] = points.geometry.apply(lambda p: p.x)
points['y'] = points.geometry.apply(lambda p: p.y)

# 定义栅格的分辨率和大小
resolution = 0.01  # 栅格分辨率
width = 100  # 栅格宽度
height = 100  # 栅格高度

# 计算栅格的左上角坐标
left = points['x'].min() - resolution / 2
top = points['y'].max() + resolution / 2

# 创建栅格数据集
transform = from_origin(left, top, resolution, resolution)
grid = gpd.GeoDataFrame(geometry=[Point(x, y) for x in range(width) for y in range(height)], crs=points.crs)
grid['x'] = grid.geometry.apply(lambda p: p.x)
grid['y'] = grid.geometry.apply(lambda p: p.y)

# 将栅格数据集与纬度/经度点进行空间连接
joined = gpd.sjoin(grid, points, how='left', op='contains')

# 将连接结果转换为栅格面
raster = joined['geometry'].notnull().values.reshape((height, width)).astype(int)