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

Van: postgis-users <[email protected]> Namens James Howe
Verzonden: dinsdag 31 januari 2023 12:16
Aan: [email protected]
Onderwerp: [postgis-users] ST_Transform Z heights between ODN and WSG

U ontvangt niet vaak e-mail van [email protected]<mailto:[email protected]>. Meer 
informatie over waarom dit belangrijk 
is<https://aka.ms/LearnAboutSenderIdentification>
I have elevation data in SRID 7405 (Ordinance Survey Great Britain 1936 with 
Ordinance Datum Newlyn) and I want to convert it to SRID 4979 (World Geodetic 
System 1984 ensemble).

Both SRs have three axes defined in their spatial_ref_sys entries.

However, ST_Transform doesn't change the Z coordinate.

SELECT ST_Transform('SRID=7405;POINT(545068 258591 8.51)'::geometry, 4979);
-> POINT(0.121585415593599 52.206415868084655 8.51)

Am I missing something?
Also asked on SE: https://gis.stackexchange.com/q/447648/183667

Tried with postgres12-postgis3.0 and postgres14-postgis3.2.

Thanks,
James

_______________________________________________
postgis-users mailing list
[email protected]
https://lists.osgeo.org/mailman/listinfo/postgis-users

Reply via email to