r - 在R地图类中查找国家经纬度

Question

我正在尝试在世界地图上使用R.

    # World map is available in the maps package
    library(maps)
    library(dplyr)
    
    # No margin
    par(mar=c(0,0,0,0))
    
    # World map
    m <- map('world',
        col="#f2f2f2", fill=TRUE, bg="white", lwd=0.05,
        mar=rep(0,4),border=0, ylim=c(-80,80) 
    )

# Get coordinates        
brazil <- c(m$x[which(m$names=="Brazil")], m$y[which(m$names=="Brazil")])
netherlands <- c(m$x[which(m$names=="Netherlands")], m$y[which(m$names=="Netherlands")])

# Data frame 
data <- rbind(brazil, netherlands) %>% 
  as.data.frame()
colnames(data) <- c("long","lat")

# No margin
    par(mar=c(0,0,0,0))
     map('world',
        col="#f2f2f2", fill=TRUE, bg="white", lwd=0.05,
        mar=rep(0,4),border=0, ylim=c(-80,80) 
    )
    points(x=data$long, y=data$lat, col="slateblue", cex=3, pch=20)

但是，这些点距离请求的国家还很远。

有什么想法可以正确指出这些国家吗？

Answer 1

县名与地图数据结构中的 x 和 y 值之间没有一对一的比率。

例如：

str(m)
List of 4
 $ x    : num [1:100375] -69.9 -69.9 -69.9 -70 -70.1 ...
 $ y    : num [1:100375] 12.5 12.4 12.4 12.5 12.5 ...
 $ range: num [1:4] -180 190 -80 80
 $ names: chr [1:1599] "Aruba" "Afghanistan" "Angola" "Angola:Cabinda" ...
 - attr(*, "class")= chr "map"

有 1599 个名称和超过 100,000 个 xy 对。

这是一种查找感兴趣的国家/地区，提取该信息，然后从范围值中找到中点的方法。在这个分辨率下足够接近。
您可能希望添加exact=TRUE到地图功能以隔离到该国的主要地理区域并移除外围岛屿和领土。美国是应用此选项的一个很好的例子。

# World map is available in the maps package
library(maps)

# World map
m <- map('world',
         col="#f2f2f2", fill=TRUE, bg="white", lwd=0.05,
         mar=rep(0,4),border=0, ylim=c(-80,80), plot=FALSE
)

# Get coordinates        
brazil <- map('world', regions = "brazil", plot=FALSE, exact=TRUE)
  brazil_long <- mean(brazil$range[1:2])
  brazil_lat <- mean(brazil$range[3:4])

netherlands <- map('world', regions = "Netherlands", plot=FALSE)    
  netherland_long <- mean(netherlands$range[1:2])
  netherland_lat <- mean(netherlands$range[3:4])

# Data frame 
data <- data.frame(long = c(brazil_long, netherland_long), lat= c(brazil_lat, netherland_lat))

# No margin
par(mar=c(0,0,0,0))
map('world',
    col="#f2f2f2", fill=TRUE, bg="white", lwd=0.05,
    mar=rep(0,4),border=0, ylim=c(-80,80) 
)
points(x=data$long, y=data$lat, col="slateblue", cex=3, pch=20)