[gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff

Joaquim Luis jluis at ualg.pt
Sat Jan 29 12:20:25 EST 2011


Anna,

I wonder how ENVI does the warping and if there isn't any hidden trick 
that must be applied. I did the warping in Mirone where I can choose to 
do a gdalwarp or warping using a totally independent code (from Matlab) 
and the results are equal and ... wrong.

Joaquim

>
> hi all,
> excuse me for not giving up. i undertook some additional tests, and i 
> got the same problems i described below. But i can't believe that 
> nobody is working with gdal and modis level 1b data. At least this 
> format is listet as supported at http://www.gdal.org/frmt_hdf4.html.
> So my question: does anybody in this list work with modis level 1b? If 
> yes, is the georeference of the resulting image ok or does somebody 
> knows an workaround?
> maybe the answer is, that it is not possible to warp modis level 1b 
> data with gdal because the modis geolocation (mod3xxx) file is missing?
> i appreciate any help.
> cheers
> anna
>
>     ------------------------------------------------------------------------
>     *Von:* "Chaitanya kumar CH" <chaitanya.ch at gmail.com>
>     *Gesendet:* 19.01.2011 10:13:49
>     *An:* "anna auge" <annaauge at web.de>
>     *Betreff:* Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to
>     geotiff
>
>     Anna,
>
>     Check a couple of GCPs to rule out bad data. If they are OK, file
>     a ticket at http://trac.osgeo.org/gdal/newticket
>
>     On Wed, Jan 19, 2011 at 1:23 PM, anna auge <annaauge at web.de
>     <mailto:annaauge at web.de>> wrote:
>
>         So it is a problem with gdalwarp and the GCPs ín the GeoTiff.
>         Do you have any hints how to go on with this problem?
>         Cheers, Anna
>
>             ------------------------------------------------------------------------
>             *Von:* "Nikolaos Hatzopoulos" <nhatzop at gmail.com
>             <mailto:nhatzop at gmail.com>>
>             *Gesendet:* 19.01.2011 01:33:43
>             *An:* "anna auge" <annaauge at web.de <mailto:annaauge at web.de>>
>             *Betreff:* Re: [gdal-dev] Modis L1B SWATH: georef problem
>             hdf to geotiff
>
>
>             You are right, I notice that there isn't any difference
>             from the
>             band_1.tiff and band_1_warp.tiff
>
>             --Nikos Hatzopoulos
>
>             On Mon, Jan 17, 2011 at 8:15 AM, anna auge
>             <annaauge at web.de <mailto:annaauge at web.de>> wrote:
>
>                 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
>
>
>                 WEB.DE <http://WEB.DE> DSL Doppel-Flat ab 19,99
>                 €/mtl.! Jetzt mit
>                 gratis Handy-Flat!
>                 *http://produkte.web.de/go/DSL_Doppel_Flatrate/2*
>
>
>                 _______________________________________________
>                 gdal-dev mailing list
>                 gdal-dev at lists.osgeo.org <mailto:gdal-dev at lists.osgeo.org>
>                 http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
>         NEU: FreePhone - kostenlos mobil telefonieren und surfen!
>         Jetzt informieren: *http://produkte.web.de/go/webdefreephone*
>
>
>         _______________________________________________
>         gdal-dev mailing list
>         gdal-dev at lists.osgeo.org <mailto:gdal-dev at lists.osgeo.org>
>         http://lists.osgeo.org/mailman/listinfo/gdal-dev
>
>
>
>
>     -- 
>     Best regards,
>     Chaitanya kumar CH.
>     /tʃaɪθənjə/ /kʊmɑr/
>     +91-9494447584
>     17.2416N 80.1426E
>
>
>
> Empfehlen Sie WEB.DE DSL Ihren Freunden und Bekannten und wir
> belohnen Sie mit bis zu 50,- Euro! *https://freundschaftswerbung.web.de*
>
>
> _______________________________________________
> gdal-dev mailing list
> gdal-dev at lists.osgeo.org
> http://lists.osgeo.org/mailman/listinfo/gdal-dev

-------------- next part --------------
An HTML attachment was scrubbed...
URL: http://lists.osgeo.org/pipermail/gdal-dev/attachments/20110129/3e017b2a/attachment.html


More information about the gdal-dev mailing list