[GRASS-dev] import hdf files with r.import?

Veronica Andreo veroandreo at gmail.com
Mon Oct 19 08:06:36 PDT 2015


Hi devs,

Would that be possible to have a simpler way to import hdf (modis
especially) files into GRASS, without the need of reprojecting and
translating into .tif the desired subdatasets first?

I'm using pymodis to download, reproject and convert into tif, and just
then r.in.gdal... The thing is that for this particular case r.import is
just the same as r.in.gdal, and IIUC, r.modis still uses MRT

There's another problem, though... We all know MODIS hdf tiles are in SIN
proj but they do not have it defined, as you can see with:

{{{
gdalinfo
/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf
Driver: HDF4/Hierarchical Data Format Release 4
Files:
/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf
Size is 512, 512
Coordinate System is `'
Metadata:
  ASSOCIATEDINSTRUMENTSHORTNAME=MODIS
  ASSOCIATEDPLATFORMSHORTNAME=Terra
  ASSOCIATEDSENSORSHORTNAME=MODIS
 ...
  GRINGPOINTLATITUDE=-40.000000, -29.855157, -29.884849, -40.035199
  GRINGPOINTLONGITUDE=-65.270364, -57.692450, -45.983544, -52.017065
 ...
Subdatasets:

SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf":MOD44W_250m_GRID:water_mask
  SUBDATASET_1_DESC=[4800x4800] water_mask MOD44W_250m_GRID (8-bit unsigned
integer)

SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf":MOD44W_250m_GRID:water_mask_QA
  SUBDATASET_2_DESC=[4800x4800] water_mask_QA MOD44W_250m_GRID (8-bit
unsigned integer)
Corner Coordinates:
Upper Left  (    0.0,    0.0)
Lower Left  (    0.0,  512.0)
Upper Right (  512.0,    0.0)
Lower Right (  512.0,  512.0)
Center      (  256.0,  256.0)
}}}

It would be nice to have, in grass, the functionality to select the desired
layers of an hdf and then get them reprojected and imported in one step as
r.import does :) (this is something more general than the MODIS issue, I
know)

Sorry for bothersome! :P

Cheers,
Vero
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/grass-dev/attachments/20151019/d00cb50e/attachment.html>


More information about the grass-dev mailing list