if (FALSE) {
# 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
}
Run the code above in your browser using DataLab