Konwersja szerokości i długości geograficznej na układ współrzędnych ECEF

Question

Studiuję pArk przykładowy kod Apple i jak to działa. ktoś wie, w jaki sposób przekonwertować szerokość i długość geograficzną na współrzędne ECEF, i Zamień ECEF na współrzędne ENU wyśrodkowane na podaną liczbę lat, funkcje lon działają? Chcę tylko zrozumieć, co się dzieje w tej funkcji!

dzięki.

void latLonToEcef(double lat, double lon, double alt, double *x, double *y, double *z) 
{ 
    double clat = cos(lat * DEGREES_TO_RADIANS); 
    double slat = sin(lat * DEGREES_TO_RADIANS); 
    double clon = cos(lon * DEGREES_TO_RADIANS); 
    double slon = sin(lon * DEGREES_TO_RADIANS); 

    double N = WGS84_A/sqrt(1.0 - WGS84_E * WGS84_E * slat * slat); 

    *x = (N + alt) * clat * clon; 
    *y = (N + alt) * clat * slon; 
    *z = (N * (1.0 - WGS84_E * WGS84_E) + alt) * slat; 
} 
// Coverts ECEF to ENU coordinates centered at given lat, lon 
void ecefToEnu(double lat, double lon, double x, double y, double z, double xr, double yr, double zr, double *e, double *n, double *u) 
{ 
    double clat = cos(lat * DEGREES_TO_RADIANS); 
    double slat = sin(lat * DEGREES_TO_RADIANS); 
    double clon = cos(lon * DEGREES_TO_RADIANS); 
    double slon = sin(lon * DEGREES_TO_RADIANS); 
    double dx = x - xr; 
    double dy = y - yr; 
    double dz = z - zr; 

    *e = -slon*dx + clon*dy; 
    *n = -slat*clon*dx - slat*slon*dy + clat*dz; 
    *u = clat*clon*dx + clat*slon*dy + slat*dz; 
} 

Answer 1

Sposób latLonToEcef stanowi realizację algorytmu przedstawionego na stronie Wikipedia:

gdzie

Φ ma szerokość, i N; jest długość, a

Podobnie metoda ecefToEnu jest implementacją algorytmu na ECEF to ENU:

Jeśli potrzebujesz dalszych odniesień, można je znaleźć na dole, że Strona w Wikipedii. Możesz także zapoznać się ze specyfikacją World Geodetic System 1984.