Title:	Comprehensive Geospatiotemporal Analysis and Multimodal Integration Toolkit
Version:	0.1.0
Maintainer:	Olatunde D. Akanbi <olatunde.akanbi@case.edu>
Description:	A comprehensive toolkit for geospatiotemporal analysis featuring 60+ vegetation indices, advanced raster visualization, universal spatial mapping, water quality analysis, CDL crop analysis, spatial interpolation, temporal analysis, and terrain analysis. Designed for agricultural research, environmental monitoring, remote sensing applications, and publication-quality mapping with support for any geographic region and robust error handling. Methods include vegetation indices calculations (Rouse et al. 1974), NDVI and enhanced vegetation indices (Huete et al. 1997) <doi:10.1016/S0034-4257(97)00104-1>, spatial interpolation techniques (Cressie 1993, ISBN:9780471002556), water quality indices (McFeeters 1996) <doi:10.1080/01431169608948714>, and crop data layer analysis (USDA NASS 2024) https://www.nass.usda.gov/Research_and_Science/Cropland/. Funding: This material is based upon financial support by the National Science Foundation, EEC Division of Engineering Education and Centers, NSF Engineering Research Center for Advancing Sustainable and Distributed Fertilizer production (CASFER), NSF 20-553 Gen-4 Engineering Research Centers award 2133576.
License:	MIT + file LICENSE
Depends:	R (≥ 3.5.0)
Imports:	dplyr, ggplot2, graphics, grDevices, htmlwidgets, leaflet, magrittr, mice, parallel, RColorBrewer, rnaturalearth, sf, stats, stringr, terra, tigris, tools, utils, viridis
Suggests:	knitr, pkgnet, rmarkdown, testthat (≥ 3.0.0)
VignetteBuilder:	knitr
Config/testthat/edition:	3
Encoding:	UTF-8
RoxygenNote:	7.3.2
NeedsCompilation:	no
Packaged:	2025-09-11 07:55:14 UTC; oda10
Author:	Olatunde D. Akanbi [aut, cre, cph], Vibha Mandayam [aut], Yinghui Wu [aut], Jeffrey Yarus [aut], Erika I. Barcelos [aut, cph], Roger H. French [aut, cph]
Repository:	CRAN
Date/Publication:	2025-09-16 06:00:02 UTC

geospatialsuite: Comprehensive Geospatiotemporal Analysis and Multimodal Integration Toolkit

Description

GeoSpatialSuite provides universal functions for geospatial analysis and reliable visualization that work with any region for multimodal data. Features include 60+ vegetation indices, efficient terra-based visualization, universal spatial mapping, dynamic crop analysis, water quality assessment, and publication-quality mapping with support for any geographic region and robust error handling.

Details

Key Features:

Universal Spatial Analysis:

• Universal region support (states, countries, CONUS, custom boundaries)

• Universal spatial join (works with ANY raster-vector combination)

• Multi-dataset integration and temporal analysis

• Spatial interpolation and terrain analysis

Advanced Vegetation Analysis:

• 60+ vegetation indices including NDVI, EVI, SAVI, ARVI, PRI, SIPI, etc.

• Specialized crop analysis with stress detection and yield assessment

• Auto band detection from multi-band satellite imagery

• Quality filtering and temporal smoothing for time series

Reliable Visualization:

• Universal mapping with auto-detection (quick_map() function)

• Terra-based plotting using reliable terra::plot() and terra::plotRGB()

• Interactive maps with leaflet integration (optional)

• RGB composites with stretching algorithms

• Comparison maps for before/after analysis

Agricultural Applications:

• Dynamic CDL crop analysis (all crop codes and categories)

• NDVI time series with classification

• Crop-specific analysis (corn, soybeans, wheat, etc.)

• Water quality assessment with threshold analysis

Performance & Reliability Features:

• Standard terra plotting - no complex dependencies required

• Robust error handling throughout all functions

• Simplified visualization for maps

• Smart fallback systems when optional packages unavailable

Quick Start Examples:

# One-line mapping (auto-detects everything!)
quick_map("mydata.shp")

# Calculate multiple vegetation indices
indices <- calculate_multiple_indices(
  red = red_band, nir = nir_band,
  indices = c("NDVI", "EVI", "SAVI", "PRI")
)

# Comprehensive crop analysis
crop_analysis <- analyze_crop_vegetation(
  spectral_data = sentinel_data,
  crop_type = "corn",
  analysis_type = "comprehensive"
)

# Enhanced NDVI calculation
ndvi_enhanced <- calculate_ndvi_enhanced(
  red_data = red_raster,
  nir_data = nir_raster,
  quality_filter = TRUE
)

# Fast, reliable RGB plotting
plot_rgb_raster(satellite_data, r = 4, g = 3, b = 2,
               stretch = "hist", title = "False Color")

Recommended Optional Packages:

For enhanced features, consider installing these optional packages:

# For interactive mapping
install.packages("leaflet")

# For enhanced colors
install.packages(c("viridis", "RColorBrewer"))

# For advanced remote sensing (optional)
install.packages("RStoolbox")

# For multi-panel plots (optional)
install.packages("patchwork")

Core Dependencies Only:

The package works reliably with just the core dependencies:

• terra (raster operations and plotting)

• sf (vector operations)

• ggplot2 (static mapping)

• dplyr (data manipulation)

No complex visualization dependencies required!

Author(s)

Olatunde D. Akanbi olatunde.akanbi@case.edu

Erika I. Barcelos erika.barcelos@case.edu

Roger H. French roger.french@case.edu

Add boundary overlay to existing plot

Description

Add boundary overlay to existing plot

Usage

add_boundary_overlay(region_boundary, verbose)

Add crop statistics overlay

Description

Add crop statistics overlay

Usage

add_crop_statistics_overlay(plot_obj, cdl_raster, crop_codes, crop_names)

Add lines to leaflet -

Description

Add lines to leaflet -

Usage

add_lines_to_leaflet(
  map_obj,
  spatial_data,
  fill_variable,
  popup_vars,
  color_scheme
)

Add points to leaflet -

Description

Add points to leaflet -

Usage

add_points_to_leaflet(
  map_obj,
  spatial_data,
  fill_variable,
  popup_vars,
  color_scheme,
  cluster_points
)

Add polygons to leaflet -

Description

Add polygons to leaflet -

Usage

add_polygons_to_leaflet(
  map_obj,
  spatial_data,
  fill_variable,
  popup_vars,
  color_scheme
)

Analyze CDL crops dynamically

Description

Perform dynamic analysis of CDL crop data including area calculations, classification, and masking for any crop combination. Now accepts file paths directly.

Usage

analyze_cdl_crops_dynamic(
  cdl_data,
  crop_selection,
  region_boundary = NULL,
  analysis_type = "classify",
  output_folder = NULL,
  year = NULL
)

Arguments

Details

Usage Tips:

Area Analysis Results:

# Access specific results
result$total_area_ha          # Total area in hectares
result$crop_areas_ha          # Area by crop code
result$crop_areas_ha[["5"]]   # Soybean area (code 5)

Visualization:

# For mask/classify results (SpatRaster):
terra::plot(result)           # Plot the raster

# For area results (list):
print(result$total_area_ha)   # Print total area
barplot(unlist(result$crop_areas_ha))  # Simple bar plot

Value

Analysis results depend on analysis_type:

• "mask": SpatRaster with binary mask (1 = crop, 0 = other)

• "classify": SpatRaster with binary classification

• "area": List with detailed area analysis:

  • crop_areas_ha: Named list of areas by crop code (hectares)

  • total_area_ha: Total crop area (hectares)

  • pixel_size_ha: Individual pixel size (hectares)

  • crop_codes: CDL codes analyzed

  • crop_name: Crop name(s)

  • region: Region analyzed

  • raster_info: Technical raster information

  • total_pixels: Total number of pixels

  • valid_pixels: Number of valid (non-NA) pixels

Examples

## Not run: 
# These examples require actual CDL data files
# Analyze soybean area in Iowa - accepts file path directly
soybean_area <- analyze_cdl_crops_dynamic(
  "/path/to/cdl_2023.tif", "soybeans", "Iowa", "area"
)

# Access area results
soybean_area$total_area_ha                    # Total hectares
soybean_area$crop_areas_ha[["5"]]            # Soybean area (code 5)
soybean_area$total_area_ha * 2.47105         # Convert to acres

# Create grain classification
grain_classes <- analyze_cdl_crops_dynamic(
  cdl_data, "grains", "CONUS", "classify"
)
terra::plot(grain_classes)  # Plot the classification

# Works with directories too
results <- analyze_cdl_crops_dynamic(
  "/path/to/cdl/files/", "corn", "Ohio", "area"
)

## End(Not run)


# Example with mock CDL data (this can run)
# Create sample CDL raster
mock_cdl <- terra::rast(nrows = 10, ncols = 10, xmin = 0, xmax = 10,
                       ymin = 0, ymax = 10, crs = "EPSG:4326")
terra::values(mock_cdl) <- sample(c(1, 5, 24), 100, replace = TRUE)  # corn, soy, wheat

# Analyze mock data
result <- analyze_cdl_crops_dynamic(mock_cdl, "corn", analysis_type = "mask")
print(class(result))  # Should be SpatRaster

Specialized crop vegetation analysis

Description

Perform comprehensive vegetation analysis specifically designed for crop monitoring including growth stage detection, stress identification, and yield prediction support. Handles test scenarios properly with better input validation.

Usage

analyze_crop_vegetation(
  spectral_data,
  crop_type = "general",
  growth_stage = "unknown",
  analysis_type = "comprehensive",
  cdl_mask = NULL,
  reference_data = NULL,
  output_folder = NULL,
  verbose = FALSE
)

Arguments

Details

Crop-Specific Index Selection:

• Corn: NDVI, EVI, GNDVI, DVI, RVI, PRI

• Soybeans: NDVI, EVI, SAVI, GNDVI, PRI

• Wheat: NDVI, EVI, SAVI, DVI

• General: NDVI, EVI, SAVI, GNDVI, DVI, RVI

Analysis Types:

• comprehensive: All analyses (stress, growth, yield)

• stress: Focus on stress detection indices

• growth: Growth stage analysis

• yield: Yield prediction support

Value

List with comprehensive vegetation analysis results:

• vegetation_indices: SpatRaster with calculated indices

• analysis_results: Detailed analysis results by type

• metadata: Analysis metadata and parameters

Examples

## Not run: 
# These examples require actual spectral data
# Comprehensive corn analysis
corn_analysis <- analyze_crop_vegetation(
  spectral_data = sentinel_data,
  crop_type = "corn",
  analysis_type = "comprehensive",
  cdl_mask = corn_mask
)

# Access results
corn_analysis$vegetation_indices        # SpatRaster with indices
corn_analysis$analysis_results$stress_analysis  # Stress detection results
corn_analysis$metadata$indices_used     # Which indices were calculated

# Stress detection in soybeans
stress_analysis <- analyze_crop_vegetation(
  spectral_data = landsat_stack,
  crop_type = "soybeans",
  analysis_type = "stress",
  growth_stage = "mid"
)

## End(Not run)


# Example with mock spectral data
# Create mock multi-band raster (simulating satellite data)
red_band <- terra::rast(nrows = 5, ncols = 5, crs = "EPSG:4326")
nir_band <- terra::rast(nrows = 5, ncols = 5, crs = "EPSG:4326")
terra::values(red_band) <- runif(25, 0.1, 0.3)   # Typical red values
terra::values(nir_band) <- runif(25, 0.4, 0.8)   # Typical NIR values
spectral_stack <- c(red_band, nir_band)
names(spectral_stack) <- c("red", "nir")

# Analyze with mock data
result <- analyze_crop_vegetation(spectral_stack, crop_type = "general")
print(names(result))  # Should show analysis components

Analyze growth stage

Description

Analyze growth stage

Usage

analyze_growth_stage(vegetation_indices, crop_type, growth_stage, verbose)

Analyze seasonal patterns

Description

Internal function to analyze seasonal patterns in temporal data.

Usage

analyze_seasonal_patterns(rasters, dates)

Arguments

Value

List with seasonal analysis results

Analyze temporal changes in geospatial data

Description

Analyze temporal changes in raster data including trend analysis, change detection, and seasonal patterns. Works with any time series data.

Usage

analyze_temporal_changes(
  data_list,
  dates = NULL,
  region_boundary = NULL,
  analysis_type = "trend",
  output_folder = NULL
)

Arguments

Value

Temporal analysis results

Examples

## Not run: 
# These examples require external data files not included with the package
# Analyze NDVI trends over time
ndvi_trend <- analyze_temporal_changes(
  data_list = c("ndvi_2020.tif", "ndvi_2021.tif", "ndvi_2022.tif"),
  dates = c("2020", "2021", "2022"),
  region_boundary = "Iowa",
  analysis_type = "trend"
)

# Detect land cover changes
land_changes <- analyze_temporal_changes(
  data_list = land_cover_files,
  dates = land_cover_dates,
  analysis_type = "change_detection"
)

## End(Not run)

Analyze correlations between multiple variables

Description

Analyze correlations between multiple raster variables and create correlation matrices and plots.

Usage

analyze_variable_correlations(
  variable_list,
  output_folder = NULL,
  region_boundary = NULL,
  method = "pearson",
  create_plots = TRUE
)

Arguments

Value

List with correlation results

Examples

## Not run: 
# These examples require directory structures with multiple data files
# Analyze correlations between multiple variables
variables <- list(
  ndvi = "ndvi.tif",
  nitrogen = "soil_nitrogen.tif",
  elevation = "dem.tif",
  precipitation = "precip.tif"
)

correlation_results <- analyze_variable_correlations(
  variables,
  output_folder = "correlations/",
  region_boundary = "Ohio"
)

## End(Not run)

Analyze water body characteristics using multiple indices

Description

Comprehensive water body analysis using multiple water indices to classify and characterize water features.

Usage

analyze_water_bodies(
  green,
  nir,
  swir1 = NULL,
  region_boundary = NULL,
  water_threshold_ndwi = 0.3,
  water_threshold_mndwi = 0.5,
  output_folder = NULL,
  verbose = FALSE
)

Arguments

Value

List with water analysis results

Examples

## Not run: 
# These examples require external data files not included with the package
# Comprehensive water analysis
water_analysis <- analyze_water_bodies(
  green = "green.tif",
  nir = "nir.tif",
  swir1 = "swir1.tif",
  region_boundary = "study_area.shp",
  verbose = TRUE
)

# Access results
water_analysis$water_indices     # All calculated indices
water_analysis$water_mask        # Binary water mask
water_analysis$statistics        # Water body statistics

## End(Not run)

Analyze water quality comprehensively with flexible data handling

Description

Complete water quality analysis with flexible data input handling, robust error checking, and comprehensive spatial integration. Supports any water quality dataset format with automatic column detection and standardized processing.

Usage

analyze_water_quality_comprehensive(
  water_data,
  variable = NULL,
  region_boundary = NULL,
  river_network = NULL,
  output_folder = tempdir(),
  thresholds = NULL,
  coord_cols = NULL,
  date_column = NULL,
  station_id_col = NULL,
  quality_filters = list(),
  verbose = FALSE
)

Arguments

Value

List with comprehensive water quality analysis results:

• water_data: Processed spatial data

• statistics: Summary statistics by variable and category

• spatial_analysis: Spatial pattern analysis

• temporal_analysis: Temporal trends (if date data available)

• threshold_analysis: Threshold exceedance analysis

• output_files: Paths to generated output files

• metadata: Analysis metadata and parameters

Examples

## Not run: 
# These examples require external data files not included with the package
# Flexible data input - auto-detects columns
results <- analyze_water_quality_comprehensive("water_stations.csv")

# Specify parameters for custom data
results <- analyze_water_quality_comprehensive(
  water_data = my_data,
  variable = "nitrate_concentration",
  region_boundary = "Ohio",
  coord_cols = c("longitude", "latitude"),
  thresholds = list(
    Normal = c(0, 2),
    Elevated = c(2, 5),
    High = c(5, 10),
    Critical = c(10, Inf)
  )
)

# Multi-dataset analysis
results <- analyze_water_quality_comprehensive(
  water_data = list(
    surface = "surface_water.csv",
    groundwater = "groundwater.csv"
  ),
  variable = "total_nitrogen"
)

## End(Not run)

Analyze spatial patterns in water quality data

Description

Analyze spatial patterns in water quality data

Usage

analyze_water_spatial_patterns(water_sf, variable, verbose)

Analyze temporal patterns in water quality data

Description

Analyze temporal patterns in water quality data

Usage

analyze_water_temporal_patterns(water_sf, variable, date_column, verbose)

Analyze yield potential

Description

Analyze yield potential

Usage

analyze_yield_potential(vegetation_indices, crop_type, indices, verbose)

Apply color scheme to ggplot

Description

Apply color scheme to ggplot

Usage

apply_color_scheme_ggplot(plot, color_scheme, aesthetic = "fill")

Apply quality filter to remove outliers

Description

Apply quality filter to remove outliers

Usage

apply_quality_filter(ndvi, verbose)

Apply temporal smoothing

Description

Apply temporal smoothing

Usage

apply_temporal_smoothing(ndvi_stack, verbose)

Apply water quality filters

Description

Apply water quality filters

Usage

apply_water_quality_filters(water_sf, variable, quality_filters, verbose)

Assess data characteristics for interpolation method selection

Description

Assess data characteristics for interpolation method selection

Usage

assess_data_for_interpolation(spatial_sf, target_variables, verbose)

Auto-detect best spatial join method

Description

Auto-detect best spatial join method

Usage

auto_detect_method(source_type, target_type, verbose = FALSE)

Arguments

Value

Character string of recommended method

Auto-detect spectral bands

Description

Auto-detect spectral bands

Usage

auto_detect_spectral_bands(spectral_data, layer_names, verbose)

Auto-detect title for quick mapping

Description

Auto-detect title for quick mapping

Usage

auto_detect_title(spatial_data, variable)

Auto-detect variable for quick mapping

Description

Auto-detect variable for quick mapping

Usage

auto_detect_variable(spatial_data)

Calculate advanced terrain metrics

Description

Calculate advanced terrain metrics from DEM including curvature, wetness index, and stream power index.

Usage

calculate_advanced_terrain_metrics(
  elevation_raster,
  metrics = c("wetness_index", "curvature", "convergence"),
  region_boundary = NULL
)

Arguments

Value

List of terrain metric rasters

Examples

## Not run: 
# These examples require external data files not included with the package
# Calculate advanced terrain metrics
terrain_metrics <- calculate_advanced_terrain_metrics(
  elevation_raster = "dem.tif",
  metrics = c("wetness_index", "curvature", "convergence"),
  region_boundary = "watershed.shp"
)

## End(Not run)

Calculate comprehensive vegetation statistics -

Description

Calculate comprehensive vegetation statistics -

Usage

calculate_comprehensive_vegetation_stats(vegetation_stack, indices)

Calculate contributing area

Description

Internal function to calculate contributing area for wetness index.

Usage

calculate_contributing_area(dem)

Arguments

Value

SpatRaster of contributing area

Calculate index by type

Description

Calculate index by type

Usage

calculate_index_by_type(
  index_type,
  red,
  nir,
  blue,
  green,
  swir1,
  swir2,
  red_edge,
  coastal,
  nir2,
  verbose
)

Calculate multiple vegetation indices at once

Description

Calculate multiple vegetation indices from the same spectral data in a single operation. Efficient for comparative analysis and comprehensive vegetation assessment. Supports directory input and automatic CRS handling.

Usage

calculate_multiple_indices(
  spectral_data = NULL,
  indices = c("NDVI", "EVI", "SAVI"),
  output_stack = TRUE,
  region_boundary = NULL,
  parallel = FALSE,
  verbose = FALSE,
  ...
)

Arguments

Value

SpatRaster stack or list of indices

Examples

## Not run: 
# These examples require satellite imagery files (Landsat/Sentinel data etc.)
# Calculate multiple basic indices from directory
multi_indices <- calculate_multiple_indices(
  spectral_data = "/path/to/sentinel/bands/",
  indices = c("NDVI", "EVI", "SAVI", "MSAVI"),
  auto_detect_bands = TRUE
)

# Comprehensive vegetation analysis from individual files
veg_analysis <- calculate_multiple_indices(
  red = red_band, nir = nir_band, blue = blue_band,
  indices = c("NDVI", "EVI", "ARVI", "GNDVI", "DVI"),
  output_stack = TRUE,
  region_boundary = "Iowa"
)

# Directory with custom band matching
stress_indices <- calculate_multiple_indices(
  spectral_data = "/path/to/bands/",
  indices = c("PRI", "SIPI", "NDRE"),
  band_names = c("red", "green", "nir", "red_edge"),
  output_stack = TRUE
)

## End(Not run)

Calculate multiple water indices at once

Description

Calculate multiple water indices from the same spectral data in a single operation. Efficient for comprehensive water and moisture analysis.

Usage

calculate_multiple_water_indices(
  green,
  nir,
  swir1 = NULL,
  indices = c("NDWI", "MNDWI", "NDMI"),
  output_stack = TRUE,
  clamp_values = TRUE,
  mask_invalid = TRUE,
  verbose = FALSE
)

Arguments

Value

SpatRaster stack or list of water indices

Examples

## Not run: 
# These examples require external data files not included with the package
# Calculate multiple water indices
water_indices <- calculate_multiple_water_indices(
  green = green_band,
  nir = nir_band,
  swir1 = swir1_band,
  indices = c("NDWI", "MNDWI", "NDMI", "MSI"),
  output_stack = TRUE,
  verbose = TRUE
)

# Access individual indices
ndwi <- water_indices[["NDWI"]]
mndwi <- water_indices[["MNDWI"]]

## End(Not run)

Calculate NDVI with time series options

Description

NDVI calculation specifically designed for time series analysis with date matching, quality filtering, temporal smoothing, and multi-temporal support. Use this for time series analysis, use calculate_vegetation_index() for single dates.

Usage

calculate_ndvi_enhanced(
  red_data,
  nir_data,
  clamp_range = c(-0.2, 1),
  match_by_date = FALSE,
  quality_filter = FALSE,
  temporal_smoothing = FALSE,
  verbose = FALSE,
  date_patterns = NULL
)

Arguments

Details

When to Use Enhanced vs Basic NDVI:

Use calculate_ndvi_enhanced() for:

• Time series analysis: Multiple dates, trend analysis

• Quality control: Remove outliers, temporal smoothing

• Date matching: Automatic pairing of red/NIR by date

• Multi-temporal studies: Seasonal analysis, change detection

Use calculate_vegetation_index(index_type="NDVI") for:

• Single date analysis: One-time calculation

• Different indices: Want to calculate EVI, SAVI, etc. too

• Quick calculations: Simple, fast NDVI

• Mixed workflows: Part of larger vegetation index analysis

Value

SpatRaster with NDVI layers (single or multi-layer for time series)

Examples

## Not run: 
# These examples require external data files not included with the package
# Time series NDVI with date matching
ndvi_series <- calculate_ndvi_enhanced(
  red_data = "/path/to/red/time_series/",
  nir_data = "/path/to/nir/time_series/",
  match_by_date = TRUE,
  quality_filter = TRUE,
  temporal_smoothing = TRUE
)

# Simple NDVI (single date with quality control)
ndvi_clean <- calculate_ndvi_enhanced(
  red_data = red_raster,
  nir_data = nir_raster,
  quality_filter = TRUE
)

## End(Not run)

Calculate skewness for data assessment

Description

Calculate skewness for data assessment

Usage

calculate_skewness(x)

Calculate spatial correlation between raster layers

Description

Calculate spatial correlation between two raster layers using various methods. Supports pixel-wise correlation and local correlation analysis.

Usage

calculate_spatial_correlation(
  raster1,
  raster2,
  method = "pearson",
  local_correlation = FALSE,
  window_size = 3
)

Arguments

Value

Correlation coefficient or SpatRaster of local correlations

Examples

## Not run: 
# These examples require external data files not included with the package
# Global correlation between NDVI and soil nitrogen
correlation <- calculate_spatial_correlation(ndvi_raster, nitrogen_raster)

# Local correlation with moving window
local_corr <- calculate_spatial_correlation(
  ndvi_raster, nitrogen_raster,
  local_correlation = TRUE,
  window_size = 5
)

## End(Not run)

Calculate temporal statistics

Description

Internal function to calculate temporal statistics.

Usage

calculate_temporal_statistics(rasters, dates)

Arguments

Value

List with temporal statistics

Calculate temporal trend using linear regression

Description

Internal function to calculate pixel-wise temporal trends.

Usage

calculate_temporal_trend(rasters, dates)

Arguments

Value

List with trend analysis results

Calculate comprehensive vegetation indices

Description

Calculate a wide range of vegetation indices from spectral bands with automatic band detection, comprehensive error handling, and validation. Supports 40+ different vegetation indices for various applications. Accepts directories, file lists, and automatic CRS handling.

Usage

calculate_vegetation_index(
  spectral_data = NULL,
  red = NULL,
  nir = NULL,
  blue = NULL,
  green = NULL,
  swir1 = NULL,
  swir2 = NULL,
  red_edge = NULL,
  coastal = NULL,
  nir2 = NULL,
  index_type = "NDVI",
  auto_detect_bands = FALSE,
  band_names = NULL,
  clamp_range = NULL,
  mask_invalid = TRUE,
  scale_factor = 1,
  auto_crs_fix = TRUE,
  verbose = FALSE
)

Arguments

Details

Input Format Support:

Single Calculation:

# Individual band files
ndvi <- calculate_vegetation_index(red = "red.tif", nir = "nir.tif", index_type = "NDVI")

# Multi-band raster
evi <- calculate_vegetation_index(spectral_data = "landsat.tif", index_type = "EVI",
                                 auto_detect_bands = TRUE)

Directory/Multiple Files:

# Directory with band files
savi <- calculate_vegetation_index(spectral_data = "/path/to/bands/",
                                  band_names = c("red", "nir"), index_type = "SAVI")

# File list
arvi <- calculate_vegetation_index(spectral_data = c("red.tif", "nir.tif", "blue.tif"),
                                  index_type = "ARVI")

Enhanced vs Basic NDVI:

Basic calculate_vegetation_index():

• Single time point calculation

• 40+ different indices

• Directory/file support

• Automatic CRS fixing

• Use for: Single-date analysis, comparing different indices

calculate_ndvi_enhanced():

• Time series support

• Quality filtering

• Temporal smoothing

• Date matching between red/NIR

• Use for: Multi-temporal analysis, time series trends

Value

SpatRaster of vegetation index

Examples

## Not run: 
# These examples require satellite imagery files (Landsat/Sentinel data etc.)
# Basic NDVI calculation
ndvi <- calculate_vegetation_index(red = red_band, nir = nir_band, index_type = "NDVI")

# Multi-band raster with auto-detection
evi <- calculate_vegetation_index(spectral_data = landsat_stack,
                                 index_type = "EVI", auto_detect_bands = TRUE)

# Directory with automatic band detection
savi <- calculate_vegetation_index(spectral_data = "/path/to/sentinel/bands/",
                                  index_type = "SAVI", auto_detect_bands = TRUE)

# Advanced index with custom parameters
pri <- calculate_vegetation_index(red = red_band, nir = nir_band, green = green_band,
                                 index_type = "PRI", clamp_range = c(-1, 1))

# Custom band names for multi-band data
ndvi <- calculate_vegetation_index(spectral_data = sentinel_data,
                                  band_names = c("B4", "B3", "B2", "B8"),
                                  index_type = "NDVI")

## End(Not run)

Calculate vegetation statistics

Description

Calculate vegetation statistics

Usage

calculate_vegetation_statistics(vegetation_indices, indices, verbose)

Calculate water indices including both NDWI variants

Description

Calculate various water indices including NDWI (McFeeters 1996), MNDWI (Xu 2006), and NDMI (Gao 1996) for water body detection and moisture content. Updated formulas based on latest research and satellite missions (2024).

Usage

calculate_water_index(
  green,
  nir,
  swir1 = NULL,
  index_type = "NDWI",
  clamp_range = NULL,
  mask_invalid = TRUE,
  verbose = FALSE
)

Arguments

Details

Available water indices with their specific applications:

Primary Water Detection Indices:

• NDWI (McFeeters 1996): (Green - NIR) / (Green + NIR) - Use: Open water body detection, flood mapping - Range: Values from -1 to 1, water bodies typically > 0.3 - Pros: Simple, effective for clear water - Cons: Sensitive to built-up areas, can overestimate water

• MNDWI (Xu 2006): (Green - SWIR1) / (Green + SWIR1) - Use: Enhanced water detection, urban water bodies - Range: Values from -1 to 1, water bodies typically > 0.5 - Pros: Better separation of water from built-up areas - Cons: Requires SWIR band, less effective with turbid water

Vegetation Moisture Indices:

• NDMI (Gao 1996): (NIR - SWIR1) / (NIR + SWIR1) - Use: Vegetation water content, drought monitoring - Range: Values from -1 to 1, higher values = more water content - Application: Agriculture, forest fire risk assessment

• MSI: SWIR1 / NIR - Moisture Stress Index - Use: Plant water stress detection - Range: [0, 5+], lower values = higher moisture

• NDII: (NIR - SWIR1) / (NIR + SWIR1) - Same as NDMI - Use: Alternative name for NDMI, vegetation moisture

Specialized Water Indices:

• WI: NIR / SWIR1 - Water Index (simple ratio)

• SRWI: NIR / SWIR1 - Simple Ratio Water Index

• LSWI: (NIR - SWIR1) / (NIR + SWIR1) - Land Surface Water Index

Band Requirements by Satellite:

• Landsat 8/9: Green=Band 3, NIR=Band 5, SWIR1=Band 6

• Sentinel-2: Green=Band 3, NIR=Band 8, SWIR1=Band 11

• MODIS: Green=Band 4, NIR=Band 2, SWIR1=Band 6

Value

SpatRaster of water index

Examples

## Not run: 
# These examples require external data files not included with the package
# Original NDWI for water body detection
ndwi <- calculate_water_index(green_band, nir_band, index_type = "NDWI")

# Modified NDWI for enhanced water detection (requires SWIR1)
mndwi <- calculate_water_index(green_band, nir_band, swir1_band, index_type = "MNDWI")

# NDMI for vegetation moisture monitoring
ndmi <- calculate_water_index(green_band, nir_band, swir1_band, index_type = "NDMI")

# With quality control
water_index <- calculate_water_index(
  green = "green.tif",
  nir = "nir.tif",
  swir1 = "swir1.tif",
  index_type = "MNDWI",
  clamp_range = c(-1, 1),
  mask_invalid = TRUE,
  verbose = TRUE
)

## End(Not run)

Calculate comprehensive water quality statistics

Description

Calculate comprehensive water quality statistics

Usage

calculate_water_quality_statistics(
  water_sf,
  variable,
  numeric_vars,
  threshold_results,
  verbose
)

Check raster compatibility

Description

Check raster compatibility

Usage

check_raster_compatibility(raster1, raster2, auto_align = TRUE)

Arguments

Value

Aligned raster2 or original if compatible

Check and load required packages

Description

Check and load required packages

Usage

check_required_packages(packages)

Classify spatial data type and load data

Description

Classify spatial data type and load data

Usage

classify_spatial_data(data_input, verbose = FALSE)

Arguments

Value

List with data, type, class, and path information

Combine multiple water quality datasets

Description

Combine multiple water quality datasets

Usage

combine_water_datasets(water_sf_list, verbose)

Compare interpolation methods

Description

Compare multiple interpolation methods using cross-validation and return performance metrics for method selection.

Usage

compare_interpolation_methods(
  spatial_data,
  target_variable,
  methods = c("NN", "simple", "spline"),
  cv_folds = 5,
  verbose = TRUE
)

Arguments

Value

Data frame with method comparison results

Examples

## Not run: 
# These examples require external data files not included with the package
# Compare interpolation methods
method_comparison <- compare_interpolation_methods(
  soil_data,
  target_variable = "nitrogen",
  methods = c("NN", "simple", "spline"),
  cv_folds = 10
)

# View results
print(method_comparison)
# Best method
best_method <- method_comparison$method[which.min(method_comparison$rmse)]

## End(Not run)

Compile interpolation results into requested format

Description

Compile interpolation results into requested format

Usage

compile_interpolation_results(
  interpolation_results,
  original_data,
  target_variables,
  output_format,
  verbose
)

Convert Julian date to standard format

Description

Internal function to convert Julian dates (YYYY-DDD) to standard YYYY-MM-DD format.

Usage

convert_julian_date(julian_str)

Arguments

Value

Standard date string "YYYY-MM-DD"

Convert data.frame to sf with flexible coordinate detection

Description

Convert data.frame to sf with flexible coordinate detection

Usage

convert_to_spatial_flexible(water_df, coord_cols, verbose)

Create base plot map as fallback

Description

Create base plot map as fallback

Usage

create_base_plot_map(spatial_sf, fill_variable, map_type, title, verbose)

Create comparison map (before/after, side-by-side)

Description

Create comparison maps showing before/after analysis or side-by-side comparisons using reliable terra plotting.

Usage

create_comparison_map(
  data1,
  data2,
  comparison_type = "side_by_side",
  titles = c("Dataset 1", "Dataset 2"),
  region_boundary = NULL,
  color_scheme = "viridis",
  output_file = NULL,
  verbose = FALSE
)

Arguments

Value

NULL (plots directly to device) or file path if saved

Examples

## Not run: 
# These examples require external data files not included with the package
# Before/after NDVI comparison
create_comparison_map("ndvi_2020.tif", "ndvi_2023.tif",
                     comparison_type = "side_by_side",
                     titles = c("2020", "2023"),
                     output_file = "ndvi_comparison.png")

## End(Not run)

Create correlation plots

Description

Internal function to create correlation plots and heatmaps.

Usage

create_correlation_plots(correlation_matrix, output_folder)

Arguments

Create crop area map -

Description

Create crop area map -

Usage

create_crop_area_map(cdl_raster, crop_codes, crop_names, boundary)

Create crop diversity map -

Description

Create crop diversity map -

Usage

create_crop_diversity_map(cdl_raster, crop_codes, boundary)

Create crop map visualization

Description

Create crop map visualization

Usage

create_crop_map(
  cdl_data,
  crop_selection,
  region_boundary = NULL,
  style = "categorical",
  output_file = NULL
)

Arguments

Value

Character string with output file path if saved, or NULL if plotted to screen. Called primarily for side effects (creating plots).

Create crop mask from CDL data

Description

Create binary or classified crop mask from USDA CDL data for specified crops. Fixed to handle terra operations properly.

Usage

create_crop_mask(
  cdl_data,
  crop_codes,
  region_boundary = NULL,
  mask_type = "binary"
)

Arguments

Value

SpatRaster with crop mask

Examples

## Not run: 
# These examples require actual CDL data files
# Create corn mask
corn_mask <- create_crop_mask("cdl_2023.tif", "corn", "Iowa")

# Create grain crops mask
grain_mask <- create_crop_mask(cdl_raster, "grains", mask_type = "preserve")

## End(Not run)


# Example with mock data (this can run)
mock_cdl <- terra::rast(nrows = 5, ncols = 5, crs = "EPSG:4326")
terra::values(mock_cdl) <- c(1, 1, 5, 5, 24, 1, 5, 5, 24, 24,
                            1, 1, 5, 24, 24, 5, 5, 24, 24, 1,
                            1, 5, 5, 24, 1)  # corn, soy, wheat

# Create corn mask from mock data
corn_mask <- create_crop_mask(mock_cdl, "corn")
print(terra::values(corn_mask))  # Should show 1s and 0s

Create dominant crop map -

Description

Create dominant crop map -

Usage

create_dominant_crop_map(cdl_raster, crop_codes, crop_names, boundary)

Create ggplot map safely

Description

Create ggplot map safely

Usage

create_ggplot_map_safe(
  spatial_sf,
  fill_variable,
  map_type,
  region_boundary,
  color_scheme,
  title,
  point_size,
  verbose
)

Create interactive map using leaflet (if available)

Description

Create interactive maps with leaflet integration when available. Falls back gracefully when leaflet is not installed.

Usage

create_interactive_map(
  spatial_data,
  fill_variable = NULL,
  popup_vars = NULL,
  basemap = "terrain",
  color_scheme = "viridis",
  title = "Interactive Map",
  verbose = FALSE
)

Arguments

Value

leaflet map object or NULL if leaflet unavailable

Examples

## Not run: 
# These examples demonstrate workflows with user's own spatial data
# Simple interactive point map
map <- create_interactive_map(study_sites, fill_variable = "ndvi_mean")

# Polygon map with custom basemap
map <- create_interactive_map(counties, fill_variable = "population",
                             basemap = "satellite")

## End(Not run)

Create interactive map safely

Description

Create interactive map safely

Usage

create_interactive_map_safe(
  spatial_sf,
  fill_variable,
  map_type,
  title,
  verbose
)

Create NDVI map visualization

Description

Create NDVI map visualization

Usage

create_ndvi_map(
  ndvi_data,
  region_boundary = NULL,
  ndvi_classes = "none",
  title = "NDVI Analysis",
  output_file = NULL
)

Arguments

Value

Character string with output file path if saved, or NULL if plotted to screen. Called primarily for side effects (creating plots).

Create overlay comparison

Description

Create overlay comparison

Usage

create_overlay_comparison(data1, data2, titles, color_scheme)

Create reliable raster map

Description

Create reliable raster map

Usage

create_raster_map_reliable(
  raster_data,
  region_boundary,
  color_scheme,
  title,
  output_file,
  verbose
)

Create raster mosaic with intelligent file selection

Description

Create mosaics from multiple raster files with various methods and intelligent file selection based on region boundaries.

Usage

create_raster_mosaic(
  input_data,
  method = "merge",
  region_boundary = NULL,
  output_file = NULL,
  parallel = FALSE
)

Arguments

Value

SpatRaster object

Examples

## Not run: 
# These examples require external data files not included with the package
# Basic mosaic
mosaic <- create_raster_mosaic("/path/to/rasters", method = "merge")

# Mosaic for specific region
ohio_mosaic <- create_raster_mosaic("/aster/files", "merge", "Ohio")

# Mean composite
mean_mosaic <- create_raster_mosaic(raster_list, method = "mean")

## End(Not run)

Create universal spatial map with reliable terra plotting

Description

Universal mapping function that works with any spatial data type using reliable terra and base R plotting. No complex dependencies required. Falls back gracefully when optional packages are unavailable.

Usage

create_spatial_map(
  spatial_data,
  fill_variable = NULL,
  coord_cols = c("lon", "lat"),
  region_boundary = NULL,
  map_type = "auto",
  color_scheme = "viridis",
  interactive = FALSE,
  title = NULL,
  point_size = 3,
  output_file = NULL,
  verbose = FALSE
)

Arguments

Value

ggplot2 object, leaflet map, or file path (depending on options)

Examples

## Not run: 
# These examples require external data files not included with the package
# Simple point map
create_spatial_map(study_sites, fill_variable = "ndvi_mean")

# Raster map with region boundary
create_spatial_map(ndvi_raster, region_boundary = "Ohio",
                  color_scheme = "ndvi")

# Interactive map (if leaflet available)
create_spatial_map(counties, fill_variable = "population",
                  interactive = TRUE)

## End(Not run)

Create water quality plot

Description

Create water quality plot

Usage

create_water_quality_plot(
  water_data,
  variable,
  region_boundary = NULL,
  river_network = NULL,
  thresholds = NULL,
  output_file = NULL
)

Arguments

Value

Character string with output file path if saved, or NULL if plotted to screen. Called primarily for side effects (creating plots).

Create robust water quality visualization

Description

Create robust water quality visualization

Usage

create_water_quality_plot_robust(
  water_data,
  variable,
  region_boundary,
  river_network,
  threshold_results,
  output_file
)

Detect coordinate columns automatically

Description

Detect coordinate columns automatically

Usage

detect_coordinate_columns(water_df, verbose)

Detect geometry type automatically

Description

Detect geometry type automatically

Usage

detect_geometry_type(spatial_sf)

Detect temporal changes between periods

Description

Internal function to detect changes between time periods.

Usage

detect_temporal_changes(rasters, dates)

Arguments

Value

List with change detection results

Detect vegetation stress

Description

Detect vegetation stress

Usage

detect_vegetation_stress(vegetation_indices, crop_type, indices, verbose)

Detect water quality data columns intelligently

Description

Detect water quality data columns intelligently

Usage

detect_water_quality_columns(
  water_sf,
  variable,
  date_column,
  station_id_col,
  verbose
)

Execute specific interpolation method

Description

Execute specific interpolation method

Usage

execute_interpolation_method(
  complete_data,
  variable,
  method,
  target_grid,
  power,
  max_distance,
  min_points,
  max_points,
  mice_method,
  mice_iterations,
  verbose
)

Execute MICE interpolation

Description

Execute MICE interpolation

Usage

execute_mice_interpolation(
  complete_data,
  variable,
  mice_method,
  mice_iterations,
  verbose
)

Execute nearest neighbor interpolation

Description

Execute nearest neighbor interpolation

Usage

execute_nn_interpolation(complete_data, variable, target_grid, verbose)

Execute simple distance weighting interpolation

Description

Execute simple distance weighting interpolation

Usage

execute_simple_interpolation(
  complete_data,
  variable,
  target_grid,
  power,
  max_distance,
  max_points,
  verbose
)

Execute spline interpolation

Description

Execute spline interpolation

Usage

execute_spline_interpolation(complete_data, variable, target_grid, verbose)

Extract coordinates from ASTER filename

Description

Internal function to extract lat/lon coordinates from ASTER filenames.

Usage

extract_aster_coordinates(filename)

Arguments

Value

List with lat and lon coordinates, or NULL if pattern not matched

Extract bands from multi-band raster

Description

Extract bands from multi-band raster

Usage

extract_bands_from_raster(spectral_data, auto_detect, band_names, verbose)

Extract dates from filenames using various patterns

Description

Universal function to extract dates from filenames or provide custom labels. Enhanced with more flexible regex patterns that work with any filename prefix.

Usage

extract_dates_universal(input_data, date_patterns = NULL, verbose = FALSE)

Arguments

Value

Character vector of extracted or inferred date labels

Examples

## Not run: 
# These examples require external data files not included with the package
# Extract dates from filenames
dates <- extract_dates_universal(c("ndvi_2023-05-15.tif", "evi_2023-06-15.tif"))

# Custom date patterns
custom_patterns <- list("MMDDYYYY" = "\\b[0-9]{2}[0-9]{2}[0-9]{4}\\b")
dates <- extract_dates_universal(files, custom_patterns)

## End(Not run)

Generate enhanced analysis summary

Description

Generate enhanced analysis summary

Usage

generate_enhanced_analysis_summary(results, config)

Generate HTML summary report

Description

Generate HTML summary report

Usage

generate_html_summary(summary_report, results, output_file)

Generate comprehensive statistics report

Description

Generate comprehensive statistics report

Usage

generate_statistics_report(
  statistics,
  threshold_results,
  spatial_analysis,
  temporal_analysis,
  variable,
  region_boundary
)

Get available indices

Description

Get available indices

Usage

get_available_indices()

Get comprehensive CDL crop codes

Description

Get USDA Cropland Data Layer (CDL) codes for specific crops or crop categories. Supports all major crops and predefined categories.

Usage

get_comprehensive_cdl_codes(crop_type = "all")

Arguments

Value

Vector of CDL codes

Examples

# Get corn code
corn_codes <- get_comprehensive_cdl_codes("corn")
print(corn_codes)  # Should be 1

# Get all grain crop codes
grain_codes <- get_comprehensive_cdl_codes("grains")
print(grain_codes)  # Should be vector of grain codes

# See available crop types (this will print to console)
get_comprehensive_cdl_codes("help")

Get index formulas

Description

Get index formulas

Usage

get_index_formulas(indices)

Get index typical ranges

Description

Get index typical ranges

Usage

get_index_ranges(indices)

Get index references

Description

Get index references

Usage

get_index_references(indices)

Get interpretation guidelines

Description

Get interpretation guidelines

Usage

get_interpretation_guidelines(indices)

Get region boundary for any specified region

Description

Universal function to get region boundaries for any geographic area including US states, countries, CONUS, counties, or custom bounding boxes with comprehensive error handling.

Usage

get_region_boundary(region_def, verbose = FALSE)

Arguments

Value

sf object with boundary geometry

Examples

# US State with error handling
ohio_boundary <- get_region_boundary("Ohio")

# Custom bounding box with validation
custom_area <- get_region_boundary(c(-84.5, 39.0, -82.0, 41.0))

Get reliable colors for terra plotting

Description

Get reliable colors for terra plotting

Usage

get_reliable_colors(color_scheme, n = 100)

Get satellite band information

Description

Get satellite band information

Usage

get_satellite_band_info(indices)

Get summary function for terra operations

Description

Get summary function for terra operations

Usage

get_summary_function(summary_function)

Arguments

Value

Function object

Get terra colors for plotting

Description

Get terra colors for plotting

Usage

get_terra_colors(scheme = "viridis", n = 100)

Get water index formulas

Description

Get water index formulas

Usage

get_water_index_formulas(indices)

Get water index requirements

Description

Get water index requirements

Usage

get_water_index_requirements()

Handle edge cases for index calculations

Description

Handle edge cases for index calculations

Usage

handle_index_edge_cases(index, index_type, verbose)

Handle outliers in spatial data

Description

Handle outliers in spatial data

Usage

handle_outliers_in_data(
  spatial_sf,
  target_variables,
  method,
  threshold,
  verbose
)

Handle edge cases for water index calculations

Description

Handle edge cases for water index calculations

Usage

handle_water_index_edge_cases(index, index_type, verbose)

Integrate multiple datasets

Description

Integrate multiple datasets

Usage

integrate_multiple_datasets(
  vector_data,
  raster_datasets,
  region_boundary = NULL,
  output_format = "sf"
)

Integrate terrain analysis with vector data

Description

Specialized function for terrain analysis integration. Calculates terrain variables from DEM and extracts values to vector data points/polygons.

Usage

integrate_terrain_analysis(
  vector_data,
  elevation_raster,
  terrain_vars = c("slope", "aspect", "TRI", "TPI", "flowdir"),
  custom_terrain_functions = NULL,
  extraction_method = "simple"
)

Arguments

Value

sf object with terrain attributes

Examples

## Not run: 
# These examples require external data files not included with the package
# Extract terrain variables for study sites
sites_with_terrain <- integrate_terrain_analysis(
  vector_data = "study_sites.shp",
  elevation_raster = "dem.tif",
  terrain_vars = c("slope", "aspect", "TRI", "TPI")
)

# Use custom terrain functions
custom_functions <- list(
  ruggedness = function(sf_data) {
    sf_data$slope * sf_data$TRI
  }
)

terrain_analysis <- integrate_terrain_analysis(
  vector_data = field_boundaries,
  elevation_raster = dem_raster,
  custom_terrain_functions = custom_functions
)

## End(Not run)

Validate date string

Description

Internal function to check if a date string is valid.

Usage

is_valid_date(date_str)

Arguments

Value

Logical indicating if date is valid

Get comprehensive list of available vegetation indices

Description

Returns detailed information about all 40+ available vegetation indices including formulas, required bands, applications, and references.

Usage

list_vegetation_indices(
  category = "all",
  application = "all",
  detailed = FALSE
)

Arguments

Value

Data frame with vegetation index information

Examples

# All available indices
all_indices <- list_vegetation_indices()

# Only stress detection indices
stress_indices <- list_vegetation_indices(category = "stress")

# Detailed information with formulas
detailed_info <- list_vegetation_indices(detailed = TRUE)

# Agricultural applications only
ag_indices <- list_vegetation_indices(application = "agriculture")

Get comprehensive list of available water indices

Description

Returns detailed information about all available water indices including formulas, required bands, applications, and interpretation guidelines.

Usage

list_water_indices(detailed = FALSE, application_filter = "all")

Arguments

Value

Data frame with water index information

Examples

# All available water indices
water_indices <- list_water_indices()

# Detailed information with formulas
detailed_info <- list_water_indices(detailed = TRUE)

# Only water detection indices
water_detection <- list_water_indices(application_filter = "water_detection")

Load and process spatial data for interpolation

Description

Load and process spatial data for interpolation

Usage

load_and_process_spatial_data(spatial_data, coord_cols, verbose)

Load and stack individual band files

Description

Load and stack individual band files

Usage

load_and_stack_bands(file_list, band_names = NULL, verbose = FALSE)

Load and validate spectral band

Description

Load and validate spectral band

Usage

load_and_validate_band(band_data, band_name, required = FALSE)

Load bands from directory

Description

Load bands from directory

Usage

load_bands_from_directory(directory, band_names = NULL, verbose = FALSE)

Load raster data from various sources

Description

Universal function to load raster data from files, directories, or raster objects with comprehensive error handling and validation.

Usage

load_raster_data(
  input_data,
  pattern = "\\.(tif|tiff)$",
  recursive = FALSE,
  verbose = FALSE
)

Arguments

Value

List of terra SpatRaster objects

Examples

## Not run: 
# These examples require directory structures with multiple data files
# Load from directory with error handling
rasters <- load_raster_data("/path/to/raster/files")

# Load from file list with validation
rasters <- load_raster_data(c("file1.tif", "file2.tif"))

## End(Not run)

Load raster with error handling

Description

Load raster with error handling

Usage

load_raster_safe(x, param_name = "raster")

Arguments

Value

SpatRaster object

Load single water quality dataset

Description

Load single water quality dataset

Usage

load_single_water_dataset(data_input, coord_cols, verbose)

Load vector data safely with coordinate detection

Description

Load vector data safely with coordinate detection

Usage

load_vector_data_safe(file_path)

Arguments

Value

sf object

Load water quality data with flexible format handling

Description

Load water quality data with flexible format handling

Usage

load_water_quality_data_flexible(water_data, coord_cols, verbose)

Map custom band names

Description

Map custom band names

Usage

map_custom_band_names(spectral_data, band_names, verbose = FALSE)

Mask invalid values based on index type

Description

Mask invalid values based on index type

Usage

mask_invalid_values(index, index_type, verbose)

Mask invalid values for water indices

Description

Mask invalid values for water indices

Usage

mask_invalid_water_values(index, index_type, verbose)

Match rasters by date

Description

Match rasters by date

Usage

match_rasters_by_date(
  red_data,
  nir_data,
  date_patterns = NULL,
  verbose = FALSE
)

Multi-scale spatial operations

Description

Handle multi-scale operations including up-scaling, down-scaling, and pyramid operations for efficient processing.

Usage

multiscale_operations(
  spatial_data,
  target_scales = c(1, 2, 4, 8),
  operation = "mean",
  pyramid = FALSE
)

Arguments

Value

List of results at different scales

Examples

## Not run: 
# These examples require external data files not included with the package
# Create multi-scale analysis
scales <- multiscale_operations("data.tif", c(1, 2, 4, 8), "mean")

## End(Not run)

Perform cross-validation for interpolation accuracy

Description

Perform cross-validation for interpolation accuracy

Usage

perform_cross_validation(
  spatial_sf,
  target_variables,
  method,
  cv_folds,
  power,
  max_distance,
  min_points,
  max_points,
  verbose
)

Perform extract join (vector to raster)

Description

Perform extract join (vector to raster)

Usage

perform_extract_join(
  source_info,
  target_info,
  summary_function,
  buffer_distance,
  crs_target,
  na_strategy,
  verbose
)

Perform other join methods (stubs for now)

Description

Perform other join methods (stubs for now)

Usage

perform_overlay_join(
  source_info,
  target_info,
  summary_function,
  crs_target,
  verbose
)

Perform resample join (raster to raster)

Description

Perform resample join (raster to raster)

Usage

perform_resample_join(
  source_info,
  target_info,
  scale_factor,
  summary_function,
  crs_target,
  verbose
)

Perform threshold analysis

Description

Perform threshold analysis

Usage

perform_threshold_analysis(water_sf, variable, thresholds, verbose)

Create fast raster plot using terra

Description

Create efficient raster plots using terra's native plotting capabilities. Fast and reliable without external dependencies.

Usage

plot_raster_fast(
  raster_data,
  title = "Raster Plot",
  color_scheme = "viridis",
  region_boundary = NULL,
  breaks = NULL,
  output_file = NULL,
  verbose = FALSE
)

Arguments

Value

NULL (plots directly to device) or file path if saved

Examples

## Not run: 
# These examples demonstrate workflows with user's own spatial data
# Simple raster plot
plot_raster_fast(ndvi_raster, "NDVI Analysis", "ndvi")

# With custom breaks and save to file
plot_raster_fast(elevation, "Elevation", "terrain",
                breaks = c(0, 500, 1000, 1500, 2000),
                output_file = "elevation_map.png")

## End(Not run)

Create multi-band raster RGB plot

Description

Create RGB plots from multi-band rasters using terra's native RGB plotting. Reliable and fast without external dependencies.

Usage

plot_rgb_raster(
  raster_data,
  r = 1,
  g = 2,
  b = 3,
  stretch = "lin",
  title = "RGB Composite",
  output_file = NULL,
  verbose = FALSE
)

Arguments

Value

NULL (plots directly to device) or file path if saved

Examples

## Not run: 
# These examples require external data files not included with the package
# True color composite
plot_rgb_raster(satellite_data, r = 3, g = 2, b = 1, title = "True Color")

# False color composite with histogram stretch
plot_rgb_raster(landsat_data, r = 4, g = 3, b = 2, stretch = "hist",
               title = "False Color Composite", output_file = "rgb_composite.png")

## End(Not run)

Process spatial data safely with error handling

Description

Process spatial data safely with error handling

Usage

process_spatial_data_safe(spatial_data, coord_cols, verbose)

Process vector data from data frame

Description

Process vector data from data frame

Process vector data from data frame

Usage

process_vector_data(data_input, coord_cols = c("lon", "lat"))

process_vector_data(data_input, coord_cols = c("lon", "lat"))

Arguments

Value

sf object

Quick diagnostic check

Description

Quick diagnostic to identify what might be wrong with the package.

Usage

quick_diagnostic()

Value

List containing diagnostic results with components:

r_version

Character string of R version

minimal_works

Logical indicating if basic functionality works

function_status

Logical indicating function availability status

Called primarily for side effects (printing diagnostic messages).

Quick map function - one-line mapping with auto-detection

Description

Ultra-simple function for quick spatial mapping. Auto-detects data type and creates appropriate map.

Usage

quick_map(spatial_data, variable = NULL, title = NULL, ...)

Arguments

Value

Map object

Examples

## Not run: 
# These examples require external data files not included with the package
quick_map("data.shp")
quick_map(my_raster)
quick_map(points_data, interactive = TRUE)

## End(Not run)

Raster to Raster Operations

Description

Specialized function for mathematical and overlay operations between rasters. Handles alignment, projection, and complex operations with comprehensive error handling and performance optimization.

Usage

raster_to_raster_ops(
  raster1,
  raster2,
  operation = "overlay",
  align_method = "resample",
  summary_function = "mean",
  handle_na = "propagate",
  mask_value = NA,
  output_file = NULL,
  verbose = FALSE
)

Arguments

Value

SpatRaster with operation results

Examples

## Not run: 
# These examples require external data files not included with the package
# Mathematical operations
sum_raster <- raster_to_raster_ops("ndvi.tif", "evi.tif", "add")
diff_raster <- raster_to_raster_ops("before.tif", "after.tif", "subtract")

# Masking operations
masked <- raster_to_raster_ops("data.tif", "mask.tif", "mask")

# Complex overlay with alignment
overlay <- raster_to_raster_ops(
  raster1 = "fine_res.tif",
  raster2 = "coarse_res.tif",
  operation = "overlay",
  align_method = "resample",
  summary_function = "mean",
  verbose = TRUE
)

## End(Not run)

Run comprehensive geospatial workflow -

Description

Execute complete geospatial analysis workflows with simplified visualization. to handle test cases and provide robust error handling without complex dependencies.

Usage

run_comprehensive_geospatial_workflow(analysis_config)

Arguments

Value

List containing analysis results, visualizations, summary, and configuration

Examples

## Not run: 
# These examples require external data files not included with the package
# Simple NDVI crop analysis workflow
config <- list(
  analysis_type = "ndvi_crop_analysis",
  input_data = list(red = red_raster, nir = nir_raster),
  region_boundary = "Ohio",
  output_folder = "results/"
)
results <- run_comprehensive_geospatial_workflow(config)

## End(Not run)

Run comprehensive vegetation analysis workflow -

Description

Complete vegetation analysis using multiple indices with simplified processing. for reliability and test compatibility.

Usage

run_comprehensive_vegetation_workflow(config, output_folder = tempdir())

Arguments

Value

Comprehensive vegetation analysis results

Run enhanced mosaic workflow

Description

Run enhanced mosaic workflow

Run enhanced mosaic workflow - IMPLEMENTED

Usage

run_enhanced_mosaic_workflow(config, output_folder)

run_enhanced_mosaic_workflow(config, output_folder)

Run enhanced NDVI crop analysis workflow -

Description

Enhanced NDVI workflow with quality filtering, temporal analysis, and visualization. to handle test scenarios and provide robust error handling.

Usage

run_enhanced_ndvi_crop_workflow(config, output_folder = tempdir())

Arguments

Value

List with enhanced NDVI results and visualizations

Run enhanced temporal workflow

Description

Run enhanced temporal workflow

Run enhanced temporal workflow - IMPLEMENTED

Usage

run_enhanced_temporal_workflow(config, output_folder)

run_enhanced_temporal_workflow(config, output_folder)

Run enhanced terrain analysis workflow

Description

Run enhanced terrain analysis workflow

Run enhanced terrain analysis workflow - IMPLEMENTED

Usage

run_enhanced_terrain_analysis_workflow(config, output_folder)

run_enhanced_terrain_analysis_workflow(config, output_folder)

Run enhanced water quality analysis workflow

Description

Enhanced water quality workflow with simplified visualization.

Usage

run_enhanced_water_quality_workflow(config, output_folder = tempdir())

run_enhanced_water_quality_workflow(config, output_folder = tempdir())

Arguments

Value

Enhanced water quality analysis results

Run interactive mapping workflow

Description

Create interactive mapping workflow with multiple data types and layers. Simplified for reliability.

Usage

run_interactive_mapping_workflow(config, output_folder = tempdir())

Arguments

Value

Interactive mapping results

Run multi-dataset workflow - IMPLEMENTED

Description

Run multi-dataset workflow - IMPLEMENTED

Usage

run_multi_dataset_workflow(config, output_folder)

Save interactive map to file

Description

Save interactive map to file

Usage

save_interactive_map(map_object, output_file, verbose)

Save interpolation results to file

Description

Save interpolation results to file

Usage

save_interpolation_results(results, output_file, output_format, verbose)

Save plot to file with error handling

Description

Save plot to file with error handling

Usage

save_plot_to_file(
  plot_function,
  output_file,
  verbose,
  width = 1200,
  height = 800
)

Save static map with ggplot2

Description

Save static map with ggplot2

Usage

save_static_map(plot_object, output_file, verbose)

Save temporal analysis results

Description

Internal function to save temporal analysis results to files.

Usage

save_temporal_results(result, output_folder, analysis_type)

Arguments

Save vegetation analysis results

Description

Save vegetation analysis results

Usage

save_vegetation_analysis_results(
  analysis_results,
  vegetation_indices,
  output_folder,
  crop_type,
  verbose = FALSE
)

Save water quality analysis results

Description

Save water quality analysis results

Usage

save_water_quality_results(
  water_sf,
  variable,
  statistics,
  threshold_results,
  spatial_analysis,
  temporal_analysis,
  region_boundary,
  river_network,
  output_folder,
  verbose
)

Select appropriate indices for crop analysis

Description

Select appropriate indices for crop analysis

Usage

select_crop_indices(crop_type, analysis_type, growth_stage)

Select optimal interpolation method based on data characteristics

Description

Select optimal interpolation method based on data characteristics

Usage

select_optimal_method(data_assessment, verbose)

Select rasters for specific region with intelligent filtering

Description

Intelligently select raster files that overlap with a specified region.

Usage

select_rasters_for_region(input_folder, region_boundary, buffer_size = 0.1)

Arguments

Value

Character vector of relevant file paths

Examples

# Select ASTER files for Michigan
michigan_files <- select_rasters_for_region("/aster/files", "Michigan")

# Select with custom buffer
nevada_files <- select_rasters_for_region("/data", "Nevada", buffer_size = 0.2)

Legacy spatial interpolation function (for backward compatibility)

Description

Simplified version of spatial interpolation maintaining backward compatibility. For new projects, use spatial_interpolation_comprehensive() instead.

Usage

spatial_interpolation(
  spatial_data,
  target_variables,
  method = "NN",
  power = 2,
  mice_method = "pmm"
)

Arguments

Value

sf object with interpolated values

Examples

## Not run: 
# These examples require external data files not included with the package
# Simple interpolation (legacy interface)
interpolated_data <- spatial_interpolation(
  soil_data,
  target_variables = c("nitrogen", "carbon"),
  method = "NN"
)

## End(Not run)

Perform spatial interpolation for missing data

Description

Perform spatial interpolation using reliable methods to fill missing values in spatial datasets. Supports nearest neighbor, spline interpolation, and multivariate imputation with comprehensive error handling.

Usage

spatial_interpolation_comprehensive(
  spatial_data,
  target_variables,
  method = "NN",
  target_grid = NULL,
  region_boundary = NULL,
  power = 2,
  max_distance = Inf,
  min_points = 3,
  max_points = 50,
  cross_validation = FALSE,
  cv_folds = 5,
  handle_outliers = "none",
  outlier_threshold = 3,
  coord_cols = c("lon", "lat"),
  mice_method = "pmm",
  mice_iterations = 10,
  output_format = "sf",
  output_file = NULL,
  verbose = FALSE
)

Arguments

Details

Supported Interpolation Methods:

Distance-Based Methods:

• NN (Nearest Neighbor): Assigns nearest known value - Best for: Categorical data or when preserving exact values - Fast and creates Voronoi-like patterns - No assumptions about data distribution

• Simple (Simple distance weighting): Basic distance-based averaging - Best for: Quick estimates with minimal computation - Uses inverse distance weighting without external dependencies

Statistical Methods:

• Spline: Smooth surface interpolation using thin plate splines - Best for: Smooth, continuous phenomena - Creates smooth surfaces without sharp changes - Good for environmental data with gradual spatial variation

Multivariate Methods:

• MICE: Multivariate imputation by chained equations - Best for: Multiple correlated variables with missing values - Handles complex missing data patterns - Preserves relationships between variables - Requires mice package

Input Data Support:

• sf objects with point geometries

• data.frame with coordinate columns

• File paths (CSV, shapefile, GeoJSON)

• Target grids for raster output

Quality Control Features:

• Cross-validation for method comparison

• Outlier detection and handling

• Performance metrics calculation

• Robust error handling

Value

Depending on output_format:

"sf"

sf object with interpolated values

"raster"

SpatRaster with interpolated surfaces

"both"

List containing both sf and raster results

Additional attributes include:

• interpolation_info: Method used, parameters, processing time

• cross_validation: CV results if performed

Method Selection Guide

Dense, regular data

Simple distance weighting for good balance

Sparse, irregular data

Nearest neighbor for stability

Environmental data

Spline for smooth surfaces

Categorical data

Nearest neighbor

Multiple correlated variables

MICE for multivariate patterns

Unknown data characteristics

Auto-selection based on data properties

Performance Optimization

• For large datasets: Set max_points=50-100 for faster processing

• For high accuracy: Use cross_validation=TRUE to compare methods

• For memory efficiency: Process variables individually

• For smooth results: Use spline method

See Also

• universal_spatial_join for spatial data integration

• calculate_spatial_correlation for spatial correlation analysis

• create_spatial_map for visualization

Examples

## Not run: 
# These examples require external data files not included with the package
# Basic nearest neighbor interpolation
soil_interpolated <- spatial_interpolation_comprehensive(
  spatial_data = "soil_samples.csv",
  target_variables = c("nitrogen", "phosphorus", "ph"),
  method = "NN",
  target_grid = study_area_grid,
  region_boundary = "Iowa"
)

# Simple distance weighting
temp_interp <- spatial_interpolation_comprehensive(
  spatial_data = weather_stations,
  target_variables = "temperature",
  method = "simple",
  power = 2,
  cross_validation = TRUE,
  verbose = TRUE
)

# Multivariate imputation for environmental data
env_imputed <- spatial_interpolation_comprehensive(
  spatial_data = env_monitoring,
  target_variables = c("temp", "humidity", "pressure", "wind_speed"),
  method = "mice",
  mice_iterations = 15,
  handle_outliers = "cap"
)

# Auto-method selection with comparison
best_interp <- spatial_interpolation_comprehensive(
  spatial_data = precipitation_data,
  target_variables = "annual_precip",
  method = "auto",
  cross_validation = TRUE,
  cv_folds = 10,
  target_grid = dem_template
)

# Access results and diagnostics
plot(best_interp)  # Plot interpolated surface
best_interp$cross_validation$rmse  # Cross-validation RMSE
best_interp$interpolation_info$method_selected  # Method chosen

## End(Not run)

Test individual function existence

Description

Helper function to test if core functions exist and are callable.

Usage

test_function_availability(verbose = FALSE)

Arguments

Value

List of function availability

Test GeoSpatialSuite with simplified, robust tests

Description

Simplified testing function that focuses on core functionality with minimal dependencies and robust error handling. Designed for 100% success rate. This replaces the complex testing function with simple, reliable tests.

Usage

test_geospatialsuite_package_simple(
  test_output_dir = tempdir(),
  verbose = FALSE
)

Arguments

Value

List of test results with success/failure status for each component

Examples

# Quick test (essential functions only)
test_results <- test_geospatialsuite_package_simple()

# Verbose test
test_results <- test_geospatialsuite_package_simple(verbose = TRUE)

Test package with minimal complexity

Description

Ultra-minimal test that only checks the most basic functionality. Designed to always pass if the package is minimally functional.

Usage

test_package_minimal(verbose = FALSE)

Arguments

Value

Logical indicating basic functionality

Universal Spatial Join - Complete Implementation

Description

Comprehensive spatial join system that handles ALL spatial data combinations: Vector to Vector, Vector to Raster, Raster to Raster with full documentation, error handling, and extensive examples. This replaces all previous spatial join functions with a unified, robust system.

Usage

universal_spatial_join(
  source_data,
  target_data,
  method = "auto",
  scale_factor = NULL,
  summary_function = "mean",
  buffer_distance = NULL,
  temporal_tolerance = NULL,
  crs_target = NULL,
  na_strategy = "remove",
  chunk_size = 1e+06,
  parallel = FALSE,
  verbose = FALSE
)

Arguments

Details

Quick Start Guide:

Most common use case - extract raster values to point locations:

result <- universal_spatial_join("my_points.csv", "my_raster.tif", method="extract")

Supported Operations:

Data Type Combinations:

• Vector → Raster: Extract raster values to points/polygons/lines

• Raster → Vector: Calculate zonal statistics for polygons

• Raster → Raster: Resample, overlay, mathematical operations

• Vector → Vector: Spatial intersections, overlays, nearest neighbor

Input Format Support:

• File paths: ".tif", ".shp", ".gpkg", ".geojson", ".nc"

• Directories: Automatically processes all spatial files

• R objects: SpatRaster, sf, data.frame with coordinates

• Lists: Multiple files or raster stacks

Scaling Operations:

• Up-scaling: Aggregate to coarser resolution (scale_factor > 1)

• Down-scaling: Interpolate to finer resolution (scale_factor < 1)

• Custom resolution: Match target raster geometry

Error Handling:

• Auto CRS reprojection: Handles coordinate system mismatches

• Geometry alignment: Auto-crops, extends, or resamples as needed

• NA handling: Multiple strategies for missing data

• Memory management: Chunked processing for large datasets

Method Selection Guide:

extract

Use when you have point/polygon locations and want to get values from a raster

zonal

Use when you have polygons and want statistics from raster data within each polygon

resample

Use when you need to change raster resolution or align two rasters

overlay

Use when joining two vector datasets based on spatial relationships

nearest

Use when you want to find the closest features between two vector datasets

auto

Let the function choose - works well for standard extract/resample operations

Value

Spatial data object with joined attributes. Return type depends on operation:

extract (vector → raster)

sf object with new columns containing extracted raster values. Original geometry preserved, new columns named "extracted_" followed by the raster layer name

zonal (raster → vector)

sf object with new columns containing zonal statistics. Original geometry preserved, new columns named "zonal_" followed by the statistic name and raster layer name

resample (raster → raster)

SpatRaster with resampled/processed data matching target resolution or scale factor

overlay (vector → vector)

sf object with intersected/overlaid features combining attributes from both datasets

nearest

sf object with attributes from nearest features joined

Returned objects include 'spatial_join_info' attribute containing:

• method: Join method used

• source_type, target_type: Data types processed

• processing_time: Time taken (if verbose=TRUE)

• timestamp: Processing timestamp

• summary_function: Aggregation function used

Common Error Solutions

CRS Mismatch

"CRS mismatch detected" - Function automatically reprojects data, but manual CRS checking recommended for precision

Memory Issues

"Large dataset processing" - Reduce chunk_size parameter (try 500000) or set parallel=FALSE

No Spatial Overlap

"No spatial overlap found" - Check that source and target data cover the same geographic area

File Not Found

"File does not exist" - Verify file paths and ensure files exist at specified locations

Missing Bands

"Required bands not found" - For raster operations, ensure expected spectral bands are present

Invalid Geometries

"Geometry errors" - Function attempts to fix automatically, but check input data quality

Performance Tips

• For large datasets (>1M cells): set chunk_size=500000 and parallel=TRUE

• Use method="resample" with scale_factor > 1 to reduce data size before complex operations

• For time series analysis: consider temporal_tolerance to balance accuracy vs processing speed

• When processing multiple datasets: ensure consistent CRS to avoid reprojection overhead

• For point extraction: use smaller buffer_distance when possible to reduce processing time

See Also

• raster_to_raster_ops for specialized raster operations

• multiscale_operations for multi-scale analysis

• process_vector_data for vector data preprocessing

Examples

## Not run: 
# These examples require satellite imagery files (Landsat/Sentinel data etc.)
# =================================================================
# MOST COMMON USE CASE: Extract raster values to CSV points
# =================================================================

# Your typical workflow: CSV file with coordinates + raster file
results <- universal_spatial_join(
  source_data = "my_field_sites.csv",    # CSV with lon, lat columns
  target_data = "satellite_image.tif",   # Any raster file
  method = "extract",                     # Extract raster values to points
  buffer_distance = 100,                  # 100m buffer around each point
  summary_function = "mean",              # Average within buffer
  verbose = TRUE                          # See what's happening
)

# Check results - original data + new columns with raster values
head(results)
#   site_id       lon       lat           geometry extracted_satellite_image
# 1       1 -83.12345  40.12345 POINT (-83.1 40.1)                   0.752
# 2       2 -83.23456  40.23456 POINT (-83.2 40.2)                   0.681
# 3       3 -83.34567  40.34567 POINT (-83.3 40.3)                   0.594

# Access the extracted values
results$extracted_satellite_image

# =================================================================
# ZONAL STATISTICS: Calculate statistics by polygon areas
# =================================================================

# Calculate average precipitation by watershed
watershed_precip <- universal_spatial_join(
  source_data = "precipitation_raster.tif",  # Raster data
  target_data = "watershed_boundaries.shp",  # Polygon boundaries
  method = "zonal",                           # Calculate zonal statistics
  summary_function = "mean",                  # Average precipitation per watershed
  verbose = TRUE
)

# Result: polygons with precipitation statistics
head(watershed_precip)
#   watershed_id                     geometry zonal_mean_precipitation_raster
# 1            1 POLYGON ((-84.2 40.1, ...))                             42.3
# 2            2 POLYGON ((-84.5 40.3, ...))                             38.7

# =================================================================
# RESAMPLE RASTER: Change resolution or align rasters
# =================================================================

# Aggregate 30m Landsat to 250m MODIS resolution
landsat_resampled <- universal_spatial_join(
  source_data = "landsat_30m.tif",      # High resolution input
  target_data = "modis_250m.tif",       # Target resolution template
  method = "resample",                   # Resample operation
  summary_function = "mean",             # Average when aggregating
  verbose = TRUE
)

# Check new resolution
terra::res(landsat_resampled)
# [1] 250 250

# Scale by factor instead of template
coarser_raster <- universal_spatial_join(
  source_data = "fine_resolution.tif",
  target_data = NULL,                    # No template needed
  method = "resample",
  scale_factor = 5,                      # 5x coarser resolution
  summary_function = "mean"
)

# =================================================================
# VECTOR OVERLAY: Join two vector datasets
# =================================================================

# Find which counties contain each field site
sites_with_counties <- universal_spatial_join(
  source_data = "field_sites.shp",      # Point data
  target_data = "county_boundaries.shp", # Polygon data
  method = "overlay",                     # Spatial intersection
  verbose = TRUE
)

# Result: points with county attributes added
head(sites_with_counties)
#   site_id           geometry county_name state_name
# 1       1 POINT (-83.1 40.1)    Franklin       Ohio
# 2       2 POINT (-83.2 40.2)    Delaware       Ohio

# =================================================================
# AUTO-DETECTION: Let function choose best method
# =================================================================

# Function automatically detects: points + raster = extract method
auto_result <- universal_spatial_join(
  source_data = my_points,               # Any point data
  target_data = my_raster,               # Any raster data
  method = "auto",                       # Automatically choose method
  verbose = TRUE                         # See what method was chosen
)
# Output: "Auto-detected method: extract for vector to raster"

# =================================================================
# ERROR HANDLING EXAMPLES
# =================================================================

# Function handles common issues automatically
robust_result <- universal_spatial_join(
  source_data = "points_wgs84.csv",     # WGS84 coordinate system
  target_data = "raster_utm.tif",       # UTM coordinate system
  method = "extract",
  na_strategy = "nearest",               # Handle missing values
  verbose = TRUE                         # See CRS handling messages
)
# Output: "CRS mismatch detected. Reprojecting to match raster CRS..."

## End(Not run)

Validate method parameter

Description

Validate method parameter

Usage

validate_method(method, valid_methods, param_name = "method")

Arguments

Value

TRUE if valid, stops with error if invalid

Validate method compatibility

Description

Validate method compatibility

Usage

validate_method_compatibility(method, source_type, target_type)

Arguments

Validate numeric range

Description

Validate numeric range

Usage

validate_numeric_range(
  x,
  min_val = NULL,
  max_val = NULL,
  param_name = "parameter",
  allow_null = FALSE
)

Arguments

Value

TRUE if valid, stops with error if invalid

Validate output quality

Description

Validate output quality

Usage

validate_output(index, index_type, verbose)

Validate raster input

Description

Validate raster input

Usage

validate_raster_input(x, param_name = "raster_data", allow_null = FALSE)

Arguments

Value

TRUE if valid, stops with error if invalid

Validate required bands for specific indices

Description

Validate required bands for specific indices

Usage

validate_required_bands(
  index_type,
  blue,
  green,
  swir1,
  swir2,
  red_edge,
  coastal,
  nir2,
  verbose = FALSE
)

Validate vector input

Description

Validate vector input

Usage

validate_vector_input(x, param_name = "vector_data", allow_null = FALSE)

Arguments

Value

TRUE if valid, stops with error if invalid

Validate vegetation analysis

Description

Validate vegetation analysis

Usage

validate_vegetation_analysis(
  vegetation_indices,
  reference_data,
  verbose = FALSE
)

Validate water index output

Description

Validate water index output

Usage

validate_water_index_output(index, index_type, verbose)

Input Validation Helpers

Description

Internal validation functions for robust parameter checking