PostGIS  2.5.0dev-r@@SVN_REVISION@@

◆ test_lwpoly_construct_circle()

static void test_lwpoly_construct_circle ( void  )
static

Definition at line 1356 of file cu_algorithm.c.

References ASSERT_INT_EQUAL, LW_TRUE, lwgeom_area(), lwgeom_count_vertices(), lwgeom_get_bbox(), lwgeom_get_srid(), lwpoly_as_lwgeom(), lwpoly_construct_circle(), lwpoly_free(), r, pixval::x, GBOX::xmax, GBOX::xmin, pixval::y, GBOX::ymax, and GBOX::ymin.

Referenced by algorithms_suite_setup().

1357 {
1358  LWPOLY* p;
1359  const GBOX* g;
1360  const int srid = 4326;
1361  const int segments_per_quad = 17;
1362  const int x = 10;
1363  const int y = 20;
1364  const int r = 5;
1365 
1366  /* With normal arguments you should get something circle-y */
1367  p = lwpoly_construct_circle(srid, x, y, r, segments_per_quad, LW_TRUE);
1368 
1369  ASSERT_INT_EQUAL(lwgeom_count_vertices(lwpoly_as_lwgeom(p)), segments_per_quad * 4 + 1);
1371 
1373  CU_ASSERT_DOUBLE_EQUAL(g->xmin, x-r, 0.1);
1374  CU_ASSERT_DOUBLE_EQUAL(g->xmax, x+r, 0.1);
1375  CU_ASSERT_DOUBLE_EQUAL(g->ymin, y-r, 0.1);
1376  CU_ASSERT_DOUBLE_EQUAL(g->ymax, y+r, 0.1);
1377 
1378  CU_ASSERT_DOUBLE_EQUAL(lwgeom_area(lwpoly_as_lwgeom(p)), M_PI*5*5, 0.1);
1379 
1380  lwpoly_free(p);
1381 
1382  /* No segments? No circle. */
1383  p = lwpoly_construct_circle(srid, x, y, r, 0, LW_TRUE);
1384  CU_ASSERT_TRUE(p == NULL);
1385 
1386  /* Negative radius? No circle. */
1387  p = lwpoly_construct_circle(srid, x, y, -1e-3, segments_per_quad, LW_TRUE);
1388  CU_ASSERT_TRUE(p == NULL);
1389 
1390  /* Zero radius? Invalid circle */
1391  p = lwpoly_construct_circle(srid, x, y, 0, segments_per_quad, LW_TRUE);
1392  CU_ASSERT_TRUE(p != NULL);
1393  lwpoly_free(p);
1394 }
char * r
Definition: cu_in_wkt.c:24
double xmax
Definition: liblwgeom.h:295
int32_t lwgeom_get_srid(const LWGEOM *geom)
Return SRID number.
Definition: lwgeom.c:916
LWGEOM * lwpoly_as_lwgeom(const LWPOLY *obj)
Definition: lwgeom.c:320
double ymin
Definition: liblwgeom.h:296
double xmin
Definition: liblwgeom.h:294
void lwpoly_free(LWPOLY *poly)
Definition: lwpoly.c:175
#define LW_TRUE
Return types for functions with status returns.
Definition: liblwgeom.h:75
#define ASSERT_INT_EQUAL(o, e)
const GBOX * lwgeom_get_bbox(const LWGEOM *lwgeom)
Get a non-empty geometry bounding box, computing and caching it if not already there.
Definition: lwgeom.c:734
double ymax
Definition: liblwgeom.h:297
uint32_t lwgeom_count_vertices(const LWGEOM *geom)
Count the total number of vertices in any LWGEOM.
Definition: lwgeom.c:1235
double lwgeom_area(const LWGEOM *geom)
Definition: lwgeom.c:1857
LWPOLY * lwpoly_construct_circle(int srid, double x, double y, double radius, uint32_t segments_per_quarter, char exterior)
Definition: lwpoly.c:120
Here is the call graph for this function:
Here is the caller graph for this function: