I want to convert my coodrinates from lon, lat
to x, y
of orthographic projection. I tried to use PyProj
to do it, but it uses sphere, not ellipsoid. Here it is:
lonlat_proj = pyproj.Proj(proj="longlat", ellps="WGS84")
lat = 90
lon = 0
alt = 0
# ECEF IS OKAY
ecef_proj = pyproj.Proj(proj="geocent", ellps="WGS84")
ortho_proj = pyproj.Proj(proj="ortho", lat_0=0, lon_0=0, ellps="WGS84")
# ORTHO IS NOT!!!
x, y, z = pyproj.transform(lonlat_proj, ecef_proj, lon, lat, alt)
print(x, y, z)
# Gives: 3.918e-10 (okay), 0.0, 6356752.314 (okay, it's earth semiminor)
x_ortho, y_ortho = pyproj.transform(lonlat_proj, ortho_proj, lon, lat)
print(x_ortho, y_ortho)
# Gives: 0.0 (okay), 6378137,0 (earth SEMIMAJOR instead of semiminor)
# WHY?!!
How can I cope with it?