In hypertidy/cartilage: Spatializing 'rayshader'

knitr::opts_chunk$set(
  collapse = TRUE,
  comment = "#>",
  fig.path = "man/figures/README-",
  out.width = "100%"
)

cartilage

The goal of cartilage is to blend rayshader image forms with absolute coordinates in rgl.

This is very much a work in progress. See functions sphere, ray, lamb and ambient, which expect a raster and provide the same interfaces as the corresponding rayshader *_shade functions. Be careful not to input very large data, keep to something like 300*300 or so unless you outlandishly keen (we'll add helpers).

The output of functions is the data used to build the scene, and currently includes the quadmesh, the hillshade and the raster - but we don't have helpers for those yet. An internal function plots and constructs this list output.

Feedback welcome!

Installation

Install from github with

devtools::install_github("hypertidy/cartilage")

Example

We use the built-in GeoTIFF file of a reduced resolution Monterey Bay elevation.

f <- system.file("extdata/monterey/monterey_dem.tif", package= "cartilage")

r <- raster::raster(f)

library(cartilage)
x <- sphere(r)
rgl::axes3d()

Now in in x$qm we have the rgl mesh3d form, so we can do some tricks:

## in x$qm we have the rgl mesh3d form, so we can do some tricks
qm <- x$quadmesh  ## convert to Geocentric coordinates
qm$vb[1:3, ] <- t(proj4::ptransform(t(qm$vb[1:3, ]), raster::projection(r), "+proj=geocent +datum=WGS84"))
rgl::rgl.clear()  
rgl::shade3d(qm)  ## we are now geocentric
library(graticule)
g <- graticule(seq(-123, -120, by = 0.25), seq(35, 38, by = .2))
g <- silicate::SC(g)  ## hypertidy/silicate
library(anglr)  ## hypertidy/anglr
plot3d(anglr::globe(g), add = TRUE)
rgl::aspect3d(1, 1, 1)

Getting raster data ...

(... at reasonable resolution)

The 'raster' package has a pretty efficient 'aggregate' function, and 'rgdal' can read data at lower resolution than native using the 'output.dim' argument of 'readGDAL' - the stars project aims to supersede both packages and provide the raster sf. Our favourite method currently is to use lazyraster which can read directly from a file via GDAL into a raster of a desired resolution, with helpers for setting an extent prior to reading any data.

The output of cartilage functions is the data used to build the scene, but we don't have helpers for those yet. A related idea is being pursued in SOmap for Antarctica maps.

Please note that the 'cartilage' project is released with a Contributor Code of Conduct. By contributing to this project, you agree to abide by its terms.

'cartilage' was originally named 'montereybayshader', as it was motivated by getting a shareable data set and settled on Monterey Bay. The original code was just to plot rayshader products in absolute coordinates, or to mesh other products into a rayshader scene. Some older notes are in data-raw/.

hypertidy/cartilage documentation built on June 30, 2021, 10:06 p.m.