CL_co_car2ell — Cartesien to elliptical coordinates
[pos_ell[,jacob]] = CL_co_car2ell(pos_car[,er[,obl]])
If no optionnal arguments are passed, defaults are the earth ellipsoid semi-major axis(equatorial radius) and oblateness. Coordinates are then called geodetic coordinates
Transformation jacobian is optionnally computed.
Warning : pos_ell contains longitude, elliptical latitude, and elliptical altitude (whereas in spherical coordinates it is longitude, spherical latitude and distance to center
positions in cartesien coordinates (in terrestrial frame if used with earth) [X;Y;Z] [m] (3xN)
(optional) reference ellipsoid semi-major axis [m] (default is earth equatorial radius : %CL_eqRad)
(optional) reference ellipsoid oblateness (default is earth oblateness : %CL_obla)
positions in elliptical (or geodetic) coordinates [longitude;latitude;alt] [rad],[m] (3xN)
(optional) transformation jacobian (3x3xN)
1 Mecanique Spatiale, Cnes - Cepadues Editions, Tome I, section 3.2.3 (Les reperes de l'espace et du temps, Relations entre les coordonnees)
2 CNES - MSLIB FORTRAN 90, Volume T (mt_car_geod)
// Example 1 pos_car = [4637885.347 ; 121344.608 ; 4362452.869]; [pos_ell] = CL_co_car2ell(pos_car); // Example 2 : reciprocity test pos_car = [4637885.347,121344.608,4362452.869]'; [pos_ell,jacob1] = CL_co_car2ell(pos_car); [pos_car2,jacob2] = CL_co_ell2car(pos_ell); M = jacob1*jacob2; // unitary matrix