Instantly share code, notes, and snippets.

# geotheory/3d_mapping

Last active Jan 30, 2018
Plotting 3D maps using OpenStreetMap and RGL.
 # Plotting 3D maps using OpenStreetMap and RGL. For info see: # http://geotheory.co.uk/blog/2013/04/26/plotting-3d-maps-with-rgl/ map3d <- function(map, ...){ if(length(map\$tiles)!=1){stop("multiple tiles not implemented") } nx = map\$tiles[]\$xres ny = map\$tiles[]\$yres xmin = map\$tiles[]\$bbox\$p1 xmax = map\$tiles[]\$bbox\$p2 ymin = map\$tiles[]\$bbox\$p1 ymax = map\$tiles[]\$bbox\$p2 xc = seq(xmin,xmax,len=ny) yc = seq(ymin,ymax,len=nx) colours = matrix(map\$tiles[]\$colorData,ny,nx) m = matrix(0,ny,nx) surface3d(xc,yc,m,col=colours, ...) } #Sys.setenv(NOAWT=1) # Mac users: fixes an OSM/X11 issue that may arise require(rgl) require(OpenStreetMap) require(ggplot2) require(maptools) # download map tile (the '8' parameter for map resolution) lat <- c(51.7, 51.3); lon <- c(-0.53, 0.3) map <- openproj(openmap(c(lat,lon),c(lat,lon), 8, 'osm')) # import London rents data (originally from London Data Store) rents <- read.csv("http://bit.ly/YZRYEC", header=T) head(rents) # create xyz matrix for line heights (row pairs for segment points) m <- mat.or.vec(66, 3) for(i in 1:66) for(j in 1:3) m[i,j] = rents[ceiling(i/2),j+2] for(i in 1:33) m[i*2,3] = 0; m[,3] = m[,3]/15000 head(m) # draw map run <- function(){ open3d(windowRect=c(100,100,800,600)) map3d(map, lit=F) segments3d(m, lwd=5, col='red', alpha=0.6) play3d(spin3d(axis=c(0,0,0.5), rpm=24), duration=2.5) } run() # adjust to custom view and save to file rgl.snapshot("output.png", fmt="png", top=TRUE)

### jasswitz commented Jan 30, 2018

 Hi there, this post is from a few years ago, but nevertheless interesting for me. however, i cannot find the used rents-file anywhere which makes it hard for me to understand the sought after structure for the matrix that will be plotted. do you still have the file and could put it somewhere to download? or can you add a screenshot of the data structure of the rents file, so one can reconstruct something in a similar structure? thank you so much