lwgeom_length_spheroid()

double lwgeom_length_spheroid	(	const LWGEOM *	geom,
		const SPHEROID *	s
	)

Calculate the geodetic length of a lwgeom on the unit sphere.

The result will have to by multiplied by the real radius to get the real length.

Definition at line 3281 of file lwgeodetic.c.

References LWCOLLECTION::geoms, LINETYPE, lwerror(), lwgeom_is_empty(), lwtype_is_collection(), MULTIPOINTTYPE, LWCOLLECTION::ngeoms, LWPOLY::nrings, POINTTYPE, POLYGONTYPE, ptarray_length_spheroid(), LWPOLY::rings, TRIANGLETYPE, ovdump::type, and LWGEOM::type.

Referenced by geography_length(), geography_perimeter(), LWGEOM_length2d_ellipsoid(), and LWGEOM_length_ellipsoid_linestring().

{
        int type;
        int i = 0;
        double length = 0.0;

        assert(geom);

        /* No area in nothing */
        if ( lwgeom_is_empty(geom) )
                return 0.0;

        type = geom->type;

        if ( type == POINTTYPE || type == MULTIPOINTTYPE )
                return 0.0;

        if ( type == LINETYPE )
                return ptarray_length_spheroid(((LWLINE*)geom)->points, s);

        if ( type == POLYGONTYPE )
        {
                LWPOLY *poly = (LWPOLY*)geom;
                for ( i = 0; i < poly->nrings; i++ )
                {
                        length += ptarray_length_spheroid(poly->rings[i], s);
                }
                return length;
        }

        if ( type == TRIANGLETYPE )
                return ptarray_length_spheroid(((LWTRIANGLE*)geom)->points, s);

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

                for ( i = 0; i < col->ngeoms; i++ )
                {
                        length += lwgeom_length_spheroid(col->geoms[i], s);
                }
                return length;
        }

        lwerror("unsupported type passed to lwgeom_length_sphere");
        return 0.0;
}

Here is the call graph for this function:

Here is the caller graph for this function: