data-raw/22-variazioni_zone_confini.R

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

# VARIAZIONI

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

fn <- max(list.files(file.path(bnd_path, 'OMI'), '20[0-9]'))
yc <- fread('./data-raw/csv/comuni.csv', select = c('catasto', 'CMN'))
yb <- qs::qread(file.path(bnd_path, 'OMI', fn))
yo <- fread('./data-raw/csv/zone.csv', select = c('OMI', 'catasto', 'termine'))
yo <- yc[yo[is.na(termine)], on = 'catasto'][, termine := NULL]

## OMI presenti in `zone` (yo) ma mancanti in `confini` (yb)
yob <- yo[!OMI %chin% yb$OMI]
leaflet() |> 
    addTiles() |> 
    addPolygons(data = yob, label = ~OMI, ) |> 
    addPolygons(data = CMN |> subset(CMN %in% comuni[catasto %chin% unique(substr(ybo$OMI, 1, 4)), CMN]), label = ~CMN, color = 'red', fillOpacity = 0)

## OMI presenti in `confini` (yb) ma mancanti in `zone` (yb)
ybo <- yb |> subset(!OMI %in% yo$OMI) |> st_transform(4326) |> dplyr::mutate(catasto = substr(OMI,1,4)) |> merge(yc, all.x = TRUE)
leaflet() |> 
    addTiles() |> 
    addPolygons(data = ybo, label = ~OMI, ) |> 
    addPolygons(data = CMN |> subset(CMN %in% comuni[catasto %chin% unique(substr(ybo$OMI, 1, 4)), CMN]), label = ~CMN, color = 'red', fillOpacity = 0)

## CMN ISTAT (masteRgeo::comuni) *totalmente* non presenti in zone (yo)
masteRgeo::comuni[!CMN %in% unique(yo$CMN), .(CMN, CMNd, PRVs)]


#=== comuni da aggiungere in zone ===#
######################################

#>>> M432 Misiliscemi 81025, da L331 TRAPANI 81021 <<<# ----
x <- 'M432'
xv <- 'L331'
masteRgeo::comuni[catasto == x]
masteRomi::zone[catasto == xv]
masteRgeo::comuni[catasto == xv]
masteRgeo::zone[codice == 81025]
yb2 <- qs::qread(file.path(bnd_path, 'OMI', 20202))

leaflet() |> 
    addTiles() |> 
    addPolygons(data = yb |> subset(substr(OMI,1,4) == x) |> st_transform(4326), label = ~OMI, color = 'red') |> 
    addPolygons(data = yb2 |> subset(substr(OMI,1,4) == xv) |> st_transform(4326), label = ~OMI, color = 'black') |> 
    addPolygons(data = CMN |> merge(comuni) |> subset(RGNd == 'SICILIA'), label = ~CMNd, fillOpacity = 0)


# aggiungere in tabella `comuni`
#   CMN catasto nazionale    CMNd PRVs    RGNd  RPTd inizio termine  bb_xmin  bb_ymin  bb_xmax  bb_ymax
# 81021    L331      U9AA TRAPANI   TP SICILIA ISOLE  20102    <NA> 12.48797 37.83693 12.80075 38.03455
data.table(81025, 'M432', NA, 'MISILISCEMI', 'TP', 'SICILIA', 'ISOLE', 20211, NA, 12.47942, 37.85658, 12.61547, 37.97099)

# aggiungere in confini  `CMN`

# cambiare in tabella `zone`


# M43200B1 da 
# M43200E1 da 
# M43200R1 da 


#>>> M404 Valle Cannobina 103079, da ? <<<# ----
x <- 'M404'
masteRgeo::comuni[catasto == x]
masteRomi::zone[catasto == 'M']
masteRgeo::comuni[catasto == 'L331']
masteRgeo::zone[codice == 81025]
yb2 <- qs::qread(file.path(bnd_path, 'OMI', 20202))



#>>> M382 Valfornace 43058, da ? <<<# ----
x <- 'M382'
masteRgeo::comuni[catasto == x]
masteRomi::zone[catasto == 'M']
masteRgeo::comuni[catasto == 'L331']
masteRgeo::zone[codice == 81025]
yb2 <- qs::qread(file.path(bnd_path, 'OMI', 20202))



#>>> M309 Fonte Nuova 58122, da ? <<<# ----
x <- 'M309'
masteRgeo::comuni[catasto == x]
leaflet() |> 
    addTiles() |> 
    addPolygons(data = yb |> subset(substr(OMI,1,4) == x) |> st_transform(4326), label = ~OMI, color = 'red') |> 
    addPolygons(data = yb2 |> subset(substr(OMI,1,4) == xv) |> st_transform(4326), label = ~OMI, color = 'black') |> 
    addPolygons(data = CMN |> merge(comuni) |> subset(RGNd == 'SICILIA'), label = ~CMNd, fillOpacity = 0)
masteRomi::zone[catasto == 'M']
masteRgeo::comuni[catasto == 'L331']
masteRgeo::zone[codice == 81025]
yb2 <- qs::qread(file.path(bnd_path, 'OMI', 20202))