Great Circle Arc
Get points along great circle between two locations
getgreatarc(lat1, lon1, lat2, lon2, num)
lat1
: Latitude, point 1 (degrees)lon1
: Longitude, point 1 (degrees)lat2
: Latitude, point 2 (degrees)lon2
: Longitude, point 2 (degrees)num
: number of points along arclat: Latitude
lon: Longitude
Jonathan M. Leesjonathan.lees@unc.edu
getgreatarc, distaz
PARIS = c(48.8666666666667, 2.33333333333333) RIODEJANEIRO =c( -22.9, -43.2333333333333) g = getgreatarc(PARIS[1],PARIS[2], RIODEJANEIRO[1], RIODEJANEIRO[2], 100) library(geomapdata) data(worldmap) plotGEOmap(worldmap, add=FALSE, shiftlon=180) lines(g$lon+180, g$lat)
Useful links