[GRASS-user] reprojection from Projection Mollweide

Hy everyone.
I have a problem with a reprojection.
I have a global raster (a grid esri) that I imported (by choosing the file
w001001.adf) in one location. I created the location with a reference system
that has these parameters (in a file prj.adf):
Projection Mollweide
Units METERS
Zunits NO
Xshift 0.0
Yshift 0.0
Parameters 6370997.0 6370997.0
   0 0 0.0 / * Longitude of projection center

So far so ok.

Then I tried to reproject the raster in a location with ETRS_1989_LAEA but
I get (at least on Italy) a strong localization error.

I tried also with WGS84 (EPSG 4326) but nothing....and even from the
command line making various adjustments to the parameters .. for example:
gdalwarp-s_srs "proj = moll ellps sphere lon_0 = 0 = x_0 y_0 = 0 = 0
a=6370997 b=6370997 no_defs units=m "-t_srs "EPSG : 4326 "map.tif
map_wgs84.tif

Does anyone have experience on?

Thanks

Gabriele
--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/reprojection-from-Projection-Mollweide-tp5797013p5797013.html
Sent from the Grass - Users mailing list archive at Nabble.com.

I found some other posts .... I'm getting close to the solution (I get good
results):

gdalwarp -s_srs '+proj=goode +ellps=sphere +lon_0=30E +towgs84=0,0,0,0,0,0,0
+x_0=3335846.22854 +y_0=20000' -t_srs '+proj=latlong +ellps=wgs84
towgs84=0,0,0,0,0,0,0' -te 7 33 34 59 map.tif map_wgs84.tif

I do not know why, but that's okay with +proj=good" and not "+proj = moll"

In addition, by removing the parameter y_0 = 20000 I get a more accurate
overlay.
So by the time I reprojected raster in the WGS84 as follows:

gdalwarp -s_srs ' proj=goode ellps=sphere lon_0=30E towgs84=0,0,0,0,0,0,0
x_0=3335846.22854' -t_srs ' proj=latlong ellps=wgs84 towgs84=0,0,0,0,0,0,0'
-te 7 33 34 59 map.tif map_wgs84.tif

Gabriele
--
View this message in context: http://osgeo-org.1803224.n2.nabble.com/reprojection-from-Projection-Mollweide-tp5797013p5797451.html
Sent from the Grass - Users mailing list archive at Nabble.com.