PostGIS  2.5.2dev-r@@SVN_REVISION@@
cu_clip_by_rect.c
Go to the documentation of this file.
1 /**********************************************************************
2  *
3  * PostGIS - Spatial Types for PostgreSQL
4  * http://postgis.net
5  *
6  * Copyright (C) 2014 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 static void test_lwgeom_clip_by_rect(void)
20 {
21  LWGEOM *in, *out;
22  const char *wkt;
23  char *tmp;
24 
25  /* Because i don't trust that much prior tests... ;) */
27 
28  wkt = "LINESTRING(0 0, 5 5, 10 0)";
30  out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
31  tmp = lwgeom_to_ewkt(out);
32  /* printf("%s\n", tmp); */
33  CU_ASSERT_STRING_EQUAL("LINESTRING(5 5,10 0)", tmp)
34  lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
35 
36  wkt = "LINESTRING EMPTY";
38  out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
39  tmp = lwgeom_to_ewkt(out);
40  /* printf("%s\n", tmp); */
41  CU_ASSERT_STRING_EQUAL(wkt, tmp)
42  lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
43 
44  wkt = "MULTIPOINT EMPTY";
46  out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
47  tmp = lwgeom_to_ewkt(out);
48  /* printf("%s\n", tmp); */
49  CU_ASSERT_STRING_EQUAL(wkt, tmp)
50  lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
51 
52  wkt = "MULTIPOINT(0 0, 6 6, 7 5)";
54  out = lwgeom_clip_by_rect(in, 5, 0, 10, 10);
55  tmp = lwgeom_to_ewkt(out);
56  /* printf("%s\n", tmp); */
57  CU_ASSERT_STRING_EQUAL("MULTIPOINT(6 6,7 5)", tmp)
58  lwfree(tmp); lwgeom_free(out); lwgeom_free(in);
59 
60  /* Disjoint polygon */
61  wkt = "POLYGON((311017 4773762,311016 4773749,311006 4773744,310990 4773748,310980 4773758,310985 4773771,311003 4773776,311017 4773762))";
63  out = lwgeom_clip_by_rect(in, -80, -80, 80, 80);
64  //tmp = lwgeom_to_ewkt(out); printf("%s\n", tmp); lwfree(tmp);
65  CU_ASSERT(lwgeom_is_empty(out));
66  lwgeom_free(out); lwgeom_free(in);
67 }
68 
69 /*
70 ** Used by test harness to register the tests in this file.
71 */
72 void clip_by_rect_suite_setup(void);
74 {
75  CU_pSuite suite = CU_add_suite("clip_by_rectangle", NULL, NULL);
77 }
static void test_lwgeom_clip_by_rect(void)
void lwfree(void *mem)
Definition: lwutil.c:244
char * lwgeom_to_ewkt(const LWGEOM *lwgeom)
Return an alloced string.
Definition: lwgeom.c:556
void lwgeom_free(LWGEOM *geom)
Definition: lwgeom.c:1144
LWGEOM * lwgeom_clip_by_rect(const LWGEOM *geom1, double x0, double y0, double x1, double y1)
LWGEOM * lwgeom_from_wkt(const char *wkt, const char check)
Definition: lwin_wkt.c:904
void clip_by_rect_suite_setup(void)
#define LW_PARSER_CHECK_NONE
Definition: liblwgeom.h:2004
void cu_error_msg_reset()
#define PG_ADD_TEST(suite, testfunc)
int lwgeom_is_empty(const LWGEOM *geom)
Return true or false depending on whether a geometry is an "empty" geometry (no vertices members) ...
Definition: lwgeom.c:1393
This library is the generic geometry handling section of PostGIS.