Aggiorna versione a 1.4.3: WMS Cartografia Catastale AdE, versione in GUI

Aggiunta opzione per caricare il WMS del Catasto (Agenzia delle Entrate)
con 7 sublayer (province, mappe, acque, strade, vestizioni, fabbricati,
particelle). Connessione WMS salvata nel profilo QGIS. Layer posizionato
sopra eventuali basemap XYZ. Versione plugin nel titolo della finestra.

Co-Authored-By: Claude Opus 4.6 <noreply@anthropic.com>

Author: pigreco

metadata.txt

@@ -6,7 +6,7 @@
 name=WFS Catasto Download Particelle BBox
 qgisMinimumVersion=3.0
 description=Download particelle catastali dal WFS dell'Agenzia delle Entrate tramite BBox, selezione poligono, linea con buffer o layer di punti.
-version=1.4.2
+version=1.4.3
 author=Salvatore Fiandaca
 email=pigrecoinfinito@gmail.com
 
@@ -30,6 +30,11 @@ repository=https://github.com/pigreco/wfs_catasto_download_particelle_bbox
 
 hasProcessingProvider=no
 changelog=
+    1.4.3 - Opzione per caricare il WMS Cartografia Catastale dell'Agenzia delle Entrate
+            (province, mappe, acque, strade, vestizioni, fabbricati, particelle).
+            Connessione WMS aggiunta automaticamente al profilo QGIS se non presente.
+            Layer posizionato sopra eventuali basemap XYZ (es. OpenStreetMap).
+            Versione plugin visualizzata nel titolo della finestra di dialogo.
     1.4.2 - Nuova funzione per il calcolatore di campi: get_particella_info($geometry)
             nel gruppo 'Catasto'. Restituisce un array con riferimento catastale,
             foglio, label, codice comune, geometria WKT, sezione e allegato

wfs_catasto_download_particelle_bbox_d.py

@@ -6,6 +6,7 @@
 Compatibile con QGIS 3 (Qt5) e QGIS 4 (Qt6).
 """
 
+import configparser
 import os
 
 from qgis.PyQt.QtCore import Qt
@@ -44,6 +45,14 @@
     _SmoothTransformation = Qt.SmoothTransformation
 
 
+def _plugin_version():
+    """Legge la versione dal file metadata.txt del plugin."""
+    meta_path = os.path.join(os.path.dirname(__file__), "metadata.txt")
+    cfg = configparser.ConfigParser()
+    cfg.read(meta_path)
+    return cfg.get("general", "version", fallback="")
+
+
 class AvvisoDialog(QDialog):
     """Finestra di avviso sull'uso responsabile del plugin (una volta per sessione QGIS)."""
 
@@ -180,7 +189,11 @@ def __init__(self, parent=None, default_buffer_m=50,
         self._init_ui()
 
     def _init_ui(self):
-        self.setWindowTitle("WFS Catasto - Scelta modalità")
+        ver = _plugin_version()
+        titolo = "WFS Catasto - Scelta modalità"
+        if ver:
+            titolo += f"  v{ver}"
+        self.setWindowTitle(titolo)
         self.setMinimumWidth(500)
         self.setWindowFlags(
             self.windowFlags()
@@ -241,6 +254,12 @@ def _init_ui(self):
         self.check_espandi_catastale.setStyleSheet("font-weight: normal;")
         self.check_espandi_catastale.setChecked(False)
         go_layout.addWidget(self.check_espandi_catastale)
+        self.check_carica_wms = QCheckBox(
+            "Carica WMS Cartografia Catastale"
+        )
+        self.check_carica_wms.setStyleSheet("font-weight: normal;")
+        self.check_carica_wms.setChecked(False)
+        go_layout.addWidget(self.check_carica_wms)
         group_opzioni.setLayout(go_layout)
         layout.addWidget(group_opzioni)
 
@@ -447,6 +466,10 @@ def _svg_label(self, filename, max_w=140, max_h=121):
     def espandi_catastale(self):
         return self.check_espandi_catastale.isChecked()
 
+    @property
+    def carica_wms(self):
+        return self.check_carica_wms.isChecked()
+
     # ---- Slot ----
 
     def _on_disegna(self):

wfs_catasto_download_particelle_bbox_p.py

@@ -21,6 +21,7 @@
     Qgis,
     QgsProject,
     QgsVectorLayer,
+    QgsRasterLayer,
     QgsCoordinateReferenceSystem,
     QgsCoordinateTransform,
     QgsRectangle,
@@ -33,7 +34,7 @@
     QgsExpression,
 )
 from qgis.gui import QgsMapTool, QgsRubberBand
-from qgis.PyQt.QtCore import Qt, QVariant, QTimer
+from qgis.PyQt.QtCore import Qt, QVariant, QTimer, QSettings
 from qgis.PyQt.QtGui import QColor, QIcon
 from qgis.PyQt.QtWidgets import (
     QAction,
@@ -136,6 +137,20 @@ def _is_point_layer(layer):
 # Distanza buffer in metri di default
 BUFFER_DISTANCE_M = 50
 
+# WMS Catasto
+WMS_BASE_URL = "https://wms.cartografia.agenziaentrate.gov.it/inspire/wms/ows01.php"
+WMS_LAYERS = [
+    "province",
+    "CP.CadastralZoning",
+    "acque",
+    "strade",
+    "vestizioni",
+    "fabbricati",
+    "CP.CadastralParcel",
+]
+WMS_CONNECTION_NAME = "Catasto AdE"
+WMS_LAYER_NAME = "Cartografia Catastale WMS AdE"
+
 
 # =============================================================================
 # FUNZIONI COMUNI
@@ -215,6 +230,68 @@ def _determina_utm_epsg(lon, lat):
     return f"EPSG:{epsg}"
 
 
+def carica_wms_catasto():
+    """
+    Aggiunge la connessione WMS del Catasto al profilo QGIS (se non presente)
+    e carica un layer WMS combinato nel progetto corrente.
+    Il layer viene posizionato in fondo al pannello Layer, prima di eventuali basemap XYZ.
+    """
+    # Controlla se il layer WMS è già nel progetto
+    for layer in QgsProject.instance().mapLayers().values():
+        if layer.providerType() == "wms" and WMS_BASE_URL in layer.source():
+            print(f"[WMS] Layer WMS Catasto già presente nel progetto: {layer.name()}")
+            return
+
+    # Aggiungi connessione WMS in QSettings se non presente
+    settings = QSettings()
+    key_prefix = f"qgis/connections-wms/{WMS_CONNECTION_NAME}"
+    existing_url = settings.value(f"{key_prefix}/url", "")
+    if not existing_url:
+        settings.setValue(f"{key_prefix}/url", WMS_BASE_URL)
+        print(f"[WMS] Connessione '{WMS_CONNECTION_NAME}' aggiunta al profilo QGIS")
+    else:
+        print(f"[WMS] Connessione '{WMS_CONNECTION_NAME}' già presente nel profilo")
+
+    # Costruisci URI con tutti i layer combinati
+    layers_params = "&".join(f"layers={l}" for l in WMS_LAYERS)
+    styles_params = "&".join("styles=" for _ in WMS_LAYERS)
+    uri = (
+        f"contextualWMSLegend=0"
+        f"&crs=EPSG:6706"
+        f"&dpiMode=7"
+        f"&featureCount=10"
+        f"&format=image/png"
+        f"&{layers_params}"
+        f"&{styles_params}"
+        f"&url={WMS_BASE_URL}"
+    )
+
+    wms_layer = QgsRasterLayer(uri, WMS_LAYER_NAME, "wms")
+    if not wms_layer.isValid():
+        print("[WMS] ERRORE: impossibile caricare il layer WMS Catasto")
+        return
+
+    QgsProject.instance().addMapLayer(wms_layer, False)
+    root = QgsProject.instance().layerTreeRoot()
+
+    # Inserisci in fondo ma prima di eventuali mappe di sfondo (XYZ tiles)
+    children = root.children()
+    pos = len(children)  # default: ultima posizione
+    for i in range(len(children) - 1, -1, -1):
+        node = children[i]
+        if hasattr(node, "layer") and node.layer():
+            src = node.layer().source()
+            if "type=xyz" in src:
+                pos = i
+            else:
+                break
+        else:
+            break
+    root.insertLayer(pos, wms_layer)
+
+    print(f"[WMS] Layer '{WMS_LAYER_NAME}' caricato nel progetto")
+
+
 def calcola_griglia_tile(min_lat, min_lon, max_lat, max_lon, max_tile_km2):
     """
     Suddivide il bbox in una griglia di tile, ciascuna con area <= max_tile_km2.
@@ -308,7 +385,8 @@ def scarica_singolo_tile(min_lat, min_lon, max_lat, max_lon):
 def esegui_download_e_caricamento(min_lat, min_lon, max_lat, max_lon, filter_geom=None,
                                   layer_name="Particelle WFS",
                                   espandi_catastale=False,
-                                  post_filter_points=None):
+                                  post_filter_points=None,
+                                  carica_wms=False):
     """
     Gestisce il download WFS: singolo o multi-tile con progress bar.
 
@@ -744,6 +822,10 @@ def esegui_download_e_caricamento(min_lat, min_lon, max_lat, max_lon, filter_geo
     canvas.setExtent(extent_proj)
     canvas.refresh()
 
+    # Carica WMS Catasto se richiesto
+    if carica_wms:
+        carica_wms_catasto()
+
     # Riepilogo finale
     print("\n" + "=" * 60)
     print("  COMPLETATO!")
@@ -776,13 +858,15 @@ class BBoxDrawTool(QgsMapTool):
     clicca il secondo angolo per confermare e avviare il download.
     """
 
-    def __init__(self, canvas, on_completed=None, espandi_catastale=False):
+    def __init__(self, canvas, on_completed=None, espandi_catastale=False,
+                 carica_wms=False):
         super().__init__(canvas)
         self.canvas = canvas
         self.first_point = None
         self.preview_rb = None
         self.on_completed = on_completed
         self.espandi_catastale = espandi_catastale
+        self.carica_wms = carica_wms
         self._create_preview_rubberband()
 
     def _create_preview_rubberband(self):
@@ -835,7 +919,8 @@ def canvasPressEvent(self, event):
             esegui_download_e_caricamento(
                 min_lat, min_lon, max_lat, max_lon,
                 layer_name="Particelle WFS (BBox)",
-                espandi_catastale=self.espandi_catastale
+                espandi_catastale=self.espandi_catastale,
+                carica_wms=self.carica_wms,
             )
 
             # Ripristina Pan e riapri dialog
@@ -875,11 +960,13 @@ class PolySelectTool(QgsMapTool):
     Estrae il bbox della geometria e avvia il download WFS.
     """
 
-    def __init__(self, canvas, on_completed=None, espandi_catastale=False):
+    def __init__(self, canvas, on_completed=None, espandi_catastale=False,
+                 carica_wms=False):
         super().__init__(canvas)
         self.canvas = canvas
         self.on_completed = on_completed
         self.espandi_catastale = espandi_catastale
+        self.carica_wms = carica_wms
 
     def canvasPressEvent(self, event):
         click_map_point = self.toMapCoordinates(event.pos())
@@ -968,7 +1055,8 @@ def canvasPressEvent(self, event):
                     min_lat, min_lon, max_lat, max_lon,
                     filter_geom=poly_geom_wfs,
                     layer_name="Particelle WFS (Poligono)",
-                    espandi_catastale=self.espandi_catastale
+                    espandi_catastale=self.espandi_catastale,
+                    carica_wms=self.carica_wms,
                 )
 
                 qgis_iface.actionPan().trigger()
@@ -1001,12 +1089,13 @@ class LineSelectTool(QgsMapTool):
     """
 
     def __init__(self, canvas, buffer_distance=BUFFER_DISTANCE_M, on_completed=None,
-                 espandi_catastale=False):
+                 espandi_catastale=False, carica_wms=False):
         super().__init__(canvas)
         self.canvas = canvas
         self.buffer_distance = buffer_distance
         self.buffer_rb = None  # Rubberband per visualizzare il buffer
         self.espandi_catastale = espandi_catastale
+        self.carica_wms = carica_wms
         self.on_completed = on_completed
         # Stato per modalità disegno polilinea
         self._draw_points = []  # Vertici della polilinea in coordinate mappa
@@ -1075,7 +1164,8 @@ def _esegui_download_da_linea(self, line_geom, geom_crs):
             min_lat, min_lon, max_lat, max_lon,
             filter_geom=buffer_geom_wfs,
             layer_name=f"Particelle WFS (Linea buffer {self.buffer_distance} m)",
-            espandi_catastale=self.espandi_catastale
+            espandi_catastale=self.espandi_catastale,
+            carica_wms=self.carica_wms,
         )
 
         qgis_iface.actionPan().trigger()
@@ -1334,13 +1424,14 @@ class PointSelectTool(QgsMapTool):
     """
 
     def __init__(self, canvas, buffer_distance=1, snap_tolerance=15,
-                 on_completed=None, espandi_catastale=False):
+                 on_completed=None, espandi_catastale=False, carica_wms=False):
         super().__init__(canvas)
         self.canvas = canvas
         self.buffer_distance = buffer_distance
         self.snap_tolerance = snap_tolerance
         self.buffer_rb = None
         self.espandi_catastale = espandi_catastale
+        self.carica_wms = carica_wms
         self.on_completed = on_completed
 
     def _visualizza_buffer(self, buffer_geom, buffer_crs):
@@ -1598,6 +1689,7 @@ def _processa_layer_punti(self, layer):
             layer_name=f"Particelle WFS (Punti buffer {self.buffer_distance} m)",
             espandi_catastale=self.espandi_catastale,
             post_filter_points=wfs_points,
+            carica_wms=self.carica_wms,
         )
 
         qgis_iface.actionPan().trigger()
@@ -1663,6 +1755,7 @@ def _processa_click_singolo(self, click_point, click_crs):
             layer_name="Particelle WFS (click punto)",
             espandi_catastale=self.espandi_catastale,
             post_filter_points=[pt_wfs],
+            carica_wms=self.carica_wms,
         )
 
         qgis_iface.actionPan().trigger()
@@ -1761,6 +1854,7 @@ def _on_modalita_scelta(self):
         dlg = self._dlg
         canvas = self.iface.mapCanvas()
         espandi = dlg.espandi_catastale
+        wms = dlg.carica_wms
 
         if dlg.scelta == "disegna":
             print("\n  MODALITÀ: Disegna BBox")
@@ -1769,7 +1863,7 @@ def _on_modalita_scelta(self):
             print("  >>> Clicca per il SECONDO angolo")
             print("  >>> Il download partirà automaticamente\n")
             tool = BBoxDrawTool(canvas, on_completed=self._reopen_dialog,
-                                espandi_catastale=espandi)
+                                espandi_catastale=espandi, carica_wms=wms)
             canvas.setMapTool(tool)
             self._active_tool = tool
 
@@ -1779,7 +1873,7 @@ def _on_modalita_scelta(self):
             print("  >>> Il bbox verrà estratto e il CRS verificato")
             print("  >>> Il download partirà automaticamente\n")
             tool = PolySelectTool(canvas, on_completed=self._reopen_dialog,
-                                  espandi_catastale=espandi)
+                                  espandi_catastale=espandi, carica_wms=wms)
             canvas.setMapTool(tool)
             self._active_tool = tool
 
@@ -1792,7 +1886,7 @@ def _on_modalita_scelta(self):
             print("  >>> ATTENZIONE: Il layer deve avere un CRS proiettato (metri)\n")
             tool = LineSelectTool(canvas, buffer_distance=buffer_m,
                                   on_completed=self._reopen_dialog,
-                                  espandi_catastale=espandi)
+                                  espandi_catastale=espandi, carica_wms=wms)
             canvas.setMapTool(tool)
             self._active_tool = tool
 
@@ -1808,7 +1902,7 @@ def _on_modalita_scelta(self):
             tool = PointSelectTool(canvas, buffer_distance=buffer_m,
                                    snap_tolerance=snap_px,
                                    on_completed=self._reopen_dialog,
-                                   espandi_catastale=espandi)
+                                   espandi_catastale=espandi, carica_wms=wms)
             canvas.setMapTool(tool)
             self._active_tool = tool