[postgis-users] ST_Transform Z heights between ODN and WSG

[Paul Ramsey]

I'm seeing something similar in my environment, even going between two
systems with explicit vertical datums.

pramsey=# SELECT ST_AsTexT(ST_Transform('SRID=7405;POINT(545068 258591
8.51)'::geometry, 3901));
                      st_astext
-----------------------------------------------------
 POINT Z (1681647.7511311213 6134585.208145529 8.51)

Time to put it in the debugger.

On Wed, Feb 1, 2023 at 4:13 AM James Howe <jmh205 at cam.ac.uk> wrote:
>
> I can confirm it works in PostGIS 2.5, but not PostGIS 3.0+.
> I've been using the official Docker images to test.
>
> Time to file a bug then?
>
> James
>
> > I had postgresql10-postgis2.5 running.
> > It didn't have 4979, so I inserted it.
> > -- ERROR: GetProj4StringSPI: Cannot find SRID (4979) in spatial_ref_sys
> > -- https://epsg.io/4979
> > INSERT into spatial_ref_sys (srid, auth_name, auth_srid, proj4text, srtext) values ( 4979, 'EPSG', > 4979, '+proj=longlat +datum=WGS84 +no_defs +type=crs', 'None');
> > SELECT ST_Transform('SRID=7405;POINT(545068 258591 8.51)'::geometry, 4979);
> > That rendered me:
> > st_transform                                                   |
> > ---------------------------------------------------------------+
> > POINT Z(0.1215563454880247 52.20645574758174 55.02789907436818)|
> >
> >
> > Met vriendelijke groet,
> > Wouter Scherphof
> _______________________________________________
> postgis-users mailing list
> postgis-users at lists.osgeo.org
> https://lists.osgeo.org/mailman/listinfo/postgis-users