Dear GDAL'ers, <div><br></div><div>I have currently some issues with reprojecting global datasets in geographic coordinates to the Mollweide projection (GDAL 1.8.0). When using gdalwarp like:</div><div><br></div><div>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</div>
<div><br></div><div>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). </div>
<div><br></div><div>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?</div>
<div><br></div><div>Happy about any feedback on this issue. Thanks, Michael</div><div><br></div><div><div>-- </div>-----------------------------------------------------------<br>Michael Schulz<br><br><br><br><br>
</div>