在 R 中创建地图就像 rworldmap 的方式一样，但针对具有省份的特定国家

Question

我将 rworldmap 软件包与 WorldBank Data 一起使用，我很喜欢。我想绘制伊朗地图，其中包含与每个省相关的数据。这样做的步骤是什么？我知道我们可以在 R 中绘制像美国这样的某些国家/地区的地图，但并非所有国家/地区都如此。

Answer 1

您可以使用 raster 包直接导入 GADM 数据。然后，您可以使用 ggplot2 绘制地图。下载数据时，您可以指定不同的级别。根据这个你会看到不同的边界。

library(raster)
library(ggplot2)

### Get data
iran<- getData("GADM", country = "Iran", level = 2)

### SPDF to DF
map <- fortify(iran)

### Draw a map
ggplot() +
geom_map(data = map, map = map, aes(x = long, y = lat, map_id = id, group = group))

编辑

看到安迪的回答，我想展示一下如何在ggplot2中添加省份名称。这是我从@hrbrmstr 那里学到的东西。 rworldmap 包允许你输入更少。 ggplot2 软件包仍然提供非常漂亮的图形。

library(raster)
library(rgdal)
library(rgeos)
library(ggplot2)
library(dplyr)

### Get data
iran<- getData("GADM", country = "Iran", level = 1)

### SPDF to DF
map <- fortify(iran)

map$id <- as.integer(map$id)

dat <- data.frame(id = 1:(length(iran@data$NAME_1)), state = iran@data$NAME_1)
map_df <- inner_join(map, dat, by = "id")

# Find a center point for each province
centers <- data.frame(gCentroid(iran, byid = TRUE))
centers$state <- dat$state


### This is hrbrmstr's own function
theme_map <- function (base_size = 12, base_family = "") {
theme_gray(base_size = base_size, base_family = base_family) %+replace% 
theme(
axis.line=element_blank(),
axis.text.x=element_blank(),
axis.text.y=element_blank(),
axis.ticks=element_blank(),
axis.ticks.length=unit(0.3, "lines"),
axis.ticks.margin=unit(0.5, "lines"),
axis.title.x=element_blank(),
axis.title.y=element_blank(),
legend.background=element_rect(fill="white", colour=NA),
legend.key=element_rect(colour="white"),
legend.key.size=unit(1.5, "lines"),
legend.position="right",
legend.text=element_text(size=rel(1.2)),
legend.title=element_text(size=rel(1.4), face="bold", hjust=0),
panel.background=element_blank(),
panel.border=element_blank(),
panel.grid.major=element_blank(),
panel.grid.minor=element_blank(),
panel.margin=unit(0, "lines"),
plot.background=element_blank(),
plot.margin=unit(c(1, 1, 0.5, 0.5), "lines"),
plot.title=element_text(size=rel(1.8), face="bold", hjust=0.5),
strip.background=element_rect(fill="grey90", colour="grey50"),
strip.text.x=element_text(size=rel(0.8)),
strip.text.y=element_text(size=rel(0.8), angle=-90) 
)   
}

ggplot() +
geom_map(data = map_df, map = map_df,
         aes(map_id = id, x = long, y = lat, group = group),
         color = "#ffffff", fill = "#bbbbbb", size = 0.25) +
geom_text(data = centers, aes(label = state, x = x, y = y), size = 3) +
coord_map() +
labs(x = "", y = "", title = "Iran Province") +
theme_map()

Answer 2

您可以将 rworldmap 与@jazzurro 提出的使用 raster 获得 GADM 边界的重要建议相结合。

我怀疑您的主要困难可能是让省份名称在您的数据和地图之间匹配。

下面的示例使用您可以更改的默认值，并为每个省提供不同的颜色。

library(raster)
library(rworldmap)
## 1 Get map of provinces (level 1 admin regions)
iranLevel1<- raster::getData("GADM", country = "Iran", level = 1)

## 2 join your [data] onto the map by specifying the join column in each
## this will report on any mis-matches between province names
#iranLevel1 <- rworldmap::joinData2Map([data],nameMap="iranLevel1",nameJoinIDMap="NAME_1",nameJoinColumnData=[insert])

## 3 plot map (change NAME_1 to the data you want to plot)
rworldmap::mapPolys(iranLevel1, nameColumnToPlot="NAME_1", addLegend=FALSE)

## 4 add text labels for provinces
text(iranLevel1, label="NAME_1", cex=0.7)

请注意，joinData2Map()、mapPolys() 是 joinCountryData2Map()、mapCountryData() 的更通用的等价物。

另一种方法是使用 choroplethr 包。

在 R 中创建地图就像 rworldmap 的方式一样，但针对具有省份的特定国家

Creating maps in R just like the way rworldmap does but for specific country with provinces

r

plot

ggplot2

maps

rworldmap