test_gbox_utils()

static void test_gbox_utils ( void  )

static

Definition at line 1422 of file cu_geodetic.c.

{
        LWGEOM *lwg;
        GBOX gbox;
        double a1, a2;
        SPHEROID s;
        POINT2D pt;
 
        /* Init to WGS84 */
        spheroid_init(&s, WGS84_MAJOR_AXIS, WGS84_MINOR_AXIS);
 
        gbox.flags = gflags(0, 0, 1);
 
        /* One-degree square by equator */
        lwg = lwgeom_from_wkt("POLYGON((1 20,1 21,2 21,2 20,1 20))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        a1 = gbox_angular_width(&gbox);
        a2 = gbox_angular_height(&gbox);
        CU_ASSERT_DOUBLE_EQUAL(a1, 0.0177951, 0.0000001);
        CU_ASSERT_DOUBLE_EQUAL(a2, 0.017764, 0.0000001);
        lwgeom_free(lwg);
 
        /* One-degree square *across* 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);
        a1 = gbox_angular_width(&gbox);
        a2 = gbox_angular_height(&gbox);
        //printf("a1=%g a2=%g\n", a1, a2);
        CU_ASSERT_DOUBLE_EQUAL(a1, 0.0174613, 0.0000001);
        CU_ASSERT_DOUBLE_EQUAL(a2, 0.0174553, 0.0000001);
        lwgeom_free(lwg);
 
        /* One-degree square *across* antimeridian */
        lwg = lwgeom_from_wkt("POLYGON((178.5 2,178.5 1,-179.5 1,-179.5 2,178.5 2))", LW_PARSER_CHECK_NONE);
        lwgeom_calculate_gbox_geodetic(lwg, &gbox);
        gbox_centroid(&gbox, &pt);
        //printf("POINT(%g %g)\n", pt.x, pt.y);
        CU_ASSERT_DOUBLE_EQUAL(pt.x, 179.5, 0.0001);
        CU_ASSERT_DOUBLE_EQUAL(pt.y, 1.50024, 0.0001);
        lwgeom_free(lwg);
 
}

References GBOX::flags, gbox_angular_height(), gbox_angular_width(), gbox_centroid(), gflags(), LW_PARSER_CHECK_NONE, lwgeom_calculate_gbox_geodetic(), lwgeom_free(), lwgeom_from_wkt(), s, spheroid_init(), WGS84_MAJOR_AXIS, WGS84_MINOR_AXIS, POINT2D::x, and POINT2D::y.

Referenced by geodetic_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: