lwgeom_length_spheroid()

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

extern

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 3297 of file lwgeodetic.c.

{
        int type;
        uint32_t 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;
}

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

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

Here is the call graph for this function:

Here is the caller graph for this function: