R: Native interface to gdal utils

gdal_utils {sf} | R Documentation

## Native interface to gdal utils

### Description

Native interface to gdal utils

### Usage

    gdal_utils(
      util = "info",
      source,
      destination,
      options = character(0),
      quiet = !(util %in% c("info", "mdiminfo")),
      processing = character(0),
      colorfilename = character(0)
    )

### Arguments

`util` | character; one of `info`, `warp`, `rasterize`, `translate`, `vectortranslate` (for ogr2ogr), `buildvrt`, `demprocessing`, `nearblack`, `grid`, `mdiminfo` and `mdimtranslate` (the last two requiring GDAL 3.1)
`source` | character; name of input layer(s); for `warp`, `buidvrt` or `mdimtranslate` this can be more than one
`destination` | character; name of output layer
`options` | character; options for the utility
`quiet` | logical; if `TRUE`, suppress printing the output for `info` and `mdiminfo`, and suppress printing progress
`processing` | character; processing options for `demprocessing`
`colorfilename` | character; name of color file for `demprocessing` (mandatory if `processing="color-relief"`)

### Value

`info` returns a character vector with the raster metadata; all other utils return (invisibly) a logical indicating success (i.e., `TRUE`); in case of failure, an error is raised.

### Examples

    if (sf_extSoftVersion()["GDAL"] > "2.1.0") {
    # info utils can be used to list information about about a raster
    # dataset. More info: https://gdal.org/programs/gdalinfo.html
    in_file <- system.file("tif/geomatrix.tif", package = "sf")
    gdal_utils("info", in_file, options = c("-mm", "-proj4"))

    # vectortranslate utils can be used to convert simple features data between
    # file formats. More info: https://gdal.org/programs/ogr2ogr.html
    in_file <- system.file("shape/storms_xyz.shp", package="sf")
    out_file <- paste0(tempfile(), ".gpkg")
    gdal_utils(
      util = "vectortranslate",
      source = in_file,
      destination = out_file, # output format must be specified for GDAL < 2.3
      options = c("-f", "GPKG")
    )
    # The parameters can be specified as c("name") or c("name", "value"). The
    # vectortranslate utils can perform also various operations during the
    # conversion process. For example we can reproject the features during the
    # translation.
    gdal_utils(
      util = "vectortranslate",
      source = in_file,
      destination = out_file,
      options = c(
      "-f", "GPKG", # output file format for GDAL < 2.3
      "-s_srs", "EPSG:4326", # input file SRS
      "-t_srs", "EPSG:2264", # output file SRS
      "-overwrite"
      )
    )
    st_read(out_file)
    # The parameter s_srs had to be specified because, in this case, the in_file
    # has no associated SRS.
    st_read(in_file)
    }

[Package sf version 1.0-7 Index]