I use a similar command, without the -a_ullr argument e.g.
gdal_translate -a_srs '+proj=sinu +R=6371007.181 +nadgrids=@null +wktext' <input HDF subdataset> <output tiff file> On Thu, Apr 18, 2013 at 6:10 AM, Rutger <kass...@gmail.com> wrote: > Hey Laura, > > Since the MxD14.A2 products are already projected in the sinusoidal > projection there is no need to warp using GCP's. Fitting a polynomial > through some GCP's will always be less accurate then using the projection > the data is already in. > > Somewhere in the metadata in the HDF there should be a boundingbox (corner > coordinates, in sinusoidal projection) for the specific tile you are using, > you can then assign these with the projection using gdal_translate. Once > assigned, you can use gdalwarp in the same way as you already described. > > Youre gdal_translate command should look somewhat like this: > gdal_translate -of "VRT" -a_srs "+proj=sinu +lon_0=0 +x_0=0 +y_0=0 > +a=6371007.181 +b=6371007.181 +units=m +no_defs" -a_ullr 2223901.039333 > -3335851.559000 3335851.559000 -4447802.078667 <input HDF subdataset> > <output vrt file> > > You can get the correct subdataset reference by running gdalinfo on you > HDF. > The coordinates i used are for tile h20v12, so make sure you use the ones > from your tile. I used a vrt as output which you can use as an input for > gdal_translate, but its also fine to use Geotiff or anything else you > prefer. > > Hope this helps. > > Regards, > Rutger > > > > > > > > > > > -- > View this message in context: > http://osgeo-org.1560.x6.nabble.com/MODIS-reprojection-tp5047700p5047887.html > Sent from the GDAL - Dev mailing list archive at Nabble.com. > _______________________________________________ > 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