Bin continuous depth values into discrete intervals
Examples
data(seg_Oulanka2023_Session01_T067)
img = terra::rast(seg_Oulanka2023_Session01_T067)
mask = img[[1]] - img[[2]]
mask[mask == 255] <- NA
img = img
depthmap = create_depthmap(img,mask,start.soil = 2.9, select.layer = 2 )
#> Processing row 100 of 4900
#> Processing row 200 of 4900
#> Processing row 300 of 4900
#> Processing row 400 of 4900
#> Processing row 500 of 4900
#> Processing row 600 of 4900
#> Processing row 700 of 4900
#> Processing row 800 of 4900
#> Processing row 900 of 4900
#> Processing row 1000 of 4900
#> Processing row 1100 of 4900
#> Processing row 1200 of 4900
#> Processing row 1300 of 4900
#> Processing row 1400 of 4900
#> Processing row 1500 of 4900
#> Processing row 1600 of 4900
#> Processing row 1700 of 4900
#> Processing row 1800 of 4900
#> Processing row 1900 of 4900
#> Processing row 2000 of 4900
#> Processing row 2100 of 4900
#> Processing row 2200 of 4900
#> Processing row 2300 of 4900
#> Processing row 2400 of 4900
#> Processing row 2500 of 4900
#> Processing row 2600 of 4900
#> Processing row 2700 of 4900
#> Processing row 2800 of 4900
#> Processing row 2900 of 4900
#> Processing row 3000 of 4900
#> Processing row 3100 of 4900
#> Processing row 3200 of 4900
#> Processing row 3300 of 4900
#> Processing row 3400 of 4900
#> Processing row 3500 of 4900
#> Processing row 3600 of 4900
#> Processing row 3700 of 4900
#> Processing row 3800 of 4900
#> Processing row 3900 of 4900
#> Processing row 4000 of 4900
#> Processing row 4100 of 4900
#> Processing row 4200 of 4900
#> Processing row 4300 of 4900
#> Processing row 4400 of 4900
#> Processing row 4500 of 4900
#> Processing row 4600 of 4900
#> Processing row 4700 of 4900
#> Processing row 4800 of 4900
#> Processing row 4900 of 4900
binned.map = binning(depthmap,nn = 5)