0

即使我按照另一篇文章中的建议加载了 maptools 库，我也无法让 ggplot 的 fortify 或 broom 的 tidy 包含区域。

我首先加载了一堆库，包括 maptools (0.8-41)、rgeos (0.3-22)、broom (0.4.1) 和 ggplot2 (2.2.1.9000)。接下来，我使用以下命令抓取一张世界地图，其中包含一系列不同的区域选择，包括我感兴趣的一个 - ISO_a3。

world <- readOGR(dsn="https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson", layer="OGRGeoJSON")

接下来我跑map <- broom::tidy(world, region = "iso_a3")

head(map) 生成以下内容。请注意，不包括区域列。

      long       lat order  hole piece group  id
1 48.93857 11.258447     1 FALSE     1 -99.1 -99
2 48.93848 10.982324     2 FALSE     1 -99.1 -99
3 48.93848 10.714209     3 FALSE     1 -99.1 -99
4 48.93838 10.433252     4 FALSE     1 -99.1 -99
5 48.93828  9.973486     5 FALSE     1 -99.1 -99
6 48.93828  9.807617     6 FALSE     1 -99.1 -99
4

1 回答 1

1

Afaik，fortify并且tidy不要附加具有相同名称的列。另见示例

if (require("maptools")) {
 sids <- system.file("shapes/sids.shp", package="maptools")
 nc1 <- readShapePoly(sids,
   proj4string = CRS("+proj=longlat +datum=NAD27"))
 nc1_df <- fortify(nc1)
}
head(nc1_df)

您可能正在寻找id包含iso_a3s 的列：

library(rgdal)
library(ggplot2)
fn <- file.path(tempdir(), "ne_50m_admin_0_countries.geojson.gz")
download.file("https://raw.githubusercontent.com/nvkelso/natural-earth-vector/master/geojson/ne_50m_admin_0_countries.geojson.gz", fn)
R.utils::gunzip(fn)
world <- readOGR(dsn = tools::file_path_sans_ext(fn), layer = "OGRGeoJSON")
map <- broom::tidy(world, region="iso_a3")

all(sort(unique(map$id))==sort(levels(world$iso_a3)))
# [1] TRUE

ggplot(map, aes(long, lat, group=group)) + 
  geom_polygon(color="white") + 
  coord_quickmap()
于 2017-01-31T21:48:24.260 回答