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
|
function geocentricforward(~, ~, ~)
%geocentricforward Convert geographic coordinates to geocentric
%
% [geocentric, rot] = geocentricforward(geodetic)
% [geocentric, rot] = geocentricforward(geodetic, a, f)
%
% geodetic is an M x 3 or M x 2 matrix of geodetic coordinates
% lat = geodetic(:,1) in degrees
% lon = geodetic(:,2) in degrees
% h = geodetic(:,3) in meters (default 0 m)
%
% geocentric is an M x 3 matrix of geocentric coordinates
% x = geocentric(:,1) in meters
% y = geocentric(:,2) in meters
% z = geocentric(:,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 (east, north, up)
% by M transforms the vector to geocentric coordinates.
%
% 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
% Geocentric::Forward
% See the documentation on this function for more information:
% http://geographiclib.sf.net/html/classGeographicLib_1_1Geocentric.html
error('Error: executing .m file instead of compiled routine');
end
% geocentricforward.m
% Matlab .m file for geographic to geocentric conversions
%
% Copyright (c) Charles Karney (2011) <charles@karney.com> and licensed
% under the MIT/X11 License. For more information, see
% http://geographiclib.sourceforge.net/
|
