scientimate.distancegc#
arclen, azimuthdir, metedir = scientimate.distancegc(lat1, lon1, lat2, lon2, CalcMethod='haversine', R=6371000, dispout='no')
Description#
Calculate distance and azimuth (bearing) between (Latitude,Longitude) points using Great Circle
Inputs#
- lat1
Latitude (y) of start point (first point) in (Degree)
- lon1
Longitude (x) of start point (first point) in (Degree)
- lat2
Latitude (y) of end point (last point) in (Degree)
- lon2
Longitude (x) of end point (last point) in (Degree)
- CalcMethod=’haversine’
- Distance calculation method‘cos’: Spherical law of cosines‘haversine’: Haversine formula‘vincenty’: Vincenty formula, Vincenty (1975)
- R=6371000
Earth radius in (m), mean earth radius=6371000 m
- dispout=’no’
Define to display outputs or not (‘yes’: display, ‘no’: not display)
Outputs#
- arclen
Total distance from start point to end point in (m)
- azimuthdir
- Azimuth (bearing or compass direction) from start point to end point in (Degree)0 (degree): toward North, 90 (degree): toward East, 180 (degree): toward South, 270 (degree): toward West
- metedir
- Meteorological direction from start point to end point in (Degree)0 (degree): from North, 90 (degree): from East, 180 (degree): from South, 270 (degree): from West
Examples#
import scientimate as sm
lat1=29.5 #First point
lon1=-89.4 #First point
lat2=29.7 #last point
lon2=-89.4 #last point
arclen,azimuthdir,metedir=sm.distancegc(lat1,lon1,lat2,lon2)
lat1=[29.5,29] #First point
lon1=[-89.4,-89] #First point
lat2=[29.7,30] #Last point
lon2=[-89.4,-90] #Last point
arclen,azimuthdir,metedir=sm.distancegc(lat1,lon1,lat2,lon2,'haversine',6371000,'yes')
References#
Vincenty, T. (1975). Direct and inverse solutions of geodesics on the ellipsoid with application of nested equations. Survey review, 23(176), 88-93.