data-raw/ParacouSubset.R

In VincyaneBadouard/TreeData: Forest Inventories Harmonization & Correction

## code to prepare `ParacouSubset` dataset goes here


#### Packages libraries ####
library(TreeData)
library(data.table)
library(EcoFoG)

#### Import data ####
Paracou6 <- EcoFoG::Guyafor2df(WHERE="Forest='Paracou' AND Plot='6'")
setDT(Paracou6) # as data.table

# Only SubPlot 1
Paracou6C1 <- Paracou6[SubPlot == 1]

# We select the 5 last censuses 2016-2020:

ParacouSubset <- Paracou6C1[CensusYear >= 2016] # 4904 obs

ind <- sample(unique(ParacouSubset$idTree), 1000)

ParacouSubset <- ParacouSubset[idTree %in% ind] # 4854 obs

# 4 sub-sub-squares creation:

MaxX <- max(ParacouSubset$Xutm)
MinX <- min(ParacouSubset$Xutm)
MaxY <- max(ParacouSubset$Yutm)
MinY <- min(ParacouSubset$Yutm)

HalfX <- MinX + (MaxX - MinX)/2
HalfY <- MinY + (MaxY - MinY)/2

ParacouSubset[, SubSubPlot := NA_real_]
ParacouSubset[Xutm < HalfX & Yutm > HalfY, SubSubPlot := "1"]
ParacouSubset[Xutm > HalfX & Yutm > HalfY, SubSubPlot := "2"]
ParacouSubset[Xutm < HalfX & Yutm < HalfY, SubSubPlot := "3"]
ParacouSubset[Xutm > HalfX & Yutm < HalfY, SubSubPlot := "4"]

ParacouSubset[, Project := NULL]

#### Save this data in the package ####
usethis::use_data(ParacouSubset, overwrite = TRUE)