地理區域的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