[GRASSLIST:6350] Re: Goode in Grass

[Markus Neteler]

On Tue, Apr 05, 2005 at 10:13:58AM +0200, Maciek Sieczka wrote:
> From: "Markus Neteler" <neteler at itc.it>
> >On Sun, Apr 03, 2005 at 11:26:52AM +0200, Maciek Sieczka wrote:
> >>Hi
> >>
> >>Could somebody please share the PROJ_INFO for Goode's projection? Or  the
> >>proj4 definition?
> >>
> >>I'd like to work on MODIS imagery of Poland, Eastern Europe. I couldn't
> >>figure out the appropriate projection definition by myself.
> 
> >Most MODIS maps are in Sinusoidal (on MODIS sphere).
> >Where did you get the map?
> 
> say ftp://ftp.glcf.umiacs.umd.edu/modis/500m/Eurasia/Goodes.EUAS.2000321/
> via the GLCF UMIACS www

Looking at
 http://xserve.flids.com/pipermail/gdal-dev/2004-March/005095.html
 "... cutting the image into pieces which I seperately reprojected using either a 
  sinusoidal or a mollweide projection ..."

which contains the link to
 http://gis.esri.com/library/userconf/proc98/PROCEED/TO850/PAP844/P844.HTM

>From this I tried:

     /*
     /*   Northern Eurasia
     /*   REGION = "02   "
 4,   30.0000000,  65.0000000
     -40.0000000,  40.7366111
     -40.0000000,  60.0000000
     -10.0000000,  60.0000000
     -10.0000000,  90.0000000
     180.0000000,  90.0000000
     180.0000000,  40.7366111
end

#####################################
dd.r02:

input
projection GEOGRAPHIC
units DD
parameters
output
projection MOLLWEIDE
units METERS
spheroid SPHERE
xshift 3335846.22854
yshift -336410.83237
parameters
30 00 00
end
#####################################

#note that we cannot cut out on the fly:
gdalwarp -s_srs '+proj=moll +ellps=sphere +lon_0=30E \
 +y_0=-336410.83237 +x_0=3335846.22854' \
 Goodes.EUAS.2000321.band4.tif d02_europe_EUAS.2000321.band4.moll.tif

#Since it's only defined for a subregion, extract the valid part
#(not completely true for part of greenland, I think):
gdalwarp -s_srs '+proj=moll +ellps=sphere +lon_0=30E \
 +y_0=-336410.83237 +x_0=3335846.22854' \
 -t_srs '+init=epsg:4326' -te -40 40 180 90 \
 d02_europe_EUAS.2000321.band4.moll.tif d02_europe_EUAS.2000321.band4.LL.tif

... I didn't fully try, was running out of disk space.

gdalinfo d02_europe_EUAS.2000321.band4.LL.tif
qgis d02_europe_EUAS.2000321.band4.LL.tif
r.in.gdal ....

Maybe this gives you some ideas.
Please post a solution if you find anything.

Markus