大円を使用してRでマップ上のポイントを接続するためにマップと地球圏を使用する方法に関するこの素晴らしいチュートリアルを見つけました:
https://flowingdata.com/2011/05/11/how-to-map-connections-with-great-circles/
線が地図の端からはみ出す場合 (国際日付変更線を越える場合)、線を逆にする方法はありますか?
私の現在のコードは、DC と北京を接続する回線の場合、次のようになります。
library(maps)
library(geosphere)
map("world")
lat_ca <- 38.89511
lon_ca <- -77.03637
lat_me <- 39.917
lon_me <- 116.383
inter <- gcIntermediate(c(lon_ca, lat_ca), c(lon_me, lat_me), n=50, addStartEnd=TRUE)# breakAtDateLine=TRUE)
lines(inter)
明確にするために、次のようなものになりたいと思います。
これは、最初のものと同じ平面にある大円のもう 1 つの弧のようですが、描き方がわかりません。