如何在一个国家/地区绘制区域，每个区域都有一些相应的值

Question

特别是我在印度有一系列不同的地区，每个地区都有一个＆＃34; markethare＆＃34;值。

regions <- data.frame(Region=c("TELANGANA", "TAMIL NADU", "JAMMU & KASHMIR"), Marketshare=c(.25, .30, .15))
regions
           Region Marketshare
1       TELANGANA        0.25
2      TAMIL NADU        0.30
3 JAMMU & KASHMIR        0.15

我可以使用

轻松地使用ggmap真实地绘制印度

library(ggmap)
map <- get_map("India", zoom=6, maptype="toner-lite")
ggmap(map)

但是，我如何才能在我的数据集中专门包含这些区域，并使用 Marketshare 列对每个区域进行着色？

注意：我已经看过this question and answer，但对我的需求来说似乎过于复杂，并且想知道是否有更简单的方法。

Answer 1

正如hrbrmstr所说，我认为如果你能表现出一些努力会很棒。在这里，我创建了两个地图。第一个是ggmap，另一个是ggplot。在示例代码中，您的地图缩放已关闭，以显示您指定的所有三个区域。所以我修改了你的代码。您需要的地图顶部是您想要的区域的多边形。要获取它们，请使用get_data()下载数据。然后，您需要对其进行子集化。如链接问题所示，您需要完成合并过程才能将marketshare添加到最终数据框。我收集了每个地区的ID，并创建了一个marketshare的数据框，即foo。 temp是包含多边形信息以及marketshare的最终数据框。您使用marketshare作为fill并绘制地图。

library(raster)
library(ggmap)
library(ggplot2)

### Get India data
india <- getData("GADM", country = "India", level = 2)

map <- get_map("India", zoom = 4, maptype = "toner-lite")
regions <- data.frame(Region=c("Telangana", "Tamil Nadu", "Jammu and Kashmir"), Marketshare=c(.25, .30, .15))

states <- subset(india, NAME_1 %in% regions$Region)

ind1 <- states$ID_2[states$NAME_1 == "Telangana"]
ind2 <- states$ID_2[states$NAME_1 == "Tamil Nadu"]
ind3 <- states$ID_2[states$NAME_1 == "Jammu and Kashmir"]

states <- fortify(states)

foo <- data.frame(id = c(ind1, ind2, ind3),
                  marketshare = rep(regions$Marketshare, times = c(length(ind1), length(ind2), length(ind3))))

temp <- merge(states, foo, by = "id")

ggmap(map) +
geom_map(data = temp, map = temp,
         aes(x = long, y = lat, map_id = id, group = group, fill = marketshare),
         colour = "grey50", size = 0.3) +
theme(legend.position = "none")

或者，您可以使用ggplot绘制地图，以便仅显示印度。在这种情况下，您可以从下载数据india创建数据框。您绘制一次印度，然后在地图上绘制特定区域。我在这里将marketshare转换为颜色因素。

india.map <- fortify(india)

ggplot() +
geom_map(data = india.map, map = india.map,
         aes(x = long, y = lat, map_id = id, group = group),
         colour = "grey50", size = 0.3) +
geom_map(data = temp, map = temp,
         aes(x = long, y = lat, map_id = id, group = group, fill = factor(marketshare)),
         colour = "grey50", size = 0.3) +
theme(legend.position = "none")