test_spheroid_area()

static void test_spheroid_area ( void  )

static

Definition at line 1341 of file cu_geodetic.c.

References GBOX::flags, gflags(), LW_PARSER_CHECK_NONE, lwgeom_area_sphere(), lwgeom_area_spheroid(), lwgeom_calculate_gbox_geodetic(), lwgeom_free(), lwgeom_from_wkt(), s, spheroid_init(), WGS84_MAJOR_AXIS, and WGS84_MINOR_AXIS.

Referenced by geodetic_suite_setup().

{
        LWGEOM *lwg;
        GBOX gbox;
        double a1, a2;
        SPHEROID s;

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

        gbox.flags = gflags(0, 0, 1);

        /* Medford lot test polygon */
        lwg = lwgeom_from_wkt("POLYGON((-122.848227067007 42.5007249610493,-122.848309475585 42.5007179884263,-122.848327688675 42.500835880696,-122.848245279942 42.5008428533324,-122.848227067007 42.5007249610493))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        /* sphere: Planimeter -E -p 20 -e $WGS84_SPHERE -r --input-string \
        "42.5007249610493 -122.848227067007;42.5007179884263 -122.848309475585;"\
        "42.500835880696 -122.848327688675;42.5008428533324 -122.848245279942" */
        a1 = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a1, 89.721147136698008, 0.1);
        /* spheroid: Planimeter -E -p 20 -r --input-string \
        "42.5007249610493 -122.848227067007;42.5007179884263 -122.848309475585;"\
        "42.500835880696 -122.848327688675;42.5008428533324 -122.848245279942" */
        a2 = lwgeom_area_spheroid(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a2, 89.868413479309585, 0.1);
        lwgeom_free(lwg);

        /* Big-ass polygon */
        lwg = lwgeom_from_wkt("POLYGON((-2 3, -2 4, -1 4, -1 3, -2 3))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        /* sphere: Planimeter -E -p 20 -e $WGS84_SPHERE -r --input-string "3 -2;4 -2;4 -1;3 -1" */
        a1 = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a1, 12341436880.106982993974659, 0.1);
        /* spheroid: Planimeter -E -p 20 -r --input-string "3 -2;4 -2;4 -1;3 -1" */
#if PROJ_GEODESIC
        // printf("XXXXX %d\n", PJ_VERSION);
        a2 = lwgeom_area_spheroid(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a2, 12286884908.946891319597874, 0.1);
#endif
        lwgeom_free(lwg);

        /* One-degree square */
        lwg = lwgeom_from_wkt("POLYGON((8.5 2,8.5 1,9.5 1,9.5 2,8.5 2))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        /* sphere: Planimeter -E -p 20 -e $WGS84_SPHERE --input-string "2 8.5;1 8.5;1 9.5;2 9.5" */
        a1 = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a1, 12360265021.368023059138681, 0.1);
        /* spheroid: Planimeter -E -p 20 --input-string "2 8.5;1 8.5;1 9.5;2 9.5" */
#if PROJ_GEODESIC
        a2 = lwgeom_area_spheroid(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a2, 12305128751.042900673161556, 0.1);
#endif
        lwgeom_free(lwg);

        /* One-degree square *near* the antimeridian */
        lwg = lwgeom_from_wkt("POLYGON((179.5 2,179.5 1,178.5 1,178.5 2,179.5 2))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        /* sphere: Planimeter -E -p 20 -e $WGS84_SPHERE -r --input-string "2 179.5;1 179.5;1 178.5;2 178.5" */
        a1 = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a1, 12360265021.368023059138681, 0.1);
        /* spheroid: Planimeter -E -p 20 -r --input-string "2 179.5;1 179.5;1 178.5;2 178.5" */
#if PROJ_GEODESIC
        a2 = lwgeom_area_spheroid(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a2, 12305128751.042900673161556, 0.1);
#endif
        lwgeom_free(lwg);

        /* One-degree square *across* the antimeridian */
        lwg = lwgeom_from_wkt("POLYGON((179.5 2,179.5 1,-179.5 1,-179.5 2,179.5 2))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        /* sphere: Planimeter -E -p 20 -e $WGS84_SPHERE --input-string "2 179.5;1 179.5;1 -179.5;2 -179.5" */
        a1 = lwgeom_area_sphere(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a1, 12360265021.368023059138681, 0.1);
        /* spheroid: Planimeter -E -p 20 --input-string "2 179.5;1 179.5;1 -179.5;2 -179.5" */
#if PROJ_GEODESIC
        a2 = lwgeom_area_spheroid(lwg, &s);
        CU_ASSERT_DOUBLE_EQUAL(a2, 12305128751.042900673161556, 0.1);
#endif
        lwgeom_free(lwg);
}

Here is the call graph for this function:

Here is the caller graph for this function: