Filling map polygons using GGplot in R

Question

For my school assignment, I am trying to make a map of the Netherlands filled with a color scheme, which depends on integer numbers (number of citizens). I have a dataset called mun_neth, which is a SpatialPolygonDataFrame and contains all the polygons of the Netherlands and all the data I want to build. I tried three different methods that I added under it. I also put error messages at the bottom. I think I misunderstand the fill requirement. Where does this go wrong and how do I fix it?

After searching here on stackoverflow, I feel like I'm pretty close to building a map. But, unfortunately, it does not work yet.

# Set workspace
getwd()
setwd("~/Wageningen/2.2 Geoscripting/data")

# Load libraries
install.packages("RCurl", dependencies=TRUE)
library(RCurl)
install.packages("ggplot2", dependencies=TRUE)
library(ggplot2)
install.packages("rgdal", dependencies=TRUE)
library(rgdal)
# Load in data
dl_from_dropbox <- function(x, key) {
  require(RCurl)
  bin <- getBinaryURL(paste0("https://dl.dropboxusercontent.com/s/", key, "/", x),           
                      ssl.verifypeer = FALSE)
  con <- file(x, open = "wb")
  writeBin(bin, con)
  close(con)
  message(noquote(paste(x, "read into", getwd())))                        
}

dl_from_dropbox("Netherlands.zip", "bocfjn1l2yhxzxe")
unzip("Netherlands.zip")
mun_neth <- readOGR('gem_2012_v1.shp', layer = 'gem_2012_v1')

# First attempt
mun_neth <- readOGR('gem_2012_v1.shp', layer = 'gem_2012_v1')
mun_neth@data$id <- rownames( mun_neth@data )
mun_neth.df <- as.data.frame( mun_neth )
mun_neth.fort <- fortify( mun_neth , region = "id" )
mun_neth.gg <- join( mun_neth.fort , mun_neth.df , by = "id" )

ggplot(data = mun_neth, aes(long, lat, group=group)) +
  geom_map(aes(fill = mun_neth$AANT_INW, color = category), map =mun_neth.gg) +
  scale_fill_gradient(high = "red", low = "white", guide = "colorbar")
  scale_colour_hue(h = c(120, 240))

Regions defined for each polygon

Error: cannot allocate 9.5 MB vector

# second attempt
ggplot(mun_neth$AANT_INW, aes(x=T_MEAN))

: ggplot2 , integer

# Third attempt
ggplot(aes(x=x,y=y,fill=AANT_INW),data=mun_neth)

,

:

Answer 1

, , - . , , , .

library(rgdal)
library(ggplot2)

work.dir <- "your_work_dir"
mun.neth <- readOGR(work.dir, layer = 'gem_2012_v1')

mun.neth.fort <- fortify(mun.neth, region = "AANT_INW")
mun.neth.fort$id <- as.numeric(mun.neth.fort$id)
mun.neth.fort$id <- mun.neth.fort$id/1000 # optionally change to thousands?
mun.neth.fort[mun.neth.fort$id <= 0, 'id'] <- NA # some areas not numerically valid,
                                                 # presumably water zones

ggplot(data = mun.neth.fort, aes(x = long, y = lat, fill = id, group = group)) +
    geom_polygon(colour = "black") +
    coord_equal() +
    theme()