test_lwgeom_area_sphere()

static void test_lwgeom_area_sphere ( void  )

static

Definition at line 1561 of file cu_geodetic.c.

{
        LWGEOM *lwg;
        double area;
        SPHEROID s;
 
        /* Init to Sphere */
        spheroid_init(&s, WGS84_RADIUS, WGS84_RADIUS);
        s.a = s.b = s.radius;
 
        /* Simple case */
        lwg = lwgeom_from_wkt("POLYGON((1 1, 1 2, 2 2, 2 1, 1 1))", LW_PARSER_CHECK_NONE);
        area = lwgeom_area_sphere(lwg, &s);
 
        CU_ASSERT_DOUBLE_EQUAL(area, 12360265021.3561, 1.0);
        lwgeom_free(lwg);
 
        /* Robustness tests, from ticket #3393 */
        lwg = lwgeom_from_wkt("POLYGON((0 78.703946026663,0 0,179.999997913235 0,179.999997913235 -33.0888306884702,0 78.703946026663))", LW_PARSER_CHECK_NONE);
        area = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(area, 127516466960671, 1.0);
        lwgeom_free(lwg);
 
        lwg = lwgeom_from_wkt("POLYGON((0 78.703946026662,0 0,179.999997913235 0,179.999997913235 -33.0888306884702,0 78.703946026662))", LW_PARSER_CHECK_NONE);
        area = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(area, 127516466960671, 1.0);
        lwgeom_free(lwg);
 
        lwg = lwgeom_from_wkt("POLYGON((0 78.703946026664,0 0,179.999997913235 0,179.999997913235 -33.0888306884702,0 78.703946026664))", LW_PARSER_CHECK_NONE);
        area = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(area, 127516466960671, 1.0);
        lwgeom_free(lwg);
        /* end #3393 */
}

References LW_PARSER_CHECK_NONE, lwgeom_area_sphere(), lwgeom_free(), lwgeom_from_wkt(), s, spheroid_init(), and WGS84_RADIUS.

Referenced by geodetic_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: