tree cover max year

Author: zivankaraman

R/calc_treecover_area.R

@@ -155,7 +155,7 @@ calc_treecover_area <- function(years = 2000:2024,
   data[data[["patches"]] %in% keep_patches, ]
 }
 
-.sum_gfw <- function(data, what = "coverage_area", max_year = 2023) {
+.sum_gfw <- function(data, what, max_year) {
   # calculate loss area by year
   df <- data.frame(years = 2000:max_year, var = 0)
   names(df)[2] <- what
@@ -244,7 +244,7 @@ calc_treecover_area <- function(years = 2000:2024,
   rcl<- matrix(c(NA, NA, 0,
                  0, cover, 0,
                  cover, 100.001, 1), ncol = 3, byrow = TRUE)
-  # need vale > 100 to include 100 with 'right = FALSE'
+  # need value > 100 to include 100 with 'right = FALSE'
   binary_treecover <- terra::classify(treecover,
                                       rcl = rcl,
                                       right = FALSE, others = NA)

R/calc_treecover_area_and_emissions.R

@@ -102,7 +102,7 @@ calc_treecover_area_and_emissions <- function(years = 2000:2024,
         data <- .prep_gfw_data(data, min_size)
         losses <- .sum_gfw(data, "coverage_area", max_year)
         names(losses)[2] <- "loss"
-        emissions <- .sum_gfw(data, "emissions")
+        emissions <- .sum_gfw(data, "emissions", max_year)
 
         org_coverage <- sum(data[["coverage_area"]])
 

tests/testthat/test-calc_treecover_area.R

@@ -63,9 +63,9 @@ test_that("test gfw utils", {
     "patches", "emissions", "coverage_area"
   ))
 
-  expect_silent(area <- .sum_gfw(df, "coverage_area"))
+  expect_silent(area <- .sum_gfw(df, "coverage_area", 2023))
   expect_equal(names(area), c("years", "coverage_area"))
-  expect_silent(emis <- .sum_gfw(df, "emissions"))
+  expect_silent(emis <- .sum_gfw(df, "emissions", 2023))
   expect_equal(names(emis), c("years", "emissions"))
 })