gdalraster is an R interface to the Raster and Vector APIs of the Geospatial Data Abstraction Library (GDAL).
API-level bindings are implemented in the exposed C++ classes
GDALRaster
and
GDALVector,
along with several stand-alone
functions.
Bindings to the GDAL Virtual Systems Interface
(VSI) are also
included to support file system operations and binary I/O on URLs, cloud
storage services, Zip/GZip/7z/RAR, and in-memory files, as well as
regular file systems. Calling signatures resemble the native C, C++ and
Python APIs provided by the GDAL project. Supports:
- manual creation of uninitialized raster and vector datasets
- vector layer creation and schema management
- read/set parameters and metadata for raster bands and vector layers
- low-level I/O
- build/read/set color tables and raster attribute tables
- virtual raster (VRT) for virtual cropping, resampling, kernel filtering, mosaicing
gdalwarpwrapper for reproject/resample/crop/mosaic- coordinate transformation
- spatial reference systems
- geometry operations on raw vectors of WKB, or WKT strings
- GDAL algorithms (
dem_proc(),polygonize(),rasterize(),...) - OGR vector utilities (
ogrinfo(),ogr2ogr(),ogr_manageinterface) - OGR facilities for vector geoprocessing
(
ogr_proc()) - raster and vector dataset management (copy files/move/rename/delete)
- create/append to Seek-Optimized ZIP (SOZip)
- abstraction of file system operations on URLs, cloud storage, in-memory files, etc.
- Standard C binary file I/O through VSI (class
VSIFile)
Additional functionality includes:
- class
RunningStatscalculates mean and variance in one pass, and tracks the min, max, sum, and count (i.e., summary statistics on a data stream). The input data values are not stored in memory, so this class can be used to compute statistics for very large data streams. - class
CmbTableidentifies and counts unique combinations of integer values using a hash table. combine()overlays multiple rasters so that a unique ID is assigned to each unique combination of input values. Pixel counts for each unique combination are obtained, and combination IDs are optionally written to an output raster.calc()evaluates an R expression for each pixel in a raster layer or stack of layers. Individual pixel coordinates are available as variables in the R expression, as either x/y in the raster projected coordinate system or inverse projected longitude/latitude.plot_raster()displays raster data using base R graphics.
gdalraster may be useful in applications that need scalable, low-level I/O, or prefer a direct GDAL API. Comprehensive documentation is provided in the package and online.
Install the released version from CRAN with:
install.packages("gdalraster")CRAN provides pre-compiled binary packages for Windows and macOS. These do not require any separate installation of external libraries for GDAL.
GDAL >= 3.1.0 built with GEOS is required, but a more recent version is recommended.
On Ubuntu, recent versions of geospatial libraries can be installed from the ubuntugis-unstable PPA with the following commands. Note that libxml2 is required for the R package xml2 which is a dependency, so we go ahead and install it here as well:
sudo add-apt-repository ppa:ubuntugis/ubuntugis-unstable
sudo apt update
sudo apt install libgdal-dev libgeos-dev libxml2-dev
The versions in ubuntugis-unstable generally work well and are more up-to-date, but less recent versions in the ubuntugis-stable PPA could be used instead.
Package sf provides helpful instructions for installing the geospatial libraries on other Linux distributions.
With the dependent libraries available on the system, install from CRAN:
install.packages("gdalraster")Or install the development version from GitHub using package remotes:
remotes::install_github("USDAForestService/gdalraster")RTools is needed to install from source on Windows. RTools since version 4.2 includes GDAL and all other dependent libraries that are needed to compile gdalraster. Note that CRAN releases periodic revisions to RTools that often include updates to the libraries as new versions become available. Release 6536 of RTools 4.5 contains GDAL 3.10.2, GEOS 3.13.1 and PROJ 9.5.1.
With RTools installed:
# Install the development version from GitHub
remotes::install_github("USDAForestService/gdalraster")GDAL can be installed with Homebrew:
brew install gdal
then
# Install the development version from GitHub
remotes::install_github("USDAForestService/gdalraster")Caution is warranted on macOS with regard to mixing a source installation with installation of binaries from CRAN. Consider installing the development version from R-universe instead.
R-universe
provides pre-compiled binary packages for Windows and macOS that track
the development version of gdalraster. New packages are built
usually within ~1 hour of the most recent commit in branch main.
# Install the development version from r-universe
install.packages("gdalraster", repos = c("https://usdaforestservice.r-universe.dev", "https://cran.r-project.org"))