Edzer Pebesma
Good point. Currently, `read_stars` can't, as GDAL doesn't give you dimension names, and also doesn't tell you the order in which they come in the source file.
OK, I'm in - use a google hangout? together with @adrfantini @dblodgett-usgs you and me?
sounds good; what time in CEST?
OK, @mdsumner please suggest a time and a medium.
Here's a suggestion: pkg `stars` will * provide three functions, `read_ncdf`, `read_gdal` and `read_stars` * the first two do what they say, * `read_stars` when given a NetCDF file will...
For large number of GTiff layers you may try to set e.g., before loading terra, ```r Sys.setenv(GDAL_MAX_BAND_COUNT = 1000000) ``` See https://github.com/r-spatial/stars/issues/467#issuecomment-961876870
It's indeed a hard limit in the geotiff code: https://github.com/OSGeo/gdal/blob/master/frmts/gtiff/geotiff.cpp#L16791
> Probably some issue in stars, but good to keep in mind. Not sure ;-) It looks like mapview tries but fails to create a title from everything between `mapview(`...
> But perhaps it's something with my software stack... maybe switch to using `stars` from github? I'm pretty close to a new CRAN release.
Could you please clean up the formatting of this issue and reduce it to a minimal reproducible example?