GRASS logo

NOME

r.in.gdal - Importa un file raster supportato da GDAL in un layer raster binario

PAROLE CHIAVE

raster, importazione

SINOSSI

r.in.gdal
r.in.gdal help
r.in.gdal [-oeflk] [input=nome] [output=nome] [band=integer] [memory=integer] [target=nome] [title="phrase"] [location=nome] [--overwrite] [--verbose] [--quiet]

Flag:

-o
Sovrascrive la proiezione (usa la proiezione della location)
-e
Estende le dimensioni della regione basandosi sul nuovo dataset
Inoltre aggiorna la regione predefinita nel mapset PERMANENT
-f
Elenca i formati supportati ed esce
-l
Forza le mappe Lat/Lon per ricadere nelle coordinate geografiche (90N,S; 180E,W)
-k
Usa i numeri delle bande invece che i nomi dei colori
--overwrite
Sovrascrive i file esistenti
--verbose
Output verboso del modulo
--quiet
Output quieto del modulo

Parametri:

input=nome
Il file raster da importare
output=nome
Nome della mappa raster di output
band=integer
Banda da selezionare (il valore predefinito è tutte)
memory=integer
Dimensione della memoria cache (MiB)
target=nome
Nome della location da cui leggere la proiezione per la trasformazione dei GPG
title="phrase"
Titolo della mappa raster risultante
location=nome
Nome della nuova location da creare

DESCRIZIONE

r.in.gdal permette all'utente di creare un tematismo raster nel formato binario di GRASS, o un gruppo di immagini, partendo da un file in un qualsiasi qualsiasi formato raster supportato dalle GDAL, con un titolo. Il file importato può essere usato anche per creare una nuova location.

Formati raster suportati da GDAL

Tutti i dettagli sui formati raster supportati dalle GDAL sono reperibili al seguente indirizzo:

http://www.gdal.org/formats_list.html

Alcuni formati, selezionati tra gli oltre 80 formati supportati:


Nome completo del formato                               Codice                  Creazione       Georeferenziazione      Dimensione massima del file

Arc/Info ASCII Grid                                     AAIGrid                 Yes             Yes                     No limits
Arc/Info Binary Grid                                    AIG                     No              Yes                     --
AIRSAR Polarimetric                                     AIRSAR                  No              No                      --
Microsoft Windows Device Independent Bitmap (.bmp)      BMP                     Yes             Yes                     4GiB
BSB Nautical Chart Format (.kap)                        BSB                     No              Yes                     --
VTP Binary Terrain Format (.bt)                         BT                      Yes             Yes                     --
CEOS (Spot for instance)                                CEOS                    No              No                      --
First Generation USGS DOQ (.doq)                        DOQ1                    No              Yes                     --
New Labelled USGS DOQ (.doq)                            DOQ2                    No              Yes                     --
Data (.dt0, .dt1)                                       DTED                    No              Yes                     --
ERMapper Compressed Wavelets (.ecw)                     ECW                     Yes             Yes			
ESRI .hdr Labelled                                      EHdr                    No              Yes                     --
ENVI .hdr Labelled Raster                               ENVI                    Yes             Yes                     No limits
Envisat Image Product (.n1)                             Envisat                 No              No                      --
EOSAT FAST Format                                       FAST                    No              Yes                     --
FITS (.fits)                                            FITS                    Yes             No			
Graphics Interchange Format (.gif)                      GIF                     Yes             No                      2GB
Arc/Info Binary Grid (.adf)                             GIO                     Yes             Yes		     
GRASS Rasters                                           GRASS                   No              Yes                     --
TIFF / GeoTIFF (.tif)                                   GTiff                   Yes             Yes                     4GiB
Hierarchical Data Format Release 4 (HDF4)               HDF4                    Yes             Yes                     2GiB
Erdas Imagine (.img)                                    HFA                     Yes             Yes                     No limits
Atlantis MFF2e                                          HKV                     Yes             Yes                     No limits
Image Display and Analysis (WinDisp)                    IDA                     Yes             Yes                     2GB
ILWIS Raster Map (.mpr,.mpl)                            ILWIS                   Yes             Yes                     --
Japanese DEM (.mem)                                     JDEM                    No              Yes                     --
JPEG JFIF (.jpg)                                        JPEG                    Yes             Yes                     4GiB (max dimentions 65500x65500)
JPEG2000 (.jp2, .j2k)                                   JPEG2000                Yes             Yes                     2GiB
JPEG2000 (.jp2, .j2k)                                   JP2KAK                  Yes             Yes                     No limits
NOAA Polar Orbiter Level 1b Data Set (AVHRR)            L1B                     No              Yes                     --
Erdas 7.x .LAN and .GIS                                 LAN                     No              Yes                     2GB
In Memory Raster                                        MEM                     Yes             Yes                     2GiB
Atlantis MFF                                            MFF                     Yes             Yes                     No limits
Multi-resolution Seamless Image Database                MrSID                   No              Yes                     --
NDF                                                     NLAPS Data Format       No              Yes                     No limits
NITF                                                    NITF                    Yes             Yes		
NetCDF                                                  netCDF                  Yes             Yes                     2GB
OGDI Bridge                                             OGDI                    No              Yes                     --
PCI .aux Labelled                                       PAux                    Yes             No                      No limits
PCI Geomatics Database File                             PCIDSK                  Yes             Yes                     No limits
Portable Network Graphics (.png)                        PNG                     Yes             No		
PCRaster (.map)                                         PCRaster                Yes             No			
Netpbm (.ppm,.pgm)                                      PNM                     Yes             No                      No limits
RadarSat2 XML (product.xml)                             RS2                     No              Yes                     4GB
USGS SDTS DEM (*CATD.DDF)                               SDTS                    No              Yes                     --
SAR CEOS                                                SAR_CEOS                No              Yes                     --
USGS ASCII DEM (.dem)                                   USGSDEM                 No              Yes                     --
X11 Pixmap (.xpm)                                       XPM                     Yes             No			

Creazione della location

r.in.gdal cerca di conservare le informazioni sulla proiezione quando importa i dataset, se il formato d'origine include informazioni sulla proiezione, e se il driver delle GDAL lo supporta. Se la proiezione del dataset d'origine non corrisponde a quella della location corrente, r.in.gdal produce un messaggio d'errore (Projection of dataset does not appear to match current location) e riporta i parametri PROJ_INFO del dataset d'origine.

Se l'utente preferisce ignorare le differenze fra il sistema di coordinate che appare nei dati d'origine e quello della location corrente, può passare il flag -o per tralasciare la verifica della proiezione.

Se l'utente preferisce importare i dati con la definizione completa della proiezione, è possibile fare in modo che r.in.gdal crei automaticamente una nuova location basata sulla proiezione e sull'estensione del file da leggere. Si può fare passando il nome da utilizzare per la nuova location attraverso il parametro location. Una volta completato il comando, verrà creata una nuova location (con solo un mapset PERMANENT), e importato il raster col nome indicato in output nel mapset PERMANENT.

Supporto per i GCP: Nel caso in cui l'immagine contenga GCP [Ground Control Points - NdT], questi verranno scritti in un file POINTS all'interno di un gruppo imagery, e possono essere usati direttamente con i.rectify. L'opzione target permette di riproiettare automaticamente i GCP dalla proiezione propria in un'altra proiezione letta dal file PROJ_INFO della location denominata target.

NOTE

In futuro sono previsti dei miglioramenti in r.in.gdal, che includeranno il supporto per riportare tutte le informazioni possibili di un dataset, se non viene impostato il parametro output.

Il comando r.in.gdal supporta le seguenti funzionalità, sempre che il driver per quel formato le supporti:

Tavola dei colori
Le tavole dei colori associate alle bande di colore vengono trasferite. Da notare che se il file d'origine non ha una mappa dei colori, r.in.gdal in GRASS 5.0 non ne crea una. Usare r.colors map=... color=grey per assegnare una mappa dei colori a scala di grigi. In una successiva versione di GRASS r.in.gdal probabilmente verrà aggiornato per creare automaticamente mappe dei colori a scala di grigi.
Tipi di dati
È supportata la maggior parte dei tipi di dati delle GDAL. I tipi di banda di colore Float32 e Float64 sono tradotti in GRASS come celle con valori a virgola mobile (ma non a precisone doppia ... questo verrà aggiunto se necessario), e la maggior parte degli altri tipi sono tradotti in GRASS come celle con valori interi, compresi i valori interi a 16bit. Bande con dati complessi (alcuni formati di dati del segnale SAR) vengono convertiti in due tematismi con celle a virgola mobile (*.real e *.imaginary).
Georeferenziazione
Se il dataset contiene informazioni sulla georeferenziazione, questa verrà usata per impostare i bordi nord, sud, est e ovest. I coefficienti rotazionali verranno ignorati, poiché risultano con un orientamento non corretto per i dataset ruotati.
Proiezione
La proiezione dei dataset è utilizzata per verificare la location corrente o per definirne una nuova. Internamente le GDAL reppresentano le proiezioni nel formato Well Known Text di OpenGIS. È supportata gran parte della serie di proiezioni di GRASS.
Valori nulli
I pixel con valore nullo delle bande di colore raster per le quali le GDAL riconoscono un valore nullo, vengono trasformati durante l'importazione in valori nulli nello stile di GRASS. Per molti formati generici (e per i formati poco supportati dalle GDAL), non c'è un modo per riconoscere i pixel nulli, nel qual caso dovrebbe essere usato r.null dopo l'importazione.
GCP
I punti di controllo a terra (Ground Control Points) dei dataset contenenti GCP vengono importati come un file POINTS associato ad un gruppo imagery. I dataset con una sola banda, che altrimenti sarebbero stati convertiti come una semplice mappa raster, se contengono punti di controllo a terra hanno anche un gruppo imagery associato. Il sistema di coordinate dei punti di controllo viene riportato da r.in.gdal, ma non conservato. È compito dell'utente assicurarsi che la location definita con i.target abbia un sistema di coordinate compatibile, prima di usare i.rectify.

Messagi d'errore

"ERROR: Input map is rotated - cannot import."
In questo caso l'immagine dev'essere prima ruotata esternamente, applicando le informaziooni sulla rotazione memorizzate nel campo metadati del file dell'immagine raster. Per esempio, il software gdalwarp può essere usato per trasformare la mappa in modo che abbia il Nord in alto (nota: ci sono diversi parametri di gdalwarp per selezionare l'algoritmo di ricampionamento):
gdalwarp rotated.tif northup.tif

"ERROR: Projection of dataset does not appear to match the current location."
Si deve creare una location la cui proiezione corrisponda ai dati da importare. Col parametro location si può creare una nuova location in base alle informazioni sulla proiezione contenute nel file. Se si desidera, si può riproiettarlo su un'altra location con r.proj. In alternativa è possibile ignorare quest'errore usando il flag -o

"WARNING: G_set_window(): Illegal latitude for North"
Le location Latitudine/Longitudine in GRASS non possono avere regioni che superano i 90° Nord o Sud. L'immagine non georeferenziata avrà coordinate basate sul numero di pixel dell'immagine: 0,0 in basso a sinistra; colonne,righe in alto a destra. Tipicamente l'immagine sarà più grande di 90 pixel e perciò il GIS si rifiuta di importarla. Se si è sicuri che i dati sono appropriati per la location Lat/Lon e si intende reimpostare i confini della mappa col modulo r.region direttamente dopo l'importazione, si può usare il flag -l per forzare le coordinate della mappa a valori validi. Mentre i confini e la risoluzione risultante probabilmente saranno errati per la mappa, i dati della mappa saranno inalterati e sicuri. Dopo aver reimpostato ai confini noti con r.region si dovrebbe fare una doppia verifica con r.info, mettendo particolare attenzione alla risoluzione della mappa. Nella maggior parte dei casi sarà necessario importare nella proiezione nativa del file di dati, o in una seplice location XY e utilizzare gli strumenti di georettificazione (i.rectify et al.) per proiettare correttamente. Il flag -l dovrebbe essere usato solo se si sa che la proiezione è corretta ma la georeferenziazione interna è andata perduta, e limiti e risoluzione sono già noti.

ESEMPI

DEM GTOPO30

Per evitare che i dati GTOPO30 vengano letti in maniera non corretta, si può aggiungere nel file .HDR una nuova riga: "PIXELTYPE SIGNEDINT", per forzare l'interpretazione del file come file di numeri relativi anziché come file di numeri naturali. Quindi può essere importato il file .DEM. Infine, alla mappa importata può essere assegnata una tavola dei colori con r.colors, p.es. la tavola "terrain".

DEM GLOBE

Per importare le tile dei DEM GLOBE (risoluzione approssimativa 1km, più elevata di quella dei dati DEM GTOPO30), l'utente deve scaricare anche il/i relativo/i file HDR). Infine, alla mappa importata può essere assegnata una tavola dei colori con r.colors, p.es. la tavola "terrain".

Worldclim.org

Per importare dati Worldclim, bisogna aggiungere le seguent righe ad ogni file .hdr:
PIXELTYPE SIGNEDINT

HDF

L'importazione di bande HDF richiede la specificazione delle singole bande come sono riconosciute dalle GDAL:
# Esempio MODIS FPAR
gdalinfo MOD15A2.A2003153.h18v04.004.2003171141042.hdf
...
Subdatasets:
  SUBDATASET_1_NOME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km
  SUBDATASET_1_DESC=[1200x1200] Fpar_1km MOD_Grid_MOD15A2 (8-bit unsigned integer)
  SUBDATASET_2_NOME=HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Lai_1km
  SUBDATASET_2_DESC=[1200x1200] Lai_1km MOD_Grid_MOD15A2 (8-bit unsigned integer)
...

# importazione della prima banda, nell'esempio FPAR 1km:
r.in.gdal HDF4_EOS:EOS_GRID:"MOD15A2.A2003153.h18v04.004.2003171141042.hdf":MOD_Grid_MOD15A2:Fpar_1km \
          out=fpar_1km_2003_06_02
# ... in modo analogo per le altre bande HDF del file.

VEDI ANCHE

r.colors, r.in.ascii, r.in.arc, r.in.bin, r.null

RIFERIMENTI

Pagine delle GDAL: http://www.gdal.org/

AUTORE

Frank Warmerdam (email).

Last changed: $Date: 2010-09-21 14:19:28 +0200(mar, 21 set 2010) $


Main index - raster index - Full index

© 2003-2011 GRASS Development Team