地理区域的R包

Question

有人告诉我一个R包，其中包含一些区域和国家边界的物体定义。我试图找到意大利地区威尼托的边界。包装的名称是什么？你知道另一种方法来找到边界（我不确定威尼托是否包含在包中）？

Answer 1

Howabout：

require(raster) 
veneto = subset(getData('GADM', country='ITA', level=1), NAME_1=="Veneto") 
plot(veneto) 

Answer 2

您也可以与那些流行与D3人群TopoJSON/GeoJSON的文件的工作：

library(rgdal) 
library(RCurl) 
library(ggplot2) 
library(rgeos) 

# for theme_map() 
devtools::source_gist("https://gist.github.com/hrbrmstr/33baa3a79c5cfef0f6df") 

# location of Italy TopoJSON file 
url <- "https://gist.githubusercontent.com/andreaangeli/5b64a102e357198780e9/raw/e076b45eb2fcb8de6a0ec612f863499181bb494b/ita.json" 

download.file(url, destfile="ita.json", method="curl") 

ogrListLayers("ita.json") # to see that the layer is "reg2011" 

ita <- readOGR("ita.json", "reg2011") 

veneto <- ita[ita$id == 5,] # region 5 is Veneto 

# could just plot(veneto) now, but I like ggplot better 

veneto_map <- fortify(veneto, region="id") 

gg <- ggplot(data=veneto_map) 
gg <- gg + geom_map(map=veneto_map, aes(map_id=id, group=group, x=long, y=lat), color="black", fill="white") 
gg <- gg + coord_map() 
gg <- gg + theme_map() 
gg 

你也可以与“正常”形状文件一起工作，这是一个（大多数可能）有什么在raster包Spacedman的帖子引用：

download.file("http://biogeo.ucdavis.edu/data/diva/adm/ITA_adm.zip", "ITA_adm.zip") 
unzip("ITA_adm.zip", exdir="ITA_adm") 

ita <- readOGR("ITA_adm/", "ITA_adm1") 

veneto <- ita[ita$NAME_1 == "Veneto",] # this shapefile has region names in it! 

veneto_map <- fortify(veneto) 

gg <- ggplot(data=veneto_map) 
gg <- gg + geom_map(map=veneto_map, aes(map_id=id, group=group, x=long, y=lat), color="black", fill="white") 
gg <- gg + coord_map() 
gg <- gg + theme_map() 
gg