{"id":4163,"date":"2017-03-06T13:30:13","date_gmt":"2017-03-06T12:30:13","guid":{"rendered":"https:\/\/oneteamgis.wordpress.com\/?p=4163"},"modified":"2018-03-09T12:39:41","modified_gmt":"2018-03-09T11:39:41","slug":"arcpy-come-verificare-le-interferenze-fra-le-geometrie-utilizzando-cursori","status":"publish","type":"post","link":"https:\/\/gis.oneteam.it\/index.php\/2017\/03\/06\/arcpy-come-verificare-le-interferenze-fra-le-geometrie-utilizzando-cursori\/","title":{"rendered":"ArcPy: come verificare le interferenze fra le geometrie"},"content":{"rendered":"

Avete mai avuto la necessit\u00e0 di eseguire un’operazione tanto semplice quando indispensabile come il dover estrarre l’area non condivisa fra due geometrie sovrapposte? In sostanza, l’opposto di un’operazione di Intersect<\/em>?<\/p>\n

\"symmetrical_difference\"
Differenza simmetrica<\/figcaption><\/figure>\n

Tale operazione pu\u00f2 essere eseguita con lo strumento Symmetrical Difference<\/a><\/em>. L’unico problema \u00e8 che si tratta di un tool che necessita di una licenza di livello Advanced<\/em>, non sempre a disposizione di tutti gli utenti. Nessun problema! In questo post vi faremo vedere come sfruttare gli oggetti Geometria<\/a> per rendere possibili questa e molte altre operazioni utilissime, superando il limite della licenza. Il tutto sfruttando Python e il modulo arcpy.Ma cominciamo subito con un esempio. Il seguente codice crea, a partire da uno shapefile poligonale contenente una sola geometria rettangolare, un oggetto Geometria:<\/p>\n

\r\nimport arcpy\r\ngeom_1 = arcpy.CopyFeatures_management("geom_1.shp", arcpy.Geometry())\r\n<\/pre>\n
Abbiamo sfruttato lo strumento Copy Features<\/a><\/em> per copiare le sole geometrie contenute nello shapefile geom_1.shp<\/code> in una variabile che abbiamo denominato geom_1<\/code>. Il tutto specificando come parametro di output della funzione la chiamata a arcpy.Geometry()<\/code>, che crea l’oggetto Geometria. Quest’oggetto ha una serie di metodi (intesi in senso pythonico<\/em>), fra cui uno chiamato symmetricDifference()<\/code>, che andremo a utilizzare per il nostro scopo.<\/p>\n
Andiamo oltre e procediamo con la creazione del secondo oggetto Geometria che sar\u00e0 sovrapposto in parte al primo precedentemente creato (proprio come nell’immagine a inizio post).<\/p>\n
\r\ngeom_2 = arcpy.CopyFeatures_management("geom_2.shp", arcpy.Geometry())\r\n<\/pre>\n
Nulla di nuovo: abbiamo replicato il codice per la creazione della prima geometria, semplicemente sostituendo il nome della variabile (geom_2<\/code> anzich\u00e9 geom_1<\/code>), e lo shapefile da copiare. Ora abbiamo due oggetti Geometria da poter usare come input per ottenere una terza Geometria che sia l’unione delle due geometrie meno l’intersezione fra di esse.<\/p>\n
Un’illustrazione render\u00e0 pi\u00f9 chiaro il funzionamento del metodo symmetricDifference<\/code>.<\/p>\n
\"symmetrical_difference\"
simmetricDifference<\/figcaption><\/figure>\n
Prima di procedere col nostro codice, vanno precisate alcune cose da tenere a mente quando si usano gli oggetti Geometria:<\/p>\n
    \n
  1. ogni poligono, inteso come entit\u00e0 geometrica individuale, che compone lo shapefile in input, corrisponder\u00e0 ad un oggetto Geometria di tipo poligonale. Questo, insieme agli altri (eventuali) oggetti, vengono raccolti in una lista di Python. In sostanza, geom_1<\/code> e geom_2<\/code> sono due liste, ciascuna avente al proprio interno un solo oggetto Geometria (poligonale in questo caso come gli shapefile da cui sono derivati). In pratica se stampassimo (con la funzione print<\/code>) il valore di una delle due variabili, otterremmo qualcosa di simile a:<\/li>\n<\/ol>\n
    [Polygon object at 0x4d34bb0[0x26e2040]]<\/code><\/p>\n
      \n
    1. gli oggetti Geometria rappresentano una posizione nello spazio e la forma ad essa associata, pertanto non contengono attributi.<\/li>\n<\/ol>\n
      Chiariti questi punti, finalizziamo ora il codice.<\/p>\n
      \r\ngeom_diff = geom_1[0].symmetricDifference(geom_2[0])\r\narcpy.CopyFeatures_management(geom_diff,"geom_diff.shp")\r\n<\/pre>\n
      Fatto! Abbiamo usato il metodo symmetricDifference()<\/code> per ottenere un nuovo oggetto Geometria che rappresenta l’unione delle due geometrie meno l’intersezione fra di esse (proprio quello che fa lo strumento Symmetrical Difference<\/em> dell’Analysis Toolbox<\/em> in ArcMap), e l’abbiamo salvato in uno shapefile geom_diff.shp<\/code> usando ancora lo strumento Copy Features<\/em> (in maniera opposta rispetto alle righe 2 e 3 nel codice). Notare che alla riga 4 abbiamo dovuto specificare un indice per la variabile geom_1<\/code> e geom_2<\/code>, dato che in realt\u00e0 si tratta di due liste contenenti ciascuna un solo elemento, l’oggetto Geometria. Ricordiamo che [0]<\/code> rappresenta un indice che punta al primo (e in questo caso unico) elemento all’interno di una lista in Python.<\/p>\n
      Riportiamo di seguito il codice completo:<\/p>\n
      \r\nimport arcpy\r\ngeom_1 = arcpy.CopyFeatures_management("geom_1.shp", arcpy.Geometry())\r\ngeom_2 = arcpy.CopyFeatures_management("geom_2.shp", arcpy.Geometry())\r\ngeom_diff = geom_1[0].difference(geom_2[0])\r\narcpy.CopyFeatures_management(geom_diff,"geom_diff.shp")\r\n<\/pre>\n
      Abbiamo visto in questo post come sfruttare gli oggetti Geometria e il loro metodo symmetricDifference()<\/code> per replicare il funzionamento dello strumento Symmetrical Difference<\/em> presente in ArcMap, ma utilizzabile solo con una licenza di livello Advanced<\/em>. Questo \u00e8 utile in svariate casistiche (clausole di contenimento, individuazione delle interferenze, ecc.). Ovviamente non \u00e8 il solo metodo degli oggetti Geometria: equals()<\/code> consente di confrontare due geometrie e restituisce true<\/code> se sono identiche, ed \u00e8 utile per verificare che non vi siano duplicati geometrici all’interno di un proprio dataset; generalize()<\/code> crea una geometria che \u00e8 la versione generalizzata (semplificata) di quella in input, utile quando si hanno troppi vertici e si vuole “snellire” il dato; l’opposto \u00e8 il metodo densify()<\/code>, che invece aggiunge vertici. Per un elenco completo, si faccia riferimento a questa pagina<\/a> dell’help online di ArcGIS.<\/p>\n
      Gli oggetti Geometria possono essere utilizzati come input o output al posto di creare nuove Feature Class temporanee, applicarvi dei cursori, e alla fine cancellarle durante una fase intermedia di un processo di geoprocessing. Possono anche essere creati da zero chiamando la classe Geometry()<\/code>, perfino “al volo”, se necessario, inserendoli come parametri durante una chiamata a una funzione di geoprocessing.<\/p>\n
      Nei prossimi post vedremo come mantenere le informazioni (gli attributi) degli shapefile di partenza agganciati agli oggetti Geometria tramite l’uso dei cursori.<\/p>\n","protected":false},"excerpt":{"rendered":"
      Avete mai avuto la necessit\u00e0 di eseguire un’operazione tanto semplice quando indispensabile come il dover estrarre l’area non condivisa fra due geometrie sovrapposte? In sostanza, l’opposto di un’operazione di Intersect? Tale operazione pu\u00f2 essere eseguita con lo strumento Symmetrical Difference. L’unico problema \u00e8 che si tratta di un tool che necessita di una licenza di…<\/p>\n","protected":false},"author":1,"featured_media":4498,"comment_status":"open","ping_status":"open","sticky":false,"template":"","format":"standard","meta":{"footnotes":""},"categories":[6,14,29],"tags":[75,93,470,507],"class_list":["post-4163","post","type-post","status-publish","format-standard","has-post-thumbnail","hentry","category-arcgis-for-desktop","category-articoli-tecnici","category-python","tag-arcgis-desktop","tag-arcpy","tag-python","tag-script"],"post_mailing_queue_ids":[],"_links":{"self":[{"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/posts\/4163","targetHints":{"allow":["GET"]}}],"collection":[{"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/posts"}],"about":[{"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/types\/post"}],"author":[{"embeddable":true,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/users\/1"}],"replies":[{"embeddable":true,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/comments?post=4163"}],"version-history":[{"count":2,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/posts\/4163\/revisions"}],"predecessor-version":[{"id":6753,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/posts\/4163\/revisions\/6753"}],"wp:featuredmedia":[{"embeddable":true,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/media\/4498"}],"wp:attachment":[{"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/media?parent=4163"}],"wp:term":[{"taxonomy":"category","embeddable":true,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/categories?post=4163"},{"taxonomy":"post_tag","embeddable":true,"href":"https:\/\/gis.oneteam.it\/index.php\/wp-json\/wp\/v2\/tags?post=4163"}],"curies":[{"name":"wp","href":"https:\/\/api.w.org\/{rel}","templated":true}]}}