绘制投影地图上的连线(mapproj,gcIntermediate)

时间:2017-10-21 20:20:39

标签: r map-projections geosphere

我试图用Mollweide投影在世界地图上绘制线条是不成功的。我还在同一张地图上绘制了点,并且效果很好。对于这些行,我尝试根据我的需要调整此示例:http://flowingdata.com/2011/05/11/how-to-map-connections-with-great-circles/。 我已经失败了预测试(前面的第4步)。在下面的代码中,该线路应该连接肯尼亚和澳大利亚。它运行没有错误,但输出中没有行。 (我也测试了没有mapproj的例子,那条线就在那里。)

library("maps")
library("mapproj")
library("geosphere")

map("world",proj="mollweide", orientation= c(90,0,0), fill=TRUE, col="white", bg="lightblue")

lon_ke <- 38
lat_ke <- 1

lon_aus <- 133
lat_aus <- -27

inter <- gcIntermediate(c(mapproject(lon_ke,lat_ke), proj="mollweide", orientation= c(90,0,0)),
                        c(mapproject(lon_aus,lat_aus), proj="mollweide", orientation= c(90,0,0)),
                        n=50, addStartEnd=TRUE)
lines(inter)

1 个答案:

答案 0 :(得分:0)

根据Thomas Rahlf的书(见评论),我找到了解决问题的方法。这是我的脚本(它有助于可视化作者发表文章的位置)。

import * as DatGui from "@types/dat-gui";
var gui = new DatGui.GUI()