如何基于Shapely MULTILINESTRING正确使用PyDeck PathLayer？

我查看了PathLayer示例，发现其输入数据是普通Pandas DataFrame中的列表嵌套结构，而非GeoPandas DataFrame。我有一个名为consline的GeoPandas DataFrame，其中包含Shapely MULTILINESTRING类型的geometry列，尝试了PathLayer与GeoJsonLayer结合get_path参数的多种组合均无效，仅通过json.loads(consline.geometry.to_json())['features'][0]['geometry']['coordinates']这种先转JSON再转字典的方式成功实现，但这种方式需要PyDeck再转回JSON给deck.gl，显然不是最优解。

consline可通过以下代码加载：

import geopandas as gpd
consline = gpd.read_file("https://github.com/deanm0000/SOexamples/raw/main/consline.gpkg")

无效的尝试代码 #

直接将GeoPandas DataFrame传入PathLayer并指定get_path="geometry"无法正常工作：

import pydeck as pdk

mylayers=[
    pdk.Layer(
        'PathLayer',
        data=consline,
        get_path="geometry",
        pickable=True,
        get_fill_color=[0,255,0],
        width_scale=20,
        width_min_pixels=2,   
    )
]
view_state = pdk.ViewState(
    longitude=-98.99,
    latitude=31.79,
    zoom=5,
    min_zoom=1,
    max_zoom=15)
r = pdk.Deck(layers=mylayers, initial_view_state=view_state, map_style='light')
r.to_html("example.html")

最优解决方案 #

不需要绕JSON转换，直接从Shapely的MULTILINESTRING中提取坐标列表，生成适配PathLayer的格式：

1. 先为GeoPandas DataFrame添加一列，存储PathLayer所需的嵌套坐标列表：

consline['path_coords'] = consline.geometry.apply(
    lambda geom: [list(line.coords) for line in geom.geoms]
)

这里利用Shapely的geoms属性遍历MULTILINESTRING中的每一条LineString，再将每条LineString的坐标转为列表，完全匹配PathLayer要求的嵌套结构。

2. 传入PathLayer时，指定get_path="path_coords"即可：

import pydeck as pdk

mylayers=[
    pdk.Layer(
        'PathLayer',
        data=consline,
        get_path="path_coords",
        pickable=True,
        get_fill_color=[0,255,0],
        width_scale=20,
        width_min_pixels=2,   
    )
]
view_state = pdk.ViewState(
    longitude=-98.99,
    latitude=31.79,
    zoom=5,
    min_zoom=1,
    max_zoom=15)
r = pdk.Deck(layers=mylayers, initial_view_state=view_state, map_style='light')
r.to_html("example.html")

这种方式直接操作Shapely对象提取坐标，避免了JSON序列化/反序列化的额外开销，是更高效、更简洁的实现方式。

内容的提问来源于stack exchange，提问作者Dean MacGregor