data-raw/23c-controllo_omi_senza_szn.R

In master-info/masteRomi: Dati, Poligoni e altro per analisi di Zone OMI

#################################
# CONTROLLO "OMI SENZA SEZIONI" #
#################################

masteRfun::load_pkgs('data.table', 'sf')
load_all()

yz <- readRDS(file.path(bnd_path, 'SZN', 's20', '0'))
yo.nz <- zone[is.na(termine)][!OMI %in% szn_omi[N<=1, OMI]]
yb.nz <- OMI |> subset(OMI %in% yo.nz$OMI) |> st_transform(4326)

yy <- yz |> subset(!SZN %in% rbindlist(list(yc1, yc2t[N==1]))[, SZN]) |> st_transform(4326)
ynac <- masteRgeo::comuni[CMN %in% unique(masteRgeo::sezioni[SZN %in% yy$SZN, CMN]), .(CMN, catasto)]

leaflet() |> 
    addTiles() |> 
    addPolygons(group = 'sezioni', data = yy, label = ~SZN) |> 
    addPolygons(
        group = 'comuni',
        data = masteRconfini::CMN |> subset(CMN %in% ynac$CMN), 
        color = 'red', 
        fillOpacity = 0, 
        label = ~CMN
    ) |> 
    addPolygons(
        group = 'zone OMI',
        data = yb |> subset(OMI %in% zone[catasto %chin% ynac$catasto, OMI]) |> st_transform(4326), 
        color = 'black', 
        fillOpacity = 0,
        label = ~OMI
    ) |> 
    addLayersControl(overlayGroups = c('sezioni', 'comuni', 'zone OMI'))

x <- 'A06000R3'
yb.nz.x <- yb.nz |> subset(OMI == x)
yzx <- yz |> st_join(yb.nz.x, join = st_intersects) |> subset(!is.na(OMI))
yozx <- OMI |> subset(OMI %in% unique(szn_omi[SZN %in% yzx$SZN & N <= 1, OMI]))
gx <- paste0('OMI', x)
leaflet() |> 
    addProviderTiles(providers$Stamen.TonerLite) |> 
    addPolygons(data = yb.nz.x, group = gx) |>
    addPolygons(data = yzx, group = 'SZN', color = 'red', weight = 2, fillOpacity = 0, label = ~SZN) |>
    addPolygons(data = yozx, group = 'Altri OMI', color = 'black', weight = 2, fillOpacity = 0, label = ~OMI) |> 
    addLayersControl(overlayGroups = c(gx, 'SZN', 'Altri OMI'))