source: trunk/liblwgeom/lwmpoint.c @ 5822

Last change on this file since 5822 was 5822, checked in by colivier, 6 years ago

Few lwtype_name add

  • Property svn:keywords set to Author Date Id Revision
File size: 3.1 KB
Line 
1/**********************************************************************
2 * $Id: lwmpoint.c 5822 2010-08-13 17:30:10Z colivier $
3 *
4 * PostGIS - Spatial Types for PostgreSQL
5 * http://postgis.refractions.net
6 * Copyright 2001-2006 Refractions Research Inc.
7 *
8 * This is free software; you can redistribute and/or modify it under
9 * the terms of the GNU General Public Licence. See the COPYING file.
10 *
11 **********************************************************************/
12
13#include <stdio.h>
14#include <stdlib.h>
15#include <string.h>
16#include "liblwgeom.h"
17
18void
19lwmpoint_release(LWMPOINT *lwmpoint)
20{
21        lwgeom_release(lwmpoint_as_lwgeom(lwmpoint));
22}
23
24LWMPOINT *
25lwmpoint_construct_empty(int srid, char hasz, char hasm)
26{
27        LWMPOINT *ret = (LWMPOINT*)lwcollection_construct_empty(srid, hasz, hasm);
28        TYPE_SETTYPE(ret->type, MULTIPOINTTYPE);
29        return ret;
30}
31
32LWMPOINT *
33lwmpoint_deserialize(uchar *srl)
34{
35        LWMPOINT *result;
36        LWGEOM_INSPECTED *insp;
37        int type = lwgeom_getType(srl[0]);
38        int i;
39
40        if ( type != MULTIPOINTTYPE )
41        {
42                lwerror("lwmpoint_deserialize called on NON multipoint: %d - %s",
43                        type, lwtype_name(type));
44                return NULL;
45        }
46
47        insp = lwgeom_inspect(srl);
48
49        result = lwalloc(sizeof(LWMPOINT));
50        result->type = insp->type;
51        result->SRID = insp->SRID;
52        result->ngeoms = insp->ngeometries;
53
54        if ( insp->ngeometries )
55        {
56                result->geoms = lwalloc(sizeof(LWPOINT *)*insp->ngeometries);
57        }
58        else
59        {
60                result->geoms = NULL;
61        }
62
63        if (lwgeom_hasBBOX(srl[0]))
64        {
65                result->bbox = lwalloc(sizeof(BOX2DFLOAT4));
66                memcpy(result->bbox, srl+1, sizeof(BOX2DFLOAT4));
67        }
68        else
69        {
70                result->bbox = NULL;
71        }
72
73        for (i=0; i<insp->ngeometries; i++)
74        {
75                result->geoms[i] = lwpoint_deserialize(insp->sub_geoms[i]);
76                if ( TYPE_NDIMS(result->geoms[i]->type) != TYPE_NDIMS(result->type) )
77                {
78                        lwerror("Mixed dimensions (multipoint:%d, point%d:%d)",
79                                TYPE_NDIMS(result->type), i,
80                                TYPE_NDIMS(result->geoms[i]->type)
81                               );
82                        return NULL;
83                }
84        }
85
86        return result;
87}
88
89LWMPOINT* lwmpoint_add_lwpoint(LWMPOINT *mobj, const LWPOINT *obj)
90{
91        LWDEBUG(4, "Called");
92        return (LWMPOINT*)lwcollection_add_lwgeom((LWCOLLECTION*)mobj, (LWGEOM*)obj);
93}
94
95void lwmpoint_free(LWMPOINT *mpt)
96{
97        int i;
98        if ( mpt->bbox )
99        {
100                lwfree(mpt->bbox);
101        }
102        for ( i = 0; i < mpt->ngeoms; i++ )
103        {
104                if ( mpt->geoms[i] )
105                {
106                        lwpoint_free(mpt->geoms[i]);
107                }
108        }
109        if ( mpt->geoms )
110        {
111                lwfree(mpt->geoms);
112        }
113        lwfree(mpt);
114
115}
116
117LWGEOM*
118lwmpoint_remove_repeated_points(LWMPOINT *mpoint)
119{
120        uint32 nnewgeoms;
121        uint32 i, j;
122        LWGEOM **newgeoms;
123
124        newgeoms = lwalloc(sizeof(LWGEOM *)*mpoint->ngeoms);
125        nnewgeoms = 0;
126        for (i=0; i<mpoint->ngeoms; ++i)
127        {
128                /* Brute force, may be optimized by building an index */
129                int seen=0;
130                for (j=0; j<nnewgeoms; ++j)
131                {
132                        if ( lwpoint_same((LWPOINT*)newgeoms[j],
133                                          (LWPOINT*)mpoint->geoms[i]) )
134                        {
135                                seen=1;
136                                break;
137                        }
138                }
139                if ( seen ) continue;
140                newgeoms[nnewgeoms++] = (LWGEOM*)lwpoint_clone(mpoint->geoms[i]);
141        }
142
143        return (LWGEOM*)lwcollection_construct(mpoint->type,
144                                               mpoint->SRID, mpoint->bbox ? box2d_clone(mpoint->bbox) : NULL,
145                                               nnewgeoms, newgeoms);
146
147}
148
Note: See TracBrowser for help on using the repository browser.