double lwcompound_length(const LWCOMPOUND *comp)
{
double length = 0.0;
LWLINE *line;
if ( lwgeom_is_empty((LWGEOM*)comp) )
return 0.0;
line = lwcompound_stroke(comp, 32);
length = lwline_length(line);
lwline_free(line);
return length;
}
double lwcircstring_length(const LWCIRCSTRING *circ)
{
double length = 0.0;
LWLINE *line;
if ( lwcircstring_is_empty(circ) )
return 0.0;
line = lwcircstring_segmentize(circ, 32);
length = lwline_length(line);
lwline_free(line);
return length;
}
double lwgeom_length(const LWGEOM *geom)
{
int type = geom->type;
if ( type == LINETYPE )
return lwline_length((LWLINE*)geom);
else if ( type == CIRCSTRINGTYPE )
return lwcircstring_length((LWCIRCSTRING*)geom);
else if ( type == COMPOUNDTYPE )
return lwcompound_length((LWCOMPOUND*)geom);
else if ( lwgeom_is_collection(geom) )
{
double length = 0.0;
int i;
LWCOLLECTION *col = (LWCOLLECTION*)geom;
for ( i = 0; i < col->ngeoms; i++ )
length += lwgeom_length(col->geoms[i]);
return length;
}
else
return 0.0;
}
