Return a raster from a GDAL dataset.
2169 uint32_t height = 0;
2170 uint32_t numBands = 0;
2172 char *authname = NULL;
2173 char *authcode = NULL;
2175 GDALRasterBandH gdband = NULL;
2176 GDALDataType gdpixtype = GDT_Unknown;
2187 uint32_t nXBlocks, nYBlocks;
2188 int nXBlockSize, nYBlockSize;
2189 uint32_t iXBlock, iYBlock;
2190 uint32_t nXValid, nYValid;
2193 uint8_t *values = NULL;
2194 uint32_t valueslen = 0;
2195 uint8_t *ptr = NULL;
2200 width = GDALGetRasterXSize(
ds);
2201 height = GDALGetRasterYSize(
ds);
2202 RASTER_DEBUGF(3,
"Raster dimensions (width x height): %d x %d", width, height);
2208 rterror(
"rt_raster_from_gdal_dataset: Out of memory allocating new raster");
2214 cplerr = GDALGetGeoTransform(
ds,
gt);
2215 if (GDALGetGeoTransform(
ds,
gt) != CE_None) {
2216 RASTER_DEBUG(4,
"Using default geotransform matrix (0, 1, 0, 0, 0, -1)");
2228 RASTER_DEBUGF(3,
"Raster geotransform (%f, %f, %f, %f, %f, %f)",
2235 strcmp(authname,
"EPSG") == 0 &&
2242 if (authname != NULL)
2244 if (authcode != NULL)
2248 numBands = GDALGetRasterCount(
ds);
2250 #if POSTGIS_DEBUG_LEVEL > 3
2251 for (i = 1; i <= numBands; i++) {
2252 GDALRasterBandH _grb = NULL;
2258 _grb = GDALGetRasterBand(
ds, i);
2259 GDALComputeRasterStatistics(_grb,
FALSE, &_min, &_max, &_mean, &_stddev, NULL, NULL);
2260 RASTER_DEBUGF(4,
"GDAL Band %d stats: %f, %f, %f, %f", i, _min, _max, _mean, _stddev);
2265 for (i = 1; i <= numBands; i++) {
2268 gdband = GDALGetRasterBand(
ds, i);
2269 if (NULL == gdband) {
2270 rterror(
"rt_raster_from_gdal_dataset: Could not get GDAL band");
2277 gdpixtype = GDALGetRasterDataType(gdband);
2278 RASTER_DEBUGF(4,
"gdpixtype, size = %s, %d", GDALGetDataTypeName(gdpixtype), GDALGetDataTypeSize(gdpixtype) / 8);
2281 rterror(
"rt_raster_from_gdal_dataset: Unknown pixel type for GDAL band");
2288 width = GDALGetRasterBandXSize(gdband);
2289 height = GDALGetRasterBandYSize(gdband);
2290 RASTER_DEBUGF(3,
"GDAL band dimensions (width x height): %d x %d", width, height);
2293 nodataval = GDALGetRasterNoDataValue(gdband, &hasnodata);
2294 RASTER_DEBUGF(3,
"(hasnodata, nodataval) = (%d, %f)", hasnodata, nodataval);
2299 (hasnodata ? nodataval : 0),
2303 rterror(
"rt_raster_from_gdal_dataset: Could not allocate memory for raster band");
2311 GDALGetBlockSize(gdband, &nXBlockSize, &nYBlockSize);
2312 nXBlocks = (width + nXBlockSize - 1) / nXBlockSize;
2313 nYBlocks = (height + nYBlockSize - 1) / nYBlockSize;
2314 RASTER_DEBUGF(4,
"(nXBlockSize, nYBlockSize) = (%d, %d)", nXBlockSize, nYBlockSize);
2315 RASTER_DEBUGF(4,
"(nXBlocks, nYBlocks) = (%d, %d)", nXBlocks, nYBlocks);
2318 valueslen = ptlen * nXBlockSize * nYBlockSize;
2320 if (values == NULL) {
2321 rterror(
"rt_raster_from_gdal_dataset: Could not allocate memory for GDAL band pixel values");
2325 RASTER_DEBUGF(3,
"values @ %p of length = %d", values, valueslen);
2327 for (iYBlock = 0; iYBlock < nYBlocks; iYBlock++) {
2328 for (iXBlock = 0; iXBlock < nXBlocks; iXBlock++) {
2329 x = iXBlock * nXBlockSize;
2330 y = iYBlock * nYBlockSize;
2331 RASTER_DEBUGF(4,
"(iXBlock, iYBlock) = (%d, %d)", iXBlock, iYBlock);
2334 memset(values, 0, valueslen);
2337 if ((iXBlock + 1) * nXBlockSize > width)
2338 nXValid = width - (iXBlock * nXBlockSize);
2340 nXValid = nXBlockSize;
2343 if ((iYBlock + 1) * nYBlockSize > height)
2344 nYValid = height - (iYBlock * nYBlockSize);
2346 nYValid = nYBlockSize;
2348 RASTER_DEBUGF(4,
"(nXValid, nYValid) = (%d, %d)", nXValid, nYValid);
2350 cplerr = GDALRasterIO(
2354 values, nXValid, nYValid,
2358 if (cplerr != CE_None) {
2359 rterror(
"rt_raster_from_gdal_dataset: Could not get data from GDAL raster");
2366 if (nXBlocks == 1 && nYBlockSize > 1 && nXValid == width) {
2368 y = nYBlockSize * iYBlock;
2370 RASTER_DEBUGF(4,
"Setting set of pixel lines at (%d, %d) for %d pixels",
x,
y, nXValid * nYValid);
2375 x = nXBlockSize * iXBlock;
2376 for (iY = 0; iY < nYValid; iY++) {
2377 y = iY + (nYBlockSize * iYBlock);
2379 RASTER_DEBUGF(4,
"Setting pixel line at (%d, %d) for %d pixels",
x,
y, nXValid);
2381 ptr += (nXValid * ptlen);
void rterror(const char *fmt,...)
Wrappers used for reporting errors and info.
void * rtalloc(size_t size)
Wrappers used for managing memory.
rt_pixtype rt_util_gdal_datatype_to_pixtype(GDALDataType gdt)
Convert GDALDataType to rt_pixtype.
#define RASTER_DEBUG(level, msg)
#define RASTER_DEBUGF(level, msg,...)
rt_errorstate rt_band_set_pixel_line(rt_band band, int x, int y, void *vals, uint32_t len)
Set values of multiple pixels.
void rtdealloc(void *mem)
rt_errorstate rt_util_gdal_sr_auth_info(GDALDatasetH hds, char **authname, char **authcode)
Get auth name and code.
int rt_pixtype_size(rt_pixtype pixtype)
Return size in bytes of a value in the given pixtype.
int rt_raster_generate_new_band(rt_raster raster, rt_pixtype pixtype, double initialvalue, uint32_t hasnodata, double nodatavalue, int index)
Generate a new inline band and add it to a raster.
void rt_raster_set_geotransform_matrix(rt_raster raster, double *gt)
Set raster's geotransform using 6-element array.
void rt_raster_destroy(rt_raster raster)
Release memory associated to a raster.
rt_raster rt_raster_new(uint32_t width, uint32_t height)
Construct a raster with given dimensions.
uint16_t rt_raster_get_num_bands(rt_raster raster)
rt_band rt_raster_get_band(rt_raster raster, int n)
Return Nth band, or NULL if unavailable.
void rt_raster_set_srid(rt_raster raster, int32_t srid)
Set raster's SRID.