Dear Michael, Following your advice I have tired :
library(terra) library(raster) library(dplyr) r = terra::rast(file.choose()) plot(r$LST) r_lst = r$LST r_lst %>% raster() plot(r_lst) Then I effectively get a raster layer format from the raster package as I wanted, but looking at the plot and the raster layer information--see bellow-- the scene appears to need both latitude and longitude for correct geolocation, how can I do this procedure? > r_lst %>% raster() class : RasterLayer dimensions : 5632, 5400, 30412800 (nrow, ncol, ncell) resolution : 1, 1 (x, y) extent : 0, 5400, 0, 5632 (xmin, xmax, ymin, ymax) crs : NA source : LST names : LST values : 0, 65535 (min, max) The extent and the resolution are not the correct ones. Maybe in the same .h5 is all the data for georeferencing it although I could not see how to get to it and apply it to the layer. There is also another file with complementary data related to geolocation but how can I query it, and get the data for geolocation and apply it to the raster layer?. The objective is to save it as a GeoTIFF file to be opened in QGIS with correct geolocation. What do you suggest for converting the .h5 file to a georeferenced Raster Layer in GeoTIFF format ? Thanks a lot. Kind regards, Gabriel On Sat, Nov 13, 2021 at 12:22 PM Gabriel Cotlier <gabikl...@gmail.com> wrote: > Dear Michael Summer, > Thanks a lot for your swift reply. > I will give a try to the procedure you mentioned. > Kind regards, > Gabriel > > > On Sat, Nov 13, 2021 at 12:14 PM Michael Sumner <mdsum...@gmail.com> > wrote: > >> try terra::rast() or stars::read_stars() on the file, both use gdal to >> interrogate and read. What you get depends on the structure of the files, >> variables interpreted as subdatasets. stars is more general but sees >> variables as 2D arrays with bands, terra calls these bands layers. >> >> Alternatively try raster::raster() or RNetCDF or tidync which might be >> more straightforward than rhdf5 but all are quite different to each other >> and to the gdal view. >> >> hth >> >> On Sat, 13 Nov 2021, 18:56 Gabriel Cotlier, <gabikl...@gmail.com> wrote: >> >>> Hello, >>> >>> I would like to be able to visualize / plot a specific data layer within >>> the structure a hierarchical data file format .h5 in R. I could observe >>> the structure as follows : >>> >>> library(rhdf5) >>> r= h5ls(file.choose()) >>> r >>> >>> and also in this way identify the Group and Name of the data to be >>> plotted, >>> however I've tried I could not access the data itself, nor plot it as a >>> raster layer with its corresponding geographic coordinates. >>> I would like to be able to select such data, plot it as a raster layer >>> and >>> if possible save it as GeoTIFF file format to be further processed with >>> the >>> raster package as a common raster layer. >>> >>> If any small guidance example, reprex, of this procedure is possible >>> it would be appreciated. >>> >>> Thanks a lot in advance >>> Kind regards, >>> Gabriel Cotlier >>> -- >>> Gabriel Cotlier, PhD >>> Haifa Research Center for Theoretical Physics and Astrophysics (HCTPA) >>> University of Haifa >>> Israel >>> >>> [[alternative HTML version deleted]] >>> >>> _______________________________________________ >>> R-sig-Geo mailing list >>> R-sig-Geo@r-project.org >>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo >>> >> [[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@r-project.org https://stat.ethz.ch/mailman/listinfo/r-sig-geo
