test_lwline_interpolate_point_3d()

static void test_lwline_interpolate_point_3d ( void  )

static

Definition at line 547 of file cu_algorithm.c.

{
        LWLINE *line;
        POINT4D point;
        LWPOINT *pt;
 
        /* Empty line -> Empty point*/
        line = lwline_construct_empty(4326, LW_TRUE, LW_FALSE);
        pt = lwline_interpolate_point_3d(line, 0.5);
        CU_ASSERT(lwpoint_is_empty(pt));
        CU_ASSERT(lwgeom_has_z(lwpoint_as_lwgeom(pt)));
        CU_ASSERT_FALSE(lwgeom_has_m(lwpoint_as_lwgeom(pt)));
        lwpoint_free(pt);
        lwline_free(line);
 
        line = lwgeom_as_lwline(lwgeom_from_wkt("LINESTRING Z (0 0 0, 1 1 1, 2 2 2)", LW_PARSER_CHECK_NONE));
 
        /* distance = 0 -> first point */
        pt = lwline_interpolate_point_3d(line, 0);
        lwpoint_getPoint4d_p(pt, &point);
        CU_ASSERT_DOUBLE_EQUAL(point.x, 0, 0.0001);
        CU_ASSERT_DOUBLE_EQUAL(point.y, 0, 0.001);
        CU_ASSERT_DOUBLE_EQUAL(point.z, 0, 0.001);
        lwpoint_free(pt);
 
        /* distance = 1 -> last point */
        pt = lwline_interpolate_point_3d(line, 1);
        lwpoint_getPoint4d_p(pt, &point);
        CU_ASSERT_DOUBLE_EQUAL(point.x, 2, 0.0001);
        CU_ASSERT_DOUBLE_EQUAL(point.y, 2, 0.001);
        CU_ASSERT_DOUBLE_EQUAL(point.z, 2, 0.001);
        lwpoint_free(pt);
 
        /* simple where distance 50% -> second point */
        pt = lwline_interpolate_point_3d(line, 0.5);
        lwpoint_getPoint4d_p(pt, &point);
        CU_ASSERT_DOUBLE_EQUAL(point.x, 1, 0.0001);
        CU_ASSERT_DOUBLE_EQUAL(point.y, 1, 0.001);
        CU_ASSERT_DOUBLE_EQUAL(point.z, 1, 0.001);
        lwpoint_free(pt);
 
        /* simple where distance 80% -> between second and last point */
        pt = lwline_interpolate_point_3d(line, 0.8);
        lwpoint_getPoint4d_p(pt, &point);
        CU_ASSERT_DOUBLE_EQUAL(point.x, 1.6, 0.0001);
        CU_ASSERT_DOUBLE_EQUAL(point.y, 1.6, 0.001);
        CU_ASSERT_DOUBLE_EQUAL(point.z, 1.6, 0.001);
        lwpoint_free(pt);
 
        lwline_free(line);
}

References LW_FALSE, LW_PARSER_CHECK_NONE, LW_TRUE, lwgeom_as_lwline(), lwgeom_from_wkt(), lwgeom_has_m(), lwgeom_has_z(), lwline_construct_empty(), lwline_free(), lwline_interpolate_point_3d(), lwpoint_as_lwgeom(), lwpoint_free(), lwpoint_getPoint4d_p(), lwpoint_is_empty(), POINT4D::x, POINT4D::y, and POINT4D::z.

Referenced by algorithms_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: