¿Cómo reproyectar sólo una banda de un hdf?

Question

Tengo este archivo que contiene:

gdalinfo product.hdf
Coordinate System is `'
Metadata:
calendar=standard
center=gsfc
Subdatasets:  
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)

como puedes ver, hay 24 bandas Lo que necesito es sólo extraer la primera banda del subconjunto de datos.

Answer 1

Puedes resolver tu problema combinando gdalwarp y gdal_translate:

gdalwarp -geoloc -t_srs epsg:3410 -tr 25000 25000 HDF4_SDS:UNKOWN:"product1.hdf":5 warped.tif
gdal_translate -b 1 warped.tif band1.tif

Si quieres las 24 bandas por separado, puedes escribir un archivo por lotes con este contenido (en Linux puede ser diferente):

FOR /L %%G IN (1,1,24) DO gdal_translate -b %%G warped.tif %%G.tif