[GRASS-user] wanted: location parameters for modis import

Hello List!

I tried to import this image to GRASS GIS:
http://rapidfire.sci.gsfc.nasa.gov/gallery/?2003239-0827/Svalbard.A2003239.1300.250m.jpg

If you follow the "more info" link you'll see that the projection is soemthing like a Lambert Azimuthal (equal area?) aver a sphere, not an ellipsoid. Additionally, it's rotated.

I didn't manage to set up a location to import the file into with these parameters. Does anyone have an advice on this?

r.in.gdal can't import the file directly, as it's rotated.
I tried gdalwarp to rotate the file by the information in its .jgw file. r.in.gdal will then import the file but the problem of location steup stays.

I want the image to end up in an UTM33-location, so the import location has to be set up correctly for r.proj to perform nicely.

thanks in advance,
\flo.

Florian Kindl <florian.kindl@uibk.ac.at> writes:

> I tried to import this image to GRASS GIS:
> http://rapidfire.sci.gsfc.nasa.gov/gallery/?2003239-0827/
> Svalbard.A2003239.1300.250m.jpg

> If you follow the "more info" link you'll see that the projection is
> soemthing like a Lambert Azimuthal (equal area?) aver a sphere, not
> an ellipsoid.

        There's a hack mentioned at [1], which sometimes helps to get
        MODIS L3 data to work with GDAL.

> Additionally, it's rotated.

> I didn't manage to set up a location to import the file into with
> these parameters. Does anyone have an advice on this?

> r.in.gdal can't import the file directly, as it's rotated. I tried
> gdalwarp to rotate the file by the information in its .jgw
> file. r.in.gdal will then import the file but the problem of
> location steup stays.

> I want the image to end up in an UTM33-location, so the import
> location has to be set up correctly for r.proj to perform nicely.

        Currently, you'll probably need to use `gdalwarp' (as suggested
        in [1]) to reproject the image, since GRASS (or PROJ.4) seems to
        lack the support for the ``WGS84 datum + spherical projection''
        combination.

        I guess, PROJ 4.6.0 will support combining any datum with any
        projection. It may cause GRASS to support the exact projections
        used by MODIS L3 data.

[1] http://trac.osgeo.org/gdal/ticket/1239