Rank: 44259 based on downloads (last 30 days) and 0 files 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.92a. 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>

 

MATLAB Central Terms of Use

NOTICE: Any content you submit to MATLAB Central, including personal information, is not subject to the protections which may be afforded information collected under other sections of The MathWorks, Inc. Web site. You are entirely responsible for all content that you upload, post, e-mail, transmit or otherwise make available via MATLAB Central. The MathWorks does not control the content posted by visitors to MATLAB Central and, does not guarantee the accuracy, integrity, or quality of such content. Under no circumstances will The MathWorks be liable in any way for any content not authored by The MathWorks, or any loss or damage of any kind incurred as a result of the use of any content posted, e-mailed, transmitted or otherwise made available via MATLAB Central. Read the complete Terms prior to use.

Contact us at files@mathworks.com