在 Python 中的地图上绘制方式（线串）

Question

这是我第二次尝试同样的问题，我真的希望有人能帮助我... 甚至认为一些非常好的人试图帮助我。尽管有帮助，但仍有很多我无法弄清楚。

从头开始：我创建了一个数据框。这个数据框很大，可以提供有关城市中旅行者的信息。数据框看起来像这样。这只是头部。

在出发地和目的地中，您有城市位置的 ID，在移动中有多少人从出发地到目的地。经度和纬度是确切点所在的位置，线串是点的组合..

我用这段代码创建了线串：

erg2['Linestring'] = erg2.apply(lambda x: LineString([(x['latitude_origin'], x['longitude_origin']), (x['latitude_destination'], x['longitude_destination'])]), axis = 1)

现在我的问题是如何在地图上标出路线。甚至以为我尝试了 geopandas 纪录片等中的所有示例。我无法自拔..

我无法向您展示我已经绘制的内容，因为它没有意义，我想从头开始绘制会更聪明。

你看到在move列中有一些0。这意味着没有人走过这条路线。所以这个我不需要情节..
我必须用旅行者开始的地方 origin 和他去的地方 destination.
的信息绘制线条
我还需要根据动作勾勒出不同的线条..

用这个绘图代码

fig = px.line_mapbox(erg2, lat="latitude_origin", lon="longitude_origin", color="move", 
hover_name= gdf["origin"] + " - " + gdf["destination"],
                         center =dict(lon=13.41053,lat=52.52437), zoom=3, height=600
                        )
    
    fig.update_layout(mapbox_style="stamen-terrain", mapbox_zoom=4, mapbox_center_lat = 52.52437,
        margin={"r":0,"t":0,"l":0,"b":0})
    
    fig.show()

也许有人有想法？？？

我用thios代码试过了：

    import requests, io, json
import geopandas as gpd
import shapely.geometry
import pandas as pd
import numpy as np
import itertools
import plotly.express as px

# get some public addressess - hospitals.  data that has GPS lat / lon
dfhos = pd.read_csv(io.StringIO(requests.get("http://media.nhschoices.nhs.uk/data/foi/Hospital.csv").text),
    sep="¬",engine="python",).loc[:, ["OrganisationName", "Latitude", "Longitude"]]

a = np.arange(len(dfhos))
np.random.shuffle(a)
# establish N links between hospitals
N = 10
df = (
    pd.DataFrame({0:a[0:N], 1:a[25:25+N]}).merge(dfhos,left_on=0,right_index=True)
    .merge(dfhos,left_on=1, right_index=True, suffixes=("_origin", "_destination"))
)

# build a geopandas data frame that has LineString between two hospitals
gdf = gpd.GeoDataFrame(
    data=df,
    geometry=df.apply(
        lambda r: shapely.geometry.LineString(
            [(r["Longitude_origin"], r["Latitude_origin"]),
             (r["Longitude_destination"], r["Latitude_destination"]) ]), axis=1)
)

# sample code https://plotly.com/python/lines-on-mapbox/#lines-on-mapbox-maps-from-geopandas
lats = []
lons = []
names = []

for feature, name in zip(gdf.geometry, gdf["OrganisationName_origin"] + " - " + gdf["OrganisationName_destination"]):
    if isinstance(feature, shapely.geometry.linestring.LineString):
        linestrings = [feature]
    elif isinstance(feature, shapely.geometry.multilinestring.MultiLineString):
        linestrings = feature.geoms
    else:
        continue
    for linestring in linestrings:
        x, y = linestring.xy
        lats = np.append(lats, y)
        lons = np.append(lons, x)
        names = np.append(names, [name]*len(y))
        lats = np.append(lats, None)
        lons = np.append(lons, None)
        names = np.append(names, None)

fig = px.line_mapbox(lat=lats, lon=lons, hover_name=names)

fig.update_layout(mapbox_style="stamen-terrain",
                  mapbox_zoom=4,
                  mapbox_center_lon=gdf.total_bounds[[0,2]].mean(),
                  mapbox_center_lat=gdf.total_bounds[[1,3]].mean(),
                  margin={"r":0,"t":0,"l":0,"b":0}
                 )

这看起来是完美的代码，但我不能真正将它用于我的数据.. 我对编码很陌生。所以请耐心等待;))

非常感谢。

祝一切顺利

Answer 1

之前回答过这个问题 How to plot visualize a Linestring over a map with Python?。我建议你更新那个问题，我仍然建议你这样做
行字符串恕我直言不是要走的路。 plotly 不使用线串，因此编码为线串以解码为 numpy 数组更加复杂。查看官方文档中的示例 https://plotly.com/python/lines-on-mapbox/。这里很清楚 geopandas 只是一个必须编码成 numpy 数组的源

数据

你的样本数据应该是一个 Dataframe，不需要 geopandas 或 line strings
几乎所有示例数据都无法使用，因为起点和终点不同的每一行都有移动为零，您注意到应该将其排除在外

import pandas as pd
import numpy as np
import plotly.express as px

df = pd.DataFrame({"origin": [88, 88, 88, 88, 88, 87], 
                   "destination": [88, 89, 110, 111, 112, 83], 
                   "move": [20, 0, 5, 0, 0, 10], 
                   "longitude_origin": [13.481016, 13.481016, 13.481016, 13.481016, 13.481016, 13.479667], 
                   "latitude_origin": [52.457055, 52.457055, 52.457055, 52.457055, 52.457055, 52.4796], 
                   "longitude_destination": [13.481016, 13.504075, 13.613772, 13.586891, 13.559341, 13.481016], 
                   "latitude_destination": [52.457055, 52.443923, 52.533194, 52.523562, 52.507418, 52.457055]})

解决方案

进一步完善了 line_array() 函数，因此它可用于对我之前提供的简化解决方案中的悬停和颜色参数进行编码

# lines in plotly are delimited by none
def line_array(data, cols=[], empty_val=None):
    if isinstance(data, pd.DataFrame):
        vals = data.loc[:, cols].values
    elif isinstance(data, pd.Series):
        a = data.values
        vals = np.pad(a.reshape(a.shape[0], -1), [(0, 0), (0, 1)], mode="edge")
    return np.pad(vals, [(0, 0), (0, 1)], constant_values=empty_val).reshape(
        1, (len(df) * 3))[0]


# only draw lines where move > 0 and destination is different to origin
df = df.loc[df["move"].gt(0) & (df["origin"]!=df["destination"])]

lons = line_array(df, ["longitude_origin", "longitude_destination"])
lats = line_array(df, ["latitude_origin", "latitude_destination"])

fig = px.line_mapbox(
    lat=lats,
    lon=lons,
    hover_name=line_array(
        df.loc[:, ["origin", "destination"]].astype(str).apply(" - ".join, axis=1)
    ),
    hover_data={
        "move": line_array(df, ["move", "move"], empty_val=-99),
        "origin": line_array(df, ["origin", "origin"], empty_val=-99),
    },
    color=line_array(df, ["origin", "origin"], empty_val=-99),
).update_traces(visible=False, selector={"name": "-99"})

fig.update_layout(
    mapbox={
        "style": "stamen-terrain",
        "zoom": 9.5,
        "center": {"lat": lats[0], "lon": lons[0]},
    },
    margin={"r": 0, "t": 0, "l": 0, "b": 0},
)

在 Python 中的地图上绘制方式（线串）

Plotting ways (linestrings) over a map in Python

numpy

pandas

shapely

geopandas

plotly-express

数据

解决方案