Acquisizione contemporanea di informazioni da mappe vettoriali e raster

Lo scopo di questa esercitazione è quello di incrociare delle informazioni da mappe vettoriali e raster e visualizzarle. Si vuole ricavare il tipo di suolo che compone le proprietà di G. Tandy. Per fare questo si deve incrociare la cartografia catastale (vettoriale) con la carta geologica (raster). Dapprima si usi il GIS Manager per visualizzare il file raster geology ed il file vettoriale fields.

Si devono impostare le proprietà di rappresentazione del file vettoriale prendendo confidenza con le proprietà di visualizzazione delle mappe vettoriali. Ad esempio, partire cercare di visualizzare i contorni di ogni area colorandoli di rosso, e via via operare sulle trasparenze e visualizzare a schermo i proprietari dei fondi.


Proprietà di visualizzazione del vettoriale fields

Visualizzazione composita

Ora si estrarrà dalla mappa catastale la proprietà di George Tandy. Dapprima si interroghi la mappa vettoriale e si identifichi la particella catastale appartenente a Tandy. Nella finestra di Output del Gis.m, che visualizza la tabella associata al file vettorial,e si noterà che l'identificalivo è riportato nella colonna label sotto la stringa G. Tandy#1 .


Risultato dell'interrogazione della mappa vettoriale

Si può quindi procedere all'estrazione. Per fare ciò si utilizza il comando v.extract (Vector -> Query by attributes) come visualizzato dalla figura seguente, creando il file vettoriale Gtandy:


Risulatato dell'interrogazione della mappa vettoriale
Attenzione: la mappa vettoriale Gtandy non deve esistere, v.extract si rifiuta di sovrascrivere una mappa esistente. Eventualmente si può cancellarla con g.remove (File -> Manage maps and volumes -> Remove maps).

Si imposti ora la regione di lavoro sulle dimensione della nuova mappa estratta.

Ora si converte la mappa vettoriale in una mappa raster: per fare questo si usi il comando v.to.rast (File -> Map type conversions -> Vector to raster), usando il valore (in Source of raster values mettere val) e null'altro.


Conversione in raster della mappa vettoriale

Per creare il report del tipo di terreno presente sulla particella del sig. George Tandy occorre utilizzare il comando r.stats (Raster -> Report and statistics -> Sum area by map and category).


Report