This function unprojects cartesian projected coordinates (in a defined cartographic projection) into geodetic coordinates using the PROJ function proj_trans_generic().
X is a column vector containing the X projected coordinates. Y is a column vector containing the Y projected coordinates. params is a text string containing the projection parameters in PROJ format (ONLY format ’+’ style).
The coordinate vectors X and Y must be both scalars or both column vectors (of the same size). Linear units are by default meters, although other can be specified in params, so X and Y must be congruent with params
lon is a column vector containing the geodetic longitude. lat is a column vector containing the geodetic latitude.
If a projection error occurs, the resultant coordinates for the affected points have both Inf value and a warning message is emitted (one for each erroneous point). Angular units are by default radians, although other can be specified in params, so lon and lat will be congruent with params
See also: _op_fwd, _op_transform.
Package: octproj