Hi all, I am trying to convert modis level 1B data (type SWATH) to geotiff following these steps. 1) Examining the HDF gdalinfo d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdf gdalinfo HDF4_EOS:EOS_SWATH:"d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB 2) Extracting the 1. band of the subdataset MODIS_SWATH_Type_L1B:EV_1KM_RefSB gdal_translate -of GTiff -b 1 HDF4_EOS:EOS_SWATH:"d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB d:\modis\1b\band_1.tiff 4) Warping the image to WGS84 gdalwarp -t_srs EPSG:4326 d:\modis\1b\band_1.tiff d:\modis\1b\band_1_warp.tiff The example dataset, which i am using, can be downloaded here ftp://ladsftp.nascom.nasa.gov/allData/5/MOD021KM/2011/008//MOD021KM.A2011008.0850.005.2011008195744.hdf And now to my problem: I get always an certain shift (http://www.flickr.com/photos/54033867@N00/5363471425/in/photostream/), which is avoided if i am using the "Georeference MODIS" option in ENVI (http://www.flickr.com/photos/54033867@N00/5364083966/in/photostream/). I also tried a few workarounds, f.g. http://thread.gmane.org/gmane.comp.gis.gdal.devel/10482, but nothing worked properly. As far as I understand gdal supports Modis L1B SWATH, so what I am doing wrong? Any help is appreciated. Regards Anna P.S. I am using FWTools 2.4.7 on Windows
|
_______________________________________________ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
