data-raw/23b-controllo_szn_senza_omi.R

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

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

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'))