Hi Jan,

Thanks for you answer.
I tried your script but unfortunatly it does not work for my case.
If someone want to make some test with the kind of image I want to use, just
try this script to download and uncompress it.

TestFolder <- choose.dir()

download.file("
http://oceandata.sci.gsfc.nasa.gov/cgi/getfile/A20072132007243.L3m_MO_CHL_chlor_a_4km.bz2",
f<-paste(TestlFolder, "/hdf4test.bz2", sep = "")
library(R.utils)
gunzip(f, f2<-paste(TestFolder, "/testhdf4.hdf", sep = "")
GDALinfo(f2) # in my case it results on this message
`C:/test/testhdf4.hdf' not recognised as a supported file format.



Here is my sessionInfo() result.

R version 2.12.0 (2010-10-15)
Platform: x86_64-pc-mingw32/x64 (64-bit)

locale:
[1] LC_COLLATE=French_Canada.1252  LC_CTYPE=French_Canada.1252
LC_MONETARY=French_Canada.1252 LC_NUMERIC=C
LC_TIME=French_Canada.1252

attached base packages:
[1] grDevices datasets  splines   graphics  stats     tcltk     utils
methods   base

other attached packages:
 [1] R.utils_1.7.8     R.oo_1.8.1        R.methodsS3_1.2.1 rgdal_0.7-1
raster_1.9-5      sp_0.9-84         svSocket_0.9-51   TinnR_1.0.3
R2HTML_2.2        Hmisc_3.8-3       survival_2.35-8

loaded via a namespace (and not attached):
[1] cluster_1.14.0  grid_2.12.0     lattice_0.19-33 svMisc_0.9-61
tools_2.12.0


Thanks for any help !

Arnaud

	[[alternative HTML version deleted]]

