This function can be used to crop or reproject module inputs from raw data.
cropInputs( x, studyArea, rasterToMatch, verbose = getOption("reproducible.verbose", 1), ... ) # S3 method for default cropInputs(x, studyArea, rasterToMatch, ...) # S3 method for spatialClasses cropInputs( x, studyArea = NULL, rasterToMatch = NULL, verbose = getOption("reproducible.verbose", 1), extentToMatch = NULL, extentCRS = NULL, useGDAL = getOption("reproducible.useGDAL", TRUE), ... ) # S3 method for sf cropInputs( x, studyArea = NULL, rasterToMatch = NULL, verbose = getOption("reproducible.verbose", 1), extentToMatch = NULL, extentCRS = NULL, ... )
x | A |
---|---|
studyArea |
|
rasterToMatch | Template |
verbose | Numeric, -1 silent (where possible), 0 being very quiet,
1 showing more messaging, 2 being more messaging, etc.
Default is 1. Above 3 will output much more information about the internals of
Caching, which may help diagnose Caching challenges. Can set globally with an
option, e.g., |
... | Passed to raster::crop |
extentToMatch | Optional. Can pass an extent here and a |
extentCRS | Optional. Can pass a |
useGDAL | Logical or |
# Add a study area to Crop and Mask to # Create a "study area" library(sp) library(raster) ow <- setwd(tempdir()) # make a SpatialPolygon coords1 <- structure(c(-123.98, -117.1, -80.2, -100, -123.98, 60.9, 67.73, 65.58, 51.79, 60.9), .Dim = c(5L, 2L)) Sr1 <- Polygon(coords1) Srs1 <- Polygons(list(Sr1), "s1") shpEcozone <- SpatialPolygons(list(Srs1), 1L) crs(shpEcozone) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" # make a "study area" that is subset of larger dataset coords <- structure(c(-118.98, -116.1, -99.2, -106, -118.98, 59.9, 65.73, 63.58, 54.79, 59.9), .Dim = c(5L, 2L)) Sr1 <- Polygon(coords) Srs1 <- Polygons(list(Sr1), "s1") StudyArea <- SpatialPolygons(list(Srs1), 1L) crs(StudyArea) <- "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0" ########## shpEcozonePostProcessed <- postProcess(shpEcozone, studyArea = StudyArea)#>#> cropping ...#> useGDAL is TRUE, but problem is small enough for RAM; skipping GDAL; useGDAL = 'force' to override#> Checking for errors in SpatialPolygon#> Found no errors.#>#> Checking for errors in SpatialPolygon#> Found no errors.#> No cacheRepo supplied and getOption('reproducible.cachePath') is inside a temporary directory; #>#> reprojecting ...#> Checking for errors in SpatialPolygon#> Found no errors.#> No cacheRepo supplied and getOption('reproducible.cachePath') is inside a temporary directory; #>#> Checking for errors in SimpleFeature#> Found no errors.#> maskInputs with sf class objects is still experimental#> intersecting ...#> Checking for errors in SimpleFeature#> Found no errors.#>#>#> Skipping writeOutputs; filename2 is NULL#> reprojecting ...shpEcozoneCropped <- cropInputs(shpEcozone, StudyArea)#> cropping ...#> useGDAL is TRUE, but problem is small enough for RAM; skipping GDAL; useGDAL = 'force' to override#> Checking for errors in SpatialPolygon#> Found no errors.#>#> Checking for errors in SpatialPolygon#> Found no errors.#> Checking for errors in SimpleFeature#> Found no errors.#> maskInputs with sf class objects is still experimental#> intersecting ...#> Checking for errors in SimpleFeature#> Found no errors.#>#>setwd(ow)