How to convert wgs84 from ellipsoidal height to mean sea level i.e. height above goid?

[sebhan2]

Hello,

I want to convert wgs84 coordinates from ellipsoidal height to wgs84 coordinates with orthometric height, i.e. height above goid. Therefore I set up a new conda environment.

This is the example code I have written:

import pyproj

# define coordinates
lat = 54.44
lon = 9.88
z = 100

# define CRSs
wgs84_ellipse = pyproj.crs.CRS.from_epsg(4979)
wgs84_msl = pyproj.crs.CRS.from_epsg(9705)

# create transformer
tform = pyproj.transformer.Transformer.from_crs(crs_from=wgs84_ellipse, crs_to=wgs84_msl, always_xy=True)

# transform coordinates
print(tform.transform(lon,lat,z))

This tells me that I am missing a file:

from pyproj.transformer import TransformerGroup
t = TransformerGroup(wgs84_ellipse, wgs84_msl, always_xy=True)

However I am not able to download that file anymore.

Is there a way I can use pyproj to get the correct height conversion?

I am open also for alternative approaches that work!

Thank you.

Here is a link to my original question: https://gis.stackexchange.com/questions/454317/how-to-convert-wgs84-from-ellipsoidal-height-to-mean-sea-level-i-e-height-above?noredirect=1#comment741734_454317

---

Installation method/steps

• Installation method: conda
• How did you install PROJ? I used conda. Where is it installed? In my virtual environment.
• Please provide all commands/steps you used to install pyproj and PROJ:

conda create -n "test_env" python=3.11.0
conda install -c conda-forge pyproj
conda install -c conda-forge proj-data

Environment Information

• pyproj version you are attempting to install:
• PROJ version: 9.1.1
• Python version 3.11.0
• Operation System Information: Windows-10-10.0.19044-SP0

conda list gives me:

# Name                    Version                   Build  Channel
bzip2                     1.0.8                h8ffe710_4    conda-forge
ca-certificates           2022.12.7            h5b45459_0    conda-forge
certifi                   2022.12.7          pyhd8ed1ab_0    conda-forge
jpeg                      9e                   hcfcfb64_3    conda-forge
krb5                      1.20.1               heb0366b_0    conda-forge
lerc                      4.0.0                h63175ca_0    conda-forge
libcurl                   7.88.1               h68f0423_0    conda-forge
libdeflate                1.17                 hcfcfb64_0    conda-forge
libffi                    3.4.2                h8ffe710_5    conda-forge
libsqlite                 3.40.0               hcfcfb64_0    conda-forge
libssh2                   1.10.0               h9a1e1f7_3    conda-forge
libtiff                   4.5.0                hf8721a0_2    conda-forge
libzlib                   1.2.13               hcfcfb64_4    conda-forge
openssl                   3.0.8                hcfcfb64_0    conda-forge
pip                       23.0.1             pyhd8ed1ab_0    conda-forge
proj                      9.1.1                heca977f_2    conda-forge
proj-data                 1.13                 hd8ed1ab_0    conda-forge
pyproj                    3.4.1           py311h36482e4_1    conda-forge
python                    3.11.0          hcf16a7b_0_cpython    conda-forge
python_abi                3.11                    3_cp311    conda-forge
setuptools                67.4.0             pyhd8ed1ab_0    conda-forge
sqlite                    3.40.0               hcfcfb64_0    conda-forge
tk                        8.6.12               h8ffe710_0    conda-forge
tzdata                    2022g                h191b570_0    conda-forge
ucrt                      10.0.22621.0         h57928b3_0    conda-forge
vc                        14.3                hb6edc58_10    conda-forge
vs2015_runtime            14.34.31931         h4c5c07a_10    conda-forge
wheel                     0.38.4             pyhd8ed1ab_0    conda-forge
xz                        5.2.6                h8d14728_0    conda-forge
zstd                      1.5.2                h12be248_6    conda-forge

[russkel]

I am also having the issue where this conversion yields no difference in altitude.

I do not even get an error.