如何光栅化矢量数据线并在 python 中沿线插入值？

如何解决如何光栅化矢量数据线并在 python 中沿线插入值？

我想光栅化一个线向量。必须在“from”和“to”值之间插入栅格值。

示例在此处。我想栅格化 het geodataframe 中的线，栅格的值应该从 5 到 2（插值）。

import geopandas as gpd
from shapely.geometry import Linestring

# example line
line = Linestring([(10.00,20.00),(20.00,50.00)])
line_gdf = gpd.GeoDataFrame([line],geometry=[line])
line_gdf['from'] = 5
line_gdf['to'] = 2

预期的结果是这样的： example raster

该线以黑色可见。栅格应该只有与像元相交的线的值。这些值应该从“起始值”插入到“目标值”，如示例栅格（用红色轮廓标记的单元格）所示。最好是 DataArray (xarray)。

解决方法

您需要使用 Bresenham 算法： https://en.wikipedia.org/wiki/Bresenham_line_algorithm

您会在此处找到 Python 实现： https://pypi.org/project/bresenham/