*/
LWGEOM* lwgeom_split(const LWGEOM* lwgeom_in, const LWGEOM* blade_in);
+/*
+ * Return the set of paths shared between two linear geometries,
+ * and their direction (same or opposite).
+ *
+ * Developed by Sandro Santilli <strk@keybit.net> for Faunalia
+ * (http://www.faunalia.it) with funding from Regione Toscana - Sistema
+ * Informativo per la Gestione del Territorio e dell' Ambiente
+ * [RT-SIGTA]". For the project: "Sviluppo strumenti software per il
+ * trattamento di dati geografici basati su QuantumGIS e Postgis (CIG
+ * 0494241492)"
+ */
+LWGEOM* lwgeom_sharedpaths(const LWGEOM* geom1, const LWGEOM* geom2);
+
#endif /* !defined _LIBLWGEOM_H */
#endif /* POSTGIS_GEOS_VERSION >= 33 */
}
+
+LWGEOM*
+lwgeom_sharedpaths(const LWGEOM* geom1, const LWGEOM* geom2)
+{
+#if POSTGIS_GEOS_VERSION < 33
+ lwerror("The GEOS version this postgis binary "
+ "was compiled against (%d) doesn't support "
+ "'SharedPaths' function (3.3.0+ required)",
+ POSTGIS_GEOS_VERSION);
+ return NULL;
+#else /* POSTGIS_GEOS_VERSION >= 33 */
+ GEOSGeometry *g1, *g2, *g3;
+ LWGEOM *out;
+ int is3d, srid;
+
+ srid = geom1->srid;
+ error_if_srid_mismatch(srid, (int)(geom2->srid));
+
+ is3d = (FLAGS_GET_Z(geom1->flags) || FLAGS_GET_Z(geom2->flags)) ;
+
+ initGEOS(lwnotice, lwgeom_geos_error);
+
+ g1 = (GEOSGeometry *)LWGEOM2GEOS(geom1);
+ if ( 0 == g1 ) /* exception thrown at construction */
+ {
+ lwerror("First argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
+ return NULL;
+ }
+
+ g2 = (GEOSGeometry *)LWGEOM2GEOS(geom2);
+ if ( 0 == g2 ) /* exception thrown at construction */
+ {
+ lwerror("Second argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
+ GEOSGeom_destroy(g1);
+ return NULL;
+ }
+
+ g3 = GEOSSharedPaths(g1,g2);
+
+ if (g3 == NULL)
+ {
+ GEOSGeom_destroy(g1);
+ GEOSGeom_destroy(g2);
+ lwerror("GEOSSharedPaths: %s", lwgeom_geos_errmsg);
+ return NULL;
+ }
+
+ GEOSGeom_destroy(g1);
+ GEOSGeom_destroy(g2);
+
+ GEOSSetSRID(g3, srid);
+ out = GEOS2LWGEOM(g3, is3d);
+
+ if (out == NULL)
+ {
+ GEOSGeom_destroy(g3);
+ lwerror("GEOS2LWGEOM threw an error");
+ return NULL;
+ }
+ GEOSGeom_destroy(g3);
+
+ return out;
+#endif /* POSTGIS_GEOS_VERSION >= 33 */
+}
lwgeom_geos.o \
lwgeom_geos_prepared.o \
lwgeom_geos_clean.o \
- lwgeom_geos_sharedpaths.o \
lwgeom_geos_relatematch.o \
lwgeom_export.o \
lwgeom_in_gml.o \
}
+/**********************************************************************
+ *
+ * ST_SharedPaths
+ *
+ * Return the set of paths shared between two linear geometries,
+ * and their direction (same or opposite).
+ *
+ * Developed by Sandro Santilli (strk@keybit.net) for Faunalia
+ * (http://www.faunalia.it) with funding from Regione Toscana - Sistema
+ * Informativo per la Gestione del Territorio e dell' Ambiente
+ * [RT-SIGTA]". For the project: "Sviluppo strumenti software per il
+ * trattamento di dati geografici basati su QuantumGIS e Postgis (CIG
+ * 0494241492)"
+ *
+ **********************************************************************/
+Datum ST_SharedPaths(PG_FUNCTION_ARGS);
+PG_FUNCTION_INFO_V1(ST_SharedPaths);
+Datum ST_SharedPaths(PG_FUNCTION_ARGS)
+{
+#if POSTGIS_GEOS_VERSION < 33
+ lwerror("The GEOS version this postgis binary "
+ "was compiled against (%d) doesn't support "
+ "'ST_SharedPaths' function (3.3.0+ required)",
+ POSTGIS_GEOS_VERSION);
+ PG_RETURN_NULL();
+#else /* POSTGIS_GEOS_VERSION >= 33 */
+ PG_LWGEOM *geom1, *geom2, *out;
+ LWGEOM *g1, *g2, *lwgeom_out;
+
+ geom1 = (PG_LWGEOM *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
+ geom2 = (PG_LWGEOM *) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
+
+ g1 = pglwgeom_deserialize(geom1);
+ g2 = pglwgeom_deserialize(geom2);
+
+ lwgeom_out = lwgeom_sharedpaths(g1, g2);
+ if ( ! lwgeom_out )
+ {
+ PG_FREE_IF_COPY(geom1, 0);
+ PG_FREE_IF_COPY(geom2, 1);
+ PG_RETURN_NULL();
+ }
+
+ out = pglwgeom_serialize(lwgeom_out);
+ PG_FREE_IF_COPY(geom1, 0);
+ PG_FREE_IF_COPY(geom2, 1);
+
+ PG_RETURN_POINTER(out);
+
+#endif /* POSTGIS_GEOS_VERSION >= 33 */
+
+}
+
+++ /dev/null
-/**********************************************************************
- * $Id: lwgeom_geos.c 5258 2010-02-17 21:02:49Z strk $
- *
- * PostGIS - Spatial Types for PostgreSQL
- * http://postgis.refractions.net
- *
- * Copyright (C) 2010 Sandro Santilli <strk@keybit.net>
- *
- * This is free software; you can redistribute and/or modify it under
- * the terms of the GNU General Public Licence. See the COPYING file.
- *
- **********************************************************************
- *
- * ST_SharedPaths
- *
- * Return the set of paths shared between two linear geometries,
- * and their direction (same or opposite).
- *
- * Developed by Sandro Santilli (strk@keybit.net) for Faunalia
- * (http://www.faunalia.it) with funding from Regione Toscana - Sistema
- * Informativo per la Gestione del Territorio e dell' Ambiente
- * [RT-SIGTA]". For the project: "Sviluppo strumenti software per il
- * trattamento di dati geografici basati su QuantumGIS e Postgis (CIG
- * 0494241492)"
- *
- **********************************************************************/
-
-#include "lwgeom_geos.h"
-#include "liblwgeom_internal.h"
-
-#include <string.h>
-#include <assert.h>
-
-/* #define POSTGIS_DEBUG_LEVEL 4 */
-
-Datum ST_SharedPaths(PG_FUNCTION_ARGS);
-PG_FUNCTION_INFO_V1(ST_SharedPaths);
-Datum ST_SharedPaths(PG_FUNCTION_ARGS)
-{
-#if POSTGIS_GEOS_VERSION < 33
- lwerror("The GEOS version this postgis binary "
- "was compiled against (%d) doesn't support "
- "'ST_SharedPaths' function (3.3.0+ required)",
- POSTGIS_GEOS_VERSION);
- PG_RETURN_NULL();
-#else /* POSTGIS_GEOS_VERSION >= 33 */
- PG_LWGEOM *geom1, *geom2, *out;
- GEOSGeometry *g1, *g2, *g3;
- int is3d;
- int srid;
-
- geom1 = (PG_LWGEOM *) PG_DETOAST_DATUM(PG_GETARG_DATUM(0));
- geom2 = (PG_LWGEOM *) PG_DETOAST_DATUM(PG_GETARG_DATUM(1));
-
- is3d = ( pglwgeom_has_z(geom1) ) ||
- ( pglwgeom_has_z(geom2) );
-
- srid = pglwgeom_get_srid(geom1);
- error_if_srid_mismatch(srid, pglwgeom_get_srid(geom2));
-
- initGEOS(lwnotice, lwgeom_geos_error);
-
- g1 = (GEOSGeometry *)POSTGIS2GEOS(geom1);
- if ( 0 == g1 ) /* exception thrown at construction */
- {
- lwerror("First argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
- PG_FREE_IF_COPY(geom1, 0);
- PG_RETURN_NULL();
- }
- PG_FREE_IF_COPY(geom1, 0);
-
- g2 = (GEOSGeometry *)POSTGIS2GEOS(geom2);
- if ( 0 == g2 ) /* exception thrown at construction */
- {
- lwerror("Second argument geometry could not be converted to GEOS: %s", lwgeom_geos_errmsg);
- GEOSGeom_destroy(g1);
- PG_FREE_IF_COPY(geom2, 1);
- PG_RETURN_NULL();
- }
- PG_FREE_IF_COPY(geom2, 1);
-
- g3 = GEOSSharedPaths(g1,g2);
- if (g3 == NULL)
- {
- GEOSGeom_destroy(g1);
- GEOSGeom_destroy(g2);
- lwerror("GEOSIntersection: %s", lwgeom_geos_errmsg);
- PG_RETURN_NULL(); /* never get here */
- }
-
- GEOSSetSRID(g3, srid);
- out = GEOS2POSTGIS(g3, is3d);
- if (out == NULL)
- {
- GEOSGeom_destroy(g1);
- GEOSGeom_destroy(g2);
- GEOSGeom_destroy(g3);
- lwerror("GEOSSharedPaths() threw an error (result postgis geometry formation)!");
- PG_RETURN_NULL(); /* never get here */
- }
-
- GEOSGeom_destroy(g1);
- GEOSGeom_destroy(g2);
- GEOSGeom_destroy(g3);
-
- PG_RETURN_POINTER(out);
-
-#endif /* POSTGIS_GEOS_VERSION >= 33 */
-
-}
-
t3|GEOMETRYCOLLECTION(MULTILINESTRING EMPTY,MULTILINESTRING((0 0,10 0)))
t4|GEOMETRYCOLLECTION(MULTILINESTRING EMPTY,MULTILINESTRING EMPTY)
t5|GEOMETRYCOLLECTION(MULTILINESTRING((20 0,30 0)),MULTILINESTRING((70 0,80 0)))
-ERROR: GEOSIntersection: IllegalArgumentException: Geometry is not lineal
+ERROR: GEOSSharedPaths: IllegalArgumentException: Geometry is not lineal
t7|GEOMETRYCOLLECTION(MULTILINESTRING((4 2,7 2)),MULTILINESTRING((14 7,11 6)))
t8|GEOMETRYCOLLECTION(MULTILINESTRING((4 2,7 2)),MULTILINESTRING((14 7,11 6)))
projects / postgis / commitdiff