Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
Nikolaos I have warped modis with swath2grid and mosaicked it with mapserver, I have not seen that issue. also you don't need to use a vrt to mosaic, this can be done with one gdalwarp call. gdalwarp file1 file2 file3 outfile Brian On Mon, 2011-07-25 at 17:20 -0700, Nikolaos Hatzopoulos wrote: It seems interesting I'll keep a note on this :) The question is how you can build a VRT mosaic directly from MODIS hdf files, because the problem is that you have to mosaic first and after that to transform them to the EPSG:4326 projection. If you transform them first and mosaic them after you might have white lines. --Nikos 2011/7/25 Brian Case r...@winkey.org anna, Nikolaos, after much experimentation I found something that seems to work. gdalwarp --config GEOL_AS_GCPS FULL \ -tps -geoloc -t_srs EPSG:4326 \ -te -119.88004906118 13.5092548294379 \ -92.6036641645677 34.7070050511582 \ HDF4_EOS:EOS_SWATH:MOD02HKM.A2011169.1750.005.NRT.hdf:MODIS_SWATH_Type_L1B:EV_500_RefSB \ MOD02HKM.A2011169.1750.005.NRT.tif note i tried to escape the newlines to make it easy to read but the email may mess that up. there is about 2.7 million gcp's in a normal 5min modis swath file, by default gdal only reads every 10th row and column of gcp's this is not enough to handle the bowtie effect. Brian On Tue, 2011-01-18 at 16:33 -0800, Nikolaos Hatzopoulos wrote: 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 annaa...@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 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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
cool! does it works with 150 files? --Nikos On Tue, Jul 26, 2011 at 9:17 AM, Brian Case r...@winkey.org wrote: Nikolaos I have warped modis with swath2grid and mosaicked it with mapserver, I have not seen that issue. also you don't need to use a vrt to mosaic, this can be done with one gdalwarp call. gdalwarp file1 file2 file3 outfile Brian On Mon, 2011-07-25 at 17:20 -0700, Nikolaos Hatzopoulos wrote: It seems interesting I'll keep a note on this :) The question is how you can build a VRT mosaic directly from MODIS hdf files, because the problem is that you have to mosaic first and after that to transform them to the EPSG:4326 projection. If you transform them first and mosaic them after you might have white lines. --Nikos 2011/7/25 Brian Case r...@winkey.org anna, Nikolaos, after much experimentation I found something that seems to work. gdalwarp --config GEOL_AS_GCPS FULL \ -tps -geoloc -t_srs EPSG:4326 \ -te -119.88004906118 13.5092548294379 \ -92.6036641645677 34.7070050511582 \ HDF4_EOS:EOS_SWATH:MOD02HKM.A2011169.1750.005.NRT.hdf:MODIS_SWATH_Type_L1B:EV_500_RefSB \ MOD02HKM.A2011169.1750.005.NRT.tif note i tried to escape the newlines to make it easy to read but the email may mess that up. there is about 2.7 million gcp's in a normal 5min modis swath file, by default gdal only reads every 10th row and column of gcp's this is not enough to handle the bowtie effect. Brian On Tue, 2011-01-18 at 16:33 -0800, Nikolaos Hatzopoulos wrote: 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 annaa...@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 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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
It seems interesting I'll keep a note on this :) The question is how you can build a VRT mosaic directly from MODIS hdf files, because the problem is that you have to mosaic first and after that to transform them to the EPSG:4326 projection. If you transform them first and mosaic them after you might have white lines. --Nikos 2011/7/25 Brian Case r...@winkey.org anna, Nikolaos, after much experimentation I found something that seems to work. gdalwarp --config GEOL_AS_GCPS FULL \ -tps -geoloc -t_srs EPSG:4326 \ -te -119.88004906118 13.5092548294379 \ -92.6036641645677 34.7070050511582 \ HDF4_EOS:EOS_SWATH:MOD02HKM.A2011169.1750.005.NRT.hdf:MODIS_SWATH_Type_L1B:EV_500_RefSB \ MOD02HKM.A2011169.1750.005.NRT.tif note i tried to escape the newlines to make it easy to read but the email may mess that up. there is about 2.7 million gcp's in a normal 5min modis swath file, by default gdal only reads every 10th row and column of gcp's this is not enough to handle the bowtie effect. Brian On Tue, 2011-01-18 at 16:33 -0800, Nikolaos Hatzopoulos wrote: 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 annaa...@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 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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
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...@gmail.com *Gesendet:* 19.01.2011 10:13:49 *An:* anna auge annaa...@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 annaa...@web.de mailto:annaa...@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 nhat...@gmail.com mailto:nhat...@gmail.com *Gesendet:* 19.01.2011 01:33:43 *An:* anna auge annaa...@web.de mailto:annaa...@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 annaa...@web.de mailto:annaa...@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@lists.osgeo.org mailto:gdal-dev@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@lists.osgeo.org mailto:gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards, Chaitanya kumar CH. /tʃaɪθənjə/ /kʊmɑr
Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
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.cheersanna Von: "Chaitanya kumar CH" chaitanya...@gmail.comGesendet: 19.01.2011 10:13:49An: "anna auge" annaa...@web.deBetreff: Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiffAnna,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 annaa...@web.de wrote: So it is a problem with gdalwarp and the GCPsn the GeoTiff. Do you have any hints how to go on with this problem? Cheers,Anna Von: "Nikolaos Hatzopoulos" nhat...@gmail.comGesendet: 19.01.2011 01:33:43An: "anna auge" annaa...@web.deBetreff: Re: [gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff You are right, I notice that there isn't any difference from theband_1.tiff and band_1_warp.tiff--Nikos Hatzopoulos On Mon, Jan 17, 2011 at 8:15 AM, anna auge annaa...@web.de wrote: Hi all,I am trying to convert modis level 1B data (type SWATH) to geotiff following these steps.1) Examining the HDFgdalinfo d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdfgdalinfo 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_RefSBgdal_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 WGS84gdalwarp -t_srs EPSG:4326 d:\modis\1b\band_1.tiff d:\modis\1b\band_1_warp.tiffThe example dataset, which i am using, can be downloaded hereftp://ladsftp.nascom.nasa.gov/allData/5/MOD021KM/2011/008//MOD021KM.A2011008.0850.005.2011008195744.hdfAnd 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 AnnaP.S. I am using FWTools 2.4.7 on Windows 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 listgdal-dev@lists.osgeo.orghttp://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 listgdal-dev@lists.osgeo.orghttp://lists.osgeo.org/mailman/listinfo/gdal-dev -- Best regards,Chaitanya kumar CH. /tʃaɪənjə/ /kʊmɑr/ +91-949444758417.2416N 80.1426E Empfehlen Sie WEB.DE DSL Ihren Freunden und Bekannten und wirbelohnen Sie mit bis zu 50,- Euro! https://freundschaftswerbung.web.de ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
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 annaa...@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 nhat...@gmail.com *Gesendet:* 19.01.2011 01:33:43 *An:* anna auge annaa...@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 annaa...@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 DSL Doppel-Flat ab 19,99 €/mtl.! Jetzt mit gratis Handy-Flat! *http://produkte.web.de/go/DSL_Doppel_Flatrate/2*http://produkte.web.de/go/DSL_Doppel_Flatrate/2 ___ gdal-dev mailing list gdal-dev@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*http://produkte.web.de/go/webdefreephone ___ gdal-dev mailing list gdal-dev@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 ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
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 annaa...@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 DSL Doppel-Flat ab 19,99 €/mtl.! Jetzt mit gratis Handy-Flat! *http://produkte.web.de/go/DSL_Doppel_Flatrate/2*http://produkte.web.de/go/DSL_Doppel_Flatrate/2 ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev ___ gdal-dev mailing list gdal-dev@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev
[gdal-dev] Modis L1B SWATH: georef problem hdf to geotiff
Hi all,I am trying to convert modis level 1B data (type SWATH) to geotiff following these steps.1) Examining the HDFgdalinfo d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdfgdalinfo HDF4_EOS:EOS_SWATH:"d:\modis\1b\MOD021KM.A2011008.0850.005.2011008195744.hdf":MODIS_SWATH_Type_L1B:EV_1KM_RefSB2) Extracting the 1. band of the subdataset MODIS_SWATH_Type_L1B:EV_1KM_RefSBgdal_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.tiff4) Warping the image to WGS84gdalwarp -t_srs EPSG:4326 d:\modis\1b\band_1.tiff d:\modis\1b\band_1_warp.tiffThe example dataset, which i am using, can be downloaded hereftp://ladsftp.nascom.nasa.gov/allData/5/MOD021KM/2011/008//MOD021KM.A2011008.0850.005.2011008195744.hdfAnd 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 AnnaP.S. I am using FWTools 2.4.7 on WindowsWEB.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@lists.osgeo.org http://lists.osgeo.org/mailman/listinfo/gdal-dev