[英]Cartogram + choropleth map in R
我最近一直在使用ggplot2来创建一串choropleths。 我想知道是否可以使用ggplot2创建与此类似的地图(来自WorldMapper ):
这是一个choropleth,其中shapefile多边形变形以表示相对人口数。 我相信这称为制图。 他们使用许多其他变量来执行此操作。 本着Choropleth R Challenge的精神,有人知道如何使用R做到这一点吗?
CRAN上提供的cartogram
包 ,具有您想要的橡胶板变形样式的地图。
这可能起作用:
您将需要预安装FFTW。 Rcartogram
和getcartr
您将需要devtools
。
不确定如何在ggplot2
执行此ggplot2
,但这是另一种选择。
在这里,我使用的是Thematic World Map中的shapefile,下载并解压缩后,将得到一个名为TM_WORLD_BORDERS-0.3
的文件夹。
对于Choropleth / Cartogram,您首先要使用大小重塑形状,然后使用特征进行阴影处理:
library(rgdal)#needed for readOGR
library(sp) #needed for spplot
library(Rcartogram)
library(getcartr)
setwd("<your_directory_with_shapefile>") #to the file that has your shapefile and your information file (in this case, a csv named datR)
#read shapefile
#here i have a folder with a shapefile and a csv with columns as ISO (IS02 for convenience) country and value
worldR <- readOGR(dsn = getwd(), layer= "TM_WORLD_BORDERS-0.3") # If reading a shapefile, the data source name (dsn= argument) is the folder (directory) where the shapefile is, and the layer is the name of the shapefile (without the .shp extension)
#names(worldR) #note how here there are columns for ISO2 (which matches a column named 'iso' in datR and LAT\LON
#[1] "FIPS" "ISO2" "ISO3" "UN" "NAME" "AREA" "POP2005" "REGION" "SUBREGION" "LON" "LAT"
proj4string(worldR)
datR <- read.csv("datR.csv") #this is a file that has one column called 'score' and one column called size':
head(datR)
# iso size score
#1 AE 323 0.9819077
#2 AR 262 0.9591067
#3 AT 7481 0.9987313
#4 AU 5425 0.9837414
#5 BA 31 0.9871938
#6 BB 99 0.9715991
##Merge SpatialPolygonsDataFrame with other info
map_dat <- merge(worldR, datR, by.x="ISO2",by.y="iso")
#remove coordinate reference system arguments
proj4string(map_dat) <- CRS(as.character(NA)) # from here https://github.com/chrisbrunsdon/getcartr/issues/1
world.carto <- quick.carto(map_dat, map_dat$size, blur = 0)
#plot(world.carto) #cartogram without anything
#spplot size, color
my.palette = c("#ff0000", "#ff8000", "#ffff00", "#bfff00","#00ff00") #red, orange, yellow, light green, dark green
spplot(world.carto, 'score', col.regions = my.palette, cuts = length(my.palette)-1,main="Choropleth of score and cartogram of size")
这应该为您提供类似于此的情节:
我赶紧做,让我知道它是否有效
声明:本站的技术帖子网页,遵循CC BY-SA 4.0协议,如果您需要转载,请注明本站网址或者原文地址。任何问题请咨询:yoyou2525@163.com.