Oefening.R

library(rgdal)
library(rgeos)

URL <- "http://www.twiav.nl/files/NL_Gemeenten2014.zip"
dir.create("./Data", showWarnings = FALSE)
download.file(URL, destfile ="./Data/NL_Gemeenten2014.zip", mode = "wb")
unzip("./Data/NL_Gemeenten2014.zip", exdir = "./Data")
unlink("./Data/NL_Gemeenten2014.zip")

gem2014 <- readOGR("Data/NL_Gemeenten2014.TAB", "NL_Gemeenten2014")

# Toon de namen van de kolommen in gem2014
names(gem2014)

# Toon de eerste 25 waarden in de kolom Gemeentenaam
head(gem2014$Gemeentenaam, n = 25)

# Toon de eerste 6 rijen in de attribuuttabel van gem2014
head(gem2014@data)

# Hoeveel gemeenten hebben minder dan 20.000 inwoners?
sum(gem2014$Inwoners < 20000)

# Waar liggen de gemeenten met minder dan 20.000 inwoners?
plot(gem2014)
plot(gem2014[gem2014$Inwoners < 20000,], add = TRUE, col = "red")

# Welke gemeente heeft het minste aantal inwoners per vierkante kilometer?
gem2014$Oppervlakte = gArea(gem2014, byid = TRUE) / 1000^2
gem2014$Dichtheid = gem2014$Inwoners /gem2014$Oppervlakte
gem2014$Gemeentenaam[gem2014$Dichtheid == min(gem2014$Dichtheid)]

# Toon het aantal gemeenten per veiligheidsregio
table(gem2014$Veiligheidsregio_naam)

# Toon een kaart en de namen van de veiligheidsregio's in Nederland
vregios <- gUnaryUnion(gem2014, id = gem2014$Veiligheidsregio_naam)
plot(vregios)
row.names(vregios)