34 lines
1.2 KiB
Matlab
34 lines
1.2 KiB
Matlab
function geodesicinverse(~, ~, ~)
|
|
%geodesicinverse Solve inverse geodesic problem
|
|
%
|
|
% [geodesic, aux] = geodesicinverse(latlong)
|
|
% [geodesic, aux] = geodesicinverse(latlong, a, f)
|
|
%
|
|
% latlong is an M x 4 matrix
|
|
% latitude of point 1 = latlong(:,1) in degrees
|
|
% longitude of point 1 = latlong(:,2) in degrees
|
|
% latitude of point 2 = latlong(:,3) in degrees
|
|
% longitude of point 2 = latlong(:,4) in degrees
|
|
%
|
|
% geodesic is an M x 3 matrix
|
|
% azimuth at point 1 = geodesic(:,1) in degrees
|
|
% azimuth at point 2 = geodesic(:,2) in degrees
|
|
% distance between points 1 and 2 = geodesic(:,3) in meters
|
|
% aux is an M x 5 matrix
|
|
% spherical arc length = aux(:,1) in degrees
|
|
% reduced length = aux(:,2) in meters
|
|
% geodesic scale 1 to 2 = aux(:,3)
|
|
% geodesic scale 2 to 1 = aux(:,4)
|
|
% area under geodesic = aux(:,5) in meters^2
|
|
%
|
|
% a = equatorial radius (meters)
|
|
% f = flattening (0 means a sphere)
|
|
% If a and f are omitted, the WGS84 values are used.
|
|
%
|
|
% A native MATLAB implementation is available as GEODDISTANCE.
|
|
%
|
|
% See also GEODDISTANCE.
|
|
|
|
error('Error: executing .m file instead of compiled routine');
|
|
end
|