geost.base.CptCollection.change_vertical_reference

CptCollection.change_vertical_reference(to_epsg: str | int | CRS)

Change the vertical reference of the collection object’s surface levels

Parameters:

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.

Some often-used vertical datums are: NAP : 5709 MSL NL depth : 9288 LAT NL depth : 9287 Ostend height : 5710

See epsg.io for more.

Examples

To change the header’s current vertical reference to NAP:

>>> self.change_horizontal_reference(5709)

This would be the same as:

>>> self.change_horizontal_reference("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 repsectively, you can use:

>>> self.change_horizontal_reference("NAP")
>>> self.change_horizontal_reference("Ostend height")