Actually tests for the intersection of the extent of each raster

test_ras_xn(a, b)

Arguments

a, b

raster as SpatRaster or path to raster

Value

Logical

Examples

  #library(envRaster)

  out_dir <- file.path(system.file(package = "envRaster"), "examples")

  base_file <- fs::path(out_dir, "cube", "base.tif")

  if(!file.exists(base_file)) {

    make_base_grid(aoi
                   , out_res = 30
                   , out_epsg = 8059
                   , out_file = base_file
                   , datatype = "INT1U"
                   , overwrite = TRUE
                   )
  }

  base <- terra::rast(base_file)

  test_ras_xn(base, base)
#> [1] TRUE

  test_ras_xn(base
              , system.file("ex/meuse.tif", package="terra")
              )
#> Warning: [intersect] no intersection
#> [1] FALSE