Geodesic path

版本 1.0 (12.5 KB) 作者: Peter Seibold
Return coordinates for the shortest path between two points on the earth (WGS84 ellipsoid)
0.0
(0)
7.0 次下载
更新时间 2024/8/21

查看许可证

Get coordinates for the geodesic path that starts at lat1, lon1 and ends at lat2, lon2.
INPUT:
  • lat1: Start latitudes, double, column vector, -90 <= lat1 <= 90
  • lon1: Start longitudes, double, column vector, same size n as lat1
  • lat2: End latitudes, double, column vector, -90 <= lat2 <= 90, same size n as lat1
  • lon2: Start longitudes, double, column vector, same size n as lat1
  • points: optional, default 51 points including start and end point. They will have equal distances.
  • a) one integer number k -> k equidistant points including start and end point
  • b) a vector with k points in meter, points may be outside start and end points
OUTPUT:
  • latPath: latitude of path points, n x k double array
  • lonPath: longitude of path points, n x k double array
PRECISION:
The distance error of a path point is less than 100 nm
PROCESS TIME:
About 1 ms for 1000 geodesic path points on my PC.
Remark:
With latPath at pole, abs(latPath)>90-3*eps(90), lonPath is undefined, any lonPath is correct.
The code is based on two files in the geographiclib of Charles Karney
The original Matlab files 'geoddistance.m' and 'geodreckon.m' are at
https://de.mathworks.com/matlabcentral/fileexchange/50605-geographiclib
with copyright (c) Charles Karney (2012-2022) <charles@karney.com>.
The difference to the original files from C. Karney is:
Only calculation of a few parameters, all required subroutines included,
some unused parts removed, no Octave compatibility.

引用格式

Peter Seibold (2024). Geodesic path (https://www.mathworks.com/matlabcentral/fileexchange/171649-geodesic-path), MATLAB Central File Exchange. 检索时间: .

MATLAB 版本兼容性
创建方式 R2020a
兼容任何版本
平台兼容性
Windows macOS Linux
致谢

参考作品: geographiclib

Community Treasure Hunt

Find the treasures in MATLAB Central and discover how the community can help you!

Start Hunting!

版本 已发布 发行说明
1.0