Konvertieren von LLA zu XYZ

Können Sie mir helfen zu konvertieren lla zu xyz coordnates.

Ich bin mit

earthRadius = 6378.137;
var x = earthRadius * Math.cos(lat)*Math.cos(lon);
var y = earthRadius * Math.cos(lat)*Math.sin(lon);
var z = earthRadius * Math.sin(lat);

diese Methode, um zu konvertieren xyz-Koordinaten. Aber es gibt nicht das richtige Ergebnis, die ich möchte.
Die Erde ist definiert als wgs84-Objekt.

also, was ist das richtige Ergebnis?
X : 192.952 km Y : -1094.284 km Z : 6259.543 km sind richtige Ergebnisse, aber meine Ergebnisse sind X : 192.324 km Y : -1090.725 km Z : 6281.238 km
Ich bin noch nicht sicher, was Sie genau möchten, KM ist ein Maß für die Distanz zwischen zwei Punkten. Zwischen dem Punkt und der Breitengrad / Längengrad / Höhe versuchen Sie, die Entfernung zu berechnen?
Ich bin nicht die Berechnung der Entfernung, ich möchte konvertieren LLA Werte auf XYZ-Koordinaten meiner lla-Eingänge sind 100, 100, 0.1, und ich schrieb meine Ergebnisse und die tatsächlichen Ergebnisse vor. Funktioniert diese Formel ein korrektes Ergebnis ?
Mögliche Duplikate von Lat, Long, X-Y-Z-position in JS .. funktioniert nicht

InformationsquelleAutor mcelik | 2013-09-12

coordinates wgs84

1

Ich mein problem gelöst

Lat Long in X-Y-Z-position in JS nicht funktioniert ..

diese Formel ist die richtige. Vielen Dank für Stephen Quan für die Beantwortung dieser.

InformationsquelleAutor mcelik

Komplette Koordinatentransformation (python-code):

geodetic -> ECEF (earth centered earth fixed), Entfernungen von der Mitte der Erde

ECEF -> ENU (East-North-Up) zu erhalten, die lokalen typischen x-y Flugbahn:

from __future__ import print_function
import math
import pyproj
R = 6378137
f_inv = 298.257224
f = 1.0 /f_inv
e2 = 1 - (1 - f) * (1 - f)

coords = [
  (0,  45,  1000),
  (45,  90,  2000),
  (48.8567,  2.3508,  80),
  (61.4140105652, 23.7281341313,149.821),
]

def gps_to_ecef_pyproj(lat, lon, alt):
    ecef = pyproj.Proj(proj='geocent', ellps='WGS84', datum='WGS84')
    lla = pyproj.Proj(proj='latlong', ellps='WGS84', datum='WGS84')
    x, y, z = pyproj.transform(lla, ecef, lon, lat, alt, radians=False)

    return x, y, z

def gps_to_ecef(latitude, longitude, altitude):
    # (lat, lon) in WSG-84 degrees
    # h in meters
    cosLat = math.cos(latitude * math.pi /180)
    sinLat = math.sin(latitude * math.pi /180)

    cosLong = math.cos(longitude * math.pi /180)
    sinLong = math.sin(longitude * math.pi /180)

    c = 1 /math.sqrt(cosLat * cosLat + (1 - f) * (1 - f) * sinLat * sinLat)
    s = (1 - f) * (1 - f) * c

    x = (R*c + altitude) * cosLat * cosLong
    y = (R*c + altitude) * cosLat * sinLong
    z = (R*s + altitude) * sinLat

    return x, y, z

# ecef2enu
def ecef_to_enu(x, y, z, latRef, longRef, altRef):

    cosLatRef = math.cos(latRef * math.pi /180)
    sinLatRef = math.sin(latRef * math.pi /180)

    cosLongRef = math.cos(longRef * math.pi /180)
    sinLongRef = math.sin(longRef * math.pi /180)

    cRef = 1 /math.sqrt(cosLatRef * cosLatRef + (1 - f) * (1 - f) * sinLatRef * sinLatRef)

    x0 = (R*cRef + altRef) * cosLatRef * cosLongRef
    y0 = (R*cRef + altRef) * cosLatRef * sinLongRef
    z0 = (R*cRef*(1-e2) + altRef) * sinLatRef

    xEast = (-(x-x0) * sinLongRef) + ((y-y0)*cosLongRef)

    yNorth = (-cosLongRef*sinLatRef*(x-x0)) - (sinLatRef*sinLongRef*(y-y0)) + (cosLatRef*(z-z0))

    zUp = (cosLatRef*cosLongRef*(x-x0)) + (cosLatRef*sinLongRef*(y-y0)) + (sinLatRef*(z-z0))

    return xEast, yNorth, zUp

def geodetic_to_enu(lat, lon, h, lat_ref, lon_ref, h_ref):
    x, y, z = gps_to_ecef(lat, lon, h)

    return ecef_to_enu(x, y, z, lat_ref, lon_ref, h_ref)


def run_test():

    for pt in coords:

    xPy,yPy,zPy = gps_to_ecef_pyproj(pt[0], pt[1], pt[2])   
    xF,yF,zF = gps_to_ecef(pt[0], pt[1], pt[2])

        print("pyproj (XYZ)\t = ", xPy, yPy, zPy)
        print("ECEF (XYZ)\t = ", xF, yF, zF)

    xE, yN, zU = ecef_to_enu(xF,yF,zF,pt[0], pt[1], pt[2])
        print('ENU (XYZ) \t = ', xE, yN, zU)

    print("-------------------------------------------------")


run_test()

Den Unterschied zwischen geodetic -> ECEF und ECEF -> ENU wird anhand der zwei folgenden Abbildungen: