导入请求、io、json 导入 geopandas 作为 gpd 导入 shapely.geometry 导入 pandas 作为 pd

我有两个 GeoPandas 数据帧。

gpd1 与 POLYGON 几何 gpd2 与 POINT 几何

示例数据：

res = requests.get(
    "https://opendata.arcgis.com/datasets/69dc11c7386943b4ad8893c45648b1e1_0.geojson"
)
gpd1 = gpd.GeoDataFrame.from_features(res.json()["features"], crs="CRS84").pipe(
    lambda d: d.rename(columns={c: c.lower() for c in d.columns})
).rename(columns={"lad20cd": "areaCode","lad20nm":"areaName"})


# get some public addressess - hospitals.  data that can be scattered
dfhos = pd.read_csv(io.StringIO(requests.get("http://media.nhschoices.nhs.uk/data/foi/Hospital.csv").text),
    sep="¬",engine="python",)


# create a geo dataframe of hospitals
gpd2 = gpd.GeoDataFrame(
    data=dfhos,
        geometry=dfhos.apply(lambda r: shapely.geometry.Point(r["Longitude"],r["Latitude"]), axis=1), crs="EPSG:4326"
    )

我空间加入他们：

gpd1.sjoin(gpd2, how='left')

我还想计算加入的每一行的距离。如何计算距离并将 distance 作为一列包含在内？请注意，gpd2 中的每个 POINT 都在 POLYGON 内。

import pandas as pd
import geopandas as gpd
import requests, io


# some polygons
gpd1 = gpd.read_file(
    requests.get("https://www.geoboundaries.org/api/current/gbOpen/GBR/ADM2").json()[
        "simplifiedGeometryGeoJSON"
    ]
)

# some points
df = pd.read_csv(
    io.StringIO(
        requests.get("https://simplemaps.com/static/data/country-cities/gb/gb.csv").text
    )
)


gpd2 = gpd.GeoDataFrame(
    df, geometry=gpd.points_from_xy(df["lng"], df["lat"]), crs="epsg:4326"
)


gpd_sj = gpd1.sjoin(gpd2, how="inner").reset_index()
gpd_sj["distance"] = gpd_sj.distance(gpd2.loc[gpd_sj["index_right"]])
gpd_sj.loc[:, ["index", "shapeName", "index_right", "city", "distance"]]