Adding terrain layer about deckgl HOT 2 CLOSED

token = Sys.getenv("MAPBOX_TOKEN")
terrain_img = "terrain-rgb"
surface_img = "satellite"

elev_map = list(
  rScaler = 6553.6,
  gScaler = 25.6,
  bScaler = 0.1,
  offset = -10000  
)

props = list(
      minZoom = 0,
      maxZoom = 23,
      strategy = "no-overlap",
      elevationDecoder = elev_map,
      elevationData = sprintf(
        "https://api.mapbox.com/v4/mapbox.%s/{z}/{x}/{y}.png?access_token=%s",
        terrain_img, token),
      texture = sprintf(
        "https://api.mapbox.com/v4/mapbox.%s/{z}/{x}/{y}.png?access_token=%s",
        surface_img, token),
      wireframe = FALSE,
      color = c(255, 255, 255)
)

deckgl(
  element_id = "mapgl",
  latitude = 13,
  longitude =-6,
  zoom = 12,
  bearing = 140,
  pitch = 60,
  maxPitch = 89,
  height = "360px"
) %>%
  add_layer("TerrainLayer",
    id = "terrain",
    properties = props)

library(marmap)

# noaa bathy
b <- getNOAA.bathy(lon1 = -74.2, lon2 =  -73.3, lat1=40.3,
                   lat2 = 40.7, resolution = 0.1)

# convert to raster and save as png?
r <- as.raster(b)
png(file = "output/bathy_grey.png")
raster::plot(r, maxpixels=ncell(r), col=grey.colors(255), 
             axes=FALSE, box=FALSE, legend=FALSE)
dev.off()

elevation_data <- normalizePath( "output/bathy_grey.png")
bounds =  c(r@extent@xmin, r@extent@ymin, r@extent@xmax, r@extent@ymax)

elevation <-  sprintf("data:image/png;base64,%s", base64enc::base64encode(elevation_data))

elev_map = list(
  rScaler = 1,
  gScaler = 0,
  bScaler = 0,
  offset = 0  
)

props = list(
  strategy = "no-overlap",
  elevationDecoder = elev_map,
  elevationData = elevation,
  bounds = bounds 
)

deckgl(
  latitude = 40.5,
  longitude = -74,
  zoom = 10,
  pitch = 45,
  height = "600px") |> 
  add_basemap() |> 
  add_layer("TerrainLayer",
            id = "terrain",
            properties = props)