1
2
3
4
5
6
7
8
9
10
11
12
13
14
15
16
17
18
19
20
21
22
23
24
25
26
27
28
29
30
31
32
33
34
35
36
37
38
39
40
41
42
43
44
45
46
47
48
49
50
51
52
53
54
55
56
57
58
59
60
61
62
63
64
65
66
67
68
69
70
71
72
73
74
75
76
77
78
79
80
81
82
83
84
85
86
87
88
89
90
91
92
93
94
95
96
97
98
99
100
101
102
103
104
105
106
107
108
109
|
/**********************************************************************
*
* PostGIS - Spatial Types for PostgreSQL
* http://postgis.refractions.net
*
* Copyright (C) 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_internal.h"
#include "lwgeom_log.h"
void
lwmpoint_release(LWMPOINT *lwmpoint)
{
lwgeom_release(lwmpoint_as_lwgeom(lwmpoint));
}
LWMPOINT *
lwmpoint_construct_empty(int srid, char hasz, char hasm)
{
LWMPOINT *ret = (LWMPOINT*)lwcollection_construct_empty(MULTIPOINTTYPE, srid, hasz, hasm);
return ret;
}
LWMPOINT* lwmpoint_add_lwpoint(LWMPOINT *mobj, const LWPOINT *obj)
{
LWDEBUG(4, "Called");
return (LWMPOINT*)lwcollection_add_lwgeom((LWCOLLECTION*)mobj, (LWGEOM*)obj);
}
LWMPOINT *
lwmpoint_construct(int srid, const POINTARRAY *pa)
{
int i;
int hasz = ptarray_has_z(pa);
int hasm = ptarray_has_m(pa);
LWMPOINT *ret = (LWMPOINT*)lwcollection_construct_empty(MULTIPOINTTYPE, srid, hasz, hasm);
for ( i = 0; i < pa->npoints; i++ )
{
LWPOINT *lwp;
POINT4D p;
getPoint4d_p(pa, i, &p);
lwp = lwpoint_make(srid, hasz, hasm, &p);
lwmpoint_add_lwpoint(ret, lwp);
}
return ret;
}
void lwmpoint_free(LWMPOINT *mpt)
{
int i;
if ( ! mpt ) return;
if ( mpt->bbox )
lwfree(mpt->bbox);
for ( i = 0; i < mpt->ngeoms; i++ )
if ( mpt->geoms && mpt->geoms[i] )
lwpoint_free(mpt->geoms[i]);
if ( mpt->geoms )
lwfree(mpt->geoms);
lwfree(mpt);
}
LWGEOM*
lwmpoint_remove_repeated_points(LWMPOINT *mpoint)
{
uint32_t nnewgeoms;
uint32_t i, j;
LWGEOM **newgeoms;
newgeoms = lwalloc(sizeof(LWGEOM *)*mpoint->ngeoms);
nnewgeoms = 0;
for (i=0; i<mpoint->ngeoms; ++i)
{
/* Brute force, may be optimized by building an index */
int seen=0;
for (j=0; j<nnewgeoms; ++j)
{
if ( lwpoint_same((LWPOINT*)newgeoms[j],
(LWPOINT*)mpoint->geoms[i]) )
{
seen=1;
break;
}
}
if ( seen ) continue;
newgeoms[nnewgeoms++] = (LWGEOM*)lwpoint_clone(mpoint->geoms[i]);
}
return (LWGEOM*)lwcollection_construct(mpoint->type,
mpoint->srid, mpoint->bbox ? gbox_copy(mpoint->bbox) : NULL,
nnewgeoms, newgeoms);
}
|
