웹 메르카토르 맵에서보기 좋은 원호 만들기?

일부 항공편을 나타내는 데이터가 포함 된지도를 만들고 그레이트 서클 호를 사용하여 소스와 목적지를 연결하려고합니다.

기본적으로 유명한 Facebook 맵과 비슷한 것을하고 싶습니다.

이 게시물에 제공된 기능을 사용했습니다 : https://gis.stackexchange.com/a/5205/442 (예 :이 블로그 기사 : http://anitagraser.com/2011/08/20/visualizing-global-connections / ) 및 줄을 얻을 수는 있지만 국제 날짜 줄을지나 극에서 묶습니다.

@underdark는 링크 된 블로그 포스트 에서이 라인을 분리해야한다고 언급했지만 PostGIS에서 자동으로 분리하는 방법을 모르겠습니다.

또한 극 근처의 전선 다발도 해결해야합니다.

이 항공편의 출발 및 도착 지점이있을 때 두 가지를 모두 수행하려면 어떻게합니까?

이 블로그 게시물의 원칙은 일반적인 PostGIS로 이전됩니다.

http://blog.cartodb.com/jets-and-datelines/

기본적으로 ST_Segmentize지리학을 사용 하고 데이트 라인 교차 선을 슬라이스하는 약간의 마술을 사용하십시오.

측지선을 계산할 수 있습니다. 측지선을 A에서 B로 표시하려면 먼저 A에서 B까지의 거리와 방위 (역 측지 문제)를 계산 한 다음 A에서 B와 A와 B 사이의 여러 지점 (직접 측지 문제)을 계산할 수 있습니다. GeographicLib을 사용하여 Python에서 간단한 스크립트를 추가하여 GeoJSON에서 물건을 출력했습니다.

from geographiclib.geodesic import Geodesic
from geojson import MultiLineString

def geodesic(lat1, lon1, lat2, lon2, steps):
    inverse = Geodesic.WGS84.Inverse(lat1, lon1, lat2, lon2)
    linestrings = []
    coordinates = []

    for i in range(0, steps + 1):
        direct = Geodesic.WGS84.Direct(inverse['lat1'], inverse['lon1'], inverse['azi1'], (i / float(steps)) * inverse['s12'])
        if len(coordinates) > 0:
            if (coordinates[-1][0] < -90 and direct['lon2'] > 90) or (coordinates[-1][0] > 90 and direct['lon2'] < -90):
                linestrings.append(coordinates)
                coordinates = []
        coordinates.append((direct['lon2'], direct['lat2']))

    linestrings.append(coordinates)
    geojson = MultiLineString(linestrings)
    return geojson

linestrings = []

# San Francisco: 37.7793, -122.4192
# Bangalore: 12.9, 77.616667
for linestring in geodesic(37.7793, -122.4192, 12.95, 77.616667, 100)['coordinates']:
    linestrings.append(linestring)

# Boston: 42.357778, -71.059444
# Bangalore: 12.9, 77.616667
for linestring in geodesic(42.357778, -71.059444, 12.95, 77.616667, 100)['coordinates']:
    linestrings.append(linestring)

print(MultiLineString(linestrings))

결과는 WGS-84의 점들 사이의 진정한 측지선입니다. 물론 좌표를 필요한 투영으로 변환 할 수 있습니다. geojson.io에 시각화 된 결과 는 다음과 같습니다.