<div dir="ltr"><div>Hi devs,</div><div><br></div><div>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? </div><div><br></div><div>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</div><div><br></div><div>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:  </div><div><br></div><div>{{{</div><div>gdalinfo /home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf</div><div>Driver: HDF4/Hierarchical Data Format Release 4</div><div>Files: /home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf</div><div>Size is 512, 512</div><div>Coordinate System is `'</div><div>Metadata:</div><div>  ASSOCIATEDINSTRUMENTSHORTNAME=MODIS</div><div>  ASSOCIATEDPLATFORMSHORTNAME=Terra</div><div>  ASSOCIATEDSENSORSHORTNAME=MODIS</div><div> ...</div><div>  GRINGPOINTLATITUDE=-40.000000, -29.855157, -29.884849, -40.035199</div><div>  GRINGPOINTLONGITUDE=-65.270364, -57.692450, -45.983544, -52.017065</div><div> ...</div><div>Subdatasets:</div><div>  SUBDATASET_1_NAME=HDF4_EOS:EOS_GRID:"/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf":MOD44W_250m_GRID:water_mask</div><div>  SUBDATASET_1_DESC=[4800x4800] water_mask MOD44W_250m_GRID (8-bit unsigned integer)</div><div>  SUBDATASET_2_NAME=HDF4_EOS:EOS_GRID:"/home/veroandreo/gisdata/MOD44W/MOD44W.A2000055.h13v12.005.2009212173504.hdf":MOD44W_250m_GRID:water_mask_QA</div><div>  SUBDATASET_2_DESC=[4800x4800] water_mask_QA MOD44W_250m_GRID (8-bit unsigned integer)</div><div>Corner Coordinates:</div><div>Upper Left  (    0.0,    0.0)</div><div>Lower Left  (    0.0,  512.0)</div><div>Upper Right (  512.0,    0.0)</div><div>Lower Right (  512.0,  512.0)</div><div>Center      (  256.0,  256.0)</div><div>}}}</div><div><br></div><div>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)</div><div><br></div><div>Sorry for bothersome! :P</div><div><br></div><div>Cheers, </div><div>Vero</div></div>