ptarray_grid()

POINTARRAY* ptarray_grid	(	const POINTARRAY *	pa,
		const gridspec *	grid
	)

Definition at line 1890 of file ptarray.c.

References POINTARRAY::flags, FLAGS_GET_M, FLAGS_GET_Z, getPoint4d_p(), gridspec_t::ipm, gridspec_t::ipx, gridspec_t::ipy, gridspec_t::ipz, LW_FALSE, LWDEBUGF, POINT4D::m, gridspec_t::msize, POINTARRAY::npoints, ptarray_append_point(), ptarray_construct_empty(), POINT4D::x, gridspec_t::xsize, POINT4D::y, gridspec_t::ysize, POINT4D::z, and gridspec_t::zsize.

Referenced by lwcircstring_grid(), lwline_grid(), lwpoint_grid(), and lwpoly_grid().

{
        POINT4D pt;
        int ipn; /* input point numbers */
        POINTARRAY *dpa;

        LWDEBUGF(2, "ptarray_grid called on %p", pa);

        dpa = ptarray_construct_empty(FLAGS_GET_Z(pa->flags),FLAGS_GET_M(pa->flags), pa->npoints);

        for (ipn=0; ipn<pa->npoints; ++ipn)
        {

                getPoint4d_p(pa, ipn, &pt);

                if ( grid->xsize )
                        pt.x = rint((pt.x - grid->ipx)/grid->xsize) *
                                 grid->xsize + grid->ipx;

                if ( grid->ysize )
                        pt.y = rint((pt.y - grid->ipy)/grid->ysize) *
                                 grid->ysize + grid->ipy;

                if ( FLAGS_GET_Z(pa->flags) && grid->zsize )
                        pt.z = rint((pt.z - grid->ipz)/grid->zsize) *
                                 grid->zsize + grid->ipz;

                if ( FLAGS_GET_M(pa->flags) && grid->msize )
                        pt.m = rint((pt.m - grid->ipm)/grid->msize) *
                                 grid->msize + grid->ipm;

                ptarray_append_point(dpa, &pt, LW_FALSE);

        }

        return dpa;
}

Here is the call graph for this function:

Here is the caller graph for this function: