I have been using Pyproj for some time now doing only simple projection and transformation. I am totally new to vertical datum and differences in GEOID among Coordinate systems. I have been trying to do vertical transformation using pyproj but I feel like I am missing something. For example I tried this :
import pyproj lat = 43.70012234 lng = -79.41629234 z = 100 wgs84 = pyproj.Proj("+init=EPSG:4326") NAD83NDV88 = pyproj.Proj("+init=EPSG:5498") results = pyproj.transform(wgs84, NAD83NDV88, lng, lat, z)
which returns :
(43.70012234, -79.41629234, 100.0).
My logic was to take a WGS84 coordinate with an ellipsoid elevation and try to convert that into NAD83 NADV88 (epsg:5498). But it seems like it is not that simple.
I found some EPSG codes for transformation related to this question like : EPSG:3858 +> WGS 84 to EGM2008 height (1) but I haven't been able to figure out how to use those codes in Pyproj. I am missing some pieces of the puzzle!
So to summarize can this be done?