mirror of
https://git.osgeo.org/gitea/postgis/postgis
synced 2024-10-26 01:52:37 +00:00
abf5606427
git-svn-id: http://svn.osgeo.org/postgis/trunk@4057 b70326c6-7e19-0410-871a-916f4a2858ee
545 lines
12 KiB
C
545 lines
12 KiB
C
/**********************************************************************
|
|
* $Id$
|
|
*
|
|
* 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"
|
|
|
|
|
|
#define CHECK_LWGEOM_ZM 1
|
|
|
|
void
|
|
lwcollection_release(LWCOLLECTION *lwcollection)
|
|
{
|
|
lwgeom_release(lwcollection_as_lwgeom(lwcollection));
|
|
}
|
|
|
|
|
|
LWCOLLECTION *
|
|
lwcollection_construct(unsigned int type, int SRID, BOX2DFLOAT4 *bbox,
|
|
unsigned int ngeoms, LWGEOM **geoms)
|
|
{
|
|
LWCOLLECTION *ret;
|
|
int hasz, hasm;
|
|
#ifdef CHECK_LWGEOM_ZM
|
|
char zm;
|
|
unsigned int i;
|
|
#endif
|
|
|
|
LWDEBUGF(2, "lwcollection_construct called with %d, %d, %p, %d, %p.", type, SRID, bbox, ngeoms, geoms);
|
|
|
|
hasz = 0;
|
|
hasm = 0;
|
|
if ( ngeoms > 0 )
|
|
{
|
|
hasz = TYPE_HASZ(geoms[0]->type);
|
|
hasm = TYPE_HASM(geoms[0]->type);
|
|
#ifdef CHECK_LWGEOM_ZM
|
|
zm = TYPE_GETZM(geoms[0]->type);
|
|
|
|
LWDEBUGF(3, "lwcollection_construct type[0]=%d", geoms[0]->type);
|
|
|
|
for (i=1; i<ngeoms; i++)
|
|
{
|
|
LWDEBUGF(3, "lwcollection_construct type=[%d]=%d", i, geoms[i]->type);
|
|
|
|
if ( zm != TYPE_GETZM(geoms[i]->type) )
|
|
lwerror("lwcollection_construct: mixed dimension geometries: %d/%d", zm, TYPE_GETZM(geoms[i]->type));
|
|
}
|
|
#endif
|
|
}
|
|
|
|
|
|
ret = lwalloc(sizeof(LWCOLLECTION));
|
|
ret->type = lwgeom_makeType_full(hasz, hasm, (SRID!=-1),
|
|
type, 0);
|
|
ret->SRID = SRID;
|
|
ret->ngeoms = ngeoms;
|
|
ret->geoms = geoms;
|
|
ret->bbox = bbox;
|
|
|
|
return ret;
|
|
}
|
|
|
|
LWCOLLECTION *
|
|
lwcollection_construct_empty(int SRID, char hasz, char hasm)
|
|
{
|
|
LWCOLLECTION *ret;
|
|
|
|
ret = lwalloc(sizeof(LWCOLLECTION));
|
|
ret->type = lwgeom_makeType_full(hasz, hasm, (SRID!=-1),
|
|
COLLECTIONTYPE, 0);
|
|
ret->SRID = SRID;
|
|
ret->ngeoms = 0;
|
|
ret->geoms = NULL;
|
|
ret->bbox = NULL;
|
|
|
|
return ret;
|
|
}
|
|
|
|
|
|
LWCOLLECTION *
|
|
lwcollection_deserialize(uchar *srl)
|
|
{
|
|
LWCOLLECTION *result;
|
|
LWGEOM_INSPECTED *insp;
|
|
char typefl = srl[0];
|
|
int type = lwgeom_getType(typefl);
|
|
int i;
|
|
|
|
if ( type != COLLECTIONTYPE )
|
|
{
|
|
lwerror("lwcollection_deserialize called on NON geometrycollection: %d", type);
|
|
return NULL;
|
|
}
|
|
|
|
insp = lwgeom_inspect(srl);
|
|
|
|
result = lwalloc(sizeof(LWCOLLECTION));
|
|
result->type = typefl;
|
|
result->SRID = insp->SRID;
|
|
result->ngeoms = insp->ngeometries;
|
|
|
|
if (lwgeom_hasBBOX(srl[0]))
|
|
{
|
|
result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
|
|
memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
|
|
}
|
|
else result->bbox = NULL;
|
|
|
|
|
|
if ( insp->ngeometries )
|
|
{
|
|
result->geoms = lwalloc(sizeof(LWGEOM *)*insp->ngeometries);
|
|
for (i=0; i<insp->ngeometries; i++)
|
|
{
|
|
result->geoms[i] =
|
|
lwgeom_deserialize(insp->sub_geoms[i]);
|
|
}
|
|
}
|
|
|
|
return result;
|
|
}
|
|
|
|
LWGEOM *
|
|
lwcollection_getsubgeom(LWCOLLECTION *col, int gnum)
|
|
{
|
|
return (LWGEOM *)col->geoms[gnum];
|
|
}
|
|
|
|
/**
|
|
* @brief find serialized size of this collection
|
|
* @param col #LWCOLLECTION to find serialized size of
|
|
*/
|
|
size_t
|
|
lwcollection_serialize_size(LWCOLLECTION *col)
|
|
{
|
|
size_t size = 5; /* type + nsubgeoms */
|
|
int i;
|
|
|
|
if ( col->SRID != -1 ) size += 4; /* SRID */
|
|
if ( col->bbox ) size += sizeof(BOX2DFLOAT4);
|
|
|
|
LWDEBUGF(2, "lwcollection_serialize_size[%p]: start size: %d", col, size);
|
|
|
|
|
|
for (i=0; i<col->ngeoms; i++)
|
|
{
|
|
size += lwgeom_serialize_size(col->geoms[i]);
|
|
|
|
LWDEBUGF(3, "lwcollection_serialize_size[%p]: with geom%d: %d", col, i, size);
|
|
}
|
|
|
|
LWDEBUGF(3, "lwcollection_serialize_size[%p]: returning %d", col, size);
|
|
|
|
return size;
|
|
}
|
|
|
|
/** @brief convert an #LWCOLLECTION into its serialized form writing it into
|
|
* the given buffer, and returning number of bytes written into
|
|
* the given int pointer.
|
|
*/
|
|
void
|
|
lwcollection_serialize_buf(LWCOLLECTION *coll, uchar *buf, size_t *retsize)
|
|
{
|
|
size_t size=1; /* type */
|
|
size_t subsize=0;
|
|
char hasSRID;
|
|
uchar *loc;
|
|
int i;
|
|
|
|
LWDEBUGF(2, "lwcollection_serialize_buf called (%s with %d elems)",
|
|
lwgeom_typename(TYPE_GETTYPE(coll->type)), coll->ngeoms);
|
|
|
|
hasSRID = (coll->SRID != -1);
|
|
|
|
buf[0] = lwgeom_makeType_full(
|
|
TYPE_HASZ(coll->type), TYPE_HASM(coll->type),
|
|
hasSRID, TYPE_GETTYPE(coll->type), coll->bbox ? 1 : 0);
|
|
loc = buf+1;
|
|
|
|
/* Add BBOX if requested */
|
|
if ( coll->bbox )
|
|
{
|
|
memcpy(loc, coll->bbox, sizeof(BOX2DFLOAT4));
|
|
size += sizeof(BOX2DFLOAT4);
|
|
loc += sizeof(BOX2DFLOAT4);
|
|
}
|
|
|
|
/* Add SRID if requested */
|
|
if (hasSRID)
|
|
{
|
|
memcpy(loc, &coll->SRID, 4);
|
|
size += 4;
|
|
loc += 4;
|
|
}
|
|
|
|
/* Write number of subgeoms */
|
|
memcpy(loc, &coll->ngeoms, 4);
|
|
size += 4;
|
|
loc += 4;
|
|
|
|
/* Serialize subgeoms */
|
|
for (i=0; i<coll->ngeoms; i++)
|
|
{
|
|
lwgeom_serialize_buf(coll->geoms[i], loc, &subsize);
|
|
size += subsize;
|
|
loc += subsize;
|
|
}
|
|
|
|
if (retsize) *retsize = size;
|
|
|
|
LWDEBUG(3, "lwcollection_serialize_buf returning");
|
|
}
|
|
|
|
int
|
|
lwcollection_compute_box2d_p(LWCOLLECTION *col, BOX2DFLOAT4 *box)
|
|
{
|
|
BOX2DFLOAT4 boxbuf;
|
|
uint32 i;
|
|
|
|
if ( ! col->ngeoms ) return 0;
|
|
if ( ! lwgeom_compute_box2d_p(col->geoms[0], box) ) return 0;
|
|
for (i=1; i<col->ngeoms; i++)
|
|
{
|
|
if ( ! lwgeom_compute_box2d_p(col->geoms[i], &boxbuf) )
|
|
return 0;
|
|
if ( ! box2d_union_p(box, &boxbuf, box) ) return 0;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
/**
|
|
* @brief Clone #LWCOLLECTION object. #POINTARRAY are not copied.
|
|
* Bbox is cloned if present in input.
|
|
*/
|
|
LWCOLLECTION *
|
|
lwcollection_clone(const LWCOLLECTION *g)
|
|
{
|
|
uint32 i;
|
|
LWCOLLECTION *ret = lwalloc(sizeof(LWCOLLECTION));
|
|
memcpy(ret, g, sizeof(LWCOLLECTION));
|
|
if ( g->ngeoms > 0 )
|
|
{
|
|
ret->geoms = lwalloc(sizeof(LWGEOM *)*g->ngeoms);
|
|
for (i=0; i<g->ngeoms; i++)
|
|
{
|
|
ret->geoms[i] = lwgeom_clone(g->geoms[i]);
|
|
}
|
|
if ( g->bbox ) ret->bbox = box2d_clone(g->bbox);
|
|
}
|
|
else
|
|
{
|
|
ret->bbox = NULL; /* empty collection */
|
|
ret->geoms = NULL;
|
|
}
|
|
return ret;
|
|
}
|
|
|
|
/** @brief Add 'what' to this to #LWCOLLECTION at position 'where'.
|
|
* where=0 then prepend
|
|
* where=-1 then append
|
|
* Returns a GEOMETRYCOLLECTION
|
|
*/
|
|
LWGEOM *
|
|
lwcollection_add(const LWCOLLECTION *to, uint32 where, const LWGEOM *what)
|
|
{
|
|
LWCOLLECTION *col;
|
|
LWGEOM **geoms;
|
|
uint32 i;
|
|
|
|
if ( where == -1 ) where = to->ngeoms;
|
|
else if ( where < -1 || where > to->ngeoms )
|
|
{
|
|
lwerror("lwcollection_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(to->geoms[i]);
|
|
lwgeom_dropSRID(geoms[i]);
|
|
lwgeom_drop_bbox(geoms[i]);
|
|
}
|
|
geoms[where] = lwgeom_clone(what);
|
|
lwgeom_dropSRID(geoms[where]);
|
|
lwgeom_drop_bbox(geoms[where]);
|
|
for (i=where; i<to->ngeoms; i++)
|
|
{
|
|
geoms[i+1] = lwgeom_clone(to->geoms[i]);
|
|
lwgeom_dropSRID(geoms[i+1]);
|
|
lwgeom_drop_bbox(geoms[i+1]);
|
|
}
|
|
|
|
col = lwcollection_construct(COLLECTIONTYPE,
|
|
to->SRID, NULL,
|
|
to->ngeoms+1, geoms);
|
|
|
|
return (LWGEOM *)col;
|
|
|
|
}
|
|
|
|
LWCOLLECTION *
|
|
lwcollection_segmentize2d(LWCOLLECTION *col, double dist)
|
|
{
|
|
unsigned int i;
|
|
LWGEOM **newgeoms;
|
|
|
|
if ( ! col->ngeoms ) return col;
|
|
|
|
newgeoms = lwalloc(sizeof(LWGEOM *)*col->ngeoms);
|
|
for (i=0; i<col->ngeoms; i++)
|
|
newgeoms[i] = lwgeom_segmentize2d(col->geoms[i], dist);
|
|
|
|
return lwcollection_construct(col->type, col->SRID, NULL,
|
|
col->ngeoms, newgeoms);
|
|
}
|
|
|
|
/** @brief check for same geometry composition
|
|
*
|
|
*/
|
|
char
|
|
lwcollection_same(const LWCOLLECTION *c1, const LWCOLLECTION *c2)
|
|
{
|
|
unsigned int i, j;
|
|
unsigned int *hit;
|
|
|
|
LWDEBUG(2, "lwcollection_same called");
|
|
|
|
if ( TYPE_GETTYPE(c1->type) != TYPE_GETTYPE(c2->type) ) return 0;
|
|
if ( c1->ngeoms != c2->ngeoms ) return 0;
|
|
|
|
hit = lwalloc(sizeof(unsigned int)*c1->ngeoms);
|
|
memset(hit, 0, sizeof(unsigned int)*c1->ngeoms);
|
|
|
|
for (i=0; i<c1->ngeoms; i++)
|
|
{
|
|
char found=0;
|
|
for (j=0; j<c2->ngeoms; j++)
|
|
{
|
|
if ( hit[j] ) continue;
|
|
if ( lwgeom_same(c1->geoms[i], c2->geoms[j]) )
|
|
{
|
|
hit[j] = 1;
|
|
found=1;
|
|
break;
|
|
}
|
|
}
|
|
if ( ! found ) return 0;
|
|
}
|
|
return 1;
|
|
}
|
|
|
|
int lwcollection_ngeoms(const LWCOLLECTION *col)
|
|
{
|
|
int i;
|
|
int ngeoms = 0;
|
|
|
|
if( ! col ) {
|
|
lwerror("Null input geometry.");
|
|
return 0;
|
|
}
|
|
|
|
for( i = 0; i < col->ngeoms; i++ )
|
|
{
|
|
if( col->geoms[i]) {
|
|
switch(TYPE_GETTYPE(col->geoms[i]->type)) {
|
|
case POINTTYPE:
|
|
case LINETYPE:
|
|
case CIRCSTRINGTYPE:
|
|
case POLYGONTYPE:
|
|
ngeoms += 1;
|
|
break;
|
|
case MULTIPOINTTYPE:
|
|
case MULTILINETYPE:
|
|
case MULTICURVETYPE:
|
|
case MULTIPOLYGONTYPE:
|
|
ngeoms += col->ngeoms;
|
|
break;
|
|
case COLLECTIONTYPE:
|
|
ngeoms += lwcollection_ngeoms((LWCOLLECTION*)col->geoms[i]);
|
|
break;
|
|
}
|
|
}
|
|
}
|
|
return ngeoms;
|
|
}
|
|
|
|
/*
|
|
** Given a generic collection, return the "simplest" form.
|
|
** eg: GEOMETRYCOLLECTION(MULTILINESTRING()) => MULTELINESTRING()
|
|
** GEOMETRYCOLLECTION(MULTILINESTRING(), MULTILINESTRING(), POINT()) => GEOMETRYCOLLECTION(MULTILINESTRING(), MULTIPOINT())
|
|
**
|
|
** In general, if the subcomponents are homogeneous, return a properly typed collection.
|
|
** Otherwise, return a generic collection, with the subtypes in minimal typed collections.
|
|
LWCOLLECTION *lwcollection_homogenize(const LWCOLLECTION *c1)
|
|
{
|
|
TODO: pramsey
|
|
}
|
|
*/
|
|
|
|
/*
|
|
** Given a generic collection, extract and return just the desired types.
|
|
LWGEOM *lwcollection_extract(const LWCOLLECTION *col, char type)
|
|
{
|
|
LWGEOM **extracted_geoms;
|
|
extracted_geoms = lwalloc(sizeof(void*)*col->ngeoms);
|
|
extracted_curgeom = 0;
|
|
char reqtype = TYPE_GETTYPE(type);
|
|
for ( i = 0; i < col->ngeoms; i++ )
|
|
{
|
|
if( col->geoms[i] )
|
|
char geomtype = TYPE_GETTYPE(col->geoms[i]->type);
|
|
if ( geomtype == reqtype ) {
|
|
extracted_geoms[extracted_curgeom] = col->geoms[i];
|
|
extracted_curgeom++;
|
|
continue;
|
|
}
|
|
else {
|
|
if ( geomtype == COLLECTIONTYPE ) {
|
|
LWGEOM *colgeom;
|
|
colgeom = lwcollection_extract(col->geoms[i], type);
|
|
extracted_geoms[extracted_curgeom] = colgeom->geoms;
|
|
extracted_curgeom++;
|
|
if( colgeom->bbox ) lwfree(colgeom->bbox);
|
|
lwfree(colgeom);
|
|
continue;
|
|
}
|
|
}
|
|
TODO: pramsey
|
|
}
|
|
*/
|
|
|
|
|
|
|
|
void lwcollection_free(LWCOLLECTION *col)
|
|
{
|
|
int i;
|
|
if( col->bbox )
|
|
{
|
|
lwfree(col->bbox);
|
|
}
|
|
for ( i = 0; i < col->ngeoms; i++ )
|
|
{
|
|
if( col->geoms[i] ) {
|
|
switch( TYPE_GETTYPE(col->geoms[i]->type) )
|
|
{
|
|
case POINTTYPE:
|
|
lwpoint_free((LWPOINT*)col->geoms[i]);
|
|
break;
|
|
case LINETYPE:
|
|
lwline_free((LWLINE*)col->geoms[i]);
|
|
break;
|
|
case POLYGONTYPE:
|
|
lwpoly_free((LWPOLY*)col->geoms[i]);
|
|
break;
|
|
case MULTIPOINTTYPE:
|
|
lwmpoint_free((LWMPOINT*)col->geoms[i]);
|
|
break;
|
|
case MULTILINETYPE:
|
|
lwmline_free((LWMLINE*)col->geoms[i]);
|
|
break;
|
|
case MULTIPOLYGONTYPE:
|
|
lwmpoly_free((LWMPOLY*)col->geoms[i]);
|
|
break;
|
|
case COLLECTIONTYPE:
|
|
lwcollection_free((LWCOLLECTION*)col->geoms[i]);
|
|
break;
|
|
}
|
|
}
|
|
}
|
|
if( col->geoms )
|
|
{
|
|
lwfree(col->geoms);
|
|
}
|
|
lwfree(col);
|
|
|
|
};
|
|
|
|
BOX3D *lwcollection_compute_box3d(LWCOLLECTION *col)
|
|
{
|
|
int i;
|
|
BOX3D *boxfinal = NULL;
|
|
BOX3D *boxtmp1 = NULL;
|
|
BOX3D *boxtmp2 = NULL;
|
|
for ( i = 0; i < col->ngeoms; i++ )
|
|
{
|
|
if( col->geoms[i] ) {
|
|
switch( TYPE_GETTYPE(col->geoms[i]->type) )
|
|
{
|
|
case POINTTYPE:
|
|
boxtmp1 = lwpoint_compute_box3d((LWPOINT*)(col->geoms[i]));
|
|
break;
|
|
case LINETYPE:
|
|
boxtmp1 = lwline_compute_box3d((LWLINE*)(col->geoms[i]));
|
|
break;
|
|
case POLYGONTYPE:
|
|
boxtmp1 = lwpoly_compute_box3d((LWPOLY*)(col->geoms[i]));
|
|
break;
|
|
case CIRCSTRINGTYPE:
|
|
boxtmp1 = lwcircstring_compute_box3d((LWCIRCSTRING *)(col->geoms[i]));
|
|
break;
|
|
case COMPOUNDTYPE:
|
|
case CURVEPOLYTYPE:
|
|
case MULTIPOINTTYPE:
|
|
case MULTILINETYPE:
|
|
case MULTIPOLYGONTYPE:
|
|
case MULTICURVETYPE:
|
|
case MULTISURFACETYPE:
|
|
case COLLECTIONTYPE:
|
|
boxtmp1 = lwcollection_compute_box3d((LWCOLLECTION*)(col->geoms[i]));
|
|
boxfinal = box3d_union(boxtmp1, boxtmp2);
|
|
break;
|
|
}
|
|
boxtmp2 = boxfinal;
|
|
boxfinal = box3d_union(boxtmp1, boxtmp2);
|
|
if( boxtmp1 && boxtmp1 != boxfinal )
|
|
{
|
|
lwfree(boxtmp1);
|
|
boxtmp1 = NULL;
|
|
}
|
|
if( boxtmp2 && boxtmp2 != boxfinal )
|
|
{
|
|
lwfree(boxtmp2);
|
|
boxtmp2 = NULL;
|
|
}
|
|
}
|
|
}
|
|
return boxfinal;
|
|
}
|