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
|
/**********************************************************************
*
* PostGIS - Spatial Types for PostgreSQL
* http://postgis.net
*
* PostGIS is free software: you can redistribute it and/or modify
* it under the terms of the GNU General Public License as published by
* the Free Software Foundation, either version 2 of the License, or
* (at your option) any later version.
*
* PostGIS is distributed in the hope that it will be useful,
* but WITHOUT ANY WARRANTY; without even the implied warranty of
* MERCHANTABILITY or FITNESS FOR A PARTICULAR PURPOSE. See the
* GNU General Public License for more details.
*
* You should have received a copy of the GNU General Public License
* along with PostGIS. If not, see <http://www.gnu.org/licenses/>.
*
**********************************************************************
*
* Copyright (C) 2001-2006 Refractions Research Inc.
*
**********************************************************************/
#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(int32_t 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(int32_t srid, const POINTARRAY *pa)
{
uint32_t 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)
{
uint32_t 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);
}
LWMPOINT*
lwmpoint_from_lwgeom(const LWGEOM *g)
{
LWPOINTITERATOR* it = lwpointiterator_create(g);
int has_z = lwgeom_has_z(g);
int has_m = lwgeom_has_m(g);
LWMPOINT* result = lwmpoint_construct_empty(g->srid, has_z, has_m);
POINT4D p;
while(lwpointiterator_next(it, &p)) {
LWPOINT* lwp = lwpoint_make(g->srid, has_z, has_m, &p);
lwmpoint_add_lwpoint(result, lwp);
}
lwpointiterator_destroy(it);
return result;
}
|
