CircTreePIP()

static int CircTreePIP ( const CIRC_NODE *  tree1, const GSERIALIZED *  g1, const POINT4D *  in_point  )

static

Definition at line 104 of file geography_measurement_trees.c.

References circ_tree_contains_point(), circ_tree_get_point_outside(), gbox_contains_point3d(), gbox_pt_outside(), geog2cart(), geographic_point_init(), gserialized_get_gbox_p(), gserialized_get_type(), LW_FAILURE, LW_FALSE, lwgeom_calculate_gbox_geodetic(), lwgeom_free(), lwgeom_from_gserialized(), MULTIPOLYGONTYPE, POLYGONTYPE, POINT2D::x, POINT4D::x, POINT2D::y, and POINT4D::y.

Referenced by geography_distance_cache_tolerance(), and geography_tree_distance().

{
        int tree1_type = gserialized_get_type(g1);
        GBOX gbox1;
        GEOGRAPHIC_POINT in_gpoint;
        POINT3D in_point3d;

        POSTGIS_DEBUGF(3, "tree1_type=%d", tree1_type);

        /* If the tree'ed argument is a polygon, do the P-i-P using the tree-based P-i-P */
        if ( tree1_type == POLYGONTYPE || tree1_type == MULTIPOLYGONTYPE )
        {
                POSTGIS_DEBUG(3, "tree is a polygon, using tree PiP");
                /* Need a gbox to calculate an outside point */
                if ( LW_FAILURE == gserialized_get_gbox_p(g1, &gbox1) )
                {
                        LWGEOM* lwgeom1 = lwgeom_from_gserialized(g1);
                        POSTGIS_DEBUG(3, "unable to read gbox from gserialized, calculating from scratch");
                        lwgeom_calculate_gbox_geodetic(lwgeom1, &gbox1);
                        lwgeom_free(lwgeom1);
                }

                /* Flip the candidate point into geographics */
                geographic_point_init(in_point->x, in_point->y, &in_gpoint);
                geog2cart(&in_gpoint, &in_point3d);

                /* If the candidate isn't in the tree box, it's not in the tree area */
                if ( ! gbox_contains_point3d(&gbox1, &in_point3d) )
                {
                        POSTGIS_DEBUG(3, "in_point3d is not inside the tree gbox, CircTreePIP returning FALSE");
                        return LW_FALSE;
                }
                /* The candidate point is in the box, so it *might* be inside the tree */
                else
                {
                        POINT2D pt2d_outside; /* latlon */
                        POINT2D pt2d_inside;
                        pt2d_inside.x = in_point->x;
                        pt2d_inside.y = in_point->y;
                        /* Calculate a definitive outside point */
                        if (gbox_pt_outside(&gbox1, &pt2d_outside) == LW_FAILURE)
                                if (circ_tree_get_point_outside(tree1, &pt2d_outside) == LW_FAILURE)
                                        lwpgerror("%s: Unable to generate outside point!", __func__);

                        POSTGIS_DEBUGF(3, "p2d_inside=POINT(%g %g) p2d_outside=POINT(%g %g)", pt2d_inside.x, pt2d_inside.y, pt2d_outside.x, pt2d_outside.y);
                        /* Test the candidate point for strict containment */
                        POSTGIS_DEBUG(3, "calling circ_tree_contains_point for PiP test");
                        return circ_tree_contains_point(tree1, &pt2d_inside, &pt2d_outside, 0, NULL);
                }
        }
        else
        {
                POSTGIS_DEBUG(3, "tree1 not polygonal, so CircTreePIP returning FALSE");
                return LW_FALSE;
        }
}

Here is the call graph for this function:

Here is the caller graph for this function: