I am using gdalwarp in python to extract and re-project reflectance bands from MODIS level 1B (MOD021km). An example of the command that I am using is the following:
gdalwarp -of GTIFF -tps -t_srs EPSG:4326 -tr 0.01 0.01 -r bilinear -te [extent] -overwrite HDF4_EOS:EOS_SWATH:F:\Data\MYD021KM.A2007094.0930.006.2012075054822.hdf:MODIS_SWATH_Type_L1B:EV_500_Aggr1km_RefSB F:\Data\Output.tifIt converts the hdf subdataset without any problem but the output file seems to have some sort of artifact, specially close to the edges. (See below)
I converted the same data using Modis Conversion Toolkit (MCTK) and it seems to be much better. (Below)
I tried different resampling method and different output projection but still the same results! Is there anyone having this issue with gdalwarp?
أكثر...
gdalwarp -of GTIFF -tps -t_srs EPSG:4326 -tr 0.01 0.01 -r bilinear -te [extent] -overwrite HDF4_EOS:EOS_SWATH:F:\Data\MYD021KM.A2007094.0930.006.2012075054822.hdf:MODIS_SWATH_Type_L1B:EV_500_Aggr1km_RefSB F:\Data\Output.tifIt converts the hdf subdataset without any problem but the output file seems to have some sort of artifact, specially close to the edges. (See below)
I converted the same data using Modis Conversion Toolkit (MCTK) and it seems to be much better. (Below)
I tried different resampling method and different output projection but still the same results! Is there anyone having this issue with gdalwarp?
أكثر...