test_raster_to_gdal()

static void test_raster_to_gdal ( )

static

Definition at line 304 of file cu_gdal.c.

                                  {
        rt_pixtype pixtype = PT_64BF;
        rt_raster raster = NULL;
        rt_band band = NULL;
        uint32_t x;
        uint32_t width = 100;
        uint32_t y;
        uint32_t height = 100;
        char 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\"]]";
 
        uint64_t gdalSize;
        uint8_t *gdal = NULL;
 
        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);
                }
        }
 
        gdal = rt_raster_to_gdal(raster, srs, "GTiff", NULL, &gdalSize);
        /*printf("gdalSize: %d\n", (int) gdalSize);*/
        CU_ASSERT(gdalSize);
 
        /*
        FILE *fh = NULL;
        fh = fopen("/tmp/out.tif", "w");
        fwrite(gdal, sizeof(uint8_t), gdalSize, fh);
        fclose(fh);
        */
 
        if (gdal) CPLFree(gdal);
 
        cu_free_raster(raster);
 
        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, x, NULL);
                }
        }
 
        /* add check that band isn't NODATA */
        CU_ASSERT_EQUAL(rt_band_check_is_nodata(band), FALSE);
 
        gdal = rt_raster_to_gdal(raster, srs, "PNG", NULL, &gdalSize);
        /*printf("gdalSize: %d\n", (int) gdalSize);*/
        CU_ASSERT(gdalSize);
 
        if (gdal) CPLFree(gdal);
 
        gdal = rt_raster_to_gdal(raster, srs, "PCIDSK", NULL, &gdalSize);
        CU_ASSERT(gdal == NULL);
 
        cu_free_raster(raster);
}

References cu_add_band(), cu_free_raster(), FALSE, PT_64BF, rt_band_check_is_nodata(), rt_band_set_pixel(), rt_raster_new(), rt_raster_set_offsets(), rt_raster_set_scale(), and rt_raster_to_gdal().

Referenced by gdal_suite_setup().

Here is the call graph for this function:

Here is the caller graph for this function: