WGS84, 두 지점 사이의 거리

남한 정도의 지형에서라면 굳이 구면 혹은 타원체 두 점 사이의 거리를 이용해야 할 정도는 아닐 것이다. 그럼에도 불구하고 처리 시간이 엄청 바뀌지 않는다면 제대로 측정해 보는 것도 좋을 것이다.

좌표계 변환에 사용했던 pyproj 의 Geod 를 이용한다. 입력값으로 scalar와 array 모두 가능하다. 그렇다면 편리하게 array 형식으로 한 번에 계산한다.

geo = pandas.read_csv("area_data.csv")

import pyproj
g = pyproj.Geod(ellps='WGS84')
_,_,dist = g.inv([target_long] * len(geo), 
                 [target_lati] * len(geo), 
                 geo['longitude'],
                 geo['latitude'])
geo['distance'] = dist
geo.sort_values(by='distance', inplace = True)

결과적으로 보면 이렇게 할 경우 거리는 m 단위로 나오는 것 같다.