if (interactive() && requireNamespace("terra", quietly = TRUE) && terra_is_working()) {
tiff_dir <- system.file("demo-geotiff",
package = "snic",
mustWork = TRUE
)
files <- file.path(
tiff_dir,
c(
"S2_20LMR_B02_20220630.tif",
"S2_20LMR_B04_20220630.tif",
"S2_20LMR_B08_20220630.tif",
"S2_20LMR_B12_20220630.tif"
)
)
s2 <- terra::aggregate(terra::rast(files), fact = 8)
seeds <- snic_grid_manual(
s2,
snic_args = list(compactness = 0.1),
snic_plot_args = list(r = 4, g = 3, b = 1)
)
seg <- snic(s2, seeds, compactness = 0.1)
snic_plot(
s2,
r = 4, g = 3, b = 1,
stretch = "lin",
seeds = seeds,
seg = seg
)
}
Run the code above in your browser using DataLab