geost.accessors.header.PointHeader.change_vertical_reference#

PointHeader.change_vertical_reference(from_epsg: str | int | CRS, to_epsg: str | int | CRS) → None[source]#

Change the vertical datum of the object’s surface levels.

Some often-used vertical datums are:

  • NAP : 5709

  • MSL NL depth : 9288

  • LAT NL depth : 9287

  • Ostend height : 5710

See epsg.io for more.

Parameters:

  • from_epsg (str | int | CRS) – EPSG of the source vertical datum. Takes anything that can be interpreted by pyproj.crs.CRS.from_user_input(). However, it must be a vertical datum.

  • to_epsg (str | int | CRS) – EPSG of the target vertical datum. Takes anything that can be interpreted by pyproj.crs.CRS.from_user_input(). However, it must be a vertical datum.

Examples

To change the current vertical reference from Ostend height to NAP:

>>> self.change_vertical_reference(5710, 5709)

This would be the same as:

>>> self.change_vertical_reference("epsg:5710", "epsg:5709")

As the Pyproj constructors are very flexible, you can even use the CRS’s full official name instead of an EPSG number. E.g. for changing to NAP and the Belgian Ostend height vertical datums respectively, you can use:

>>> self.change_vertical_reference("Ostend height", "NAP")