This function transforms X/Y/Z/t, lon/lat/h/t points between two coordinate systems 1 and 2 using the PROJ function proj_trans_generic().
X1 is a column vector containing the first coordinates in the source coordinate system, geodetic longitude or coordinate X. Y1 is a column vector containing the second coordinates in the source coordinate system, geodetic latitude or coordinate Y. Z1 is a column vector containing the third first coordinates in the source coordinate system, ellipsoidal height or coordinate Z. t1 is a column vector containing the time coordinates in the source coordinate system. par1 is a text string containing the projection parameters for the source system, in PROJ ’+’ format, as EPSG code or as WKT2 code. par2 is a text string containing the projection parameters for the destination system, in PROJ ’+’ format, as EPSG code or as WKT2 code.
The coordinate vectors X1, Y1, Z1, and t1 must be all scalars or all column vectors (of the same size).
X2 is a column vector containing the first coordinates in the destination coordinate system, geodetic longitude or coordinate X. Y2 is a column vector containing the second coordinates in the destination coordinate system, geodetic latitude or coordinate Y. Z2 is a column vector containing the third coordinates in the destination coordinate system, ellipsoidal height or coordinate Z. t2 is a column vector containing the time coordinates in the destination coordinate system.
Angular units are by default radians, and linear meters, although other can be specified in par1, and par2, so the input data must be congruent, and output data will be congruent with the definitions.
See also: _op_fwd, _op_inv.
Package: octproj