subset.layer.Rd
Function to subset and clip a GIS vector or raster layer by a bounding box or clip/subset layer (e.g. shapefile/KML)
# S3 method for layer subset(file, coords = NULL, sub.layer = NULL, clip = FALSE, out.dir = NULL, out.name = NULL)
file | input file to be subset |
---|---|
coords | vector with xmin,ymin,xmax,ymax defing the bounding box for subset |
sub.layer | Vector layer defining the subset region |
clip | clip geometries to bounding box/subset layer? TRUE/FALSE |
out.dir | output directory for subset layer. Defaults to location of input file. Can also set to 'pwd' |
out.name | filename for subset layer. Defaults to original filename with the suffix *.sub |
# NOT RUN { # Test dataset file <- Sys.glob(file.path(R.home(), 'library', 'PEcAn.data.land','data','*.shp')) out.dir <- path.expand('~/temp') # with clipping enabled subset.layer(file=file,coords=c(-95,42,-84,47),clip=TRUE,out.dir=out.dir) # without clipping enables subset.layer(file=file,coords=c(-95,42,-84,47),out.dir=out.dir) system(paste('rm -r',out.dir,sep='')) # }# NOT RUN { # }