void object::test<1>()
{
geom1_ = GEOSGeomFromWKT("LINESTRING(0 0, 10 10, 10 0, 0 10)");
geom2_ = GEOSNode(geom1_);
ensure(0 != geom2_);
GEOSNormalize(geom2_);
char* wkt_c = GEOSWKTWriter_write(w_, geom2_);
std::string out(wkt_c);
free(wkt_c);
ensure_equals(out,
"MULTILINESTRING ((5 5, 10 0, 10 10, 5 5), (0 10, 5 5), (0 0, 5 5))"
);
}
void object::test<2>()
{
geom1_ = GEOSGeomFromWKT("MULTILINESTRING((0 0, 2 0, 4 0),(5 0, 3 0, 1 0))");
geom2_ = GEOSNode(geom1_);
ensure(0 != geom2_);
GEOSNormalize(geom2_);
char* wkt_c = GEOSWKTWriter_write(w_, geom2_);
std::string out(wkt_c);
free(wkt_c);
ensure_equals(out,
"MULTILINESTRING ((4 0, 5 0), (3 0, 4 0), (2 0, 3 0), (1 0, 2 0), (0 0, 1 0))"
);
}
LWGEOM *
lwgeom_normalize(const LWGEOM *geom1)
{
LWGEOM *result ;
GEOSGeometry *g1;
int is3d ;
int srid ;
srid = (int)(geom1->srid);
is3d = FLAGS_GET_Z(geom1->flags);
initGEOS(lwnotice, lwgeom_geos_error);
g1 = LWGEOM2GEOS(geom1);
if ( 0 == g1 ) /* exception thrown at construction */
{
lwerror("First argument geometry could not be converted to GEOS.");
return NULL ;
}
if ( -1 == GEOSNormalize(g1) )
{
lwerror("Error in GEOSNormalize: %s", lwgeom_geos_errmsg);
return NULL; /* never get here */
}
GEOSSetSRID(g1, srid); /* needed ? */
result = GEOS2LWGEOM(g1, is3d);
GEOSGeom_destroy(g1);
if (result == NULL)
{
lwerror("Error performing intersection: GEOS2LWGEOM: %s",
lwgeom_geos_errmsg);
return NULL ; /* never get here */
}
return result ;
}
