rt_raster_to_gdal()

uint8_t * rt_raster_to_gdal	(	rt_raster	raster,
		const char *	srs,
		char *	format,
		char **	options,
		uint64_t *	gdalsize
	)

Return formatted GDAL raster from raster.

Parameters

raster	: the raster to convert
srs	: the raster's coordinate system in OGC WKT
format	: format to convert to. GDAL driver short name
options	: list of format creation options. array of strings
gdalsize	: will be set to the size of returned bytea

Returns

formatted GDAL raster. the calling function is responsible for freeing the returned data using CPLFree()

Definition at line 1593 of file rt_raster.c.

  {
        const char *cc;
        const char *vio;
 
        GDALDriverH src_drv = NULL;
        int destroy_src_drv = 0;
        GDALDatasetH src_ds = NULL;
 
        vsi_l_offset rtn_lenvsi;
        uint64_t rtn_len = 0;
 
        GDALDriverH rtn_drv = NULL;
        GDALDatasetH rtn_ds = NULL;
        uint8_t *rtn = NULL;
 
        assert(NULL != raster);
        assert(NULL != gdalsize);
 
        /* any supported format is possible */
        rt_util_gdal_register_all(0);
        RASTER_DEBUG(3, "loaded all supported GDAL formats");
 
        /* output format not specified */
        if (format == NULL || !strlen(format))
                format = "GTiff";
        RASTER_DEBUGF(3, "output format is %s", format);
 
        /* load raster into a GDAL MEM raster */
        src_ds = rt_raster_to_gdal_mem(raster, srs, NULL, NULL, 0, &src_drv, &destroy_src_drv);
        if (NULL == src_ds) {
                rterror("rt_raster_to_gdal: Could not convert raster to GDAL MEM format");
                return 0;
        }
 
        /* load driver */
        rtn_drv = GDALGetDriverByName(format);
        if (NULL == rtn_drv) {
                rterror("rt_raster_to_gdal: Could not load the output GDAL driver");
                GDALClose(src_ds);
                if (destroy_src_drv) GDALDestroyDriver(src_drv);
                return 0;
        }
        RASTER_DEBUG(3, "Output driver loaded");
 
        /* CreateCopy support */
        cc = GDALGetMetadataItem(rtn_drv, GDAL_DCAP_CREATECOPY, NULL);
        /* VirtualIO support */
        vio = GDALGetMetadataItem(rtn_drv, GDAL_DCAP_VIRTUALIO, NULL);
 
        if (cc == NULL || vio == NULL) {
                rterror("rt_raster_to_gdal: Output GDAL driver does not support CreateCopy and/or VirtualIO");
                GDALClose(src_ds);
                if (destroy_src_drv) GDALDestroyDriver(src_drv);
                return 0;
        }
 
        /* convert GDAL MEM raster to output format */
        RASTER_DEBUG(3, "Copying GDAL MEM raster to memory file in output format");
        rtn_ds = GDALCreateCopy(
                rtn_drv,
                "/vsimem/out.dat", /* should be fine assuming this is in a process */
                src_ds,
                FALSE, /* should copy be strictly equivelent? */
                options, /* format options */
                NULL, /* progress function */
                NULL /* progress data */
        );
 
        /* close source dataset */
        GDALClose(src_ds);
        if (destroy_src_drv) GDALDestroyDriver(src_drv);
        RASTER_DEBUG(3, "Closed GDAL MEM raster");
 
        if (NULL == rtn_ds) {
                rterror("rt_raster_to_gdal: Could not create the output GDAL dataset");
                return 0;
        }
 
        RASTER_DEBUGF(4, "dataset SRS: %s", GDALGetProjectionRef(rtn_ds));
 
        /* close dataset, this also flushes any pending writes */
        GDALClose(rtn_ds);
        RASTER_DEBUG(3, "Closed GDAL output raster");
 
        RASTER_DEBUG(3, "Done copying GDAL MEM raster to memory file in output format");
 
        /* from memory file to buffer */
        RASTER_DEBUG(3, "Copying GDAL memory file to buffer");
        rtn = VSIGetMemFileBuffer("/vsimem/out.dat", &rtn_lenvsi, TRUE);
        RASTER_DEBUG(3, "Done copying GDAL memory file to buffer");
        if (NULL == rtn) {
                rterror("rt_raster_to_gdal: Could not create the output GDAL raster");
                return 0;
        }
 
        rtn_len = (uint64_t) rtn_lenvsi;
        *gdalsize = rtn_len;
 
        return rtn;
}

References FALSE, RASTER_DEBUG, RASTER_DEBUGF, rt_raster_to_gdal_mem(), rt_util_gdal_register_all(), rterror(), and TRUE.

Referenced by RASTER_asGDALRaster(), and test_raster_to_gdal().

Here is the call graph for this function:

Here is the caller graph for this function: