Python中点和MultiPolygon Geoseries对象之间的距离

Question

我想计算一个点与使用 GeoPandas 的 cascaded_union 函数创建的 Multiploygon geoseries 对象的距离。以下是我的代码：

import numpy as np
import pandas as pd
import geopandas as gpd
from shapely.ops import cascaded_union
from matplotlib import pyplot as plt
from shapely.geometry import LineString, Point, MultiPoint

fp = "\\TM_WORLD_BORDERS-0.3\\TM_WORLD_BORDERS-0.3.shp"
data = gpd.read_file(fp)

mycountries = ['Canada','Russia']

req_data = data[data.NAME.isin(mycountries)]
req_data = req_data.reset_index()
req_data = req_data.drop('index',axis=1)

polygons = [req_data['geometry'][0],req_data['geometry'][1]]
boundary = gpd.GeoSeries(cascaded_union(polygons))
boundary.plot(color = 'black')
plt.show()

我要求当我计算一个点与多面体的距离时，它应该考虑最近的国家/多边形，并且距离应该以英里/公里为单位。对此的任何帮助将不胜感激。

我尝试使用距离函数，但得到的结果不正确。此外，如果我考虑印度洋某处的某个点，它给出的结果为 0，这又不是我所需要的。

形状文件可以从 http://thematicmapping.org/downloads/world_borders.php

上的 TM_WORLD_BORDERS-0.3.zip 文件下载

Answer 1

这里有一个使用 Point.distance 的解决方案：

from shapely.geometry import Point
point = Point(34.23, 86.44)
dist={}
for index, row in req_data.iterrows():
    if row.geometry.geom_type == 'MultiPolygon':
        for  r in row.geometry:
            print(row.NAME,'  dist', point.distance(r))
            if row.NAME in dist:
                l = dist[row.NAME];l.append(point.distance(r))
                dist[row.NAME] = l
            else:
                dist[row.NAME] = [point.distance(r)]
    else: #Polygon
        print(row.NAME,'  dist', point.distance(row.geometry))
        dist[row.NAME] = [point.distance(row.geometry)]


for x in dist:
    print(x, min(dist[x]))

输出：与每个国家/地区的最小距离

 Canada 95.0619464401356
 Russia 12.118509615115519