[gdal-dev] gdalwarp incomplete output with geostationary satellite projection

[Marcel Derrien]

Hi all,

I am using gdalwarp (1.9.0, Linux x86_64) to reproject data from Plate 
carre projection to the target Geostationary Satellite projection 
(proj.4 definition) in an image defined by its georeferenced extent and 
resolution (target georeferenced units, metres).

The output image is as expected  but the earth disk is not completely 
filled when approaching the poles.
 
During execution the following error messages are output:
...
ERROR 1: tolerance condition error
...
ERROR 1: Reprojection failed, err = -20, further errors will be 
supressed on the transform object.

Here is the command I am using:

gdalwarp -t_srs "+proj=geos +a=6378169 +b=6356583.8 +lon_0=0 
+h=35785831" -tr 3000.403357 3000.403357 -te -5570248.832537 
-5567248.428647 5567248.428647 5570248.832537  source.tif dest.tif
 
Has anyone an idea to obtain a full coverage of the earth disk ?

Thanks in advance,
---
Marcel Derrien

-------------- next part --------------
A non-text attachment was scrubbed...
Name: marcel_derrien.vcf
Type: text/x-vcard
Size: 375 bytes
Desc: not available
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20120611/2bd5b924/attachment.vcf>