R Under development (unstable) (2025-08-17 r88631 ucrt) -- "Unsuffered Consequences" Copyright (C) 2025 The R Foundation for Statistical Computing Platform: x86_64-w64-mingw32/x64 R is free software and comes with ABSOLUTELY NO WARRANTY. You are welcome to redistribute it under certain conditions. Type 'license()' or 'licence()' for distribution details. R is a collaborative project with many contributors. Type 'contributors()' for more information and 'citation()' on how to cite R or R packages in publications. Type 'demo()' for some demos, 'help()' for on-line help, or 'help.start()' for an HTML browser interface to help. Type 'q()' to quit R. > # This file is part of the standard setup for testthat. > # It is recommended that you do not modify it. > # > # Where should you do additional test configuration? > # Learn more about the roles of various files in: > # * https://r-pkgs.org/testing-design.html#sec-tests-files-overview > # * https://testthat.r-lib.org/articles/special-files.html > > library(testthat) > library(geospatialsuite) ================================================== GeoSpatialSuite v0.1.0 - Reliable Geospatiotemporal Analysis ================================================== Universal - 60+ Vegetation Indices - Robust & Simple Reliable terra plotting (no complex dependencies) 60+ vegetation indices with auto band detection Universal mapping with quick_map() Enhanced NDVI with quality filtering Robust error handling throughout Quick Start: help(package='geospatialsuite') Test Package: test_geospatialsuite_package() ================================================== All optional packages available! > > test_check("geospatialsuite") | | | 0% | |= | 1% | |=== | 4% | |==== | 5% | |========== | 14% | |============ | 17% | |============= | 18% | |============= | 19% | |================ | 23% | |================= | 25% | |================== | 25% | |================== | 26% | |=================== | 27% | |===================== | 31% | |======================== | 35% | |=========================== | 39% | |================================= | 48% | |=================================== | 50% | |====================================== | 55% | |======================================== | 57% | |============================================== | 65% | |=============================================== | 68% | |=================================================== | 72% | |==================================================== | 74% | |========================================================== | 83% | |================================================================ | 91% | |================================================================== | 94% | |======================================================================| 100% The rnaturalearthdata package needs to be installed. Installing the rnaturalearthdata package. Installing package into 'D:/temp/2025_08_18_22_00_17_18652/RtmpqQOSoy/RLIBS_899c3f1f45d5' (as 'lib' is unspecified) Starting dynamic CDL crop analysis... Analyzing 2 crop codes for custom_selection Raster info: 25 cells, resolution: 0.200000 x 0.200000 CRS: GEOGCRS["WGS 84", ENSEMBLE["World Geodetic System 1984 ensemble", MEMBER["World Geodetic System 1984 (Transit)"], MEMBER["World Geodetic System 1984 (G730)"], MEMBER["World Geodetic System 1984 (G873)"], MEMBER["World Geodetic System 1984 (G1150)"], MEMBER["World Geodetic System 1984 (G1674)"], MEMBER["World Geodetic System 1984 (G1762)"], MEMBER["World Geodetic System 1984 (G2139)"], MEMBER["World Geodetic System 1984 (G2296)"], ELLIPSOID["WGS 84",6378137,298.257223563, LENGTHUNIT["metre",1]], ENSEMBLEACCURACY[2.0]], PRIMEM["Greenwich",0, ANGLEUNIT["degree",0.0174532925199433]], CS[ellipsoidal,2], AXIS["geodetic latitude (Lat)",north, ORDER[1], ANGLEUNIT["degree",0.0174532925199433]], AXIS["geodetic longitude (Lon)",east, ORDER[2], ANGLEUNIT["degree",0.0174532925199433]], USAGE[ SCOPE["Horizontal component of 3D system."], AREA["World."], BBOX[-90,-180,90,180]], ID["EPSG",4326]] Is geographic (lon/lat): TRUE Geographic coordinates detected - calculating pixel area using cellSize() Calculated pixel size: 37651.396360 hectares Code 1: 9 pixels = 338862.57 hectares Code 5: 8 pixels = 301211.17 hectares Area Analysis Summary: Crop: custom_selection (codes: 1, 5) Region: Full extent Total crop area: 640073.74 hectares (1581654.21 acres) Pixel size: 37651.396360 hectares Raster: 25 total pixels, 25 valid pixels CDL analysis completed! Rasters have different geometries. Attempting to align... Starting comprehensive geospatial workflow with simplified visualization... Starting enhanced NDVI crop analysis workflow... Step 1: Calculating enhanced NDVI... Starting NDVI calculation for time series analysis... Using provided SpatRaster objects directly NDVI calculation completed successfully! - 1 layer(s) created - Value range: [0.188, 0.762] Step 2: Calculating additional vegetation indices... Step 3: Applying region boundary... Step 5: Creating visualizations... Step 6: Saving raster results... Step 7: Generating statistics... Enhanced NDVI crop workflow completed successfully! Comprehensive workflow completed! Results saved to: D:\temp\2025_08_18_22_00_17_18652\RtmpwtAa1T [ FAIL 0 | WARN 0 | SKIP 0 | PASS 68 ] > > proc.time() user system elapsed 5.95 0.46 7.59