Function file: [lon,lat] = _op_inv(X,Y,params)

This function unprojects cartesian projected coordinates (in a defined cartographic projection) into geodetic coordinates using the PROJ function proj_trans_generic().

INPUT ARGUMENTS:

  • X is a column vector containing by default the X projected coordinates.
  • Y is a column vector containing by default the Y projected coordinates.
  • params is a text string containing the projection parameters in PROJ format (ONLY format ’+’ style is allowed, in any other case results are not guaranteed).

OUTPUT ARGUMENTS:

  • X is a column vector containing by default the X projected coordinates.
  • Y is a column vector containing by default the Y projected coordinates.

The coordinate vectors X and Y must be both scalars or both column vectors (of the same size). lon and lat will be according to the input dimensions. Angular units are by default radians and linear meters, although other can be specified in params, so X and Y must be congruent with params. The same applies to the coordinate order at input and output.

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).

See also: _op_fwd, _op_transform.

Package: octproj