raster-particle.R

library(mapgl)

mapboxgl(
  zoom = 2,
  center = c(-28, 47),
  style = mapbox_style("dark")
) |> 
  add_layer(
    id = "wind-layer", 
    type = "raster-particle",
    source = list(
      type = "raster-array",
      url = "mapbox://rasterarrayexamples.gfs-winds",
      tileSize = 512
    ),
    source_layer = "10winds", 
    paint = list(
      `raster-particle-speed-factor` = 0.4,
      `raster-particle-fade-opacity-factor` = 0.9,
      `raster-particle-reset-rate-factor` = 0.4,
      `raster-particle-count` = 4000,
      `raster-particle-max-speed` = 40,
      `raster-particle-color` = interpolate(
        property = "raster-particle-speed",
        type = "linear",
        values = c(1.5, 2.5, 4.12, 6.17, 9.26, 11.83, 14.92, 18.0, 21.6, 25.21, 29.32, 33.44, 43.72, 50.41, 59.16, 69.44),
        stops = c(
          'rgba(134,163,171,256)',
          'rgba(126,152,188,256)',
          'rgba(110,143,208,256)',
          'rgba(15,147,167,256)',
          'rgba(57,163,57,256)',
          'rgba(194,134,62,256)',
          'rgba(200,66,13,256)',
          'rgba(210,0,50,256)',
          'rgba(175,80,136,256)',
          'rgba(117,74,147,256)',
          'rgba(68,105,141,256)',
          'rgba(194,251,119,256)',
          'rgba(241,255,109,256)',
          'rgba(256,256,256,256)',
          'rgba(0,256,256,256)',
          'rgba(256,37,256,256)'
        )
      )
    ),
  )