[gdal-dev] Transformation of a vertical datum

Previous Topic Next Topic
 
classic Classic list List threaded Threaded
3 messages Options
Reply | Threaded
Open this post in threaded view
|  
Report Content as Inappropriate

[gdal-dev] Transformation of a vertical datum

openDEM

Hi,

 

I would like to transform the vertical datum EPSG:5730 to EPSG:5621

 

For more information about the European height reference systems have a look at: http://www.crs-geo.eu/

 

I used gdaltransform, but the z value is not changing:

 

The CRS is identical EPSG:3035: ETRS89 European Lambert Azimuthal Equal Area.

 

gdaltransform -s_srs EPSG:3035+5730 -t_srs EPSG:3035+5621 < xxx.xyz > xxx3.xyz

 

xxx.xyz:

4650012.500 1649987.500 332.763

4650037.500 1649987.500 337.635

4650062.500 1649987.500 342.153

 

testepsg "EPSG:4326+5773" à success in both cases
 
 
Best,
 
Martin
 
 

 

 

 

 

 


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

Re: Transformation of a vertical datum

Even Rouault-2

On jeudi 16 mars 2017 16:57:11 CET OpenDEM wrote:

> Hi,

>

>

>

> I would like to transform the vertical datum EPSG:5730 to EPSG:5621

>

>

>

> For more information about the European height reference systems have a look

> at: http://www.crs-geo.eu/

>

>

>

> I used gdaltransform, but the z value is not changing:

 

Yes, I you look at the output of gdalsrsinfo on both SRS, you'll see that there is no explicit vertical shift grid mentionned (+geoidgrids=something paramet in proj.4 string). Thus no vertical shift is done. So, you should look if there are vertical shift grids betwen EVRF2005 &/ EVRF2000 and WGS84 ellipsoidal height, so that they can be used with proj.4

 

$ gdalsrsinfo EPSG:3035+5621

 

PROJ.4 : '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +vunits=m +no_defs '

 

OGC WKT :

COMPD_CS["ETRS89 / LAEA Europe + EVRF2007 height",

PROJCS["ETRS89 / LAEA Europe",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Lambert_Azimuthal_Equal_Area"],

PARAMETER["latitude_of_center",52],

PARAMETER["longitude_of_center",10],

PARAMETER["false_easting",4321000],

PARAMETER["false_northing",3210000],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AUTHORITY["EPSG","3035"]],

VERT_CS["EVRF2007 height",

VERT_DATUM["European Vertical Reference Frame 2007",2005,

AUTHORITY["EPSG","5215"]],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Up",UP],

AUTHORITY["EPSG","5621"]]]

 

On the contrary :

 

$ gdalsrsinfo EPSG:3035+5773

 

PROJ.4 : '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=egm96_15.gtx +vunits=m +no_defs '

 

OGC WKT :

COMPD_CS["ETRS89 / LAEA Europe + EGM96 geoid height",

PROJCS["ETRS89 / LAEA Europe",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Lambert_Azimuthal_Equal_Area"],

PARAMETER["latitude_of_center",52],

PARAMETER["longitude_of_center",10],

PARAMETER["false_easting",4321000],

PARAMETER["false_northing",3210000],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AUTHORITY["EPSG","3035"]],

VERT_CS["EGM96 geoid height",

VERT_DATUM["EGM96 geoid",2005,

EXTENSION["PROJ4_GRIDS","egm96_15.gtx"],

AUTHORITY["EPSG","5171"]],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Up",UP],

AUTHORITY["EPSG","5773"]]]

 

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
|  
Report Content as Inappropriate

Re: Transformation of a vertical datum

openDEM

Hi Even,

 

Thanks for your reply.

 

Unfortunately there are no shift grids available yet, but I contacted CRS-EU.

They want to discuss this topic at the next EUREF meeting.

 

Best,

 

Martin

 

From: Even Rouault [mailto:[hidden email]]
Sent: Donnerstag, 16. März 2017 17:22
To: [hidden email]
Cc: OpenDEM
Subject: Re: [gdal-dev] Transformation of a vertical datum

 

On jeudi 16 mars 2017 16:57:11 CET OpenDEM wrote:

> Hi,

>

>

>

> I would like to transform the vertical datum EPSG:5730 to EPSG:5621

>

>

>

> For more information about the European height reference systems have a look

> at: http://www.crs-geo.eu/

>

>

>

> I used gdaltransform, but the z value is not changing:

 

Yes, I you look at the output of gdalsrsinfo on both SRS, you'll see that there is no explicit vertical shift grid mentionned (+geoidgrids=something paramet in proj.4 string). Thus no vertical shift is done. So, you should look if there are vertical shift grids betwen EVRF2005 &/ EVRF2000 and WGS84 ellipsoidal height, so that they can be used with proj.4

 

$ gdalsrsinfo EPSG:3035+5621

 

PROJ.4 : '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +vunits=m +no_defs '

 

OGC WKT :

COMPD_CS["ETRS89 / LAEA Europe + EVRF2007 height",

PROJCS["ETRS89 / LAEA Europe",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Lambert_Azimuthal_Equal_Area"],

PARAMETER["latitude_of_center",52],

PARAMETER["longitude_of_center",10],

PARAMETER["false_easting",4321000],

PARAMETER["false_northing",3210000],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AUTHORITY["EPSG","3035"]],

VERT_CS["EVRF2007 height",

VERT_DATUM["European Vertical Reference Frame 2007",2005,

AUTHORITY["EPSG","5215"]],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Up",UP],

AUTHORITY["EPSG","5621"]]]

 

On the contrary :

 

$ gdalsrsinfo EPSG:3035+5773

 

PROJ.4 : '+proj=laea +lat_0=52 +lon_0=10 +x_0=4321000 +y_0=3210000 +ellps=GRS80 +towgs84=0,0,0,0,0,0,0 +units=m +geoidgrids=egm96_15.gtx +vunits=m +no_defs '

 

OGC WKT :

COMPD_CS["ETRS89 / LAEA Europe + EGM96 geoid height",

PROJCS["ETRS89 / LAEA Europe",

GEOGCS["ETRS89",

DATUM["European_Terrestrial_Reference_System_1989",

SPHEROID["GRS 1980",6378137,298.257222101,

AUTHORITY["EPSG","7019"]],

TOWGS84[0,0,0,0,0,0,0],

AUTHORITY["EPSG","6258"]],

PRIMEM["Greenwich",0,

AUTHORITY["EPSG","8901"]],

UNIT["degree",0.0174532925199433,

AUTHORITY["EPSG","9122"]],

AUTHORITY["EPSG","4258"]],

PROJECTION["Lambert_Azimuthal_Equal_Area"],

PARAMETER["latitude_of_center",52],

PARAMETER["longitude_of_center",10],

PARAMETER["false_easting",4321000],

PARAMETER["false_northing",3210000],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AUTHORITY["EPSG","3035"]],

VERT_CS["EGM96 geoid height",

VERT_DATUM["EGM96 geoid",2005,

EXTENSION["PROJ4_GRIDS","egm96_15.gtx"],

AUTHORITY["EPSG","5171"]],

UNIT["metre",1,

AUTHORITY["EPSG","9001"]],

AXIS["Up",UP],

AUTHORITY["EPSG","5773"]]]

 

Even

 

 

--

Spatialys - Geospatial professional services

http://www.spatialys.com


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