PostGIS 3.7.0dev-r@@SVN_REVISION@@
Loading...
Searching...
No Matches
cu_buildarea.c
Go to the documentation of this file.
1/**********************************************************************
2 *
3 * PostGIS - Spatial Types for PostgreSQL
4 * http://postgis.net
5 *
6 * Copyright (C) 2012 Sandro Santilli <strk@kbt.io>
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 "CUnit/Basic.h"
14#include "cu_tester.h"
15
16#include "liblwgeom.h"
17#include "liblwgeom_internal.h"
18
19/*
20 * TODO: change lwgeom_same to lwgeom_equals
21 * (requires porting predicates to liblwgeom)
22 */
23#define check_geom_equal(gobt, gexp) do { \
24 char *obt, *exp; \
25 LWGEOM *ngobt, *ngexp; \
26 ngobt = lwgeom_normalize(gobt); \
27 ngexp = lwgeom_normalize(gexp); \
28 if ( ! lwgeom_same((ngobt), (ngexp)) ) { \
29 obt = lwgeom_to_wkt((ngobt), WKT_ISO, 8, NULL); \
30 exp = lwgeom_to_wkt((ngexp), WKT_ISO, 8, NULL); \
31 printf(" Failure at %s:%d\n", __FILE__, __LINE__); \
32 printf(" Exp: %s\n", exp); \
33 printf(" Obt: %s\n", obt); \
34 free(obt); free(exp); \
35 lwgeom_free(ngobt); lwgeom_free(ngexp); \
36 CU_ASSERT(0); \
37 } else { \
38 lwgeom_free(ngobt); lwgeom_free(ngexp); \
39 CU_ASSERT(1); \
40 } \
41} while (0)
42
43/*
44 +-----+
45 | |
46 +-----+-----+
47 | |
48 +-----+
49*/
50static void buildarea1(void)
51{
52 LWGEOM *gin, *gout, *gexp;
53
54 cu_error_msg_reset();
55
56 gin = lwgeom_from_wkt(
57"MULTILINESTRING((0 0, 10 0, 10 10, 0 10, 0 0),(10 10, 20 10, 20 20, 10 20, 10 10))",
58 LW_PARSER_CHECK_NONE);
59 CU_ASSERT( gin != NULL );
60
61 gexp = lwgeom_from_wkt(
62"MULTIPOLYGON(((0 0,0 10,10 10,10 0,0 0)),((10 10,10 20,20 20,20 10,10 10)))",
63 LW_PARSER_CHECK_NONE);
64 CU_ASSERT( gexp != NULL );
65
66 gout = lwgeom_buildarea(gin);
67 CU_ASSERT(gout != NULL);
68
69 check_geom_equal(gout, gexp);
70
71 lwgeom_free(gout);
72 lwgeom_free(gexp);
73 lwgeom_free(gin);
74}
75
76/*
77 +-----+-----+
78 | | |
79 +-----+-----+
80*/
81static void buildarea2(void)
82{
83 LWGEOM *gin, *gout, *gexp;
84
85 /* because i don't trust that much prior tests... ;) */
86 cu_error_msg_reset();
87
88 gin = lwgeom_from_wkt(
89"MULTILINESTRING((0 0, 10 0, 10 10, 0 10, 0 0),(10 10, 20 10, 20 0, 10 0, 10 10))",
90 LW_PARSER_CHECK_NONE);
91 CU_ASSERT(gin != NULL);
92
93 gexp = lwgeom_from_wkt(
94"POLYGON((0 0,0 10,10 10,20 10,20 0,10 0,0 0))",
95 LW_PARSER_CHECK_NONE);
96 CU_ASSERT(gexp != NULL);
97
98 gout = lwgeom_buildarea(gin);
99 CU_ASSERT(gout != NULL);
100
101 check_geom_equal(gout, gexp);
102
103 lwgeom_free(gout);
104 lwgeom_free(gexp);
105 lwgeom_free(gin);
106
107}
108
109/*
110 +-----------+
111 | +-----+ |
112 | | | |
113 | +-----+ |
114 +-----------+
115*/
116static void buildarea3(void)
117{
118 LWGEOM *gin, *gout, *gexp;
119
120 cu_error_msg_reset();
121
122 gin = lwgeom_from_wkt(
123"MULTILINESTRING((0 0, 20 0, 20 20, 0 20, 0 0),(2 2, 18 2, 18 18, 2 18, 2 2))",
124 LW_PARSER_CHECK_NONE);
125 CU_ASSERT(gin != NULL);
126
127 gexp = lwgeom_from_wkt(
128"POLYGON((0 0,0 20,20 20,20 0,0 0),(2 2,18 2,18 18,2 18,2 2))",
129 LW_PARSER_CHECK_NONE);
130 CU_ASSERT(gexp != NULL);
131
132 gout = lwgeom_buildarea(gin);
133 CU_ASSERT(gout != NULL);
134
135 check_geom_equal(gout, gexp);
136
137 lwgeom_free(gout);
138 lwgeom_free(gexp);
139 lwgeom_free(gin);
140
141}
142
143/*
144 +-----------+
145 | +-----+ |
146 | | +-+ | |
147 | | | | | |
148 | | +-+ | |
149 | +-----+ |
150 +-----------+
151*/
152static void buildarea4(void)
153{
154 LWGEOM *gin, *gout, *gexp;
155
156 cu_error_msg_reset();
157
158 gin = lwgeom_from_wkt(
159"MULTILINESTRING((0 0, 20 0, 20 20, 0 20, 0 0),(2 2, 18 2, 18 18, 2 18, 2 2),(8 8, 8 12, 12 12, 12 8, 8 8))",
160 LW_PARSER_CHECK_NONE);
161 CU_ASSERT(gin != NULL);
162
163 gexp = lwgeom_from_wkt(
164"MULTIPOLYGON(((0 0,0 20,20 20,20 0,0 0),(2 2,18 2,18 18,2 18,2 2)),((8 8,8 12,12 12,12 8,8 8)))",
165 LW_PARSER_CHECK_NONE);
166 CU_ASSERT(gexp != NULL);
167
168 gout = lwgeom_buildarea(gin);
169 CU_ASSERT(gout != NULL);
170
171 check_geom_equal(gout, gexp);
172
173 lwgeom_free(gout);
174 lwgeom_free(gexp);
175 lwgeom_free(gin);
176
177}
178
179/*
180 +-----------+
181 | +-----+ | This time the innermost ring has
182 | | +-+ | | more points than the other (outer) two.
183 | | | | | |
184 | | +-+ | |
185 | +-----+ |
186 +-----------+
187*/
188static void buildarea4b(void)
189{
190 LWGEOM *gin, *gout, *gexp;
191
192 cu_error_msg_reset();
193
194 gin = lwgeom_from_wkt(
195"MULTILINESTRING((0 0, 20 0, 20 20, 0 20, 0 0),(2 2, 18 2, 18 18, 2 18, 2 2), (8 8, 8 9, 8 10, 8 11, 8 12, 9 12, 10 12, 11 12, 12 12, 12 11, 12 10, 12 9, 12 8, 11 8, 10 8, 9 8, 8 8))",
196 LW_PARSER_CHECK_NONE);
197 CU_ASSERT(gin != NULL);
198
199 gexp = lwgeom_from_wkt(
200"MULTIPOLYGON(((0 0,0 20,20 20,20 0,0 0),(2 2,18 2,18 18,2 18,2 2)),((8 8,8 9,8 10,8 11,8 12,9 12,10 12,11 12,12 12,12 11,12 10,12 9,12 8,11 8,10 8,9 8,8 8)))",
201 LW_PARSER_CHECK_NONE);
202 CU_ASSERT(gexp != NULL);
203
204 gout = lwgeom_buildarea(gin);
205 CU_ASSERT(gout != NULL);
206
207 check_geom_equal(gout, gexp);
208
209 lwgeom_free(gout);
210 lwgeom_free(gexp);
211 lwgeom_free(gin);
212
213}
214
215/*
216 +---------------+
217 | +---------+ |
218 | | +--+--+ | |
219 | | | | | | |
220 | | +--+--+ | |
221 | +---------+ |
222 +---------------+
223*/
224static void buildarea5(void)
225{
226 LWGEOM *gin, *gout, *gexp;
227
228 cu_error_msg_reset();
229
230 gin = lwgeom_from_wkt(
231"MULTILINESTRING((0 0, 20 0, 20 20, 0 20, 0 0),(2 2, 18 2, 18 18, 2 18, 2 2),(8 8, 8 12, 12 12, 12 8, 8 8),(10 8, 10 12))",
232 LW_PARSER_CHECK_NONE);
233 CU_ASSERT(gin != NULL);
234
235 gexp = lwgeom_from_wkt(
236"MULTIPOLYGON(((0 0,0 20,20 20,20 0,0 0),(2 2,18 2,18 18,2 18,2 2)),((8 8,8 12,12 12,12 8,8 8)))",
237 LW_PARSER_CHECK_NONE);
238 CU_ASSERT(gexp != NULL);
239
240 gout = lwgeom_buildarea(gin);
241 CU_ASSERT(gout != NULL);
242
243 check_geom_equal(gout, gexp);
244
245 lwgeom_free(gout);
246 lwgeom_free(gexp);
247 lwgeom_free(gin);
248
249}
250
251/*
252 +---------------+
253 | +----+----+ |
254 | | | | |
255 | | | | |
256 | | | | |
257 | +----+----+ |
258 +---------------+
259*/
260static void buildarea6(void)
261{
262 LWGEOM *gin, *gout, *gexp;
263
264 cu_error_msg_reset();
265
266 gin = lwgeom_from_wkt(
267"MULTILINESTRING((0 0, 20 0, 20 20, 0 20, 0 0),(2 2, 18 2, 18 18, 2 18, 2 2),(10 2, 10 18))",
268 LW_PARSER_CHECK_NONE);
269 CU_ASSERT(gin != NULL);
270
271 gexp = lwgeom_from_wkt(
272"POLYGON((0 0,0 20,20 20,20 0,0 0),(2 2,18 2,18 18,2 18,2 2))",
273 LW_PARSER_CHECK_NONE);
274 CU_ASSERT(gexp != NULL);
275
276 gout = lwgeom_buildarea(gin);
277 CU_ASSERT(gout != NULL);
278
279 check_geom_equal(gout, gexp);
280
281 lwgeom_free(gout);
282 lwgeom_free(gexp);
283 lwgeom_free(gin);
284
285}
286
287/*
288 +--------------------+ +-------+
289 | +-----+ +----+ | | +---+ |
290 | | +-+ | | | | | | | |
291 | | | | | +----+ | | +---+ |
292 | | +-+ | | | | | |
293 | | | | | | | | | |
294 | | +-+ | | | | | |
295 | +-----+ +----+ | | |
296 +--------------------+ +-------+
297*/
298static void buildarea7(void)
299{
300 LWGEOM *gin, *gout, *gexp;
301
302 cu_error_msg_reset();
303
304 gin = lwgeom_from_wkt(
305"MULTILINESTRING( (0 0, 70 0, 70 70, 0 70, 0 0), (10 10, 10 60, 40 60, 40 10, 10 10), (20 20, 20 30, 30 30, 30 20, 20 20), (20 30, 30 30, 30 50, 20 50, 20 30), (50 20, 60 20, 60 40, 50 40, 50 20), (50 40, 60 40, 60 60, 50 60, 50 40), (80 0, 110 0, 110 70, 80 70, 80 0), (90 60, 100 60, 100 50, 90 50, 90 60))",
306 LW_PARSER_CHECK_NONE);
307 CU_ASSERT(gin != NULL);
308
309 gexp = lwgeom_from_wkt(
310"MULTIPOLYGON(((80 0,80 70,110 70,110 0,80 0),(90 60,90 50,100 50,100 60,90 60)),((20 20,20 30,20 50,30 50,30 30,30 20,20 20)),((0 0,0 70,70 70,70 0,0 0),(10 10,40 10,40 60,10 60,10 10),(50 20,60 20,60 40,60 60,50 60,50 40,50 20)))",
311 LW_PARSER_CHECK_NONE);
312 CU_ASSERT(gexp != NULL);
313
314 gout = lwgeom_buildarea(gin);
315 CU_ASSERT(gout != NULL);
316
317 check_geom_equal(gout, gexp);
318
319 lwgeom_free(gout);
320 lwgeom_free(gexp);
321 lwgeom_free(gin);
322
323}
324
325
326void buildarea_suite_setup(void);
327void buildarea_suite_setup(void)
328{
329 CU_pSuite suite = CU_add_suite("buildarea", NULL, NULL);
330 PG_ADD_TEST(suite,buildarea1);
331 PG_ADD_TEST(suite,buildarea2);
332 PG_ADD_TEST(suite,buildarea3);
333 PG_ADD_TEST(suite,buildarea4);
334 PG_ADD_TEST(suite,buildarea4b);
335 PG_ADD_TEST(suite,buildarea5);
336 PG_ADD_TEST(suite,buildarea6);
337 PG_ADD_TEST(suite,buildarea7);
338}
buildarea4b
static void buildarea4b(void)
Definition cu_buildarea.c:188
buildarea2
static void buildarea2(void)
Definition cu_buildarea.c:81
buildarea7
static void buildarea7(void)
Definition cu_buildarea.c:298
check_geom_equal
#define check_geom_equal(gobt, gexp)
Definition cu_buildarea.c:23
buildarea4
static void buildarea4(void)
Definition cu_buildarea.c:152
buildarea1
static void buildarea1(void)
Definition cu_buildarea.c:50
buildarea6
static void buildarea6(void)
Definition cu_buildarea.c:260
buildarea5
static void buildarea5(void)
Definition cu_buildarea.c:224
buildarea3
static void buildarea3(void)
Definition cu_buildarea.c:116
buildarea_suite_setup
void buildarea_suite_setup(void)
Definition cu_buildarea.c:327
cu_error_msg_reset
void cu_error_msg_reset()
Definition liblwgeom/cunit/cu_tester.c:307
PG_ADD_TEST
#define PG_ADD_TEST(suite, testfunc)
Definition liblwgeom/cunit/cu_tester.h:21
lwgeom_free
void lwgeom_free(LWGEOM *geom)
Definition lwgeom.c:1246
LW_PARSER_CHECK_NONE
#define LW_PARSER_CHECK_NONE
Definition liblwgeom.h:2149
lwgeom_buildarea
LWGEOM * lwgeom_buildarea(const LWGEOM *geom)
Take a geometry and return an areal geometry (Polygon or MultiPolygon).
Definition liblwgeom/lwgeom_geos.c:1124
lwgeom_from_wkt
LWGEOM * lwgeom_from_wkt(const char *wkt, const char check)
Definition lwin_wkt.c:940
liblwgeom.h
This library is the generic geometry handling section of PostGIS.
liblwgeom_internal.h
cu_tester.h