Sorry, this was the problem! readGDAL('HDF4_EOS:EOS_GRID:"output/h12v12/h12v12_A2000M01.hdf":MOD_Grid_MOD16A2:PET_1km') -- View this message in context: http://r-sig-geo.2731867.n2.nabble.com/readGDAL-no-bands-in-dataset-MOD16A2-hdf4-tp7587868p7587869.html Sent from the R-sig-geo mailing list archive at Nabble.com.