Special bindings for stars data structures

A good start would be to replicate the plot.stars() method (from https://github.com/r-spatial/stars)

library(stars)
etm <- st_stars(system.file("tif/L7_ETMs.tif", package = "stars"))
plot(etm)

So plot_ly(etm) would do something similar to this (automatically)

library(sf)
etm <- st_transform(etm, 3857)
band <- st_dimensions(etm)[["band"]]
n <- band$to - band$from + 1
ann <- list(x = 0.5, y = 1, yanchor = "bottom", xanchor = "middle", xref = "paper", yref = "paper", showarrow = FALSE)
plots <- lapply(seq_len(n), function(i) {
  z <- t(etm[[1]][,,i])
  ann$text <- paste(names(etm), i)
  plot_ly(z = z[seq.int(nrow(z), 1), ]) %>%
    layout(annotations = ann) %>%
    colorbar(title = i)
})
subplot(plots, nrows = 2, titleX = TRUE, shareX = TRUE, shareY = TRUE)

cc @jackparmer

まず、参照されている stars::plot.stars() の動作と、示されている plot_ly(etm) および subplot() の例を調べます。構造化データを扱う既存の plot_ly および subplot() のエントリポイントを特定し、stars オブジェクトをバンドパネルと対応するアノテーション、カラーバー、共有軸付きで自動的にプロットできることを完了条件とします。