mapmisc (version 1.1.1)

netherlands: 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.

source

See examples.

Details

The inclusion of these datasets is intended to allow the package to build when an internet connection is not present.

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)
legendBreaks('topleft', soilCol)
insetMap(meuse, "bottomright",map=world)

# location won't be marked on the inset map unless rgdal is available

# 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")