포인트와 폴리곤 간의 가장 긴 "직선"경로 찾기

Question

Point (위도와 경도)와 Polygon 사이의 가장 긴 숫자는 straight입니다. Shapely에 있습니까? 나는 가장 가까운 길을 찾는 것이 가능하다고 읽었지 만, 가장 길지는 모르겠다.

Answer 1

는 g1.hausdorff_distance(g2) fuction를 반환하는 Hausdorff distance을 시도해보십시오

from shapely.geometry import Polygon, Point 
poly = Polygon([(-1, -1), (-2, 2), (4, 4), (4, -1), (-1, -1)]) 
p = Point(0, 0) 
poly.hausdorff_distance(p) # 5.656854249492381 

매끈한는 직교 좌표 공간에서 작동, 명심하십시오. 귀하의 질문은 "위도와 경도"에 관해 묻습니다. 따라서이 거리 단위는도 단위입니다. 보다 일반적인 단위 길이를 얻으려면 이것을 적절한 좌표 참조 시스템 (CRS)에 투영해야합니다. 또한 "직선 경로"의 정의는 CRS의 선택에 따라 달라집니다.

Answer 2

한 가지 방법은 점의 집합으로 다각형의 외부를 이산화 포인트의 세트의 점과 각 점 사이의 거리를 계산한다 : 물론

def longest(poly, p, num): 

    lr = LinearRing(poly.exterior.coords) 
    dist = 0 
    for i in np.linspace(0, lr.length, num): 
     d = p.distance(lr.interpolate(i)) 
     if d > dist: 
      dist = d 
    return dist 


poly = Polygon([(-1, -1), (-2,2), (4,4), (4, -1), (-1, -1)]) 
p = Point(0,0) 
longest(poly, p, 20) 
# out: 5.428886519426354 
longest(poly, p, 100) 
# out: 5.622291976018042 

가 정확한 아니지만, 수 많은 상황에서 합리적인 근사치가되어야합니다.

비고 : 거리에 매끈함을 지닌 lon/lat를 사용하지 않아야합니다. 매끄럽게는 데카르트 좌표계를 기반으로합니다. 일반적으로 거리를 측정하기 위해 매끄럽게 사용할 수 있도록 먼저 (pyproj를 사용하여) 지오메트리를 투영해야합니다.