Nothing
library(testthat)
library(stars)
context("gdal utils")
test_that('gdal_utils work', {
skip_on_appveyor()
fname = system.file("nc/tos_O1_2001-2002.nc", package = "stars")
if (fname != "") {
info = gdal_utils("info", fname, quiet = TRUE)
sd2 = gdal_subdatasets(fname)[[4]]
info = gdal_utils("info", sd2, quiet = TRUE)
tf = tempfile()
tf2 = tempfile()
tf3 = tempfile()
#tf = "foo"
#gdal_utils("rasterize", points, tif) -> need a good example
expect_true(gdal_utils("warp", sd2, tf, c("-t_srs", "+proj=utm +zone=11 +datum=WGS84")))
# expect_true(gdal_utils("rasterize", sd2, tf)) # breaks on fedora-clang
expect_true(gdal_utils("translate", sd2, tf))
expect_true(gdal_utils("vectortranslate", sd2, tf2))
expect_warning(gdal_utils("nearblack", sd2, tf))
# create point geom:
expect_true(gdal_utils("buildvrt", sd2, tf3))
# expect_true(gdal_utils("demprocessing", sd2, tf, processing = "hillshade")) # breaks on fedora-clang
}
tf = tempfile()
points = system.file("gpkg/nc.gpkg", package="sf")
expect_true(gdal_utils("grid", points, tf))
})
# gdalwarp -t_srs '+proj=utm +zone=11 +datum=WGS84' -overwrite NETCDF:avhrr-only-v2.19810901.nc:anom utm11.tif
# becomes:
# gdal_utils("warp" "NETCDF:avhrr-only-v2.19810901.nc:anom", "utm11.tif", c("-t_srs", "+proj=utm +zone=11 +datum=WGS84"))
Any scripts or data that you put into this service are public.
Add the following code to your website.
For more information on customizing the embed code, read Embedding Snippets.