flat2lla matlab,Convert from geodetic latitude, longitude, and altitude to flat Earth position - MAT...

Algorithms

The function begins by finding the small changes in latitude and longitude from the output

latitude and longitude minus the initial latitude and longitude:

dμ=μ−μ0dι=ι−ι0.

To convert geodetic latitude and longitude to the north and east coordinates, the function

uses the radius of curvature in the prime vertical

(RN) and the radius of curvature in the meridian

(RM).

RN and

RM are defined by the following relationships:

RN=R1−(2f−f2)sin2μ0,

where (R) is the equatorial radius of the planet and f is the flattening of the planet.

Small changes in the north (dN) and east (dE) positions are approximated from small

changes in the north and east positions by

dN=dμatan(1RM),

and

dE=dιatan(1RNcosμ0).

With the conversion of the North and East coordinates to the flat Earth

x and y coordinates, the transformation has the form of

[pxpy]=[cosψsinψ−sinψcosψ][NE],

where

(ψ)

is the angle in degrees clockwise between the x-axis

and north.

The flat Earth z-axis value is the negative altitude minus the

reference height (href):

pz=−h−href.