Dear list, I am trying to convert a data matrix corresponding to a geostationary satellite data to raster format or a format supported by GDAL. I have tried using HDF5 library with the following code: n.col <- 1701 n.row <- 651 n.rec <- n.col*n.row zz <- bzfile("file_to_convert.bz2","rb") DSSF <- readBin(zz,numeric(),n=n.rec,size=4,endian="big")/86400 nDSSF <- readBin(zz,integer(),n=n.rec,size=4,endian="big") ET0 <- readBin(zz,numeric(),n=n.rec,size=4,endian="big") T2 <- readBin(zz,numeric(),n=n.rec,size=4,endian="big") Rs <- array(DSSF, dim=c(n.col,n.row)) library(hdf5) hdf5save("ET0","ET0") I am trying to convert it to HDF5 so that I can extract the raster info with gdalingo and gdaltranslate, but the output generated by hdf5save is not recognize by gdalinfo as a supported format. I don't know if the problem is the structure of the data (I have tried with all the fields available in the input binary) I use to generate it but it is just a matrix of data, so it should be able to store it in a HDF5 file. I have tried also with rgdal library and the following code: library(rgdal) writeGDAL(ET0,"file.tif",drivername = "GTiff",type="Float32") But I am newby with R and I am not familiar with rgdal library so I am sure I am missing a step before trying to write the matrix data to the file. I think it would be easier to pass the data to HDF5 due to the matrix not having spatial reference on it. If anyone could help me so that I could open HDF5 output files with gdalinfo it would be great. Thanks. Best regards, Juan Manuel Moreno Rivera Remote Sensing and G.I.S. Group Regional Development Institute University of Castilla-La Mancha 02071 Campus of Albacete, Spain Tlf. 967599200+2623 P Before printing this mail make sure it is completely necessary. The environment is everyone's business
[[alternative HTML version deleted]] _______________________________________________ R-sig-Geo mailing list R-sig-Geo@stat.math.ethz.ch https://stat.ethz.ch/mailman/listinfo/r-sig-geo