Generate error statistics from validation point datasets overlaid on a raster surface
errorGen(
finalraster,
validation.sf_ob,
validation.data,
plot = FALSE,
title = ""
)Arguments
Value
List of error statistics
Examples
library(sf)
#> Linking to GEOS 3.12.1, GDAL 3.8.4, PROJ 9.4.0; sf_use_s2() is TRUE
validation.data <- data.frame(rnorm(10, mean = 0.2, sd = 1))
names(validation.data) <- c("validation")
validation.sf_ob <- validation.data
validation.data <- as.numeric(unlist(validation.data))
xy <- data.frame(x = c(0:9), y = rep(1, 10))
validation.sf_ob <- st_as_sf(cbind(validation.sf_ob, xy), coords = c("x", "y"))
m <- matrix(1, 1, 10)
out.ras <- terra::rast(nrows = nrow(m), ncols = ncol(m),
xmin = 0, xmax = ncol(m), ymin = 0, ymax = nrow(m))
out.ras[] <- validation.data + rnorm(ncell(out.ras), mean = 0.01, sd = 0.2)
valid.stats <- errorGen(as(out.ras, "Raster"), validation.sf_ob, validation.data, plot = TRUE,
title = "Validation Plot")
valid.stats
#> [[1]]
#> r2 rmse pe logmse
#> 1 0.96 0.22 52.4 -2.99
#>
#> [[2]]
#> validation.data predicted
#> 1 -1.20004352 -1.30078339
#> 2 0.45531705 0.59111346
#> 3 -2.23726361 -1.81425863
#> 4 0.19442871 -0.12176917
#> 5 0.82155272 0.93403811
#> 6 1.34841161 0.98580931
#> 7 -1.62181766 -1.71622016
#> 8 -0.04732530 -0.04784568
#> 9 -0.04419961 0.07439966
#> 10 -0.08270545 -0.25552041
#>