Come si prende un livello vettoriale con poligoni sovrapposti e lo si trasforma in un raster in cui ogni cella conta il numero di poligoni che contengono quella cella?
Stiamo osservando i poligoni delle aree bruciate dal fuoco tra il 1976 e il 2000 e vogliamo vedere quanti incendi sarebbero contenuti in ciascuna cella raster.
Risposte:
Fallo in tre passaggi: spezza i poligoni nelle loro parti componenti, conta le sovrapposizioni e converti in raster. Questo evita il costo computazionale potenzialmente enorme di convertire separatamente ogni poligono in un raster e combinare quei raster.
Union(nel Geoprocessingmenu) suddivide i poligoni nelle loro parti.
Sfortunatamente, ogni sovrapposizione è duplicata nell'output: ha una copia identica per ogni poligono originale che la copre. Perciò
Dissolve(di nuovo nel Geoprocessingmenu) uniranno le parti sovrapposte, a condizione che sia possibile trovare un modo per identificarle in modo univoco. Leggi la finestra di dialogo: verso la fine, avrai la possibilità di calcolare "statistiche". Scegli un campo che potrebbe aver identificato l' originale poligoni e richiedi un conteggio.
In molti casi la combinazione dell'area poligonale e del perimetro identificherà in modo univoco le parti. In caso contrario, puoi aggiungere più proprietà geometriche in campi aggiuntivi, come le coordinate del centroide, fino a quando non hai accumulato abbastanza informazioni per distinguere ogni caratteristica.
Il livello risultante ha una funzione per ogni sovrapposizione di poligoni e una sorta di campo "conta" che conta il numero di sovrapposizioni.
Converti quello in un raster, usando il campo "count" per gli attributi.
Ad esempio, ecco alcuni poligoni sovrapposti e i loro identificatori con la tabella degli attributi mostrata:
Dopo il secondo passaggio abbiamo un record per ogni regione sovrapposta insieme a un conteggio che può già essere usato per simboleggiare la quantità di sovrapposizione:
Il resto è facile - ed è solo una singola operazione di rasterizzazione.
Unionè che con quasi lo stesso flusso di lavoro può affrontare la situazione in cui i poligoni risiedono in più set di dati (che di solito è una cattiva progettazione del database ma purtroppo comune): si uniscono semplicemente tutti i set di dati di input contemporaneamente.
union) seguite dalla rasterizzazione? Le operazioni vettoriali si impantaneranno quando le caratteristiche saranno digitalizzate con più dettagli del necessario, dotandole di troppi vertici. In queste circostanze estreme, l'approccio raster potrebbe essere migliore (sebbene prima semplificare i poligoni potrebbe essere un'opzione superiore). In tutte le altre situazioni, tuttavia, rasterizzare ogni poligono separatamente è un enorme spreco di tempo per computer e umani.
Il seguente post ha una domanda in qualche modo simile che potresti trovare una soluzione correlata: Creare una superficie raster di sovrapposizioni da shapefile poligonali vettoriali? .
Per questo approccio raster semplice e computazionalmente rapido, dovrai (1) separare i poligoni sovrapposti in strati separati (presumibilmente per anno per i poligoni del fuoco) usando Seleziona per attributi in ModelBuilder con un iteratore o uno strumento di script, (2) convertirli ciascuno Poligono a raster con l'assegnazione di cella MAXIMUM_COMBINED_AREA (garantendo la stessa dimensione di cella , snap raster e che l' estensione rimanga uguale all'intero set di poligoni) - utilizzando un valore di campo costante (ad esempio, utilizzare il campo o le righe dell'anno che hanno tutti 1) per convertire (utilizzare nuovamente ModelBuilder con un iteratore o uno script Python per automatizzare), quindi (3) applicare il seguente strumento Analista spaziale: Statistiche cella - usa il parametro statistics_type VARIETY se ogni raster ha un valore univoco come l'anno o SUM se tutti i valori delle celle raster sono 1s - assicurati di controllare per ignorare NoData.
I raster intermedi (dalla conversione precedente) possono essere eliminati o sono pronti per essere utilizzati nelle successive analisi raster.