[R-sig-Geo] problem reading hdf image

[Agustin Lobo]

Yes, no problem from that side. I've used
gdal_translate -sds

and get all subdatasets saved as individual tif layers.
But we should know what's wrong with the hdf flle as there are many
thousands of files like this on the
modis distribution data servers.

Agus

On Tue, Oct 8, 2013 at 1:58 PM, José Miguel Barrios
<jmbarriosg at gmail.com> wrote:
> Have you tried converting the layer(s) of interest to something like ".tif"
> by using gdal_translate?
>
>

>>
>> On Tue, Oct 8, 2013 at 1:34 PM, Matteo Mattiuzzi
>> <matteo.mattiuzzi at boku.ac.at> wrote:
>> > Dear Agus,
>> >
>> > yes please send the output of gdalinfo (close to the end is the string
>> you have to use).
>> >
>> >
>> > Matteo
>> >
>> >
>> >
>> >
>> >
>> >>>> Agustin Lobo <alobolistas at gmail.com> 10/07/13 5:42 PM >>>
>> > Hi!
>> >
>> > I get the following error when try to read an hdf image (release 4):
>> >
>> >> modisima <-
>> stack(file.path(modisdir,"MOD09A1.A2011289.h11v11.005.2011298101214.hdf"))
>> > rgdal: version: 0.8-11, (SVN revision 479M)
>> > Geospatial Data Abstraction Library extensions to R successfully loaded
>> > Loaded GDAL runtime: GDAL 1.10.0, released 2013/04/24
>> > Path to GDAL shared files: /usr/share/gdal/1.10
>> > Loaded PROJ.4 runtime: Rel. 4.8.0, 6 March 2012, [PJ_VERSION: 480]
>> > Path to PROJ.4 shared files: (autodetected)
>> > Error in GDALinfo(filename, silent = silent, returnRAT = RAT,
>> > returnCategoryNames = RAT) :
>> >   object 'RATlist' not found
>> > Error in .rasterObjectFromFile(x, objecttype = "RasterBrick", ...) :
>> >   Cannot create a RasterLayer object from this file.
>> > In addition: Warning message:
>> > In dim(x) : no bands in dataset
>> >
>> > Also:
>> >
>> >>
>> GDALinfo(file.path(modisdir,"MOD09A1.A2011289.h11v11.005.2011298101214.hdf"))
>> > Error in GDALinfo(file.path(modisdir,
>> > "MOD09A1.A2011289.h11v11.005.2011298101214.hdf")) :
>> >   object 'scaleOffset' not found
>> > In addition: Warning messages:
>> > 1: In dim(x) : no bands in dataset
>> > 2: In GDALinfo(file.path(modisdir,
>> > "MOD09A1.A2011289.h11v11.005.2011298101214.hdf")) :
>> >   GeoTransform values not available
>> >
>> > Nevertheless, hdf4 and 5 are listed by gdalDrivers() and, on the linux
>> > terminal, I do get
>> > output from gdalinfo (quite long, let me know if I should send it).
>> > Also, I can display the
>> > image in Qgis with no problem.
>> >
>> > Any idea on where the problem could be?
>> > Thanks,
>> > Agus
>> >
>> >> sessionInfo()
>> > R version 3.0.2 (2013-09-25)
>> > Platform: x86_64-pc-linux-gnu (64-bit)
>> >
>> > locale:
>> >  [1] LC_CTYPE=en_US.UTF-8       LC_NUMERIC=C
>> > LC_TIME=es_ES.UTF-8        LC_COLLATE=en_US.UTF-8
>> >  [5] LC_MONETARY=es_ES.UTF-8    LC_MESSAGES=en_US.UTF-8
>> > LC_PAPER=es_ES.UTF-8       LC_NAME=C
>> >  [9] LC_ADDRESS=C               LC_TELEPHONE=C
>> > LC_MEASUREMENT=es_ES.UTF-8 LC_IDENTIFICATION=C
>> >
>> > attached base packages:
>> > [1] stats     graphics  grDevices utils     datasets  methods   base
>> >
>> > other attached packages:
>> > [1] rgdal_0.8-11  MODIS_0.10-8  raster_2.1-49 sp_1.0-13
>> >
>> > loaded via a namespace (and not attached):
>> > [1] grid_3.0.2      lattice_0.20-23 tools_3.0.2
>> >
>> > _______________________________________________
>> > R-sig-Geo mailing list
>> > R-sig-Geo at r-project.org
>> > https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>> >
>> >
>> >
>>
>> _______________________________________________
>> R-sig-Geo mailing list
>> R-sig-Geo at r-project.org
>> https://stat.ethz.ch/mailman/listinfo/r-sig-geo
>>
>
>         [[alternative HTML version deleted]]
>
> _______________________________________________
> R-sig-Geo mailing list
> R-sig-Geo at r-project.org
> https://stat.ethz.ch/mailman/listinfo/r-sig-geo