From: Sandro Santilli <strk@keybit.net>
Date: Mon, 26 Sep 2005 12:30:27 +0000 (+0000)
Subject: Added LineMerge interface
X-Git-Tag: pgis_1_1_0~241
X-Git-Url: https://granicus.if.org/sourcecode?a=commitdiff_plain;h=838a595c497b81b781c2bb7a2a4c0dd7b0cbed40;p=postgis

Added LineMerge interface


git-svn-id: http://svn.osgeo.org/postgis/trunk@1936 b70326c6-7e19-0410-871a-916f4a2858ee
---

diff --git a/lwgeom/lwgeom_geos_c.c b/lwgeom/lwgeom_geos_c.c
index a231b9606..a4226a60c 100644
--- a/lwgeom/lwgeom_geos_c.c
+++ b/lwgeom/lwgeom_geos_c.c
@@ -49,6 +49,7 @@ Datum postgis_geos_version(PG_FUNCTION_ARGS);
 Datum postgis_jts_version(PG_FUNCTION_ARGS);
 Datum centroid(PG_FUNCTION_ARGS);
 Datum polygonize_garray(PG_FUNCTION_ARGS);
+Datum linemerge(PG_FUNCTION_ARGS);
 
 
 
@@ -2335,6 +2336,81 @@ Datum polygonize_garray(PG_FUNCTION_ARGS)
 
 }
 
+PG_FUNCTION_INFO_V1(linemerge);
+Datum linemerge(PG_FUNCTION_ARGS)
+{
+	PG_LWGEOM	*geom1;
+	double	size;
+	GEOSGeom g1,g3;
+	PG_LWGEOM *result;
+	int quadsegs = 8; // the default
+
+#ifdef PROFILE
+	profstart(PROF_QRUN);
+#endif
+
+	geom1 = (PG_LWGEOM *)  PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
+
+	initGEOS(lwnotice, lwerror);
+
+#ifdef PROFILE
+	profstart(PROF_P2G1);
+#endif
+	g1 = POSTGIS2GEOS(geom1);
+#ifdef PROFILE
+	profstop(PROF_P2G1);
+#endif
+
+#ifdef PROFILE
+	profstart(PROF_GRUN);
+#endif
+	g3 = GEOSLineMerge(g1);
+#ifdef PROFILE
+	profstop(PROF_GRUN);
+#endif
+
+	if (g3 == NULL)
+	{
+		elog(ERROR,"GEOS LineMerge() threw an error!");
+		GEOSGeom_destroy(g1);
+		PG_RETURN_NULL(); //never get here
+	}
+
+
+//	elog(NOTICE,"result: %s", GEOSasText(g3) ) ;
+
+	GEOSSetSRID(g3, pglwgeom_getSRID(geom1));
+
+#ifdef PROFILE
+	profstart(PROF_G2P);
+#endif
+	result = GEOS2POSTGIS(g3, TYPE_NDIMS(geom1->type) > 2);
+#ifdef PROFILE
+	profstop(PROF_G2P);
+#endif
+	if (result == NULL)
+	{
+		GEOSGeom_destroy(g1);
+		GEOSGeom_destroy(g3);
+		elog(ERROR,"GEOS LineMerge() threw an error (result postgis geometry formation)!");
+		PG_RETURN_NULL(); //never get here
+	}
+	GEOSGeom_destroy(g1);
+	GEOSGeom_destroy(g3);
+
+
+	//compressType(result);  // convert multi* to single item if appropriate
+
+#ifdef PROFILE
+	profstop(PROF_QRUN);
+	profreport("geos",geom1, NULL, result);
+#endif
+
+	PG_FREE_IF_COPY(geom1, 0);
+
+	PG_RETURN_POINTER(result);
+}
+
 Datum JTSnoop(PG_FUNCTION_ARGS);
 PG_FUNCTION_INFO_V1(JTSnoop);
 Datum JTSnoop(PG_FUNCTION_ARGS)