ptarray_transform()

int ptarray_transform	(	POINTARRAY *	pa,
		LWPROJ *	pj
	)

Definition at line 257 of file liblwgeom/lwgeom_transform.c.

{
        uint32_t i;
        POINT4D p;
        size_t n_converted;
        size_t n_points = pa->npoints;
        size_t point_size = ptarray_point_size(pa);
        int has_z = ptarray_has_z(pa);
        double *pa_double = (double*)(pa->serialized_pointlist);
 
        /* Convert to radians if necessary */
        if (proj_angular_input(pj->pj, PJ_FWD))
        {
                for (i = 0; i < pa->npoints; i++)
                {
                        getPoint4d_p(pa, i, &p);
                        to_rad(&p);
                }
        }
 
        if (n_points == 1)
        {
                /* For single points it's faster to call proj_trans */
                PJ_XYZT v = {pa_double[0], pa_double[1], has_z ? pa_double[2] : 0.0, 0.0};
                PJ_COORD c;
                c.xyzt = v;
                PJ_COORD t = proj_trans(pj->pj, PJ_FWD, c);
 
                int pj_errno_val = proj_errno_reset(pj->pj);
                if (pj_errno_val)
                {
                        lwerror("transform: %s (%d)", proj_errno_string(pj_errno_val), pj_errno_val);
                        return LW_FAILURE;
                }
                pa_double[0] = (t.xyzt).x;
                pa_double[1] = (t.xyzt).y;
                if (has_z)
                        pa_double[2] = (t.xyzt).z;
        }
        else
        {
                /*
                 * size_t proj_trans_generic(PJ *P, PJ_DIRECTION direction,
                 * double *x, size_t sx, size_t nx,
                 * double *y, size_t sy, size_t ny,
                 * double *z, size_t sz, size_t nz,
                 * double *t, size_t st, size_t nt)
                 */
 
                n_converted = proj_trans_generic(pj->pj,
                                                 PJ_FWD,
                                                 pa_double,
                                                 point_size,
                                                 n_points, /* X */
                                                 pa_double + 1,
                                                 point_size,
                                                 n_points, /* Y */
                                                 has_z ? pa_double + 2 : NULL,
                                                 has_z ? point_size : 0,
                                                 has_z ? n_points : 0, /* Z */
                                                 NULL,
                                                 0,
                                                 0 /* M */
                );
 
                if (n_converted != n_points)
                {
                        lwerror("ptarray_transform: converted (%d) != input (%d)", n_converted, n_points);
                        return LW_FAILURE;
                }
 
                int pj_errno_val = proj_errno_reset(pj->pj);
                if (pj_errno_val)
                {
                        lwerror("transform: %s (%d)", proj_errno_string(pj_errno_val), pj_errno_val);
                        return LW_FAILURE;
                }
        }
 
        /* Convert radians to degrees if necessary */
        if (proj_angular_output(pj->pj, PJ_FWD))
        {
                for (i = 0; i < pa->npoints; i++)
                {
                        getPoint4d_p(pa, i, &p);
                        to_dec(&p);
                }
        }
 
        return LW_SUCCESS;
}

References getPoint4d_p(), LW_FAILURE, LW_SUCCESS, lwerror(), POINTARRAY::npoints, LWPROJ::pj, ptarray_has_z(), ptarray_point_size(), POINTARRAY::serialized_pointlist, to_dec(), to_rad(), pixval::x, and pixval::y.

Referenced by gml_reproject_pa(), and lwgeom_transform().

Here is the call graph for this function:

Here is the caller graph for this function: