Nothing
test_that("prepInputs correctly unzips large files", {
skip_on_cran()
skip_on_ci()
skip_if_not(isInteractive(), "tests extracting large files should be run manually devtools::test()")
skip_if_not(getOption("reproducible.runLargeFileTests"))
## based on #145. extracted file is ~30 GB so this takes a long time to test!
testInit("terra")
# tmpdir <- "/mnt/d/temp" # need a drive that is large enough
# if (!"emcintir" %in% Sys.info()["user"] || (!dir.exists(tmpdir)))
# skip("This requires a lot of drive space")
targFile <- "CA_harvest_year_1985_2015.tif"
on.exit(
{
unlink(file.path(tmpdir, targFile), recursive = TRUE)
},
add = TRUE
)
url <- "https://opendata.nfis.org/downloads/forest_change/CA_forest_harvest_mask_year_1985_2015.zip"
options(reproducible.tempPath = file.path(tmpdir, "ttt"))
ff <- prepInputs(
url = url,
targetFile = targFile,
destinationPath = asPath(tmpdir),
fun = "terra::rast",
userTags = c("objectName:forestHarvestMask", "goal:posthocGIS")
)
fout <- file.path(tmpdir, "CA_harvest_year_1985_2015.tif")
expect_true(identical(normPath(Filenames(ff)), fout))
expect_true(file.info(fout)[["size"]] > 28 * 1024^3)
})
test_that("Issue 181 geodatabase file", {
skip_on_cran()
skip_on_ci()
skip_if_not(isInteractive(), "test #2: extracting large files should be run manually devtools::test()")
skip_if_not(getOption("reproducible.runLargeFileTests"))
## based on #145. extracted file is ~30 GB so this takes a long time to test!
testInit(c("terra", "googledrive"), needGoogleDriveAuth = TRUE)
rstLCC <- Cache(prepInputs,
targetFile = "EOSD_Mosaic.gdb",
archive = "EOSD_2000_2007_combined.zip",
alsoExtract = "similar",
url = "https://drive.google.com/file/d/1p66_P6dNdlrvAF3Mp99Xz9Bdz2lvfaQ7",
destinationPath = tmpdir,
filename2 = NULL,
fun = NA,
userTags = c(
"outFun:Cache",
"step:prepEOSD"
)
)
expect_true(is(sf::st_read(rstLCC$targetFilePath, layer = "EOSD_Mosaic_BWC_range_clip", quiet = TRUE), "sf"))
})
test_that("Issue 242 masking fail", {
skip_on_cran()
skip_on_ci()
testInit("terra", needInternet = TRUE)
# skip_if_not(isInteractive(), "test #3: extracting large files should be run manually with devtools::test()")
testInit(c("terra", "googledrive"), needGoogleDriveAuth = FALSE)
studyArea <- vect(structure(
c(
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1,
1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, 1, -1418427,
-1415931, -1398140, -1393331, -1382857, -1371358, -1365393, -1350548,
-1325832, -1323489, -1322689, -1321448, -1319605, -1318493, -1309216,
-1307387, -1305766, -1305267, -1303130, -1295515, -1296061, -1296486,
-1298156, -1299017, -1300951, -1300981, -1301682, -1302736, -1317683,
-1320782, -1322070, -1322341, -1324137, -1335839, -1336905, -1346232,
-1401188, -1422500, -1423243, -1423419, -1418427, 6990085, 6999544,
7008667, 7021910, 7028787, 7029407, 7029780, 7035153, 7036928,
7039810, 7050226, 7044431, 7035561, 7027977, 7027582, 7026342,
7021796, 7006494, 7006125, 6994606, 6994507, 6984394, 6975963,
6975848, 6968927, 6961588, 6957530, 6933733, 6933100, 6926385,
6925727, 6929217, 6932967, 6947379, 6948395, 6954427, 6980171,
6980300, 6987758, 6989296, 6990085, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0,
0, 0, 0, 0, 0, 0, 0, 0, 0, 0, 0
),
dim = c(41L, 5L),
dimnames = list(NULL, c("geom", "part", "x", "y", "hole"))
))
crs(studyArea) <- "PROJCRS[\"unknown\",\nBASEGEOGCRS[\"unknown\",\nDATUM[\"North American Datum 1983\",\nELLIPSOID[\"GRS 1980\",6378137,298.257222101,\nLENGTHUNIT[\"metre\",1]],\nID[\"EPSG\",6269]],\nPRIMEM[\"Greenwich\",0,\nANGLEUNIT[\"degree\",0.0174532925199433],\nID[\"EPSG\",8901]]],\nCONVERSION[\"unknown\",\nMETHOD[\"Lambert Conic Conformal (2SP)\",\nID[\"EPSG\",9802]],\nPARAMETER[\"Latitude of false origin\",0,\nANGLEUNIT[\"degree\",0.0174532925199433],\nID[\"EPSG\",8821]],\nPARAMETER[\"Longitude of false origin\",-95,\nANGLEUNIT[\"degree\",0.0174532925199433],\nID[\"EPSG\",8822]],\nPARAMETER[\"Latitude of 1st standard parallel\",49,\nANGLEUNIT[\"degree\",0.0174532925199433],\nID[\"EPSG\",8823]],\nPARAMETER[\"Latitude of 2nd standard parallel\",77,\nANGLEUNIT[\"degree\",0.0174532925199433],\nID[\"EPSG\",8824]],\nPARAMETER[\"Easting at false origin\",0,\nLENGTHUNIT[\"metre\",1],\nID[\"EPSG\",8826]],\nPARAMETER[\"Northing at false origin\",0,\nLENGTHUNIT[\"metre\",1],\nID[\"EPSG\",8827]]],\nCS[Cartesian,2],\nAXIS[\"(E)\",east,\nORDER[1],\nLENGTHUNIT[\"metre\",1,\nID[\"EPSG\",9001]]],\nAXIS[\"(N)\",north,\nORDER[2],\nLENGTHUNIT[\"metre\",1,\nID[\"EPSG\",9001]]]]"
studyAreaRas <- rasterize(
studyArea,
rast(extent = ext(studyArea), crs = crs(studyArea), resolution = 250)
)
sppAbundURL <- "https://ftp.maps.canada.ca/pub/nrcan_rncan/Forests_Foret/canada-forests-attributes_attributs-forests-canada/2001-attributes_attributs-2001/NFI_MODIS250m_2001_kNN_Species_Pinu_Alb_v1.tif"
a <- capture.output(
sppAbundance <- prepInputs(
targetFile = "NFI_MODIS250m_2001_kNN_Species_Pinu_Alb_v1.tif",
url = sppAbundURL,
destinationPath = tempdir2(),
fun = "terra::rast",
to = studyAreaRas,
overwrite = TRUE
)
)
expect_true(is(sppAbundance, "SpatRaster"))
expect_true(all.equal(ext(studyAreaRas), ext(sppAbundance)))
})
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.