Leaflet muddling RGB rasters from terra

[jebyrnes]

I'm working with an RGB raster in leaflet (some Landsat data). I'm stretching and sampling it down so that it's small enough. The following works great.

# stretch and then resample like geom_spatraster
landsat_rgb_stretched <- stretch(landsat_colors_1230[[1:3]],
                                 minq = 0.02, maxq = 0.98) |>
  spatSample(size = 5e5, as.raster = TRUE, method = "regular")

# set RGB channels
RGB(landsat_rgb_stretched) <- 3:1

# plot
plotRGB(landsat_rgb_stretched)

But, when I try and put it into leaflet, I get a muddled mess. It's close, but, very weird. I'm trying to figure out why.

# interactive plot
leaflet() |>
  addRasterImage(landsat_rgb_stretched)

[aloboa]

I have a similar problem using an slightly different approach:

r <- subset(rast("enmapL2A_199-114-037_gmap.tif"),1:3)
RGB(r) <- 1:3
r[r==0] <- NA
a <- colorize(r,"col")
plot(a)

Then I use
addRasterImage(a, group="a")

And the colors and grain are much worse than in the R plot() (left):

Instead, the display with 1 single band is ok:

r1 <- r[[1]]

addRasterImage(r1, group="r1")

(do not manage to get the same color table, but the leaflet display on the right is not grainy as in the previous case)