[GRASSLIST:6355] Re: Goode in Grass

[Maciek Sieczka]

Thanks Markus

I'll look into this in some spare time soon. Will let know for sure.

Maciek

----- Original Message ----- 
From: "Markus Neteler" <neteler at itc.it>
To: "Maciek Sieczka" <werchowyna at epf.pl>
Cc: <grasslist at baylor.edu>
Sent: Wednesday, April 06, 2005 6:17 PM
Subject: Re: [GRASSLIST:6332] Re: Goode in Grass


> 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
>