August 2, 2026: testing the nature of the modis object after finding NA in extraction.

Author: rdjmacklin

R/landcover_extract.R

@@ -505,6 +505,8 @@ landcover_extract <- function(
       # Open the requested layer in file j.
       modis <- terra::rast(j)[index]
 
+      modis
+
       # Loop through each site matched to file j and extract.
       for (k in unique(pts_to_fill$SurveyAreaIdentifier)) {
         # If buffered, extract using landscapemetrics::calculate_lsm(). If not,
@@ -586,8 +588,6 @@ landcover_extract <- function(
           # so handle that possibility below.
           extr_table <- terra::extract(modis, tmp, fun = unique)[, index]
 
-          print(extr_table)
-
           # Whether only a single value was extracted (class == "integer") or
           # multiple values (else) prepare to pass to input data.
           if (inherits(extr_table, "integer")) {