test_lwgeom_area_sphere()

static void test_lwgeom_area_sphere ( void  )

static

Definition at line 1544 of file cu_geodetic.c.

References area(), LW_PARSER_CHECK_NONE, lwgeom_area_sphere(), lwgeom_free(), lwgeom_from_wkt(), s, spheroid_init(), WGS84_MAJOR_AXIS, and WGS84_MINOR_AXIS.

Referenced by geodetic_suite_setup().

{
        LWGEOM *lwg;
        double area;
        SPHEROID s;

        /* Init to WGS84 */
        spheroid_init(&s, WGS84_MAJOR_AXIS, WGS84_MINOR_AXIS);

        /* 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, 127516467322130, 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, 127516467322130, 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, 127516467322130, 1.0);
        lwgeom_free(lwg);
        /* end #3393 */
}

Here is the call graph for this function:

Here is the caller graph for this function: