[gdal-dev] reprojection of global rasters to mollweide

Michael Schulz mandschulz at googlemail.com
Tue Jul 31 02:32:35 PDT 2012


Dear GDAL'ers,

I have currently some issues with reprojecting global datasets in
geographic coordinates to the Mollweide projection (GDAL 1.8.0). When using
gdalwarp like:

gdalwarp -s_srs 'EPSG:4326' -t_srs '+proj=moll +lon_0=0 +x_0=0 +y_0=0
+ellps=WGS84 +datum=WGS84 +units=m +no_defs' LatLon.2001.tree.nd.tif
Mollweide.2001.tree.nd.tif

unpredictably, it sometimes works but other times I get the "Too many
points (441 out of 441) failed to transform" error. I was wondering whether
this could be a memory size issues (although the machine I'm working has
24GB RAM).

To overcome this adding the config option CHECK_WITH_INVERT_PROJ FALSE
works, but now the projected raster shows wrapped values outside the actual
valid coordinate range for Mollweide, e.g. the Aleutian Islands appear
twice (once inside the outer ellipse. Is this due to
the CHECK_WITH_INVERT_PROJ? I haven't worked with the cutline option
before, but should/could this be used here to solve this?

Happy about any feedback on this issue. Thanks, Michael

-- 
-----------------------------------------------------------
Michael Schulz
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120731/9624dd59/attachment.html>


More information about the gdal-dev mailing list