lw_dist3d_poly_poly()

int lw_dist3d_poly_poly	(	LWPOLY *	poly1,
		LWPOLY *	poly2,
		DISTPTS3D *	dl
	)

polygon to polygon calculation

Definition at line 691 of file measures3d.c.

References define_plane(), DIST_MAX, DISTPTS3D::distance, lw_dist3d_ptarray_poly(), lw_dist3d_ptarray_ptarray(), LW_FALSE, LW_TRUE, LWDEBUG, DISTPTS3D::mode, LWPOLY::rings, DISTPTS3D::tolerance, and DISTPTS3D::twisted.

Referenced by lw_dist3d_distribute_bruteforce().

{
        PLANE3D plane1, plane2;
        int planedef1, planedef2;
        LWDEBUG(2, "lw_dist3d_poly_poly is called");
        if (dl->mode == DIST_MAX)
        {
                return lw_dist3d_ptarray_ptarray(poly1->rings[0], poly2->rings[0], dl);
        }

        planedef1 = define_plane(poly1->rings[0], &plane1);
        planedef2 = define_plane(poly2->rings[0], &plane2);

        if (!planedef1 || !planedef2)
        {
                if (!planedef1 && !planedef2)
                {
                        /* Neither polygon define a plane: Return distance line to line */
                        return lw_dist3d_ptarray_ptarray(poly1->rings[0], poly2->rings[0], dl);
                }

                if (!planedef1)
                {
                        /* Only poly2 defines a plane: Return distance from line (poly1) to poly2 */
                        return lw_dist3d_ptarray_poly(poly1->rings[0], poly2, &plane2, dl);
                }

                /* Only poly1 defines a plane: Return distance from line (poly2) to poly1 */
                return lw_dist3d_ptarray_poly(poly2->rings[0], poly1, &plane1, dl);
        }

        /*What we do here is to compare the boundary of one polygon with the other polygon
        and then take the second boundary comparing with the first polygon*/
        dl->twisted=1;
        if (!lw_dist3d_ptarray_poly(poly1->rings[0], poly2, &plane2, dl))
                return LW_FALSE;
        if (dl->distance < dl->tolerance) /*Just check if the answer already is given*/
                return LW_TRUE;

        dl->twisted=-1; /*because we switch the order of geometries we switch "twisted" to -1 which will give the right order of points in shortest line.*/
        return lw_dist3d_ptarray_poly(poly2->rings[0], poly1, &plane1, dl);
}

Here is the call graph for this function:

Here is the caller graph for this function: