How to get the points of [lat, lon] of a full great circle?

the matlab function gcwaypts generates points on the short arc of the great circle across two given points.

How to get the full great circle?

Answer 1

[Kelly Kearney]

I would use a combination of azimuth, gc2sc, and scircle1:

az = azimuth(lat1, lon1, lat2, lon2)
[lat, lon, r] = gc2sc(lat1, lon1, az)
[latc, lonc] = scircle1(lat, lon, r)

[Trung Ngo]

Hi Miss Kearney,

Between gcwaypts, track2 and your method. Which one is more preferable. I am doing this which take account of the ellipsoid.

[GClat_check,GClon_check]=track2('gc',startLat,startLon,endLat,endLon,referenceEllipsoid(7030),'degree',51);