This function takes a .csv file with gridded lat, long data and aggregates the data by spatial boundaries given different shapefiles.

metis.boundaries(
  boundaryRegShape = NULL,
  boundaryRegShpFolder = NULL,
  boundaryRegShpFile = NULL,
  boundaryRegCol = NULL,
  boundaryRegionsSelect = NULL,
  subRegShape = NULL,
  subRegShpFolder = NULL,
  subRegShpFile = NULL,
  subRegCol = NULL,
  subRegCol1 = NULL,
  subRegionsSelect = NULL,
  subRegType = "subRegType",
  dirOutputs = paste(getwd(), "/outputs", sep = ""),
  regionName = NULL,
  folderName = NULL,
  nameAppend = "",
  expandPercentWidth = 2,
  expandPercentHeight = 2,
  expandPercentWidthOV = 2,
  expandPercentHeightOV = 2,
  overlapShape = NULL,
  overlapShpFolder = NULL,
  overlapShpFile = NULL,
  labelsSize = 1.2,
  fillcolorNA = NULL,
  projX = "+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0",
  extendedFillColor = "grey75",
  extendedBGColor = "lightblue1",
  extendedHighLightColor = "cornsilk1",
  extendedLabelsColor = "grey30",
  extdendedLabelSize = 0.7,
  extension = T,
  fillPalette = "Spectral",
  cropSubShape2Bound = T,
  grids = NULL,
  innerMargins = c(0.1, 0.2, 0.1, 0.2),
  outerMargins = c(0.01, 0.01, 0.01, 0.01),
  compassScale = T,
  scalePos = c("right", "bottom"),
  compassPos = c("left", "bottom"),
  extendedBoundary = NULL,
  extendedBoundaryOV = NULL
)

Arguments

boundaryRegShape

Default=NULL. Boundary region shape if already read into R.

boundaryRegShpFolder

Default= NULL. Folder containing boundary region shapefile. Suggested: paste(getwd(),"/dataFiles/gis/metis/naturalEarth",sep Default=""),

boundaryRegShpFile

Default=NULL. Name of shapefile. Suggested: paste("ne_10m_admin_0_countries",sep Default=""),

boundaryRegCol

Default=NULL. Column name with region names. Suggested "NAME_0",

boundaryRegionsSelect

Default=NULL. The region to choose from the given shapefile.

subRegShape

Default=NULL. Sub-region shape if already read into R.

subRegShpFolder

Default=NULL. Folder containing boundary region shapefile. Suggested paste(getwd(),"/dataFiles/gis/metis/naturalEarth",sep Default=""),

subRegShpFile

Default=NULL. Name of sub-region shapefile. Suggested paste("ne_10m_admin_1_states_provinces",sep Default=""),

subRegCol

Default= NULL. Suggested for states "name",

subRegCol1

Default= NULL. Suggested for states "name",

subRegionsSelect

Default=NULL. The region to choose from the given sub-region shapefile.

subRegType

Default="subRegType". Type of subregion. Eg. "states", "basins" etc.

dirOutputs

Default=paste(getwd(),"/outputs",sep Default=""). Location for outputs.

regionName

Default = NULL,

folderName

Default = NULL,

nameAppend

Default="". Name to append to saved files.

expandPercentWidth

Default=2. Percentage to expand boundary region beyond chosen region.

expandPercentHeight

Default=2. Percentage to expand boundary region beyond chosen region.

expandPercentWidthOV

Default=2. Percentage to expand boundary region beyond chosen region.

expandPercentHeightOV

Default=2. Percentage to expand boundary region beyond chosen region.

overlapShape

Default = NULL. If boundary lines of another shapefile are desired specify the shape here.

overlapShpFolder

Default = NULL. For GCAM basins use paste(getwd(),"/dataFiles/gis/metis/gcam",sep="").

overlapShpFile

Default = NULL. For GCAM basins use ="Global235_CLM_final_5arcmin_multipart"

labelsSize

Default =1.2.

fillcolorNA

Default =NULL. Fill color for NA values.

projX

Default ="+proj=longlat +datum=WGS84 +no_defs +ellps=WGS84 +towgs84=0,0,0".

extendedFillColor

Default = "grey75". Color used to fill extended land areas.

extendedBGColor

Default = "lightblue1". Color used to fill background/water bodies.

extendedHighLightColor

Default = "cornsilk1". Color used to highlight region of analysis.

extendedLabelsColor

Default = "grey30". Color for extended country name labels.

extdendedLabelSize

Default =0.7. Size of extended country name labels.

extension

Default = T. Should the map be extended beyond chosen shapefile boudnaries.

fillPalette

Default ="Spectral". Palette to use to fill subregions.

cropSubShape2Bound

Default = T. If subregion shape file is larger than boundary file.

grids

Default = NULL. Metis comes with 0.5 and 0.25 grids in c(paste(getwd(),"/dataFiles/grids/emptyGrids/grid_025.csv",sep=""),

innerMargins

Default =c(0,0.1,0,0.1), # bottom, left, top, right

outerMargins

Default =c(0.01,0.01,0.01,0.01) # bottom, left, top, right

compassScale

Default =T

scalePos

Default = c("right","bottom"),

compassPos

Default = c("left","bottom")

extendedBoundary

Default=NULL, Shapefile

extendedBoundaryOV

Default=NULL, Shapefile paste(getwd(),"/dataFiles/grids/emptyGrids/grid_050.csv",sep="")) This may happen in the case of disputed boundaries.

Value

A table with data by polygon ID for each shapefile provided