Salve
mi trovo nella necessità di clippare un raster mediante una maschera di uno shape.
Ho provato ad utilizzare il plugin di qgis (stò su qgis 1.8.0 lisboa winXP) ma mi restituisce il seguente "ERROR: Cannot compute bounding box of cutline."
Da una breve ricerca ho visto che c'era un ticket aperto proprio su questo, che poi è stato chiuso perchè l'extent della maschera non era contenuto nella extent del raster di input.
Nel mio caso il vettore è interamente contenuto nel raster (copioincollo dalle proprietà dei layer di qgis)
ext raster:
316500.0000000000000000,3878500.0000000000000000 : 1438500.0000000000000000,5216500.0000000000000000
ext maschera:
Unità del sistema di riferimento spaziale usato nel vettore: xMin, yMin 752340.68,4432311.01 : xMax,yMax 926926.87,4551063.10
Non ci sono riproiezioni al volo, sistemi di proiezione uguali (epsg 23032) e medesimo sistema di riferimento del progetto qgis.
Da terminale ottengo lo stesso errore.
GDAL 1.9.1
GDAL Tools 1.2.29
Ogni suggerimento è prezioso, grazie
Benedetto Porfidia
Il 12/07/2012 11.51, Benedetto Porfidia ha scritto:
Salve
mi trovo nella necessità di clippare un raster mediante una maschera di
uno shape.
ext maschera:
Unità del sistema di riferimento spaziale usato nel vettore: xMin, yMin
752340.68,4432311.01 : xMax,yMax 926926.87,4551063.10
Aggiornamento:
facendo il clip sull'estensione anzichè sulla maschera, funziona correttamente.
in pratica è come se l'opzione "maschera" legge in maniera scorreta l'extent del poligono di clip
./b
Salve,
2012/7/12 Benedetto Porfidia <benedetto.porfidia@isprambiente.it>:
facendo il clip sull'estensione anzichè sulla maschera, funziona
correttamente.
sono due strumenti diversi (uno usa gdalwarp e l'altro gdal_translate
mi pare).
in pratica è come se l'opzione "maschera" legge in maniera scorreta l'extent
del poligono di clip
uhm... ora che ci rifletto questo potrebbe essere un bug di GdalTools.
In pratica non viene considerato il sistema di riferimento che imposti nella
legenda.
Quando viene fatta la chiamata al tool gdalwarp non viene passato
il CRS del vettore della maschera.
Quindi il CRS della maschera viene prelevato direttamente dalla sorgente
dati, per cui se il tuo vettore non ha un sistema di riferimento o ne ha uno
diverso potrebbe capitare che anche se in qgis lo vedi bene l'operazione
non va a buon fine.
Prova a salvare il tuo vettore nel CRS corretto, quindi usare quello nuovo
come maschera.
Potresti inoltre riaprire il ticket sul tracker di QGis? c'è in ogni
caso da vedere
se in gdalwarp ci sia un'opzione per ridefinire il sistema di riferimento della
maschera.
Saluti.
--
Giuseppe Sucameli - Faunalia
sono due strumenti diversi (uno usa gdalwarp e l'altro gdal_translate
mi pare).
Vero, non me ne ero accorto, sorry
uhm... ora che ci rifletto questo potrebbe essere un bug di GdalTools.
In pratica non viene considerato il sistema di riferimento che imposti nella
legenda.
Quando viene fatta la chiamata al tool gdalwarp non viene passato
il CRS del vettore della maschera.
Quindi il CRS della maschera viene prelevato direttamente dalla sorgente
dati, per cui se il tuo vettore non ha un sistema di riferimento o ne ha uno
diverso potrebbe capitare che anche se in qgis lo vedi bene l'operazione
non va a buon fine.
lo shape del vettore ha il suo prj impostato nello stesso sistema di riferimento
appena ho un po' di tempo vedo di riaprire il ticket 
ciao
./b