model_o3 <- terra::rast(paste0(system.file("extdata",package="eva3dm"),
"/camx_no2.Rds"))
omi_o3 <- terra::rast(paste0(system.file("extdata",package="eva3dm"),
"/omi_no2.Rds"))
# generate the statistical indexes
sat(mo = model_o3,ob = omi_o3,rname = 'NO2_statistical')
# generate categorical evaluation using 3.0 as threshold
sat(mo = model_o3,ob = omi_o3,rname = 'NO2_categorical',
eval_function = cate, threshold = 3.0)
Run the code above in your browser using DataLab