geo_coder 0.1.0

data/src/liblwgeom/lwmcurve.c

@@ -0,0 +1,125 @@
+/**********************************************************************
+ * $Id: lwmcurve.c 3639 2009-02-04 00:28:37Z pramsey $
+ *
+ * PostGIS - Spatial Types for PostgreSQL
+ * http://postgis.refractions.net
+ * Copyright 2001-2006 Refractions Research Inc.
+ *
+ * This is free software; you can redistribute and/or modify it under
+ * the terms of the GNU General Public Licence. See the COPYING file.
+ * 
+ **********************************************************************/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include "liblwgeom.h"
+
+LWMCURVE *
+lwmcurve_deserialize(uchar *srl)
+{
+        LWMCURVE *result;
+        LWGEOM_INSPECTED *insp;
+        int stype;
+        int type = lwgeom_getType(srl[0]);
+        int i;
+
+        if(type != MULTICURVETYPE)
+        {
+                lwerror("lwmcurve_deserialize called on NON multicurve: %d", type);
+                return NULL;
+        }
+
+        insp = lwgeom_inspect(srl);
+
+        result = lwalloc(sizeof(LWMCURVE));
+        result->type = insp->type;
+        result->SRID = insp->SRID;
+        result->ngeoms = insp->ngeometries;
+        result->geoms = lwalloc(sizeof(LWGEOM *)*insp->ngeometries);
+
+        if(lwgeom_hasBBOX(srl[0]))
+        {
+                result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
+                memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
+        }
+        else result->bbox = NULL;
+
+        for(i = 0; i < insp->ngeometries; i++)
+        {
+                stype = lwgeom_getType(insp->sub_geoms[i][0]);
+                if(stype == CIRCSTRINGTYPE)
+                {
+                        result->geoms[i] = (LWGEOM *)lwcircstring_deserialize(insp->sub_geoms[i]);
+                }
+                else if(stype == LINETYPE)
+                {
+                        result->geoms[i] = (LWGEOM *)lwline_deserialize(insp->sub_geoms[i]);
+                }
+                else
+                {
+                        lwerror("Only Circular and Line strings are currenly permitted in a MultiCurve.");
+                        lwfree(result);
+                        lwfree(insp);
+                        return NULL;
+                }
+                        
+                if(TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type))
+                {
+                        lwerror("Mixed dimensions (multicurve: %d, curve %d:%d)",
+                                TYPE_NDIMS(result->type), i,
+                                TYPE_NDIMS(result->geoms[i]->type));
+                        lwfree(result);
+                        lwfree(insp);
+                        return NULL;
+                }
+        }
+        return result;
+}
+
+/*
+ * Add 'what' to this multicurve at position 'where'.
+ * where=0 == prepend
+ * where=-1 == append
+ * Returns a MULTICURVE or a COLLECTION
+ */
+LWGEOM *
+lwmcurve_add(const LWMCURVE *to, uint32 where, const LWGEOM *what)
+{
+        LWCOLLECTION *col;
+        LWGEOM **geoms;
+        int newtype;
+        uint32 i;
+
+        if(where == -1) where = to->ngeoms;
+        else if(where < -1 || where > to->ngeoms)
+        {
+                lwerror("lwmcurve_add: add position out of range %d..%d",
+                        -1, to->ngeoms);
+                return NULL;
+        }
+
+        /* dimensions compatibility are checked by caller */
+
+        /* Construct geoms array */
+        geoms = lwalloc(sizeof(LWGEOM *)*(to->ngeoms+1));
+        for(i = 0; i < where; i++)
+        {
+                geoms[i] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+        }
+        geoms[where] = lwgeom_clone(what);
+        for(i = where; i < to->ngeoms; i++) 
+        {
+                geoms[i+1] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+        }
+
+        if(TYPE_GETTYPE(what->type) == CIRCSTRINGTYPE) newtype = MULTICURVETYPE;
+        else newtype = COLLECTIONTYPE;
+
+        col = lwcollection_construct(newtype,
+                to->SRID, NULL,
+                to->ngeoms + 1, geoms);
+
+        return (LWGEOM *)col;
+}
+

data/src/liblwgeom/lwmline.c

@@ -0,0 +1,137 @@
+/**********************************************************************
+ * $Id: lwmline.c 3639 2009-02-04 00:28:37Z pramsey $
+ *
+ * PostGIS - Spatial Types for PostgreSQL
+ * http://postgis.refractions.net
+ * Copyright 2001-2006 Refractions Research Inc.
+ *
+ * This is free software; you can redistribute and/or modify it under
+ * the terms of the GNU General Public Licence. See the COPYING file.
+ * 
+ **********************************************************************/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include "liblwgeom.h"
+
+void
+lwmline_release(LWMLINE *lwmline)
+{
+  lwgeom_release(lwmline_as_lwgeom(lwmline));
+}
+
+LWMLINE *
+lwmline_deserialize(uchar *srl)
+{
+	LWMLINE *result;
+	LWGEOM_INSPECTED *insp;
+	int type = lwgeom_getType(srl[0]);
+	int i;
+
+	if ( type != MULTILINETYPE ) 
+	{
+		lwerror("lwmline_deserialize called on NON multiline: %d",
+			type);
+		return NULL;
+	}
+
+	insp = lwgeom_inspect(srl);
+
+	result = lwalloc(sizeof(LWMLINE));
+	result->type = insp->type;
+	result->SRID = insp->SRID;
+	result->ngeoms = insp->ngeometries;
+	result->geoms = lwalloc(sizeof(LWLINE *)*insp->ngeometries);
+
+	if (lwgeom_hasBBOX(srl[0]))
+	{
+		result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
+		memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
+	}
+	else result->bbox = NULL;
+
+
+	for (i=0; i<insp->ngeometries; i++)
+	{
+		result->geoms[i] = lwline_deserialize(insp->sub_geoms[i]);
+		if ( TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type) )
+		{
+			lwerror("Mixed dimensions (multiline:%d, line%d:%d)",
+				TYPE_NDIMS(result->type), i,
+				TYPE_NDIMS(result->geoms[i]->type)
+			);
+			return NULL;
+		}
+	}
+
+	return result;
+}
+
+/*
+ * Add 'what' to this multiline at position 'where'.
+ * where=0 == prepend
+ * where=-1 == append
+ * Returns a MULTILINE or a COLLECTION
+ */
+LWGEOM *
+lwmline_add(const LWMLINE *to, uint32 where, const LWGEOM *what)
+{
+	LWCOLLECTION *col;
+	LWGEOM **geoms;
+	int newtype;
+	uint32 i;
+
+	if ( where == -1 ) where = to->ngeoms;
+	else if ( where < -1 || where > to->ngeoms )
+	{
+		lwerror("lwmline_add: add position out of range %d..%d",
+			-1, to->ngeoms);
+		return NULL;
+	}
+
+	/* dimensions compatibility are checked by caller */
+
+	/* Construct geoms array */
+	geoms = lwalloc(sizeof(LWGEOM *)*(to->ngeoms+1));
+	for (i=0; i<where; i++)
+	{
+		geoms[i] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+	geoms[where] = lwgeom_clone(what);
+	for (i=where; i<to->ngeoms; i++)
+	{
+		geoms[i+1] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+
+	if ( TYPE_GETTYPE(what->type) == LINETYPE ) newtype = MULTILINETYPE;
+	else newtype = COLLECTIONTYPE;
+
+	col = lwcollection_construct(newtype,
+		to->SRID, NULL,
+		to->ngeoms+1, geoms);
+	
+	return (LWGEOM *)col;
+
+}
+
+void lwmline_free(LWMLINE *mline) 
+{
+	int i;
+	if( mline->bbox ) 
+	{
+		lwfree(mline->bbox);
+	}
+	for ( i = 0; i < mline->ngeoms; i++ ) 
+	{
+		if( mline->geoms[i] ) {
+			lwline_free(mline->geoms[i]);
+		}
+	}
+	if( mline->geoms ) 
+	{
+		lwfree(mline->geoms);
+	}
+	lwfree(mline);
+	
+};

data/src/liblwgeom/lwmpoint.c

@@ -0,0 +1,138 @@
+/**********************************************************************
+ * $Id: lwmpoint.c 3639 2009-02-04 00:28:37Z pramsey $
+ *
+ * PostGIS - Spatial Types for PostgreSQL
+ * http://postgis.refractions.net
+ * Copyright 2001-2006 Refractions Research Inc.
+ *
+ * This is free software; you can redistribute and/or modify it under
+ * the terms of the GNU General Public Licence. See the COPYING file.
+ * 
+ **********************************************************************/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include "liblwgeom.h"
+
+void
+lwmpoint_release(LWMPOINT *lwmpoint)
+{
+  lwgeom_release(lwmpoint_as_lwgeom(lwmpoint));
+}
+
+
+LWMPOINT *
+lwmpoint_deserialize(uchar *srl)
+{
+	LWMPOINT *result;
+	LWGEOM_INSPECTED *insp;
+	int type = lwgeom_getType(srl[0]);
+	int i;
+
+	if ( type != MULTIPOINTTYPE ) 
+	{
+		lwerror("lwmpoint_deserialize called on NON multipoint: %d",
+			type);
+		return NULL;
+	}
+
+	insp = lwgeom_inspect(srl);
+
+	result = lwalloc(sizeof(LWMPOINT));
+	result->type = insp->type;
+	result->SRID = insp->SRID;
+	result->ngeoms = insp->ngeometries;
+	result->geoms = lwalloc(sizeof(LWPOINT *)*result->ngeoms);
+
+	if (lwgeom_hasBBOX(srl[0]))
+	{
+		result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
+		memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
+	}
+	else result->bbox = NULL;
+
+	for (i=0; i<insp->ngeometries; i++)
+	{
+		result->geoms[i] = lwpoint_deserialize(insp->sub_geoms[i]);
+		if ( TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type) )
+		{
+			lwerror("Mixed dimensions (multipoint:%d, point%d:%d)",
+				TYPE_NDIMS(result->type), i,
+				TYPE_NDIMS(result->geoms[i]->type)
+			);
+			return NULL;
+		}
+	}
+
+	return result;
+}
+
+/*
+ * Add 'what' to this multipoint at position 'where'.
+ * where=0 == prepend
+ * where=-1 == append
+ * Returns a MULTIPOINT or a COLLECTION
+ */
+LWGEOM *
+lwmpoint_add(const LWMPOINT *to, uint32 where, const LWGEOM *what)
+{
+	LWCOLLECTION *col;
+	LWGEOM **geoms;
+	int newtype;
+	uint32 i;
+
+	if ( where == -1 ) where = to->ngeoms;
+	else if ( where < -1 || where > to->ngeoms )
+	{
+		lwerror("lwmpoint_add: add position out of range %d..%d",
+			-1, to->ngeoms);
+		return NULL;
+	}
+
+	/* dimensions compatibility are checked by caller */
+
+	/* Construct geoms array */
+	geoms = lwalloc(sizeof(LWGEOM *)*(to->ngeoms+1));
+	for (i=0; i<where; i++)
+	{
+		geoms[i] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+	geoms[where] = lwgeom_clone(what);
+	for (i=where; i<to->ngeoms; i++)
+	{
+		geoms[i+1] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+
+	if ( TYPE_GETTYPE(what->type) == POINTTYPE ) newtype = MULTIPOINTTYPE;
+	else newtype = COLLECTIONTYPE;
+
+	col = lwcollection_construct(newtype,
+		to->SRID, NULL,
+		to->ngeoms+1, geoms);
+	
+	return (LWGEOM *)col;
+
+}
+
+void lwmpoint_free(LWMPOINT *mpt) 
+{
+	int i;
+	if( mpt->bbox ) 
+	{
+		lwfree(mpt->bbox);
+	}
+	for ( i = 0; i < mpt->ngeoms; i++ ) 
+	{
+		if( mpt->geoms[i] ) {
+			lwpoint_free(mpt->geoms[i]);
+		}
+	}
+	if( mpt->geoms ) 
+	{
+		lwfree(mpt->geoms);
+	}
+	lwfree(mpt);
+	
+};
+

data/src/liblwgeom/lwmpoly.c

@@ -0,0 +1,141 @@
+/**********************************************************************
+ * $Id: lwmpoly.c 3639 2009-02-04 00:28:37Z pramsey $
+ *
+ * PostGIS - Spatial Types for PostgreSQL
+ * http://postgis.refractions.net
+ * Copyright 2001-2006 Refractions Research Inc.
+ *
+ * This is free software; you can redistribute and/or modify it under
+ * the terms of the GNU General Public Licence. See the COPYING file.
+ * 
+ **********************************************************************/
+
+#include <stdio.h>
+#include <stdlib.h>
+#include <string.h>
+#include "liblwgeom.h"
+
+
+void
+lwmpoly_release(LWMPOLY *lwmpoly)
+{
+  lwgeom_release(lwmpoly_as_lwgeom(lwmpoly));
+}
+
+
+LWMPOLY *
+lwmpoly_deserialize(uchar *srl)
+{
+	LWMPOLY *result;
+	LWGEOM_INSPECTED *insp;
+	int type = lwgeom_getType(srl[0]);
+	int i;
+
+	LWDEBUG(2, "lwmpoly_deserialize called");
+
+	if ( type != MULTIPOLYGONTYPE ) 
+	{
+		lwerror("lwmpoly_deserialize called on NON multipoly: %d",
+			type);
+		return NULL;
+	}
+
+	insp = lwgeom_inspect(srl);
+
+	result = lwalloc(sizeof(LWMPOLY));
+	result->type = insp->type;
+	result->SRID = insp->SRID;
+	result->ngeoms = insp->ngeometries;
+	result->geoms = lwalloc(sizeof(LWPOLY *)*insp->ngeometries);
+
+	if (lwgeom_hasBBOX(srl[0]))
+	{
+		result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
+		memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
+	}
+	else result->bbox = NULL;
+
+	for (i=0; i<insp->ngeometries; i++)
+	{
+		result->geoms[i] = lwpoly_deserialize(insp->sub_geoms[i]);
+		if ( TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type) )
+		{
+			lwerror("Mixed dimensions (multipoly:%d, poly%d:%d)",
+				TYPE_NDIMS(result->type), i,
+				TYPE_NDIMS(result->geoms[i]->type)
+			);
+			return NULL;
+		}
+	}
+
+	return result;
+}
+
+/*
+ * Add 'what' to this multiline at position 'where'.
+ * where=0 == prepend
+ * where=-1 == append
+ * Returns a MULTIPOLY or a COLLECTION
+ */
+LWGEOM *
+lwmpoly_add(const LWMPOLY *to, uint32 where, const LWGEOM *what)
+{
+	LWCOLLECTION *col;
+	LWGEOM **geoms;
+	int newtype;
+	uint32 i;
+
+	if ( where == -1 ) where = to->ngeoms;
+	else if ( where < -1 || where > to->ngeoms )
+	{
+		lwerror("lwmline_add: add position out of range %d..%d",
+			-1, to->ngeoms);
+		return NULL;
+	}
+
+	/* dimensions compatibility are checked by caller */
+
+	/* Construct geoms array */
+	geoms = lwalloc(sizeof(LWGEOM *)*(to->ngeoms+1));
+	for (i=0; i<where; i++)
+	{
+		geoms[i] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+	geoms[where] = lwgeom_clone(what);
+	for (i=where; i<to->ngeoms; i++)
+	{
+		geoms[i+1] = lwgeom_clone((LWGEOM *)to->geoms[i]);
+	}
+
+	if ( TYPE_GETTYPE(what->type) == POLYGONTYPE ) newtype = MULTIPOLYGONTYPE;
+	else newtype = COLLECTIONTYPE;
+
+	col = lwcollection_construct(newtype,
+		to->SRID, NULL,
+		to->ngeoms+1, geoms);
+	
+	return (LWGEOM *)col;
+
+}
+
+void lwmpoly_free(LWMPOLY *mpoly) 
+{
+	int i;
+	if( mpoly->bbox ) 
+	{
+		lwfree(mpoly->bbox);
+	}
+	for ( i = 0; i < mpoly->ngeoms; i++ ) 
+	{
+		if( mpoly->geoms[i] ) {
+			lwpoly_free(mpoly->geoms[i]);
+		}
+	}
+	if( mpoly->geoms ) 
+	{
+		lwfree(mpoly->geoms);
+	}
+	lwfree(mpoly);
+	
+};
+