Dear Even, <br><br>thanks for your response, finally I got it working. The main issue was that the global dataset had an extent e.g. on the most western border of -180.0044643 which caused the error (although strangely not always, because I tried the same command with the same dataset and sometime it worked...)<div>
<br></div><div>Anyhow, clipping the data before and then specifying as you suggested the target extent (-te) works now.</div><div><br></div><div>Cheers, Michael</div><div><br><br><div class="gmail_quote">2012/8/1 Even Rouault <span dir="ltr"><<a href="mailto:even.rouault@mines-paris.org" target="_blank">even.rouault@mines-paris.org</a>></span><br>
<blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">Le mardi 31 juillet 2012 11:32:35, Michael Schulz a écrit :<br>
<div class="im">> Dear GDAL'ers,<br>
><br>
> I have currently some issues with reprojecting global datasets in<br>
> geographic coordinates to the Mollweide projection (GDAL 1.8.0). When using<br>
> gdalwarp like:<br>
><br>
> gdalwarp -s_srs 'EPSG:4326' -t_srs '+proj=moll +lon_0=0 +x_0=0 +y_0=0<br>
> +ellps=WGS84 +datum=WGS84 +units=m +no_defs' LatLon.2001.tree.nd.tif<br>
> Mollweide.2001.tree.nd.tif<br>
><br>
> unpredictably, it sometimes works but other times I get the "Too many<br>
> points (441 out of 441) failed to transform" error. I was wondering whether<br>
> this could be a memory size issues (although the machine I'm working has<br>
> 24GB RAM).<br>
<br>
</div>No, memory isn't the reason. The reason is that the algorithm used by gdalwarp<br>
doesn't manage to guess the extent of the raster once reprojected. It may or<br>
may not succeed depending on the extent of the source image and the<br>
reprojection involved. If you know the extent in the target coordinate system,<br>
specifying it with -te will solve that issue.<br></blockquote><div><br></div><div> </div><blockquote class="gmail_quote" style="margin:0 0 0 .8ex;border-left:1px #ccc solid;padding-left:1ex">
<div class="im"><br>
><br>
> To overcome this adding the config option CHECK_WITH_INVERT_PROJ FALSE<br>
> works, but now the projected raster shows wrapped values outside the actual<br>
> valid coordinate range for Mollweide, e.g. the Aleutian Islands appear<br>
> twice (once inside the outer ellipse. Is this due to<br>
> the CHECK_WITH_INVERT_PROJ ?<br>
<br>
</div>The default value of CHECK_WITH_INVERT_PROJ (YES) checks that when a<br>
reprojection of a coordinate is done, the invert reprojection gives a result<br>
close to the original coordinate, which validates the fact that we are in the<br>
validity area of the reprojection. If you override with<br>
CHECK_WITH_INVERT_PROJ=FALSE, you can have indeed strange artifacts due to non<br>
invertible transformations.<br>
<div class="im"><br>
> I haven't worked with the cutline option<br>
> before, but should/could this be used here to solve this?<br>
<br>
</div>I don't think this will help.<br>
<div class="HOEnZb"><div class="h5"><br>
><br>
> Happy about any feedback on this issue. Thanks, Michael<br>
</div></div></blockquote></div><br><br clear="all"><div><br></div>-- <br>-----------------------------------------------------------<br>Michael Schulz<br><br>Via Cime Bianche 10<br>IT-21023 Besozzo (VA)<br><br><br>
</div>