lwgeom_nudge_geodetic()

int lwgeom_nudge_geodetic

(

LWGEOM *

geom

)

When features are snapped or sometimes they are just this way, they are very close to the geodetic bounds but slightly over.

Gently move coordinates of LWGEOM if they are close enough into geodetic range.

This routine nudges those points, and only those points, back over to the bounds. http://trac.osgeo.org/postgis/ticket/1292

Definition at line 3388 of file lwgeodetic.c.

References LWCOLLECTION::geoms, LINETYPE, LW_FALSE, LW_TRUE, lwerror(), lwgeom_is_empty(), lwtype_is_collection(), lwtype_name(), LWCOLLECTION::ngeoms, LWPOLY::nrings, POINTTYPE, POLYGONTYPE, ptarray_nudge_geodetic(), LWPOLY::rings, TRIANGLETYPE, ovdump::type, and LWGEOM::type.

Referenced by geography_from_geometry(), and gserialized_geography_from_lwgeom().

{
        int type;
        int i = 0;
        int rv = LW_FALSE;

        assert(geom);

        /* No points in nothing */
        if ( lwgeom_is_empty(geom) )
                return LW_FALSE;

        type = geom->type;

        if ( type == POINTTYPE )
                return ptarray_nudge_geodetic(((LWPOINT*)geom)->point);

        if ( type == LINETYPE )
                return ptarray_nudge_geodetic(((LWLINE*)geom)->points);

        if ( type == POLYGONTYPE )
        {
                LWPOLY *poly = (LWPOLY*)geom;
                for ( i = 0; i < poly->nrings; i++ )
                {
                        int n = ptarray_nudge_geodetic(poly->rings[i]);
                        rv = (rv == LW_TRUE ? rv : n);
                }
                return rv;
        }

        if ( type == TRIANGLETYPE )
                return ptarray_nudge_geodetic(((LWTRIANGLE*)geom)->points);

        if ( lwtype_is_collection( type ) )
        {
                LWCOLLECTION *col = (LWCOLLECTION*)geom;

                for ( i = 0; i < col->ngeoms; i++ )
                {
                        int n = lwgeom_nudge_geodetic(col->geoms[i]);
                        rv = (rv == LW_TRUE ? rv : n);
                }
                return rv;
        }

        lwerror("unsupported type (%s) passed to lwgeom_nudge_geodetic", lwtype_name(type));
        return rv;
}

Here is the call graph for this function:

Here is the caller graph for this function: