Robert Hijmans

This is what I see with CRAN-stars: ``` library(stars) #Loading required package: abind #Loading required package: sf #Linking to GEOS 3.9.1, GDAL 3.3.2, PROJ 7.2.1; sf_use_s2() is TRUE (r =...

And with the development version as well. So perhaps this was fixed somewhere in GDAL > 3.3.2?

It works for me on Ubuntu with GDAL 3.4.0: ``` library(stars) #Loading required package: abind #Loading required package: sf #Linking to GEOS 3.10.1, GDAL 3.4.0, PROJ 8.2.0; sf_use_s2() is TRUE...

It could also possibly be related to how R GDAL HDF5 is built on windows: https://github.com/OSGeo/gdal/issues/1428

@edzer: terra assigns an extent of (0, ncol, 0, nrow) when creating a SpatRaster from a matrix (#1); but it is (0,1,0,1) when a file does not supply an extent....

After a bit more digging, I now see that: - This also works on OSX (in my case with GDAL 3.2.2 ) - RTools4.2 (which I was not using) has...

I also see this in R 4.1 and 4.3 --- both on windows ``` # R version 4.1.3 (2022-03-10) -- "One Push-Up" library(terra) #terra 1.5.48 gdal() #[1] "3.4.1" r =...

Perhaps @roualt can comment?

It appears that the cause of the problem is that the windows HDF5 library was built without the "szip" library. That library is needed to read this file (presumably because...

Thank you very much @kalibera, @jeroen and @rsbivand!