LWGEOM_snaptogrid_pointoff()

Datum LWGEOM_snaptogrid_pointoff

(

PG_FUNCTION_ARGS

)

Definition at line 383 of file lwgeom_functions_analytic.c.

References LWGEOM::bbox, LWPOINT::flags, FLAGS_GET_M, FLAGS_GET_Z, geometry_serialize(), getPoint4d_p(), gserialized_is_empty(), gridspec_t::ipm, gridspec_t::ipx, gridspec_t::ipy, gridspec_t::ipz, lwgeom_add_bbox(), lwgeom_as_lwpoint(), lwgeom_drop_bbox(), lwgeom_from_gserialized(), lwgeom_grid(), lwtype_name(), POINT4D::m, gridspec_t::msize, PG_FUNCTION_INFO_V1(), LWPOINT::point, ST_LineCrossingDirection(), LWGEOM::type, POINT4D::x, gridspec_t::xsize, POINT4D::y, gridspec_t::ysize, POINT4D::z, and gridspec_t::zsize.

Referenced by LWGEOM_snaptogrid().

{
        GSERIALIZED *in_geom, *in_point;
        LWGEOM *in_lwgeom;
        LWPOINT *in_lwpoint;
        GSERIALIZED *out_geom = NULL;
        LWGEOM *out_lwgeom;
        gridspec grid;
        /* BOX3D box3d; */
        POINT4D offsetpoint;

        in_geom = PG_GETARG_GSERIALIZED_P(0);

        /* Return input geometry if input geometry is empty */
        if ( gserialized_is_empty(in_geom) )
        {
                PG_RETURN_POINTER(in_geom);
        }

        in_point = PG_GETARG_GSERIALIZED_P(1);
        in_lwpoint = lwgeom_as_lwpoint(lwgeom_from_gserialized(in_point));
        if ( in_lwpoint == NULL )
        {
                lwpgerror("Offset geometry must be a point");
        }

        grid.xsize = PG_GETARG_FLOAT8(2);
        grid.ysize = PG_GETARG_FLOAT8(3);
        grid.zsize = PG_GETARG_FLOAT8(4);
        grid.msize = PG_GETARG_FLOAT8(5);

        /* Take offsets from point geometry */
        getPoint4d_p(in_lwpoint->point, 0, &offsetpoint);
        grid.ipx = offsetpoint.x;
        grid.ipy = offsetpoint.y;
        if (FLAGS_GET_Z(in_lwpoint->flags) ) grid.ipz = offsetpoint.z;
        else grid.ipz=0;
        if (FLAGS_GET_M(in_lwpoint->flags) ) grid.ipm = offsetpoint.m;
        else grid.ipm=0;

#if POSTGIS_DEBUG_LEVEL >= 4
        grid_print(&grid);
#endif

        /* Return input geometry if input grid is meaningless */
        if ( grid.xsize==0 && grid.ysize==0 && grid.zsize==0 && grid.msize==0 )
        {
                PG_RETURN_POINTER(in_geom);
        }

        in_lwgeom = lwgeom_from_gserialized(in_geom);

        POSTGIS_DEBUGF(3, "SnapToGrid got a %s", lwtype_name(in_lwgeom->type));

        out_lwgeom = lwgeom_grid(in_lwgeom, &grid);
        if ( out_lwgeom == NULL ) PG_RETURN_NULL();

        /* COMPUTE_BBOX TAINTING */
        if (in_lwgeom->bbox)
        {
                lwgeom_drop_bbox(out_lwgeom);
                lwgeom_add_bbox(out_lwgeom);
        }

        POSTGIS_DEBUGF(3, "SnapToGrid made a %s", lwtype_name(out_lwgeom->type));

        out_geom = geometry_serialize(out_lwgeom);

        PG_RETURN_POINTER(out_geom);
}

Here is the call graph for this function:

Here is the caller graph for this function: