[gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

[Even Rouault]

On mardi 24 avril 2018 13:40:24 CEST Menno van Scheers - HUSS wrote:
> Just verified this.
> 
> Warping the original dataset with the following commandline : 
> gdalwarp -t_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs " -te -65
> -65 65 65 -te_srs "+proj=longlat +ellps=WGS84 +datum=WGS84 +no_defs"
> "C:\Temp\nl\newtiff.tiff" "C:\temp\nl\oldgdalwarp.tiff"
 
> Succeeds in v2.1.0 with a lot of tolerance condition errors.
> Fails in v2.2.4 with  : ERROR 1: Too many points (5625 out of 5625) failed
> to transform, unable to compute output bounds.

Can you share this original image ? I get no error with master on the grib file with the above 
commandline.

 
> Best regards,
> 
> Menno van Scheers
> 
> 
> -----Original Message-----
> From: gdal-dev [mailto:gdal-dev-bounces at lists.osgeo.org] On Behalf Of Menno
> van Scheers - HUSS
 Sent: dinsdag 24 april 2018 15:11
> To: Andre Joost <andre+joost at nurfuerspam.de>; gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator
> projection
 
> Yeah, I already tried warping with smaller target extends.
> Im quiete certain earlier versions of gdal would throw tollerance condition
> errors, but would still warp.
 With the current version it just wont.
> 
> Menno van Scheers
> 
> -----Original Message-----
> From: gdal-dev [mailto:gdal-dev-bounces at lists.osgeo.org] On Behalf Of Andre
> Joost
 Sent: zaterdag 21 april 2018 19:28
> To: gdal-dev at lists.osgeo.org
> Subject: Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator
> projection
 
> Am 20.04.18 um 22:10 schrieb Even Rouault:
> 
> > On jeudi 19 avril 2018 09:17:35 CEST Menno van Scheers - HUSS wrote:
> > 
> >> Hi,
> >>
> >>
> >>
> >> Built version 2.2.4 of gdal with no additional libraries.
> >> Currently trying to warp a grib dataset to webmercator projection 
> >> with the following arguments : gdalwarp  -t_srs "+proj=geos +lon_0=0
> >> +h=35785831
> >> +x_0=0 +y_0=0 +a=6378140 +b=6356755 +units=m +no_defs"
> >> 201804181115_L-000-MSG4__-MPEF________-MPEG_____-000000___-201804181115-_
> >> _.
 grib webm.tiff
> >>
> >>
> >>
> >> This task fails with : ERROR 1: Too many points (5625 out of 5625) 
> >> failed to transform, unable to compute output bounds.
> >>
> >>
> >>
> >> Source dataset has a geostationary projection Basic info on the 
> >> dataset :
> >>
> >>
> >>
> >> PROJ.4 string is:
> >> '+proj=geos +lon_0=0 +h=35785831 +x_0=0 +y_0=0 +a=6378140 +b=6356755
> >> +units=m +no_defs ' Origin =
> >> (-5570248.477582973428071,5570248.477582973428071)
> >> Pixel Size = (3000.403165948275728,-3000.403165948275728)
> >> Corner Coordinates:
> >> Upper Left  (-5570248.478, 5570248.478) ERROR 1: tolerance condition 
> >> error
> >>
> >>
> >>
> >> Lower Left  (-5570248.478,-5567248.074) ERROR 1: tolerance condition 
> >> error
> >>
> >>
> >>
> >> Upper Right ( 5567248.074, 5570248.478) ERROR 1: tolerance condition 
> >> error
> >>
> >>
> >>
> >> Lower Right ( 5567248.074,-5567248.074) ERROR 1: tolerance condition 
> >> error
> >>
> >>
> >>
> >>
> >> Example of the dataset :
> >> https://hussbv-my.sharepoint.com/:u:/g/personal/menno_huss_nl/EQwTxBK
> >> wPQlOjP Z_DZcuZWcBgBs0y96_9Vev3tIqnne3KA?e=k6hmYU
> >
> >
> >
> > Menno,
> >
> >
> >
> > Yes gdalwarp has a hard-time with such images since a lot of points in 
> > them correspond to points not on the Earth surface...
> > The workaround is to explicitly specify the target extent (in the 
> > target SRS), so for full world
> > coverage:
> > -te -20037508.34 -20037508.34 20037508.34 20037508.34
> >
> >
-------------- next part --------------
An HTML attachment was scrubbed...
URL: <http://lists.osgeo.org/pipermail/gdal-dev/attachments/20180424/57314512/attachment-0001.html>