Last data update: 2014.03.03
|
R: Data from the Netherlands
netherlands | R Documentation |
Data from the Netherlands
Description
Elevation data and map tiles for the Netherlands
Usage
data("netherlands")
Format
nldElev is a raster of elevation
nltTiles is a background map
nldCities is a SpatialPointsDataFrame of city locations.
Details
The inclusion of these datasets is intended to allow the package to build
when an internet connection is not present.
Source
See examples.
See Also
meuse, getData, openmap
Examples
# soil data
library("sp")
data("meuse")
coordinates(meuse) <- ~x+y
if(require('rgdal', quietly=TRUE)) {
proj4string(meuse) <- CRS("+init=epsg:28992")
} else {
proj4string(meuse) <- CRS(
paste("+proj=sterea +lat_0=52.15616055555555 +lon_0=5.38763888888889",
"+k=0.9999079 +x_0=155000 +y_0=463000 +ellps=bessel +units=m +no_defs"
)
)
}
meuse$soilFac = factor(meuse$soil, levels=c(1,2,3),
labels=c("Calcareous","Non-Calc's","Red Brick"))
soilCol = colourScale(meuse$soilFac)
data("netherlands")
map.new(meuse)
plot(nldTiles,add=TRUE)
points(nldCities)
text(nldCities,label=nldCities$name, pos=2)
points(meuse, pch=16, col=soilCol$plot)
legend('topleft', fill=soilCol$col,legend=soilCol$legend)
insetMap(meuse, "bottomright",map=world)
# location won't be marked on the inset map unless rgdal is available
## Not run:
# this is how the data were obtained
# map tiles
nldTiles = openmap(meuse, zoom=12)
# cities
nldCities = GNcities(nldTiles, maxRows=25)
# world
world = openmap(extent(-10,30,40,60))
# elevation data
require('rgdal')
meuseLL = spTransform(meuse, CRS("+init=epsg:4326"))
getData("SRTM", lon=xmin(extent(meuseLL)),
lat=ymin(extent(meuseLL)),path=tempdir())
nldElev = raster(paste(tempdir(), "/", "srtm_38_02.tif", sep=""))
nldElev = crop(nldElev, extend(extent(meuseLL), 0.1))
nldElev = projectRaster(nldElev, crs=proj4string(meuse))
nldElev = crop(nldElev, extent(nldTiles))
# save the files where the package builder wants them
# save(nldElev, nldTiles, nldCities,world,
# file="~/workspace/diseasemapping/pkg/mapmisc/data/netherlands.RData",
# compress="xz")
## End(Not run)
Results
|