印度地图 R 中的印度每个州的五个指标的不同值答案

【问题标题】：Map of India in R with Different Values of Five Indicators for Every State of India印度地图 R 中的印度每个州的五个指标的不同值
【发布时间】：2015-04-25 07:26:05
【问题描述】：

我想在 R 中生成印度地图。我有五个指标，每个州都有不同的值。我想用五种不同的颜色绘制气泡，它们的大小应该代表它们在每个状态下的强度。例如：

State A B C D E Kerala - 39, 5, 34, 29, 11 Bihar - 6, 54, 13, 63, 81 Assam - 55, 498, 89, 15, 48, Chandigarh - 66, 11, 44, 33, 71

我浏览了一些与我的问题相关的链接：

[1]http://www.r-bloggers.com/nrega-and-indian-maps-in-r/

[2]An R package for India?

但这些链接无法达到我的目的。非常感谢这方面的任何帮助。

我也试过了

library(mapproj)
map(database= "world", regions  = "India", exact=T, col="grey80", fill=TRUE, projection="gilbert", orientation= c(90,0,90))

lat <- c(23.30, 28.38)

lon <- c(80, 77.12) # Lon and Lat for two cities Bhopal and Delhi

coord <- mapproject(lon, lat, proj="gilbert", orientation=c(90, 0, 90))

points(coord, pch=20, cex=1.2, col="red")

简而言之，问题是：（1）它没有给我地区一级的情节。甚至没有国家的边界。 (2) 如果我只有位置名称和要绘制的相应值，如何在此图中创建数据的气泡或点？ (3) 这可以在library(RgoogleMaps) 或library(ggplot2) 中轻松完成吗？（只是猜测，我对这些包了解不多）

【问题讨论】：

真的需要在每个状态下添加五个气泡吗？这似乎占用了相当多的空间；您可能没有足够的空间为某些小州添加五个气泡。我想知道您是否想要五张分开的印度地图。
你能提供你创建印度地图的代码吗？
亲爱的 lawyeR & jazzurro，我遵循了上面链接中的代码，因为我对 R 中的地理地图一无所知。但我找不到这些代码对我有用。您说得对，地图的空间会更小，特别是对于较小的州。我将尝试使用三个指标进行管理。
您链接的答案已过期，因为 GADM 链接已更改。我已经编辑了它们（编辑应该很快就会生效），这应该可以帮助您开始绘制印度地图：r-fiddle.org/#/fiddle?id=zJd1RfjS
@Molx：感谢您所做的所有工作以及 r-fiddle 上的代码。最后一行返回错误：> IND

标签： r plot ggplot2 gis rgooglemaps

【解决方案1】：

正如@lawyeR 所说，等值线（或专题）地图更常用于表示地图上的变量。这将要求您为每个变量生成一张地图。让我举个例子：

require("rgdal")  # needed to load shapefiles

# obtain India administrative shapefiles and unzip
download.file("http://biogeo.ucdavis.edu/data/diva/adm/IND_adm.zip", 
              destfile = "IND_adm.zip")
unzip("IND_adm.zip", overwrite = TRUE)

# load shapefiles
india <- readOGR(dsn = "shapes/", "IND_adm1")

# check they've loaded correctly with a plot
plot(india)

# all fine. Let's plot an example variable using ggplot2
require("ggplot2")
require("rgeos")  # for fortify() with SpatialPolygonsDataFrame types

india@data$test <- sample(65000:200000000, size = nrow(india@data),
                          replace = TRUE)

# breaks the shapefile down to points for compatibility with ggplot2
indiaF <- fortify(india, region = "ID_1")
indiaF <- merge(indiaF, india, by.x = "id", by.y = "ID_1")

# plots the polygon and fills them with the value of 'test'
ggplot() +
  geom_polygon(data = indiaF, aes(x = long, y = lat, group = group,
                                  fill = test)) +
  coord_equal()

最后，我注意到您在 GIS SE 上提出了同样的问题。这被认为是不好的做法，并且通常不受欢迎，因此我已将该问题标记为已关闭作为此问题的副本。作为一般经验法则，尽量不要创建重复项。

祝你好运！

【讨论】：

对不起！我不知道这个。我将来会处理这个问题。谢谢菲尔。
没关系，这需要一段时间才能掌握，所以我希望能有所帮助而不是评判:)
感谢您的出色回答。请注意，如果您遇到 readOGR 问题，india <- readOGR(dsn = "shapes/", "IND_adm1") 需要目录shapes 中的IND_adm1.shp 文件。我的默认解压缩位置不是shapes/，我手动将提取的文件放在shapes 文件夹中。我还必须将dsn = "shapes/" 更改为dsn = "shapes"。
如果您只想要状态数据（最有用的情况 ig），您可以删除以india@..（无论如何仅用于测试）和indiaF <- merge(.. 开头的行。现在检查indiaF，您将拥有名为id 的列，它们是各州的ID 号。您可以查看IND_adm1.csv 以查看哪个 id 属于哪个状态。如果您删除了test，只需将fill = test 更改为fill = id，即可看到每个状态的颜色不同。现在您可以通过更改fill = <your vairable> 根据您的数值变量为它们着色。使用RColorBrewer 将是等值线图的最佳选择。

【解决方案2】：

获得印度的 shapefile 后，您需要创建一个 choropleth。这将采用 shapefule 地图，并在反映您的数据的渐变上为印度的每个州着色。您可能想要创建一个包含五个图的面板，每个图根据您的五个变量之一显示印度及其各州。

对于其他可以将这个答案推得更远的人，经过一番清理，这是数据框的dput。

dput(df)
structure(list(State = c("Kerala", "Bihar", "Assam", "Chandigarh"
), A = c("39", "6", "55", "66"), B = c("5", "54", "498", "11"
), C = c("34", "13", "89", "44"), D = c("29", "63", "15", "33"
), E = c("11", "81", "48", "71")), .Names = c("State", "A", "B", 
"C", "D", "E"), row.names = c("Kerala", "Bihar", "Assam", "Chandigarh"
), class = "data.frame")

【讨论】：