test_gdal_warp()

static void test_gdal_warp ( )

static

Definition at line 528 of file cu_gdal.c.

References ovdump::band, cu_add_band(), cu_free_raster(), ES_NONE, gdal_suite_setup(), PT_64BF, rtpixdump::rast, rtrowdump::raster, rt_band_get_hasnodata_flag(), rt_band_get_nodata(), rt_band_get_pixel(), rt_band_set_pixel(), rt_raster_gdal_warp(), rt_raster_get_band(), rt_raster_get_height(), rt_raster_get_num_bands(), rt_raster_get_width(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), genraster::value, pixval::x, and pixval::y.

Referenced by gdal_suite_setup().

                             {
        rt_pixtype pixtype = PT_64BF;
        rt_band band = NULL;

        rt_raster raster;
        rt_raster rast;
        uint32_t x;
        uint32_t width = 100;
        uint32_t y;
        uint32_t height = 100;
        double value = 0;

        char src_srs[] = "PROJCS[\"unnamed\",GEOGCS[\"unnamed ellipse\",DATUM[\"unknown\",SPHEROID[\"unnamed\",6370997,0]],PRIMEM[\"Greenwich\",0],UNIT[\"degree\",0.0174532925199433]],PROJECTION[\"Lambert_Azimuthal_Equal_Area\"],PARAMETER[\"latitude_of_center\",45],PARAMETER[\"longitude_of_center\",-100],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",0],UNIT[\"Meter\",1],AUTHORITY[\"EPSG\",\"2163\"]]";

        char dst_srs[] = "PROJCS[\"NAD83 / California Albers\",GEOGCS[\"NAD83\",DATUM[\"North_American_Datum_1983\",SPHEROID[\"GRS 1980\",6378137,298.257222101,AUTHORITY[\"EPSG\",\"7019\"]],AUTHORITY[\"EPSG\",\"6269\"]],PRIMEM[\"Greenwich\",0,AUTHORITY[\"EPSG\",\"8901\"]],UNIT[\"degree\",0.01745329251994328,AUTHORITY[\"EPSG\",\"9122\"]],AUTHORITY[\"EPSG\",\"4269\"]],UNIT[\"metre\",1,AUTHORITY[\"EPSG\",\"9001\"]],PROJECTION[\"Albers_Conic_Equal_Area\"],PARAMETER[\"standard_parallel_1\",34],PARAMETER[\"standard_parallel_2\",40.5],PARAMETER[\"latitude_of_center\",0],PARAMETER[\"longitude_of_center\",-120],PARAMETER[\"false_easting\",0],PARAMETER[\"false_northing\",-4000000],AUTHORITY[\"EPSG\",\"3310\"],AXIS[\"X\",EAST],AXIS[\"Y\",NORTH]]";

        raster = rt_raster_new(width, height);
        CU_ASSERT(raster != NULL); /* or we're out of virtual memory */

        band = cu_add_band(raster, pixtype, 1, 0);
        CU_ASSERT(band != NULL);

        rt_raster_set_offsets(raster, -500000, 600000);
        rt_raster_set_scale(raster, 1000, -1000);

        for (x = 0; x < width; x++) {
                for (y = 0; y < height; y++) {
                        rt_band_set_pixel(band, x, y, (((double) x * y) + (x + y) + (x + y * x)) / (x + y + 1), NULL);
                }
        }

        rast = rt_raster_gdal_warp(
                raster,
                src_srs, dst_srs,
                NULL, NULL,
                NULL, NULL,
                NULL, NULL,
                NULL, NULL,
                NULL, NULL,
                GRA_NearestNeighbour, -1
        );
        CU_ASSERT(rast != NULL);
        CU_ASSERT_EQUAL(rt_raster_get_width(rast), 122);
        CU_ASSERT_EQUAL(rt_raster_get_height(rast), 116);
        CU_ASSERT_NOT_EQUAL(rt_raster_get_num_bands(rast), 0);

        band = rt_raster_get_band(rast, 0);
        CU_ASSERT(band != NULL);

        CU_ASSERT(rt_band_get_hasnodata_flag(band));
        rt_band_get_nodata(band, &value);
        CU_ASSERT_DOUBLE_EQUAL(value, 0., DBL_EPSILON);

        CU_ASSERT_EQUAL(rt_band_get_pixel(band, 0, 0, &value, NULL), ES_NONE);
        CU_ASSERT_DOUBLE_EQUAL(value, 0., DBL_EPSILON);

        cu_free_raster(rast);
        cu_free_raster(raster);
}

Here is the call graph for this function:

Here is the caller graph for this function: