Difference between revisions of "Population of Finland"

library(OpasnetUtils)

data <- tidy(op_baseGetData("opasnet_base", "Op_en2949"), objname = "population")
population <- new("ovariable", name = "population", data = data)
print(population@data)

objects.put(population)
cat("Object population initiated. Save the page key for further use.\n")

# (C) 2011 Leo Lahti <leo.lahti[at]iki.fi> All rights reserved. 
# License: FreeBSD, http://en.wikipedia.org/wiki/BSD_licenses

# Tama esimerkki on testattu sorvi-paketin versiolla 0.1.23

# Esimerkki Suomen kuntatason vaestonkasvutilastojen (Tilastokeskus)
# visualisoinnista Maanmittauslaitoksen karttadatalla (vuonna 2010)

# Lataa kirjastoja
library(sorvi)
library(OpasnetBaseUtils)

###############################################

# Lue Suomen kuntarajat SpatialPolygon-muodossa
# (C) Maanmittauslaitos 2011
# http://www.maanmittauslaitos.fi/aineistot-palvelut/digitaaliset-tuotteet/ilmaiset-aineistot/hankinta
data(MML)
sp <- MML[["1_milj_Shape_etrs_shape"]][["kunta1_p"]]

#################################################

# Lue kuntatason vaestonkasvutiedot tilastokeskuksen StatFin-tietokannasta
# http://www.stat.fi/tup/statfin/index.html
# PC Axis-muodossa
px <- statfin.px("vrm/synt/080_synt_tau_203_fi.px")
# statfin.px-funktio käyttää soRvin read.px-funktiota hakemaan dataa URLista, joka alkaa "http://pxweb2.stat.fi/database/StatFin/"
# URLin loppuosa pitää antaa statfin.px:lle parametrina.
# Poimi taulukosta halutut tiedot
# pxs <- subset(as.data.frame(px), Väestönmuutos.ja.väkiluku == "Luonnollinen väestönlisäys" & Vuosi == year)
pxs <- subset(as.data.frame(px), Väestönmuutos.ja.väkiluku == parameter & Vuosi == year)
# Putsaa data
vaestonkasvu <- preprocess.px(pxs)

################################################

# Lisaa tiedot karttaobjektiin
sp@data$vaestonkasvu <- vaestonkasvu$dat[match(sp$Kunta.FI, vaestonkasvu$Alue)]
# Korvaa puuttuvat arvot nollalla
sp[["vaestonkasvu"]][is.na(sp[["vaestonkasvu"]])] <- 0

################################################

# Maarittele varipaletti
my.palette <- colorRampPalette(c("blue", "white", "red"), space = "rgb")
ncol <- 40 # Number of colors

#################################################

# Piirra kuva
varname <- "vaestonkasvu"
int <- max(abs(sp[[varname]]))
q <- spplot(sp, varname,
       col.regions = my.palette(ncol),
       main = paste(parameter, year, sep=" "),
       colorkey = TRUE, 
       lwd = .4,
       col = "black", 
       at = seq(0 - int, 0 + int, length = ncol)
)

print(q)