Details
-
Improvement
-
Status: Open
-
Major
-
Resolution: Unresolved
-
None
-
None
-
None
Description
When a map projection has exact formulas for the forward projection but approximate formulas (e.g. using series expansion) for inverse projection, the accuracy of inverse projection with celestial bodies having high eccentricities may be improved as below:
Take a first estimation of (λ,φ), compute the forward projection (E′,N′) for that estimation, compute the errors compared to specified (E,N) values, convert that (ΔE,ΔN) error into a (Δλ,Δφ) error using the inverse of Jacobian matrix, correct (λ,φ) and continue iteratively until the error is small enough. This algorithm described in EPSG guidance note for the Orthographic projection could be applied to any map projection, not only Orthographic.
See org.apache.sis.referencing.operation.projection.Orthographic.inverseTransform for implementation.