ptarray_nudge_geodetic()

static int ptarray_nudge_geodetic ( POINTARRAY *  pa )

static

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

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

Definition at line 3336 of file lwgeodetic.c.

References getPoint4d_p(), LW_FALSE, LW_TRUE, lwerror(), POINTARRAY::npoints, ptarray_set_point4d(), POINT4D::x, and POINT4D::y.

Referenced by lwgeom_nudge_geodetic().

{

        int i;
        POINT4D p;
        int altered = LW_FALSE;
        int rv = LW_FALSE;
        static double tolerance = 1e-10;

        if ( ! pa )
                lwerror("ptarray_nudge_geodetic called with null input");

        for(i = 0; i < pa->npoints; i++ )
        {
                getPoint4d_p(pa, i, &p);
                if ( p.x < -180.0 && (-180.0 - p.x < tolerance) )
                {
                        p.x = -180.0;
                        altered = LW_TRUE;
                }
                if ( p.x > 180.0 && (p.x - 180.0 < tolerance) )
                {
                        p.x = 180.0;
                        altered = LW_TRUE;
                }
                if ( p.y < -90.0 && (-90.0 - p.y < tolerance) )
                {
                        p.y = -90.0;
                        altered = LW_TRUE;
                }
                if ( p.y > 90.0 && (p.y - 90.0 < tolerance) )
                {
                        p.y = 90.0;
                        altered = LW_TRUE;
                }
                if ( altered == LW_TRUE )
                {
                        ptarray_set_point4d(pa, i, &p);
                        altered = LW_FALSE;
                        rv = LW_TRUE;
                }
        }
        return rv;
}

Here is the call graph for this function:

Here is the caller graph for this function: