[gdal-dev] again .hdf to .tif

[Yann Chemin]

Hi there,

I am using this one right now on Level 2 and 3 MODIS products

-------------------------------------------
#!/bin/bash
mkdir processed
for file in *.hdf
do SDS_list=$(gdalinfo $file | grep "SUBDATASET_.*_NAME.*")
	for subdataset in $SDS_list
	do echo ${subdataset#*=} > tempfile.txt
	out=$(sed 's/\(.*\):\(.*\):"\(.*\).A\(.*\)\.\h\(.*\)\.\(.*\)\.\(.*\).hdf":\(.*\):\(.*\)/\3\_\4\_h\5\_\8\_\9\.tif/'
tempfile.txt)
	gdalwarp -of GTiff -t_srs EPSG:4326 ${subdataset#*=} processed/$out
	done
done
-------------------------------------

maybe of help,
Yann

2009/3/26 Anja Roesel <anja.roesel at zmaw.de>:
> sorry i forgot to mention the error message:
>>Segmentation fault
>
>
>
>
>
> Hi list,
>
> I am trying to convert MODIS .hdf files to .tif files (it worked with
> ASTER files earlier, thats why I am now really curious).
>
> i typed
>
>>gdalinfo MOD021KM.A2008134.1935.005.2009071201911.hdf
>
> result:
> ......
> SUBDATASET_1_NAME=HDF4_EOS:EOS_SWATH:"MOD021KM.A2008134.1935.005.2009071201911.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB
> SUBDATASET_1_DESC=[15x2030x1354] EV_1KM_RefSB (16-bit unsigned integer)
> SUBDATASET_2_NAME=HDF4_EOS:EOS_SWATH:"MOD021KM.A2008134.1935.005.2009071201911.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB_Uncert_Indexes
> SUBDATASET_2_DESC=
> ........
>
> then:
>
>>gdalwarp
> HDF4_EOS:EOS_SWATH:"MOD021KM.A2008134.1935.005.2009071201911.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB
>  modis.tif
>
> Have you any ideas, what is the problem?
> Thanks,
> Anja
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev
>



-- 
Yann Chemin
Mobile: +33 (06) 10 11 39 26
Home: +33 (02) 35 27 08 20,
Address: Gite de Mortagne,
16 rue de la chenaie,
76110 Bec de Mortagne,
France

Perso: http://www.freewebs.com/ychemin
YiKingDo: http://yikingdo.unblog.fr/