[mapping] 경도 \ 위도에서 데카르트 좌표로 변환

데카르트 좌표 로 변환 GPS（WGS84)하는 이론 https://en.wikipedia.org/wiki/Geographic_coordinate_conversion#From_geodetic_to_ECEF_coordinates

다음은 내가 사용하는 것입니다.

• GPS (WGS84)의 경도와 데카르트 좌표는 동일합니다.
• 위도는 WGS 84 타원체 매개 변수로 변환해야합니다. 반장 축은 6378137m이고
• 평탄화의 역수는 298.257223563입니다.

내가 작성한 VB 코드를 첨부했습니다 .

Imports System.Math

'Input GPSLatitude is WGS84 Latitude,h is altitude above the WGS 84 ellipsoid

Public Function GetSphericalLatitude(ByVal GPSLatitude As Double, ByVal h As Double) As Double

        Dim A As Double = 6378137 'semi-major axis
        Dim f As Double = 1 / 298.257223563  '1/f Reciprocal of flattening
        Dim e2 As Double = f * (2 - f)
        Dim Rc As Double = A / (Sqrt(1 - e2 * (Sin(GPSLatitude * PI / 180) ^ 2)))
        Dim p As Double = (Rc + h) * Cos(GPSLatitude * PI / 180)
        Dim z As Double = (Rc * (1 - e2) + h) * Sin(GPSLatitude * PI / 180)
        Dim r As Double = Sqrt(p ^ 2 + z ^ 2)
        Dim SphericalLatitude As Double =  Asin(z / r) * 180 / PI
        Return SphericalLatitude
End Function

의 h고도가 WGS 84 ellipsoid.

보통 GPS우리에게 H위의 MSL높이를 줄 것 입니다 . MSL높이는 높이로 변환하는 h상기 WGS 84 ellipsoid하여 중력 포텐샬 모델 EGM96( Lemoine은 등, 1998 ).
이것은 공간 해상도가 15 arc-minutes 인 지오이드 높이 파일의 그리드를 보간하여 수행됩니다.

당신은 몇 가지 레벨이있는 경우 또는 전문가를 GPS 고도를 가지고 H( MSL, 평균 해발 어이 )과 UNDULATION의 사이의 관계 geoid와 ellipsoid (m)선택된의 자료 출력 내부 테이블을. 당신은 얻을 수 있습니다h = H(msl) + undulation

데카르트 좌표로 XYZ하려면 :

x = R * cos(lat) * cos(lon)

y = R * cos(lat) * sin(lon)

z = R *sin(lat)

답변

LAT=40
LON=-110
echo $LON $LAT | cs2cs +proj=latlong +datum=WGS84 +to +proj=geocent +datum=WGS84

# Converting lat/long to cartesian
import numpy as np

def get_cartesian(lat=None,lon=None):
    lat, lon = np.deg2rad(lat), np.deg2rad(lon)
    R = 6371 # radius of the earth
    x = R * np.cos(lat) * np.cos(lon)
    y = R * np.cos(lat) * np.sin(lon)
    z = R *np.sin(lat)
    return x,y,z