# NOT RUN {
data(reunion)
# Smoothing all variables for Reunion (Lissage de toutes les variables pour la Reunion)
# Call mode 1: classic smoothing - automatic grid
reunionSmooth <- kernelSmoothing(dfObservations = reunion, cellSize = 200
, bandwidth = 400)
# Call mode 2: median smoothing - automatic grid
reunionSmooth <- kernelSmoothing(dfObservations = reunion, cellSize = 200
, bandwidth = 400, vQuantiles = c(0.1, 0.5, 0.9))
# Call mode 3: classic smoothing - user grid
dfCentroidsUser <- merge( x = seq(from = 314400, to = 378800, by = 200)
, y = seq(from = 7634000, to = 7691200, by = 200))
reunionSmooth <- kernelSmoothing(dfObservations = reunion, cellSize = 200
, bandwidth = 400, dfCentroids = dfCentroidsUser)
# Call mode 4: median smoothing - user grid
reunionSmooth <- kernelSmoothing(dfObservations = reunion, cellSize = 200
, bandwidth = 400, vQuantiles = c(0.1, 0.5, 0.9)
, dfCentroids = dfCentroidsUser)
# Building of the associated basemap (Creation du fond de carte associe)
grid <- smoothingToGrid(reunionSmooth, "32740")
# preview (Apercu)
library(sp)
library(cartography)
cartography::choroLayer(grid, var = "houhold", nclass = 5
, method = "fisher-jenks", border = NA)
# Export of the basemap in shapefile format (Export du fond de carte au format shapefile)
library(rgdal)
rgdal::writeOGR(grid, "reunion.shp", "reunion", driver = "ESRI Shapefile")
# }
Run the code above in your browser using DataLab