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
|
/**********************************************************************
*
* PostGIS - Spatial Types for PostgreSQL
* http://postgis.net
*
* Copyright (C) 2014 Sandro Santilli <strk@kbt.io>
*
* 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 "CUnit/Basic.h"
#include "cu_tester.h"
#include "liblwgeom.h"
#include "liblwgeom_internal.h"
static void test_lwgeom_clip_by_rect(void)
{
#if POSTGIS_GEOS_VERSION >= 35
LWGEOM *in, *out;
const char *wkt;
char *tmp;
/* Because i don't trust that much prior tests... ;) */
cu_error_msg_reset();
wkt = "LINESTRING(0 0, 5 5, 10 0)";
in = lwgeom_from_wkt(wkt, LW_PARSER_CHECK_NONE);
out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
tmp = lwgeom_to_ewkt(out);
/* printf("%s\n", tmp); */
CU_ASSERT_STRING_EQUAL("LINESTRING(5 5,10 0)", tmp)
lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
wkt = "LINESTRING EMPTY";
in = lwgeom_from_wkt(wkt, LW_PARSER_CHECK_NONE);
out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
tmp = lwgeom_to_ewkt(out);
/* printf("%s\n", tmp); */
CU_ASSERT_STRING_EQUAL(wkt, tmp)
lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
wkt = "MULTIPOINT EMPTY";
in = lwgeom_from_wkt(wkt, LW_PARSER_CHECK_NONE);
out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
tmp = lwgeom_to_ewkt(out);
/* printf("%s\n", tmp); */
CU_ASSERT_STRING_EQUAL(wkt, tmp)
lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
wkt = "MULTIPOINT(0 0, 6 6, 7 5)";
in = lwgeom_from_wkt(wkt, LW_PARSER_CHECK_NONE);
out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
tmp = lwgeom_to_ewkt(out);
/* printf("%s\n", tmp); */
CU_ASSERT_STRING_EQUAL("MULTIPOINT(6 6,7 5)", tmp)
lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
/* Disjoint polygon */
wkt = "POLYGON((311017 4773762,311016 4773749,311006 4773744,310990 4773748,310980 4773758,310985 4773771,311003 4773776,311017 4773762))";
in = lwgeom_from_wkt(wkt, LW_PARSER_CHECK_NONE);
out = lwgeom_clip_by_rect(in, -80, -80, 80, 80);
//tmp = lwgeom_to_ewkt(out); printf("%s\n", tmp); lwfree(tmp);
CU_ASSERT(lwgeom_is_empty(out));
lwgeom_free(out); lwgeom_free(in);
#endif /* POSTGIS_GEOS_VERSION >= 35 */
}
/*
** Used by test harness to register the tests in this file.
*/
void clip_by_rect_suite_setup(void);
void clip_by_rect_suite_setup(void)
{
CU_pSuite suite = CU_add_suite("clip_by_rectangle", NULL, NULL);
PG_ADD_TEST(suite, test_lwgeom_clip_by_rect);
}
|
