/src/PROJ/src/conversions/axisswap.cpp

Source
/***********************************************************************

        Axis order operation for use with transformation pipelines.

                Kristian Evers, kreve@sdfe.dk, 2017-10-31

************************************************************************

Change the order and sign of 2,3 or 4 axes. Each of the possible four
axes are numbered with 1-4, such that the first input axis is 1, the
second is 2 and so on. The output ordering is controlled by a list of the
input axes re-ordered to the new mapping. Examples:

Reversing the order of the axes:

    +proj=axisswap +order=4,3,2,1

Swapping the first two axes (x and y):

    +proj=axisswap +order=2,1,3,4

The direction, or sign, of an axis can be changed by adding a minus in
front of the axis-number:

    +proj=axisswap +order=1,-2,3,4

It is only necessary to specify the axes that are affected by the swap
operation:

    +proj=axisswap +order=2,1

************************************************************************
* Copyright (c) 2017, Kristian Evers / SDFE
*
* Permission is hereby granted, free of charge, to any person obtaining a
* copy of this software and associated documentation files (the "Software"),
* to deal in the Software without restriction, including without limitation
* the rights to use, copy, modify, merge, publish, distribute, sublicense,
* and/or sell copies of the Software, and to permit persons to whom the
* Software is furnished to do so, subject to the following conditions:
*
* The above copyright notice and this permission notice shall be included
* in all copies or substantial portions of the Software.
*
* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
* DEALINGS IN THE SOFTWARE.
*
***********************************************************************/

#include <errno.h>
#include <stdlib.h>
#include <string.h>

#include <algorithm>

#include "proj.h"
#include "proj_internal.h"

PROJ_HEAD(axisswap, "Axis ordering");

namespace { // anonymous namespace
struct pj_axisswap_data {
    unsigned int axis[4];
    int sign[4];
};
} // anonymous namespace

static int sign(int x) { return (x > 0) - (x < 0); }

static PJ_XY pj_axisswap_forward_2d(PJ_LP lp, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    PJ_XY xy;

    double in[2] = {lp.lam, lp.phi};
    xy.x = in[Q->axis[0]] * Q->sign[0];
    xy.y = in[Q->axis[1]] * Q->sign[1];
    return xy;
}

static PJ_LP pj_axisswap_reverse_2d(PJ_XY xy, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    unsigned int i;
    PJ_COORD out, in;

    in.v[0] = xy.x;
    in.v[1] = xy.y;
    out = proj_coord_error();

    for (i = 0; i < 2; i++)
        out.v[Q->axis[i]] = in.v[i] * Q->sign[i];

    return out.lp;
}

static PJ_XYZ pj_axisswap_forward_3d(PJ_LPZ lpz, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    unsigned int i;
    PJ_COORD out, in;

    in.v[0] = lpz.lam;
    in.v[1] = lpz.phi;
    in.v[2] = lpz.z;
    out = proj_coord_error();

    for (i = 0; i < 3; i++)
        out.v[i] = in.v[Q->axis[i]] * Q->sign[i];

    return out.xyz;
}

static PJ_LPZ pj_axisswap_reverse_3d(PJ_XYZ xyz, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    unsigned int i;
    PJ_COORD in, out;

    out = proj_coord_error();
    in.v[0] = xyz.x;
    in.v[1] = xyz.y;
    in.v[2] = xyz.z;

    for (i = 0; i < 3; i++)
        out.v[Q->axis[i]] = in.v[i] * Q->sign[i];

    return out.lpz;
}

static void swap_xy_4d(PJ_COORD &coo, PJ *) {
    std::swap(coo.xyzt.x, coo.xyzt.y);
}

static void pj_axisswap_forward_4d(PJ_COORD &coo, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    unsigned int i;
    PJ_COORD out;

    for (i = 0; i < 4; i++)
        out.v[i] = coo.v[Q->axis[i]] * Q->sign[i];
    coo = out;
}

static void pj_axisswap_reverse_4d(PJ_COORD &coo, PJ *P) {
    struct pj_axisswap_data *Q = (struct pj_axisswap_data *)P->opaque;
    unsigned int i;
    PJ_COORD out;

    for (i = 0; i < 4; i++)
        out.v[Q->axis[i]] = coo.v[i] * Q->sign[i];

    coo = out;
}

/***********************************************************************/
PJ *PJ_CONVERSION(axisswap, 0) {
    /***********************************************************************/
    struct pj_axisswap_data *Q = static_cast<struct pj_axisswap_data *>(
        calloc(1, sizeof(struct pj_axisswap_data)));
    char *s;
    unsigned int i, j, n = 0;

    if (nullptr == Q)
        return pj_default_destructor(P, PROJ_ERR_OTHER /*ENOMEM*/);
    P->opaque = (void *)Q;

    /* +order and +axis are mutually exclusive */
    if (!pj_param_exists(P->params, "order") ==
        !pj_param_exists(P->params, "axis")) {
        proj_log_error(P,
                       _("must provide EITHER 'order' OR 'axis' parameter."));
        return pj_default_destructor(
            P, PROJ_ERR_INVALID_OP_MUTUALLY_EXCLUSIVE_ARGS);
    }

    /* fill axis list with indices from 4-7 to simplify duplicate search further
     * down */
    for (i = 0; i < 4; i++) {
        Q->axis[i] = i + 4;
        Q->sign[i] = 1;
    }

    /* if the "order" parameter is used */
    if (pj_param_exists(P->params, "order")) {
        /* read axis order */
        char *order = pj_param(P->ctx, P->params, "sorder").s;

        /* check that all characters are valid */
        for (i = 0; i < strlen(order); i++)
            if (strchr("1234-,", order[i]) == nullptr) {
                proj_log_error(P, _("unknown axis '%c'"), order[i]);
                return pj_default_destructor(
                    P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
            }

        /* read axes numbers and signs */
        s = order;
        n = 0;
        while (*s != '\0' && n < 4) {
            Q->axis[n] = abs(atoi(s)) - 1;
            if (Q->axis[n] > 3) {
                proj_log_error(P, _("invalid axis '%d'"), Q->axis[n]);
                return pj_default_destructor(
                    P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
            }
            Q->sign[n++] = sign(atoi(s));
            while (*s != '\0' && *s != ',')
                s++;
            if (*s == ',')
                s++;
        }
    }

    /* if the "axis" parameter is used */
    if (pj_param_exists(P->params, "axis")) {
        /* parse the classic PROJ.4 enu axis specification */
        for (i = 0; i < 3; i++) {
            switch (P->axis[i]) {
            case 'w':
                Q->sign[i] = -1;
                Q->axis[i] = 0;
                break;
            case 'e':
                Q->sign[i] = 1;
                Q->axis[i] = 0;
                break;
            case 's':
                Q->sign[i] = -1;
                Q->axis[i] = 1;
                break;
            case 'n':
                Q->sign[i] = 1;
                Q->axis[i] = 1;
                break;
            case 'd':
                Q->sign[i] = -1;
                Q->axis[i] = 2;
                break;
            case 'u':
                Q->sign[i] = 1;
                Q->axis[i] = 2;
                break;
            default:
                proj_log_error(P, _("unknown axis '%c'"), P->axis[i]);
                return pj_default_destructor(
                    P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
            }
        }
        n = 3;
    }

    /* check for duplicate axes */
    for (i = 0; i < 4; i++)
        for (j = 0; j < 4; j++) {
            if (i == j)
                continue;
            if (Q->axis[i] == Q->axis[j]) {
                proj_log_error(P, _("axisswap: duplicate axes specified"));
                return pj_default_destructor(
                    P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
            }
        }

    /* only map fwd/inv functions that are possible with the given axis setup */
    if (n == 4) {
        P->fwd4d = pj_axisswap_forward_4d;
        P->inv4d = pj_axisswap_reverse_4d;
    }
    if (n == 3 && Q->axis[0] < 3 && Q->axis[1] < 3 && Q->axis[2] < 3) {
        P->fwd3d = pj_axisswap_forward_3d;
        P->inv3d = pj_axisswap_reverse_3d;
    }
    if (n == 2) {
        if (Q->axis[0] == 1 && Q->sign[0] == 1 && Q->axis[1] == 0 &&
            Q->sign[1] == 1) {
            P->fwd4d = swap_xy_4d;
            P->inv4d = swap_xy_4d;
        } else if (Q->axis[0] < 2 && Q->axis[1] < 2) {
            P->fwd = pj_axisswap_forward_2d;
            P->inv = pj_axisswap_reverse_2d;
        }
    }

    if (P->fwd4d == nullptr && P->fwd3d == nullptr && P->fwd == nullptr) {
        proj_log_error(P, _("axisswap: bad axis order"));
        return pj_default_destructor(P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
    }

    if (pj_param(P->ctx, P->params, "tangularunits").i) {
        P->left = PJ_IO_UNITS_RADIANS;
        P->right = PJ_IO_UNITS_RADIANS;
    } else {
        P->left = PJ_IO_UNITS_WHATEVER;
        P->right = PJ_IO_UNITS_WHATEVER;
    }

    /* Preparation and finalization steps are skipped, since the reason   */
    /* d'etre of axisswap is to bring input coordinates in line with the  */
    /* the internally expected order (ENU), such that handling of offsets */
    /* etc. can be done correctly in a later step of a pipeline */
    P->skip_fwd_prepare = 1;
    P->skip_fwd_finalize = 1;
    P->skip_inv_prepare = 1;
    P->skip_inv_finalize = 1;

    return P;
}

Coverage Report

Created: 2026-06-30 06:29

Line	Count	Source
1		/***********************************************************************
2
3		Axis order operation for use with transformation pipelines.
4
5		Kristian Evers, kreve@sdfe.dk, 2017-10-31
6
7		************************************************************************
8
9		Change the order and sign of 2,3 or 4 axes. Each of the possible four
10		axes are numbered with 1-4, such that the first input axis is 1, the
11		second is 2 and so on. The output ordering is controlled by a list of the
12		input axes re-ordered to the new mapping. Examples:
13
14		Reversing the order of the axes:
15
16		+proj=axisswap +order=4,3,2,1
17
18		Swapping the first two axes (x and y):
19
20		+proj=axisswap +order=2,1,3,4
21
22		The direction, or sign, of an axis can be changed by adding a minus in
23		front of the axis-number:
24
25		+proj=axisswap +order=1,-2,3,4
26
27		It is only necessary to specify the axes that are affected by the swap
28		operation:
29
30		+proj=axisswap +order=2,1
31
32		************************************************************************
33		* Copyright (c) 2017, Kristian Evers / SDFE
34		*
35		* Permission is hereby granted, free of charge, to any person obtaining a
36		* copy of this software and associated documentation files (the "Software"),
37		* to deal in the Software without restriction, including without limitation
38		* the rights to use, copy, modify, merge, publish, distribute, sublicense,
39		* and/or sell copies of the Software, and to permit persons to whom the
40		* Software is furnished to do so, subject to the following conditions:
41		*
42		* The above copyright notice and this permission notice shall be included
43		* in all copies or substantial portions of the Software.
44		*
45		* THE SOFTWARE IS PROVIDED "AS IS", WITHOUT WARRANTY OF ANY KIND, EXPRESS
46		* OR IMPLIED, INCLUDING BUT NOT LIMITED TO THE WARRANTIES OF MERCHANTABILITY,
47		* FITNESS FOR A PARTICULAR PURPOSE AND NONINFRINGEMENT. IN NO EVENT SHALL
48		* THE AUTHORS OR COPYRIGHT HOLDERS BE LIABLE FOR ANY CLAIM, DAMAGES OR OTHER
49		* LIABILITY, WHETHER IN AN ACTION OF CONTRACT, TORT OR OTHERWISE, ARISING
50		* FROM, OUT OF OR IN CONNECTION WITH THE SOFTWARE OR THE USE OR OTHER
51		* DEALINGS IN THE SOFTWARE.
52		*
53		***********************************************************************/
54
55		#include <errno.h>
56		#include <stdlib.h>
57		#include <string.h>
58
59		#include <algorithm>
60
61		#include "proj.h"
62		#include "proj_internal.h"
63
64		PROJ_HEAD(axisswap, "Axis ordering");
65
66		namespace { // anonymous namespace
67		struct pj_axisswap_data {
68		unsigned int axis[4];
69		int sign[4];
70		};
71		} // anonymous namespace
72
73	10.4k	static int sign(int x) { return (x > 0) - (x < 0); }
74
75	61.9k	static PJ_XY pj_axisswap_forward_2d(PJ_LP lp, PJ *P) {
76	61.9k	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
77	61.9k	PJ_XY xy;
78
79	61.9k	double in[2] = {lp.lam, lp.phi};
80	61.9k	xy.x = in[Q->axis[0]] * Q->sign[0];
81	61.9k	xy.y = in[Q->axis[1]] * Q->sign[1];
82	61.9k	return xy;
83	61.9k	}
84
85	0	static PJ_LP pj_axisswap_reverse_2d(PJ_XY xy, PJ *P) {
86	0	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
87	0	unsigned int i;
88	0	PJ_COORD out, in;
89
90	0	in.v[0] = xy.x;
91	0	in.v[1] = xy.y;
92	0	out = proj_coord_error();
93
94	0	for (i = 0; i < 2; i++)
95	0	out.v[Q->axis[i]] = in.v[i] * Q->sign[i];
96
97	0	return out.lp;
98	0	}
99
100	15.2k	static PJ_XYZ pj_axisswap_forward_3d(PJ_LPZ lpz, PJ *P) {
101	15.2k	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
102	15.2k	unsigned int i;
103	15.2k	PJ_COORD out, in;
104
105	15.2k	in.v[0] = lpz.lam;
106	15.2k	in.v[1] = lpz.phi;
107	15.2k	in.v[2] = lpz.z;
108	15.2k	out = proj_coord_error();
109
110	60.8k	for (i = 0; i < 3; i++)
111	45.6k	out.v[i] = in.v[Q->axis[i]] * Q->sign[i];
112
113	15.2k	return out.xyz;
114	15.2k	}
115
116	0	static PJ_LPZ pj_axisswap_reverse_3d(PJ_XYZ xyz, PJ *P) {
117	0	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
118	0	unsigned int i;
119	0	PJ_COORD in, out;
120
121	0	out = proj_coord_error();
122	0	in.v[0] = xyz.x;
123	0	in.v[1] = xyz.y;
124	0	in.v[2] = xyz.z;
125
126	0	for (i = 0; i < 3; i++)
127	0	out.v[Q->axis[i]] = in.v[i] * Q->sign[i];
128
129	0	return out.lpz;
130	0	}
131
132	244k	static void swap_xy_4d(PJ_COORD &coo, PJ *) {
133	244k	std::swap(coo.xyzt.x, coo.xyzt.y);
134	244k	}
135
136	0	static void pj_axisswap_forward_4d(PJ_COORD &coo, PJ *P) {
137	0	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
138	0	unsigned int i;
139	0	PJ_COORD out;
140
141	0	for (i = 0; i < 4; i++)
142	0	out.v[i] = coo.v[Q->axis[i]] * Q->sign[i];
143	0	coo = out;
144	0	}
145
146	0	static void pj_axisswap_reverse_4d(PJ_COORD &coo, PJ *P) {
147	0	struct pj_axisswap_data Q = (struct pj_axisswap_data )P->opaque;
148	0	unsigned int i;
149	0	PJ_COORD out;
150
151	0	for (i = 0; i < 4; i++)
152	0	out.v[Q->axis[i]] = coo.v[i] * Q->sign[i];
153
154	0	coo = out;
155	0	}
156
157		/***********************************************************************/
158	7.14k	PJ *PJ_CONVERSION(axisswap, 0) {
159		/***********************************************************************/
160	7.14k	struct pj_axisswap_data Q = static_cast<struct pj_axisswap_data >(
161	7.14k	calloc(1, sizeof(struct pj_axisswap_data)));
162	7.14k	char *s;
163	7.14k	unsigned int i, j, n = 0;
164
165	7.14k	if (nullptr == Q)
166	0	return pj_default_destructor(P, PROJ_ERR_OTHER /ENOMEM/);
167	7.14k	P->opaque = (void *)Q;
168
169		/* +order and +axis are mutually exclusive */
170	7.14k	if (!pj_param_exists(P->params, "order") ==
171	7.14k	!pj_param_exists(P->params, "axis")) {
172	9	proj_log_error(P,
173	9	_("must provide EITHER 'order' OR 'axis' parameter."));
174	9	return pj_default_destructor(
175	9	P, PROJ_ERR_INVALID_OP_MUTUALLY_EXCLUSIVE_ARGS);
176	9	}
177
178		/* fill axis list with indices from 4-7 to simplify duplicate search further
179		* down */
180	35.6k	for (i = 0; i < 4; i++) {
181	28.5k	Q->axis[i] = i + 4;
182	28.5k	Q->sign[i] = 1;
183	28.5k	}
184
185		/* if the "order" parameter is used */
186	7.13k	if (pj_param_exists(P->params, "order")) {
187		/* read axis order */
188	5.24k	char *order = pj_param(P->ctx, P->params, "sorder").s;
189
190		/* check that all characters are valid */
191	24.0k	for (i = 0; i < strlen(order); i++)
192	18.8k	if (strchr("1234-,", order[i]) == nullptr) {
193	7	proj_log_error(P, _("unknown axis '%c'"), order[i]);
194	7	return pj_default_destructor(
195	7	P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
196	7	}
197
198		/* read axes numbers and signs */
199	5.24k	s = order;
200	5.24k	n = 0;
201	15.7k	while (*s != '\0' && n < 4) {
202	10.5k	Q->axis[n] = abs(atoi(s)) - 1;
203	10.5k	if (Q->axis[n] > 3) {
204	7	proj_log_error(P, _("invalid axis '%d'"), Q->axis[n]);
205	7	return pj_default_destructor(
206	7	P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
207	7	}
208	10.4k	Q->sign[n++] = sign(atoi(s));
209	23.9k	while (s != '\0' && s != ',')
210	13.4k	s++;
211	10.4k	if (*s == ',')
212	5.26k	s++;
213	10.4k	}
214	5.24k	}
215
216		/* if the "axis" parameter is used */
217	7.12k	if (pj_param_exists(P->params, "axis")) {
218		/* parse the classic PROJ.4 enu axis specification */
219	7.55k	for (i = 0; i < 3; i++) {
220	5.66k	switch (P->axis[i]) {
221	1.23k	case 'w':
222	1.23k	Q->sign[i] = -1;
223	1.23k	Q->axis[i] = 0;
224	1.23k	break;
225	655	case 'e':
226	655	Q->sign[i] = 1;
227	655	Q->axis[i] = 0;
228	655	break;
229	1.45k	case 's':
230	1.45k	Q->sign[i] = -1;
231	1.45k	Q->axis[i] = 1;
232	1.45k	break;
233	434	case 'n':
234	434	Q->sign[i] = 1;
235	434	Q->axis[i] = 1;
236	434	break;
237	973	case 'd':
238	973	Q->sign[i] = -1;
239	973	Q->axis[i] = 2;
240	973	break;
241	913	case 'u':
242	913	Q->sign[i] = 1;
243	913	Q->axis[i] = 2;
244	913	break;
245	0	default:
246	0	proj_log_error(P, _("unknown axis '%c'"), P->axis[i]);
247	0	return pj_default_destructor(
248	0	P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
249	5.66k	}
250	5.66k	}
251	1.88k	n = 3;
252	1.88k	}
253
254		/* check for duplicate axes */
255	35.5k	for (i = 0; i < 4; i++)
256	142k	for (j = 0; j < 4; j++) {
257	113k	if (i == j)
258	28.4k	continue;
259	85.3k	if (Q->axis[i] == Q->axis[j]) {
260	16	proj_log_error(P, _("axisswap: duplicate axes specified"));
261	16	return pj_default_destructor(
262	16	P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
263	16	}
264	85.3k	}
265
266		/* only map fwd/inv functions that are possible with the given axis setup */
267	7.10k	if (n == 4) {
268	1	P->fwd4d = pj_axisswap_forward_4d;
269	1	P->inv4d = pj_axisswap_reverse_4d;
270	1	}
271	7.10k	if (n == 3 && Q->axis[0] < 3 && Q->axis[1] < 3 && Q->axis[2] < 3) {
272	1.88k	P->fwd3d = pj_axisswap_forward_3d;
273	1.88k	P->inv3d = pj_axisswap_reverse_3d;
274	1.88k	}
275	7.10k	if (n == 2) {
276	5.21k	if (Q->axis[0] == 1 && Q->sign[0] == 1 && Q->axis[1] == 0 &&
277	3.63k	Q->sign[1] == 1) {
278	3.57k	P->fwd4d = swap_xy_4d;
279	3.57k	P->inv4d = swap_xy_4d;
280	3.57k	} else if (Q->axis[0] < 2 && Q->axis[1] < 2) {
281	1.63k	P->fwd = pj_axisswap_forward_2d;
282	1.63k	P->inv = pj_axisswap_reverse_2d;
283	1.63k	}
284	5.21k	}
285
286	7.10k	if (P->fwd4d == nullptr && P->fwd3d == nullptr && P->fwd == nullptr) {
287	11	proj_log_error(P, _("axisswap: bad axis order"));
288	11	return pj_default_destructor(P, PROJ_ERR_INVALID_OP_ILLEGAL_ARG_VALUE);
289	11	}
290
291	7.09k	if (pj_param(P->ctx, P->params, "tangularunits").i) {
292	0	P->left = PJ_IO_UNITS_RADIANS;
293	0	P->right = PJ_IO_UNITS_RADIANS;
294	7.09k	} else {
295	7.09k	P->left = PJ_IO_UNITS_WHATEVER;
296	7.09k	P->right = PJ_IO_UNITS_WHATEVER;
297	7.09k	}
298
299		/* Preparation and finalization steps are skipped, since the reason */
300		/* d'etre of axisswap is to bring input coordinates in line with the */
301		/* the internally expected order (ENU), such that handling of offsets */
302		/* etc. can be done correctly in a later step of a pipeline */
303	7.09k	P->skip_fwd_prepare = 1;
304	7.09k	P->skip_fwd_finalize = 1;
305	7.09k	P->skip_inv_prepare = 1;
306	7.09k	P->skip_inv_finalize = 1;
307
308	7.09k	return P;
309	7.10k	}