Aggiorna versione a 1.4.2: funzione get_particella_info per calcolatore campi

Nuova funzione personalizzata get_particella_info($geometry) registrata
nel gruppo 'Catasto' del calcolatore di campi. Interroga il WFS e
restituisce un array con: riferimento catastale, foglio, label,
codice comune, geometria WKT, sezione e allegato.

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

Author: pigreco

get_particella_wfs.py

@@ -0,0 +1,126 @@
+# -*- coding: utf-8 -*-
+"""
+/***************************************************************************
+ WFS Catasto Agenzia delle Entrate CC BY 4.0
+                              -------------------
+        copyright            : (C) 2025 by Totò Fiandaca
+        email                : pigrecoinfinito@gmail.com
+ ***************************************************************************/
+"""
+
+from qgis.core import *
+from qgis.utils import qgsfunction
+import urllib.request
+import urllib.parse
+from xml.etree import ElementTree as ET
+
+def format_wkt(wkt, decimals=6):
+    """
+    Formatta una stringa WKT con il numero specificato di decimali.
+    """
+    import re
+    
+    def format_number(match):
+        num = float(match.group(0))
+        return f"{num:.{decimals}f}"
+    
+    formatted = re.sub(r'\d+\.\d+', format_number, wkt)
+    formatted = formatted.replace(",", ", ")
+    formatted = formatted.replace("), ", "),\n")
+    
+    return formatted
+
+@qgsfunction(args='auto', group='Catasto', usesgeometry=True)
+def get_particella_info(geom, feature, parent):
+    """
+    <h1>Catasto Agenzia delle Entrate CC BY 4.0:</h1>    
+    La funzione restituisce le informazioni WFS Catasto disponibili nella particella sottostante.
+
+    <h2>Parametri</h2>
+    <ul>
+      <li>geometry: geometria del punto (viene passata automaticamente)</li>
+    </ul>
+    
+    <h2>Returns</h2>
+    <ul>
+      <li>ARRAY: informazioni della particella</li>
+    </ul>
+    
+    <h2>Esempio</h2>
+        <pre>get_particella_info($geometry)[0]--> M011_0019C0.131</pre>
+        <pre>get_particella_info($geometry)[1]--> 0019</pre>
+        <pre>get_particella_info($geometry)[2]--> 131</pre>
+        <pre>get_particella_info($geometry)[3]--> M011</pre>
+        <pre>get_particella_info($geometry)[4]--> geometria WKT</pre>
+        <pre>get_particella_info($geometry)[5]--> _ (sezione censuaria)</pre>
+        <pre>get_particella_info($geometry)[6]--> C (allegato)</pre>
+    """
+    try:
+        # Verifica che la geometria sia un punto
+        if geom.type() != QgsWkbTypes.PointGeometry:
+            return ['ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR']
+            
+        # Prendi le coordinate del punto
+        point = geom.asPoint()
+        x = point.x()
+        y = point.y()
+        
+        # Base URL con i parametri base che sappiamo funzionare
+        uri = (f"pagingEnabled='true' "
+               f"preferCoordinatesForWfsT11='false' "
+               f"restrictToRequestBBOX='1' "
+               f"srsname='EPSG:6706' "
+               f"typename='CP:CadastralParcel' "
+               f"url='https://wfs.cartografia.agenziaentrate.gov.it/inspire/wfs/owfs01.php' "
+               f"version='2.0.0' "
+               f"language='ita'")
+               
+        # Crea un layer temporaneo per la richiesta
+        layer = QgsVectorLayer(uri, "catasto_query", "WFS")
+        
+        if not layer.isValid():
+            return ['ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR']
+            
+        # Crea il punto per il filtro spaziale
+        point_geom = QgsGeometry.fromPointXY(QgsPointXY(x, y))
+        request = QgsFeatureRequest().setFilterRect(point_geom.boundingBox())
+        
+        # Recupera le features
+        features = list(layer.getFeatures(request))
+        
+        if features:
+            # Prendi la prima particella trovata
+            feat = features[0]
+            ref = feat['NATIONALCADASTRALREFERENCE']
+            admin = feat['ADMINISTRATIVEUNIT']
+            label = feat['LABEL']
+            foglio = ref[5:9] if len(ref) > 9 else 'N/D'
+            
+            # Parsing NATIONALCADASTRALREFERENCE (formato CCCCZFFFFAS.particella)
+            sezione = 'N/D'
+            allegato = 'N/D'
+            if ref and isinstance(ref, str):
+                codice = ref.split(".")[0]  # parte prima del punto
+                if len(codice) == 11:  # CCCCZFFFFAS = 11 caratteri
+                    sez = codice[4]  # Z: sezione censuaria
+                    sezione = "" if sez == "_" else sez
+                    allegato = codice[9]  # A: allegato
+            
+            # Ottieni la geometria WKT e formattala
+            geom_wkt = 'N/D'
+            if feat.hasGeometry():
+                geom_wkt = format_wkt(feat.geometry().asWkt())
+            
+            # Restituisci la lista
+            return [ref, foglio, label, admin, geom_wkt, sezione, allegato]
+        else:
+            return ['N/D', 'N/D', 'N/D', 'N/D', 'N/D', 'N/D', 'N/D']
+                
+    except Exception as e:
+        return ['ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR', 'ERROR']
+
+# Esempio di utilizzo nel calcolatore di campi:
+# get_particella_info($geometry)[0]  # per il riferimento catastale
+# get_particella_info($geometry)[4]  # per la geometria WKT
+# get_particella_info($geometry)[5]  # per la sezione censuaria
+# get_particella_info($geometry)[6]  # per l'allegato

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.1
+version=1.4.2
 author=Salvatore Fiandaca
 email=pigrecoinfinito@gmail.com
 
@@ -30,6 +30,10 @@ repository=https://github.com/pigreco/wfs_catasto_download_particelle_bbox
 
 hasProcessingProvider=no
 changelog=
+    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
+            dalla particella WFS sottostante al punto.
     1.4.1 - Modalita' Seleziona Linea: possibilita' di disegnare una polilinea
             direttamente sulla mappa (click sinistro per aggiungere vertici,
             click destro per terminare, ESC per annullare).

wfs_catasto_download_particelle_bbox_p.py

@@ -30,6 +30,7 @@
     QgsGeometry,
     QgsFeature,
     QgsField,
+    QgsExpression,
 )
 from qgis.gui import QgsMapTool, QgsRubberBand
 from qgis.PyQt.QtCore import Qt, QVariant, QTimer
@@ -44,6 +45,7 @@
 from qgis.utils import iface as qgis_iface
 
 from .wfs_catasto_download_particelle_bbox_d import AvvisoDialog, SceltaModalitaDialog
+from .get_particella_wfs import get_particella_info
 
 
 # =============================================================================
@@ -1720,8 +1722,16 @@ def initGui(self):
         self.iface.addPluginToVectorMenu(self.menu, action)
         self.actions.append(action)
 
+        # Registra la funzione personalizzata nel calcolatore di campi
+        QgsExpression.registerFunction(get_particella_info)
+        print("[OK] Funzione personalizzata 'get_particella_info' registrata")
+
     def unload(self):
         """Rimuove azione dalla toolbar e dal menu."""
+        # Deregistra la funzione personalizzata
+        QgsExpression.unregisterFunction('get_particella_info')
+        print("[OK] Funzione personalizzata 'get_particella_info' deregistrata")
+        
         for action in self.actions:
             self.iface.removePluginVectorMenu(self.menu, action)
             self.iface.removeToolBarIcon(action)