/**
* This should be rewritten to make use of the curve itself.
*/
double
lwcurvepoly_area(const LWCURVEPOLY *curvepoly)
{
double area = 0.0;
LWPOLY *poly;
if( lwgeom_is_empty((LWGEOM*)curvepoly) )
return 0.0;
poly = lwcurvepoly_segmentize(curvepoly, 32);
area = lwpoly_area(poly);
lwpoly_free(poly);
return area;
}
double lwgeom_area(const LWGEOM *geom)
{
int type = geom->type;
if ( type == POLYGONTYPE )
return lwpoly_area((LWPOLY*)geom);
else if ( type == CURVEPOLYTYPE )
return lwcurvepoly_area((LWCURVEPOLY*)geom);
else if (type == TRIANGLETYPE )
return lwtriangle_area((LWTRIANGLE*)geom);
else if ( lwgeom_is_collection(geom) )
{
double area = 0.0;
int i;
LWCOLLECTION *col = (LWCOLLECTION*)geom;
for ( i = 0; i < col->ngeoms; i++ )
area += lwgeom_area(col->geoms[i]);
return area;
}
else
return 0.0;
}