1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
|
function localcartesianreverse(~, ~, ~, ~)
%localcartesianreverse Convert local cartesian coordinates to geographic
%
% [geodetic, rot] = localcartesianreverse(origin, cartesian)
% [geodetic, rot] = localcartesianreverse(origin, cartesian, a, f)
%
% origin is a 1 x 3 or 1 x 2 matrix
% lat0 = origin(1,1) in degrees
% lon0 = origin(1,2) in degrees
% h0 = origin(1,3) in meters (default 0 m)
% cartesian is an M x 3 matrix of local cartesian coordinates
% x = cartesian(:,1) in meters
% y = cartesian(:,2) in meters
% z = cartesian(:,3) in meters
%
% geodetic is an M x 3 matrix of geodetic coordinates
% lat = geodetic(:,1) in degrees
% lon = geodetic(:,2) in degrees
% h = geodetic(:,3) in meters
% rot is an M x 9 matrix
% M = rot(:,1:9) rotation matrix in row major order. Pre-multiplying
% a unit vector in local cartesian coordinates at (lat0, lon0, h0)
% by the transpose of M transforms the vector to local cartesian
% coordinates at (lat, lon, h).
%
% a = major radius (meters)
% f = flattening (0 means a sphere)
% If a and f are omitted, the WGS84 values are used.
%
% This is an interface to the GeographicLib C++ routine
% LocalCartesian::Reverse
% See the documentation on this function for more information:
% http://geographiclib.sf.net/html/classGeographicLib_1_1LocalCartesian.html
error('Error: executing .m file instead of compiled routine');
end
% localcartesianreverse.m
% Matlab .m file for local cartesian to geographic conversions
%
% Copyright (c) Charles Karney (2011) <charles@karney.com> and licensed under
% the MIT/X11 License. For more information, see
% http://geographiclib.sourceforge.net/
|
