Clip latitude to [90;-90]

[JanJakes]

It seems that Imposm3 produces invalid coordinates from (0 -90) for South pole:

Query SELECT ST_AsText(geometry) FROM osm_destination_points WHERE osm_id = 436012592; results in POINT(-0.00460936092720608 nan).

[ImreSamu]

It seems that Imposm3 produces invalid coordinates from (0 -90) for South pole:

imho :

• it is normal with . -srid 3857
  • http://epsg.io/3857 : Area of use: World between 85.06°S and 85.06°N.
• use : -srid 4326 for South pole
  • http://epsg.io/4326 : Area of use: World.

https://en.wikipedia.org/wiki/Web_Mercator

"Because the Mercator projects the poles at infinity, Google Maps cannot show the poles. Instead it cuts off coverage at 85.051129° north and south. This is not considered a limitation, given the purpose of the service. The value 85.051129° is the latitude at which the full map becomes a square, and is computed as φ given y = 0"

imposm3 ... -srid 3857

+ psql -d imposm3dev -c 'SELECT id,  ST_AsText(geometry) , ST_SRID( geometry)   FROM osm_spoint ;'
    id     |            st_astext            | st_srid 
-----------+---------------------------------+---------
 436012592 | POINT(-0.00460936092720608 nan) |    3857
(1 row)

imposm3 ... -srid 4326

+ psql -d imposm3dev -c 'SELECT id,  ST_AsText(geometry) , ST_SRID( geometry)   FROM osm_spoint ;'
    id     |                   st_astext                    | st_srid 
-----------+------------------------------------------------+---------
 436012592 | POINT(-4.14065937093255e-08 -90.0000000626128) |    4326
(1 row)

[sgrossberndt]

Is this the same issue as https://github.com/omniscale/imposm3/issues/3 ?