ptarray_to_GEOSCoordSeq()

GEOSCoordSeq ptarray_to_GEOSCoordSeq	(	const POINTARRAY *	pa,
		int	fix_ring
	)

Definition at line 215 of file liblwgeom/lwgeom_geos.c.

References POINTARRAY::flags, FLAGS_GET_Z, getPoint2d_cp(), getPoint3dz_cp(), LWDEBUGF, lwerror(), POINTARRAY::npoints, ptarray_is_closed_2d(), POINT2D::x, POINT3DZ::x, POINT2D::y, POINT3DZ::y, and POINT3DZ::z.

Referenced by GEOS2LWGEOM(), LWGEOM2GEOS(), and ptarray_to_GEOSLinearRing().

{
        uint32_t dims = 2;
        uint32_t i;
        int append_points = 0;
        const POINT3DZ *p3d;
        const POINT2D *p2d;
        GEOSCoordSeq sq;

        if ( FLAGS_GET_Z(pa->flags) )
                dims = 3;

        if ( fix_ring )
        {
                if (pa->npoints < 1)
                {
                        lwerror("ptarray_to_GEOSCoordSeq called with fix_ring and 0 vertices in ring, cannot fix");
                        return NULL;
                }
                else
                {
                        if ( pa->npoints < 4 )
                        {
                                append_points = 4 - pa->npoints;
                        }
                        if ( ! ptarray_is_closed_2d(pa) && append_points == 0 )
                        {
                                append_points = 1;
                        }
                }
        }

        if ( ! (sq = GEOSCoordSeq_create(pa->npoints + append_points, dims)) )
        {
                lwerror("Error creating GEOS Coordinate Sequence");
                return NULL;
        }

        for ( i=0; i < pa->npoints; i++ )
        {
                if ( dims == 3 )
                {
                        p3d = getPoint3dz_cp(pa, i);
                        p2d = (const POINT2D *)p3d;
                        LWDEBUGF(4, "Point: %g,%g,%g", p3d->x, p3d->y, p3d->z);
                }
                else
                {
                        p2d = getPoint2d_cp(pa, i);
                        LWDEBUGF(4, "Point: %g,%g", p2d->x, p2d->y);
                }

                GEOSCoordSeq_setX(sq, i, p2d->x);
                GEOSCoordSeq_setY(sq, i, p2d->y);

                if ( dims == 3 )
                {
                        GEOSCoordSeq_setZ(sq, i, p3d->z);
                }
        }

        if ( append_points )
        {
                if ( dims == 3 )
                {
                        p3d = getPoint3dz_cp(pa, 0);
                        p2d = (const POINT2D *)p3d;
                }
                else
                {
                        p2d = getPoint2d_cp(pa, 0);
                }
                for ( i = pa->npoints; i < pa->npoints + append_points; i++ )
                {
                        GEOSCoordSeq_setX(sq, i, p2d->x);
                        GEOSCoordSeq_setY(sq, i, p2d->y);

                        if ( dims == 3 )
                        {
                                GEOSCoordSeq_setZ(sq, i, p3d->z);
                        }
                }
        }

        return sq;
}

Here is the call graph for this function:

Here is the caller graph for this function: