ptarray_to_GEOSCoordSeq()

GEOSCoordSeq ptarray_to_GEOSCoordSeq	(	const POINTARRAY *	pa,
		uint8_t	fix_ring
	)

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

{
        uint32_t dims = 2;
        uint32_t i;
        int append_points = 0;
        const POINT3D *p3d = NULL;
        const POINT2D* p2d = NULL;
        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 POSTGIS_GEOS_VERSION >= 31000
        if (append_points == 0) {
                sq = GEOSCoordSeq_copyFromBuffer((const double*) pa->serialized_pointlist, pa->npoints, FLAGS_GET_Z(pa->flags), FLAGS_GET_M(pa->flags));
                if (!sq)
                {
                        lwerror("Error creating GEOS Coordinate Sequence");
                        return NULL;
                }
                return sq;
        }
#endif
        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 = getPoint3d_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);
                }
 
#if POSTGIS_GEOS_VERSION < 30800
                GEOSCoordSeq_setX(sq, i, p2d->x);
                GEOSCoordSeq_setY(sq, i, p2d->y);
                if (dims == 3) GEOSCoordSeq_setZ(sq, i, p3d->z);
#else
                if (dims == 3)
                        GEOSCoordSeq_setXYZ(sq, i, p2d->x, p2d->y, p3d->z);
                else
                        GEOSCoordSeq_setXY(sq, i, p2d->x, p2d->y);
#endif
 
        }
 
        if (append_points)
        {
                if (dims == 3)
                {
                        p3d = getPoint3d_cp(pa, 0);
                        p2d = (const POINT2D*)p3d;
                }
                else
                        p2d = getPoint2d_cp(pa, 0);
                for (i = pa->npoints; i < pa->npoints + append_points; i++)
                {
#if POSTGIS_GEOS_VERSION < 30800
                        GEOSCoordSeq_setX(sq, i, p2d->x);
                        GEOSCoordSeq_setY(sq, i, p2d->y);
#else
                        GEOSCoordSeq_setXY(sq, i, p2d->x, p2d->y);
#endif
 
                        if (dims == 3) GEOSCoordSeq_setZ(sq, i, p3d->z);
                }
        }
 
        return sq;
 
}

References POINTARRAY::flags, FLAGS_GET_M, FLAGS_GET_Z, getPoint2d_cp(), getPoint3d_cp(), LWDEBUGF, lwerror(), POINTARRAY::npoints, ptarray_is_closed_2d(), POINTARRAY::serialized_pointlist, POINT2D::x, POINT3D::x, POINT2D::y, POINT3D::y, and POINT3D::z.

Referenced by LWGEOM2GEOS(), and ptarray_to_GEOSLinearRing().

Here is the call graph for this function:

Here is the caller graph for this function: