data-raw/mapdeck/mapdeck_examples.R

In hypertidy/rangl: Mesh Topology and Visualization for Spatial Data

library(mapdeck)
mapdeck() %>% add_mesh(data = melbourne_mesh)
library(anglr)
library(ceramic)
cc <- cc_elevation(raster::extent(147.5, 147.6, -42.1, -42.0),
                   zoom = 9)

mesh <- reproj::reproj(quadmesh::triangmesh(cc), 4326)
mesh$vb[3, ] <- mesh$vb[3, ] * 35
mesh$material$col <- colourvalues::color_values(mesh$vb[3, mesh$it])
mapdeck() %>% add_mesh(data = mesh)


## from dev anglr
topo_quad <- as.mesh3d(raster::crop(quadmesh::etopo, raster::extent(100, 180, -50, -20)))
#rgl::clear3d();rgl::shade3d(topo_quad);rgl::aspect3d(1, 1, 0.01);  rgl::rglwidget()
topo_triangle <- as.mesh3d(raster::crop(quadmesh::etopo, raster::extent(100, 180, -50, -20)),
                           triangles = TRUE)
#rgl::clear3d();rgl::shade3d(topo_triangle);rgl::aspect3d(1, 1, 0.01);  rgl::rglwidget()
library(sf)
poly <- st_transform(subset(silicate::inlandwaters, Province %in% c("Victoria", "South Australia")),
                     "+proj=longlat") ## reproj via anglr doesn't work atm because WKT2
poly_triangle <- as.mesh3d(silicate::TRI0(poly))
poly_triangle_dense <- as.mesh3d(
  copy_down(anglr::DEL(poly, max_area =.005), quadmesh::etopo)
)
#rgl::clear3d();rgl::shade3d(poly_triangle_dense);rgl::aspect3d(1, 1, 0.3)

dodeca <- rgl::scale3d(rgl::translate3d(rgl::dodecahedron3d(), 147, -42, 1),
                     1, 1, 1e5)
x <- silicate::inlandwaters %>%
  dplyr::filter(Province %in%
                  c("Victoria", "South Australia", "New South Wales"))
## we need longlat
x <- sf::st_transform(x, "+proj=longlat")


## polygons + elevation + colours
x
z <- quadmesh::etopo
library(sf)
library(anglr)
## put colours on these polygons
x$color_ <- c("firebrick", "dodgerblue", "darkorange1")
poly_triangle_colours <- as.mesh3d(copy_down(
                  anglr::DEL(x, max_area = 1e9), z))
plot3d(poly_triangle_colours);


library(mapdeck)
mapdeck() %>% add_mesh(data = poly_triangle_colours)
## we preserve the colours
rgl::clear3d(); plot3d(poly_triangle_colours); rgl::aspect3d(1, 1, 0.1)
mesh_plot(poly_triangle_colours)


saveRDS(list(topo_quad = topo_quad,
             topo_triangle = topo_triangle,
             poly_triangle = poly_triangle,
             poly_triangle_dense = poly_triangle_dense,
             poly_triangle_colours = poly_triangle_colours,
             dodeca_triangle = dodeca),
        file = "data-raw/mapdeck/mapdeck_mesh3d_examples.rds", compress = "xz")


#
# md <- readRDS("data-raw/mapdeck_mesh3d_examples.rds")
# library(mapdeck)
# ptd <- md$poly_triangle_dense
# ptd$vb[3, ] <- ptd$vb[3, ] * 95
# dodeca$vb[3, ]
# mapdeck() %>% add_mesh(data = dodeca)