Rank: 221629 based on 0 downloads (last 30 days) and 0 file submitted
Personal Profile:
Professional Interests:

 

Watch this Author's files

 

Comments and Ratings by Jonathan
Updated File Comments Rating
16 Oct 2009 Geodetic Toolbox Toolbox for angle, coordinate and date conversions and transformations. Version 2.97. Author: Mike Craymer

Hello,

Thank you for such a useful toolbox! I have a question, however, regarding transforming lat,lon on the WGS84 ellipsoid to local North East Up. If I use the same lat,lon coordinate as both my point to convert and my origin, why would I get a -21.5 km local North value? Shouldn't North and East both be zero?

Thanks again,
Jon

<SNIP>

>> [a,b,e2,finv] = refell('WGS84'); %get ellipsoid params
rlat = deg2rad(40); rlon = deg2rad(-108); %lat,lon
[Xec Yec Zec] = ell2xyz(rlat,rlon,0,a,e2) %convert to ECEF xyz
[N E Up] = ct2lg(Xec,Yec,Zec,rlat,rlon) %convert to local NEU

Xec =
-1.5119e+06
Yec =
-4.6532e+06
Zec =
4.0780e+06

N =
-2.1054e+04
E =
0
Up =
6.3693e+06

</SNIP>

Contact us