How to create an interactive map of stars/raster object with arbitrary crs?

Question

Problem

I would like to create interactive maps of a regular grid that comes in a CRS which is not WGS84. In my case the CRS is EPSG:25833 which has unit metres.

Given the comment here, it does not seem to be possible to "display rasters with arbitrary crs".

Is this still the case? Or is there some kind of hack/solution now using packages like plainview, leaflet, mapview or tmap (or some other package I'm not aware of)?

Moreover, I do not understand what exactly mapview::mapview shows when I feed it a non-lat/lon stars object. See the following example, where I suppose that the sf object is rendered correctly whereas the stars object is not?

Would tmap::tmap_leaflet do a better job here?

Example

library(dplyr)
library(mapview)
library(sf)
library(raster)
library(stars)
library(tmap)

stars_obj <- raster(
  # Rough bbox of the city of Berlin, Germany
  xmn=360000, xmx=430000, 
  ymn=5790000, ymx=5850000,
  # CRS with units metres, NOT lat/lon
  crs = st_crs(25833)$wkt,
  # Grid cell size: 10km x 10km
  resolution = 10000) %>% 
  st_as_stars()
stars_obj <- stars_obj %>% 
  mutate(layer = 1:ncell(stars_obj))

sf_obj <- st_as_sf(stars_obj)

# mapview
mapview(stars_obj) + 
  mapview(sf_obj)

# tmap
tmap_leaflet(
  tm_shape(stars_obj) +
  tm_raster(alpha = 0.8)
)

This is how mapview::mapview renders it:

This is how tmap::tmap_leaflet renders it: