gdalraster 1.11.1

• fix test in test-ogr_manage.R: the test for GeoJSON layer did not need to check existence using with_update = TRUE on a file in extdata (#410)

• add apply_geotransform(): convert raster column/row to geospatial x/y coordinates, wrapper of GDALApplyGeoTransform() in the GDAL API, operating on a matrix of input col/row coordinates (the internal wrapper .apply_geotransform() is unchanged)

• add GDALRaster$apply_geotransform(): class method alternative to calling the stand-alone function apply_geotransform() on an object of class GDALRaster

• vsi_curl_clear_cache(): add parameter quiet to wrap the API call in a quiet error handler, TRUE by default

• Documentation: document the w+ access flag for class VSIFile; add CPL_VSIL_USE_TEMP_FILE_FOR_RANDOM_WRITE configuration option in vignette GDAL Config Quick Reference; replace paste0() with file.path() in the examples throughout

• code linting

gdalraster 1.11.0

System requirements

• GDAL >= 3.1.0 is now required (previously >= 2.4.0)
• package bit64 has been added to Imports
• package RcppInt64 has been added in LinkingTo

New utility functions for managing vector data sources

• initial bindings to the GDAL/OGR Vector API supporting data source management: ogr_ds_exists(), ogr_ds_format(), ogr_ds_test_cap(), ogr_ds_create(), ogr_ds_layer_count(), ogr_ds_layer_names(), ogr_layer_exists(), ogr_layer_test_cap(), ogr_layer_create(), ogr_layer_field_names(), ogr_layer_delete(), ogr_field_index(), ogr_field_create(), ogr_geom_field_create(), ogr_field_rename(), ogr_field_delete(), ogr_execute_sql()
• documentation and helper functions for feature class definition: ogr_def_field(), ogr_def_geom_field(), ogr_def_layer()

New bindings to the GDAL VSIVirtualHandle API

• class VSIFile wraps VSIVirtualHandle for Standard C binary file I/O on regular file systems, URLs, cloud storage services, Zip/GZip/7z/RAR, and in-memory files

GDAL VSI for operations on virtual file systems (bug fix / enhancements)

• bug fix in vsi_mkdir(): the file mode was set incorrectly because mode was not passed correctly as octal literal. mode is now passed as a character string containing the file mode as octal.
• add vsi_get_file_metadata(): returns metadata for network filesystem objects (/vsicurl/, /vsis3/, /vsiaz/, etc.), and with GDAL >= 3.7, /vsizip/ SOZip metadata
• add vsi_set_path_option(): set a path specific option for a given path prefix, e.g., credential setting for a virtual file system (GDAL >= 3.6)
• add vsi_clear_path_options(): clear path specific configuration options previously set with vsi_set_path_option() (GDAL >= 3.6)
• vsi_rmdir(): add argument recursive, TRUE to delete the directory and its content
• vsi_mkdir(): add argument recursive, TRUE to create the directory and its ancestors
• several VSI functions returned 0 or -1 invisibly indicating success/failure, consistent with GDAL return values. Those return values are now visible to be consistent with return values from VSIFile class methods.
• vsi_stat() with info = "size", and vsi_get_disk_free_space() now return bit64::integer64 type

GDALRaster-class

• behavior change: the class methods $info() and $infoAsJSON() now use the default command-line arguments for the underlying gdalinfo utility. Arguments are configurable in the new read/write $infoOptions field, which is an empty vector by default (character(0)).
• add support for I/O of Byte raster as R raw type, and add the setting $readByteAsRaw as a class field (#314, thanks to @mdsumner)
• add read/write fields $infoOptions and $quiet for applying per-object settings
• add an optional constructor to allow specifying whether the dataset is opened in shared mode, TRUE by default
• add method $getActualBlockSize(): retrieve the actual block size for a given block offset
• add method $get_pixel_line(): class method alternative to calling the stand-alone function get_pixel_line() on an object of class GDALRaster(#339)
• add method $getProjection(): equivalent to $getProjectionRef() (consistent with osgeo.gdal.Dataset.getProjection() / osgeo.gdal.Dataset.getProjectionRef() in the GDAL Python API)
• method $getDefaultRAT(): add progress bar since retrieving large raster attribute tables could take >30 sec

Stand-alone processing functions

• calc(): add support for multiband output (#319)
• calc(): add input validation for var.names, must be in expr
• get_pixel_line(): an object of class GDALRaster can now be passed for the gt parameter, in which case the geotransform will be obtained from the object and bounds checking on the raster extent will be done (original behavior for gt as numeric vector is unchanged) (#339)
• ogr2ogr(): add parameter open_options to support options on the source dataset
• read_ds(): add parameter as_raw to read a Byte raster as R raw type (#314, thanks to @mdsumner)

GDAL configuration

• add dump_open_datasets(): dump a list of all open datasets (shared or not) to the console
• add get_num_cpus(): get the number of processors detected by GDAL
• add get_usable_physical_ram(): get usable physical RAM reported by GDAL
• add has_spatialite(): returns TRUE if GDAL was built with SpatiaLite support
• add http_enabled(): returns TRUE if GDAL was built with libcurl support
• add .cpl_http_cleanup(): wrapper of CPLHTTPCleanup() for internal use (2024-05-29)

Geometry functions

• new additional geometry functions operating on WKT (GEOS via GDAL headers): g_is_empty(), g_is_valid(), g_name(), g_intersects(), g_equals(), g_disjoint(), g_touches(), g_contains(), g_within(), g_crosses(),g_overlaps(), g_intersection(), g_union(), g_difference(), g_sym_difference(), g_distance(), g_length(), g_area(), g_centroid()
• add bbox_transform(): transform a bounding box to a different projection
• g_transform(): now uses OGR_GeomTransformer_Create() and OGR_GeomTransformer_Transform() in the GDAL API, enhanced version of OGR_G_Transform(); add arguments wrap_date_line and date_line_offset

Documentation

• add Discussions on the GitHub repository
• add a section for HTTP/HTTPS (/vsicurl/) in the vignette GDAL Config Quick Reference
• DESCRIPTION file: add Michael D. Sumner in Authors@R

Other internal changes and fixes

• fix memory leaks detected by Valgrind in GDALRaster class methods $info(), $infoAsJSON() and $getDefaultRAT()
• register a finalizer to call CPLHTTPCleanup() upon R session exit
• add GDALRaster class method $setFilename(): set the filename of an uninitialized GDALRaster object, currently undocumented / for internal use
• add GDALRaster class method _getGDALDatasetH(): get the GDAL dataset handle for internal use
• buildRAT(): if the input raster is an object of class GDALRaster, use it by reference rather than instantiating another GDALRaster object internally
• calc(): close input raster dataset before exit when a differing extent is detected
• add some missing null checks, and object destruction on error conditions, in src/geos_wkt.cpp
• improve the check for "-json" as a cl_arg to ogrinfo()
• code linting

gdalraster 1.10.0

System requirements

• GDAL built against GEOS is now required

GDAL VSI for operations on virtual file systems

• add vsi_get_fs_prefixes(): get the list of prefixes for virtual file system handlers currently registered
• add vsi_get_fs_options(): get the list of options associated with a virtual file system handler (for setting with set_config_option())
• add vsi_supports_rnd_write() and vsi_supports_seq_write(): test whether the filesystem supports random write or sequential write, conditional on whether a local temp file is allowed before uploading to the target location
• add vsi_get_disk_free_space(): return the free disk space available on the filesystem
• fixed misspelled argument in vsi_copy_file() and vsi_sync() (#233)

Other stand-alone functions

• add ogrinfo(): wrapper of the ogrinfo command-line utility, retrieve information about a vector data source and potentially edit data with SQL statements (GDAL >= 3.7)
• add ogr2ogr(): wrapper of the ogr2ogr command-line utility, convert vector data between different formats
• add g_transform(): apply a coordinate transformation to a WKT geometry
• add geos_version(): get version information for the GEOS library in use by GDAL
• add push_error_handler(): wrapper for CPLPushErrorHandler() in the GDAL Common Portability Library
• add pop_error_handler(): wrapper for CPLPopErrorHandler() in the GDAL Common Portability Library
• calc(): the argument usePixelLonLat is deprecated as unnecessary, variables pixelLon / pixelLat are now auto-detected if used in the calc expression; small performance improvement from computing pixelY only when needed
• add optional argument quiet in several functions to configure progress reporting (#237)
• make the dataset management functions quieter (#282)
• gdal_formats() now returns a data frame with the supported raster and vector formats, and information about the capabilities of each format driver

CmbTable-class

• new(): assign default variable names in the constructor if names are not given

Internal

• src/geos_wkt.cpp, src/transform.cpp, src/wkt_conv.cpp: deallocate some OGR geometry and OSR spatial ref objects to fix memory leaks
• add more unit tests for geometry operations using GEOS via GDAL headers
• GDALRaster::getMetadataDomainList(): deallocate the returned string list to avoid memory leak
• GDALRaster::close(): clear cache if needed, and check the return values of GDALClose() and GDALFlushCache() if GDAL >= 3.7
• configure.ac: add back proj-include and proj-lib, the latter needed in some cases for source install on macOS; rework for the system requirement of GDAL built against GEOS
• remove internal has_geos() checks and update the documentation, since GDAL with GEOS is now required
• add .editorconfig file and bulk reformat code style
• fix up R code for lintr and add .lintr file
• mass replace NULL -> nullptr in C++ code
• format diagnostic messages throughout for consistency and follow guidelines given in "Writing R Extensions"
• clean up temp files in the examples throughout

gdalraster 1.9.0

Behavior change

• remove OSR_DEFAULT_AXIS_MAPPING_STRATEGY=TRADITIONAL_GIS_ORDER from _gdal_init() (#209), since this could give a different SRS comparison result depending on GDAL version in cases where axis mapping matters and the default options in OSRIsSameEx() are used

GDAL VSI for operations on virtual file systems

• add vsi_copy_file(): wrapper for VSICopyFile()
• add vsi_curl_clear_cache(): wrapper for VSICurlClearCache() and VSICurlPartialClearCache()
• add vsi_mkdir(): wrapper for VSIMkdir()
• add vsi_read_dir(): wrapper for VSIReadDirEx()
• add vsi_rename(): wrapper for VSIRename()
• add vsi_rmdir(): wrapper for VSIRmdir()
• add vsi_stat(): wrapper for VSIStatExL()
• add vsi_sync(): wrapper for VSISync()
• add vsi_unlink(): wrapper for VSIUnlink()
• add vsi_unlink_batch(): wrapper for VSIUnlinkBatch()

GDALRaster-class

• GDALRaster$new(): add a constructor for passing dataset open options

Stand-alone functions

• add footprint(): wrapper of the gdal_footprint command-line utility, compute footprint of a raster (GDAL >= 3.8)
• read_ds(): add as_list argument for option to return multi-band output in list form; attach attribute gis to the output, a list containing bbox, dimension and spatial reference (thanks to input from @mdsumner #185)
• plot_raster(): accept pixel data in list form (band vectors as list elements), and make use of gis attribute if present (thanks to input from @mdsumner #185)
• srs_is_same(): add arguments for criterion, ignore_axis_mapping and ignore_coord_epoch

Documentation

• organize the website reference index
• add Microsoft Azure information in GDAL Config Quick Reference
• update DESCRIPTION

Internal

• add _check_gdal_filename(): minimal filename check and UTF-8
• add _get_physical_RAM(): wrapper for CPLGetUsablePhysicalRAM() for internal use
• set test coverage target minimum to 90%, codecov range: 70..90 (previously 80%, 60..80)

gdalraster 1.8.0

Dependencies

• libxml2 is listed in SystemRequirements since it is required by the imported package xml2

GDALRaster-class

• read_only=TRUE is now an optional default in the class constructor
• add clearStatistics(): clear statistics on PAM supported datasets (GDAL >= 3.2)

Stand-alone functions

• add addFilesInZip(): create/append to Seek-Optimized ZIP files (GDAL >= 3.7)
• plot_raster(): now uses georeferenced coordinates by default instead of pixel/line (#184 thanks to @mdsumner)

Documentation

• example code for calc() using Landsat bands should have applied scale/offset as given in the .json metadata - this has been corrected
• add AWS_REGION, and a section for SOZip to GDAL Config Quick Ref
• update documentation and examples throughout for default read_only=TRUE in the constructor for GDALRaster
• add translate() and warp() in Raster API Tutorial

Internal

• configure.ac: remove parts of PROJ config that were unneeded since gdalraster only uses PROJ via GDAL headers; use autoconf macros for compiler invocations; get the PROJ data directory from pkg-config if possible; add action-if-cross-compiling argument (#190 and #197; thanks to Simon Urbanek for debugging an initial issue and providing helpful feedback)
• get GDAL libs from pkg-config on Windows when possible (#125 thanks to @kalibera)
• translate() and warp(): close the output dataset before closing source dataset(s) in case the output is VRT
• switch to using the R convention for x.x.x.9000 development versioning

gdalraster 1.7.0

GDAL API stand-alone functions

• add polygonize(): create a polygon feature layer from raster data, wrapper for GDALPolygonize in the GDAL Algorithms API.
• add rasterize(): burn vector geometries (points, lines, or polygons) into a raster, wrapper for the gdal_rasterize command-line utility
• add buildVRT(): build a GDAL virtual raster mosaic from a list of datasets, wrapper for the gdalbuildvrt command-line utility
• add translate(): convert raster data between different formats, wrapper for the gdal_translate command-line utility
• make t_srs optional in warp()
• the GDAL configuration option TRADITIONAL_GIS_ORDER=OSR_DEFAULT_AXIS_MAPPING_STRATEGY is now set on package load

Documentation

• add GDAL Config Quick Reference to vignettes
• warp() has additional documentation covering several processing options
• add the COMPRESS_OVERVIEW configuration option to the documentation for GDALRaster::buildOverviews()
• add str() of the GDALRaster object in the Raster API Tutorial
• add display of the EVT raster itself along with its attribute table in Raster Attribute Tables

Internal

• add src/ogr_util.cpp and src/ogr_util.h: OGR utility functions for vector data sources
• src/gdal_exp.cpp: additional error handling in several wrapper functions, increased test coverage

gdalraster 1.6.0

Dependencies

• GDAL >= 2.4.0 (previously >= 2.3.0)
• package xml2 is now required (previously in Suggests)
• package gt has been added to Suggests (required for displayRAT())

GDALRaster-class

• add methods for Raster Attribute Tables: getDefaultRAT(), setDefaultRAT()
• add getDefaultHistogram(): fetch default raster histogram for a band
• add getHistogram(): compute raster histogram for a band
• add getMinMax(): compute min/max for a raster band
• add getMetadataDomainList(): get a list of metadata domains for a dataset or raster band
• fix getMetadataItem() for a specific domain at dataset level (#109)

GDAL API stand-alone functions

• add buildRAT(): compute for a raster band the set of unique pixel values and their counts, and build a GDAL Raster Attribute Table as data frame
• add displayRAT(): generate a presentation Raster Attribute Table, showing colors if the table contains RGB columns
• add gdal_formats(): report the supported raster formats
• add getCreationOptions(): get the list of creation options of a raster format
• add copyDatasetFiles(): copy all the files associated with a dataset
• add deleteDataset(): delete a dataset in a format-specific way
• add renameDataset(): rename a dataset in a format-specific way
• add some missing error checks in src/gdal_exp.cpp (#104)

CmbTable-class

• argument incr for the count increment in CmbTable::update() can be zero

RunningStats-class

• use uint64_t for the count accumulator (previously long long) and make explicit the return cast in get_count() (no user-visible changes)
• slightly faster update

Other miscellaneous

• plot_raster(): normalize legend correctly for minmax_def and minmax_pct_cut (#131)

Documentation

• add vignette Raster Attribute Tables
• add notes for RunningStats-class
• update vignette Raster API Tutorial with gdal_formats() and getCreationOptions()
• update installation instructions in README

gdalraster 1.5.0

GDALRaster-class

• add methods for color tables: getColorTable(), getPaletteInterp(), setColorTable()
• add getRasterColorInterp(), setRasterColorInterp(): get/set color interpretation for raster bands
• add getDescription(), setDescription(): get/set description for raster band objects
• add flushCache(): flush all write cached data to disk
• add getFileList(): returns a list of files forming the dataset
• add infoAsJSON(): returns output of the gdalinfo command-line utility as a JSON-formatted string
• new(): add a warning in the class constructor if the raster has an int64 data type (would be handled as double for now)

GDAL API stand-alone functions

• add bandCopyWholeRaster(): wrapper for GDALRasterBandCopyWholeRaster(), efficiently copy a whole raster band
• add createColorRamp(): wrapper for GDALCreateColorRamp(), automatically create a ramp from one color to another
• add sieveFilter(): wrapper for GDALSieveFilter() in the Algorithms API, remove small raster polygons
• add PROJ utility functions: proj_version(), proj_search_paths(), proj_networking() (via GDAL headers)
• add g_buffer(): compute buffer of a WKT geometry (GEOS convenience function via GDAL headers)

CmbTable-class

• add updateFromMatrixByRow(): update the hash table from a matrix having integer combinations arranged in rows
• add asMatrix(): return the combinations table as a numeric matrix (alternative to asDataFrame())

Bug fix

• warp() caused segfault if proj.db could not be found (#96)

Other miscellaneous

• plot_raster(): default value of the legend argument has been changed to FALSE; legend can now use a color table for continuous data; add argument maxColorValue (e.g., to use RGB 0:255 instead of 0:1 in col_tbl)
• bbox_from_wkt(), bbox_to_wkt(): add arguments extend_x, extend_y
• on Windows, reset GDAL environment variables on package unload if they were previously set on load
• add inst/extdata/storml_tcc.tif: example dataset of NLCD Tree Canopy Cover
• update the package vignette for color tables and raster display

gdalraster 1.4.0

• add dem_proc(): wrapper for the gdaldem command-line utility to generate DEM derivatives

• add the following set methods in class GDALRaster: setMetadataItem(), setUnitType(), setScale(), setOffset()

• add GDALRaster$buildOverviews(): build raster overviews

• add GDALRaster$dim(): returns a vector of xsize, ysize, nbands

• transform_xy() and inv_project(): pts can be a data frame or matrix

• plot_raster() now accepts a GDALRaster object for the data argument

• plot_raster(): make the legend narrower and add argument digits to format legend labels when raster data are floating point

• add test suite and code coverage report

gdalraster 1.3.0

• GDALRaster::read(): data are now read as R integer type when possible for the raster data type (#23)

• add fillNodata(): wrapper for GDALFillNodata() in the GDAL Algorithms API

• add read_ds(): convenience wrapper for GDALRaster$read()

• add plot_raster(): display raster data using base R graphics

• add get_cache_used(): wrapper for GDALGetCacheUsed64() with return value in MB

• add GDALRaster$getOverviewCount(): return the number of overview layers available

• GDALRaster$info(): drop -nomd argument from the internal call

• bbox_from_wkt(): return NA if creation of the geometry object fails (#27)

• fix GDALRaster$getMetadata(): requesting band-level metadata for a domain other than the default metadata domain was returning dataset-level metadata instead

• add vignette containing an R port of the GDAL Raster API tutorial

• add description of the GDAL_RASTERIO_RESAMPLING configuration option in the documentation for GDALRaster$read()

• add web article on the GDAL block cache and configuration of GDAL_CACHEMAX