OpenStreetMap (version 0.3.2)

openproj: Projects the open street map to an alternate coordinate system

Description

Projects the open street map to an alternate coordinate system

Usage

openproj(x, projection = "+proj=longlat", ...)

Arguments

x

an OpenStreetMap object

projection

a proj4 character string or CRS object

...

additional parameters for projectRaster

Examples

## Not run: 
# library(maps)
# 
# #plot mapquest map in native mercator coords
# map <- openmap(c(70,-179),
# 		c(-70,179),zoom=1,type='mapquest-aerial')
# plot(map)
# 
# #using longlat projection lets us combine with the maps library
# map_longlat <- openproj(map)
# plot(map_longlat)
# map("world",col="red",add=TRUE)
# 
# #robinson projection. good for whole globe viewing.
# map_robinson <- openproj(map_longlat, projection=
#   "+proj=robin +lon_0=0 +x_0=0 +y_0=0 +ellps=WGS84 +datum=WGS84 +units=m +no_defs")
# plot(map_robinson)
# 
# #national parks service images
# upperMap <- openmap(c(70,-179),
# 		c(10,50),zoom=2,type='nps')
# #Lambert Conic Conformal
# map_llc <- openproj(upperMap, projection=
# 				"+proj=lcc +lat_1=33 +lat_2=45 +lat_0=39 +lon_0=-96")
# plot(map_llc,removeMargin=TRUE)
# #add choropleth
# data(states)
# library(sp)
# st_llc <- spTransform(states,CRS("+proj=lcc +lat_1=33 +lat_2=45 +lat_0=39 +lon_0=-96"))
# plot(st_llc,add=T,col=heat.colors(48,.4)[slot(st_llc,"data")[["ORDER_ADM"]]])
# 
# ## End(Not run)