GEOSCoordSeq ptarray_to_GEOSCoordSeq

(

const POINTARRAY *

pa

)

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

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

Referenced by LWGEOM2GEOS().

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

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

        if ( ! (sq = GEOSCoordSeq_create(pa->npoints, dims)) ) 
                lwerror("Error creating GEOS Coordinate Sequence");

        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);
                }

#if POSTGIS_GEOS_VERSION < 33
                /* Make sure we don't pass any infinite values down into GEOS */
                /* GEOS 3.3+ is supposed to  handle this stuff OK */
                if ( isinf(p2d->x) || isinf(p2d->y) || (dims == 3 && isinf(p3d->z)) )
                        lwerror("Infinite coordinate value found in geometry.");
                if ( isnan(p2d->x) || isnan(p2d->y) || (dims == 3 && isnan(p3d->z)) )
                        lwerror("NaN coordinate value found in geometry.");
#endif

                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: