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

classic Classic list List threaded Threaded
13 messages Options
Reply | Threaded
Open this post in threaded view
|

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

mvscheers

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/EQwTxBKwPQlOjPZ_DZcuZWcBgBs0y96_9Vev3tIqnne3KA?e=k6hmYU

 


Any help would be appreciated.

 

Menno van Scheers

 


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Even Rouault-2

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/EQwTxBKwPQlOjP

> 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

 

(note that the source image you mention appears to be fully black)

 

Even

 

--

Spatialys - Geospatial professional services

http://www.spatialys.com


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Andre Joost
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/EQwTxBKwPQlOjP
>> 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
>


I tried:
gdal_translate -a_srs "+proj=geos +h=35785831 +a=6378169 +b=6356583.8
+no_defs" -b 1
201804181115_L-000-MSG4__-MPEF________-MPEG_____-000000___-201804181115-__.grib
temp.tif
gdalwarp -t_srs EPSG:3857 -te -20037508.34 -20037508.34 20037508.34
20037508.34 temp.tif output.tif

but still get the same too much points. Even with selecting a smaller
target extent that should be within the visible sector.

Greetings,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

mvscheers
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:[hidden email]] On Behalf Of Andre Joost
Sent: zaterdag 21 april 2018 19:28
To: [hidden email]
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
>


I tried:
gdal_translate -a_srs "+proj=geos +h=35785831 +a=6378169 +b=6356583.8
+no_defs" -b 1
201804181115_L-000-MSG4__-MPEF________-MPEG_____-000000___-201804181115-__.grib
temp.tif
gdalwarp -t_srs EPSG:3857 -te -20037508.34 -20037508.34 20037508.34
20037508.34 temp.tif output.tif

but still get the same too much points. Even with selecting a smaller target extent that should be within the visible sector.

Greetings,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

mvscheers
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.

Best regards,

Menno van Scheers


-----Original Message-----
From: gdal-dev [mailto:[hidden email]] On Behalf Of Menno van Scheers - HUSS
Sent: dinsdag 24 april 2018 15:11
To: Andre Joost <[hidden email]>; [hidden email]
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:[hidden email]] On Behalf Of Andre Joost
Sent: zaterdag 21 april 2018 19:28
To: [hidden email]
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
>


I tried:
gdal_translate -a_srs "+proj=geos +h=35785831 +a=6378169 +b=6356583.8
+no_defs" -b 1
201804181115_L-000-MSG4__-MPEF________-MPEG_____-000000___-201804181115-__.grib
temp.tif
gdalwarp -t_srs EPSG:3857 -te -20037508.34 -20037508.34 20037508.34
20037508.34 temp.tif output.tif

but still get the same too much points. Even with selecting a smaller target extent that should be within the visible sector.

Greetings,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Even Rouault-2

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:[hidden email]] On Behalf Of Menno

> van Scheers - HUSS

Sent: dinsdag 24 april 2018 15:11

> To: Andre Joost <[hidden email]>; [hidden email]

> 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:[hidden email]] On Behalf Of Andre

> Joost

Sent: zaterdag 21 april 2018 19:28

> To: [hidden email]

> 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

> >

> >

>

>

>

> I tried:

> gdal_translate -a_srs "+proj=geos +h=35785831 +a=6378169 +b=6356583.8

> +no_defs" -b 1

> 201804181115_L-000-MSG4__-MPEF________-MPEG_____-000000___-201804181115-__.g

> rib

temp.tif

> gdalwarp -t_srs EPSG:3857 -te -20037508.34 -20037508.34 20037508.34

> 20037508.34 temp.tif output.tif

>

> but still get the same too much points. Even with selecting a smaller target

> extent that should be within the visible sector.

> Greetings,

> Andre Joost

>

>

> _______________________________________________

> gdal-dev mailing list

> [hidden email]

> https://lists.osgeo.org/mailman/listinfo/gdal-dev

> _______________________________________________

> gdal-dev mailing list

> [hidden email]

> https://lists.osgeo.org/mailman/listinfo/gdal-dev

> _______________________________________________

> gdal-dev mailing list

> [hidden email]

> https://lists.osgeo.org/mailman/listinfo/gdal-dev

 

 

--

Spatialys - Geospatial professional services

http://www.spatialys.com


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Andre Joost
In reply to this post by mvscheers
Am 24.04.18 um 15:40 schrieb Menno van Scheers - HUSS:

> 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.
>


Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the
full target extent, and fails with the current 2.2.4.

And the result was not all black.

HTH,
Andre Joost

_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Andre Joost
Am 25.04.18 um 14:29 schrieb Andre Joost:

> Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the
> full target extent, and fails with the current 2.2.4.
>

... and gdal 2.4.0dev from Gisinternals works the same as 2.2.3.

So only the OSGEO4W version 2.2.4 fails. Maybe because of PROJ 5.0.1?

My Linux with Ubuntugis is at v 2.2.2 and works as expected too.

HTH,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Andre Joost
Am 25.04.18 um 14:45 schrieb Andre Joost:

> Am 25.04.18 um 14:29 schrieb Andre Joost:
>
>> Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the
>> full target extent, and fails with the current 2.2.4.
>>
>
> ... and gdal 2.4.0dev from Gisinternals works the same as 2.2.3.
>
> So only the OSGEO4W version 2.2.4 fails. Maybe because of PROJ 5.0.1?
>

Seems to be. After downgrading to PROJ 4.9.3 it works again, GDAL 2.2.4
and gdal-dev.

HTH,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Kristian Evers-2
Andre,

The PROJ 5.0.1 OSGeo4W installation is faulty so that may be the problem. Specifically the null grid is not being installed in the 5.0.1 package as reported on the mailing list [0]. This has not been resolved yet. Try copying in the null grid to C:\OSGeo4W64\share\proj and see if that changes anything. You can use the file from your current 4.9.3 installation.

If that doesn’t work try to reproduce the problem using cs2cs and I’ll take a look it. 

/Kristian


On 25 Apr 2018, at 15:09, Andre Joost <[hidden email]> wrote:

Am 25.04.18 um 14:45 schrieb Andre Joost:
Am 25.04.18 um 14:29 schrieb Andre Joost:

Just checked here: It worked up to GDAL 2.2.3 (from Gisinternals) on the
full target extent, and fails with the current 2.2.4.


... and gdal 2.4.0dev from Gisinternals works the same as 2.2.3.

So only the OSGEO4W version 2.2.4 fails. Maybe because of PROJ 5.0.1?


Seems to be. After downgrading to PROJ 4.9.3 it works again, GDAL 2.2.4 and gdal-dev.

HTH,
Andre Joost


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Andre Joost
Am 25.04.18 um 15:50 schrieb Kristian Evers:
> Andre,
>
> The PROJ 5.0.1 OSGeo4W installation is faulty so that may be the
> problem. Specifically the null grid is not being installed in the
> 5.0.1 package as reported on the mailing list [0]. This has not been
> resolved yet. Try copying in the null grid to C:\OSGeo4W64\share\proj
> and see if that changes anything. You can use the file from your
> current 4.9.3 installation.
>

That does not improve the situation with gdalwarp. I will do some tatest
with cs2cs later.


Greetings,
Andre Joost

_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

Even Rouault-2

On mercredi 25 avril 2018 19:20:17 CEST Andre Joost wrote:

> Am 25.04.18 um 15:50 schrieb Kristian Evers:

> > Andre,

> >

> > The PROJ 5.0.1 OSGeo4W installation is faulty so that may be the

> > problem. Specifically the null grid is not being installed in the

> > 5.0.1 package as reported on the mailing list [0]. This has not been

> > resolved yet. Try copying in the null grid to C:\OSGeo4W64\share\proj

> > and see if that changes anything. You can use the file from your

> > current 4.9.3 installation.

>

> That does not improve the situation with gdalwarp. I will do some tatest

> with cs2cs later.

>

 

I've tracked the error down to pj_transform() implementation in proj 5. A potential fix is in https://github.com/OSGeo/proj.4/pull/955

 

This is specific to using the proj.4 API of proj5.

 

Even

 

--

Spatialys - Geospatial professional services

http://www.spatialys.com


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev
Reply | Threaded
Open this post in threaded view
|

Re: Warping EUMetsat MSG-4 MPE dataset to webmercator projection

mvscheers

Potential fix seems to be merged with master of proj4.
Can now successfully warp the dataset.

 

Thanks !

 

Menno van Scheers

 

From: gdal-dev [mailto:[hidden email]] On Behalf Of Even Rouault
Sent: woensdag 25 april 2018 22:12
To: [hidden email]
Cc: Andre Joost <andre+[hidden email]>
Subject: Re: [gdal-dev] Warping EUMetsat MSG-4 MPE dataset to webmercator projection

 

On mercredi 25 avril 2018 19:20:17 CEST Andre Joost wrote:

> Am 25.04.18 um 15:50 schrieb Kristian Evers:

> > Andre,

> >

> > The PROJ 5.0.1 OSGeo4W installation is faulty so that may be the

> > problem. Specifically the null grid is not being installed in the

> > 5.0.1 package as reported on the mailing list [0]. This has not been

> > resolved yet. Try copying in the null grid to C:\OSGeo4W64\share\proj

> > and see if that changes anything. You can use the file from your

> > current 4.9.3 installation.

>

> That does not improve the situation with gdalwarp. I will do some tatest

> with cs2cs later.

>

 

I've tracked the error down to pj_transform() implementation in proj 5. A potential fix is in https://github.com/OSGeo/proj.4/pull/955

 

This is specific to using the proj.4 API of proj5.

 

Even

 

--

Spatialys - Geospatial professional services

http://www.spatialys.com


_______________________________________________
gdal-dev mailing list
[hidden email]
https://lists.osgeo.org/mailman/listinfo/gdal-dev