crop: Crop the dataset obtained from 'crest.get_modern_data'

In crestr: A Probabilistic Approach to Reconstruct Past Climates Using Palaeoecological Datasets

Crop the dataset obtained from crest.get_modern_data

Description

Crop the dataset obtained from crest.get_modern_data according to an object of the class SpatialPolygonsDataFrame.

Usage

crop(x, shp)

Arguments

x	A crestObj produced by the crest.get_modern_data function.
shp	A shapefile to crop the data. Data points will be kept if their centroid is within the shape.

Value

An updated version of the crest.get_modern_data.

Examples

## Not run: 
  data(M1)
  ## We want only the data covering Nigeria
  M2 <- M1[M1$COUNTRY == 'Nigeria', ]
  data(reconstr)
  reconstr.cropped <- crop(reconstr, M2)
  data1 <- raster::rasterFromXYZ(reconstr$modelling$climate_space[, 1:3],
                                 crs=raster::crs(M1))
  data2 <- raster::rasterFromXYZ(reconstr.cropped$modelling$climate_space[, 1:3],
                                 crs=raster::crs(M1))
  layout(matrix(c(1,2,3,4), byrow=FALSE, ncol=2), width=1, height=c(0.2, 0.8))
  plot_map_eqearth(data1, brks.pos=seq(13,29,2), colour_scale=TRUE,
                   title='Full dataset', zlim=c(13, 29))
  plot_map_eqearth(data2, brks.pos=seq(13,29,2), colour_scale=TRUE,
                   title='Cropped dataset', zlim=c(13, 29))

## End(Not run)

crestr documentation built on Jan. 6, 2023, 5:23 p.m.