[Gfoss] Raster batch clipping

[G. Allegri]

Il giorno 20 ottobre 2011 15:56, Sandro Santilli <strk a keybit.net> ha
scritto:

> Con l'upcoming PostGIS-2.0 potrai fare tutto all'interno
> del database, a colpi di SQL.
>

Non vedo l'ora di provare le novità, prima fra tutte il supporto Raster!
E' già prevista una data di uscita?

giovanni


> --strk;
>
> On Thu, Oct 20, 2011 at 03:44:11PM +0200, Paolo Corti wrote:
> > 2011/10/20 Michele Beneventi <mbeneventi a gmail.com>:
> > > Grazie per le risposte!!!
> > > Per ora ho trovato una soluzione intermedia, considerando che sono un
> umile
> > > " utente finale", che raggiunge il suo scopo:
> > > la sintassi che ho usato è la seguente:
> > > gdalwarp -q -cutline /path/dello/shape
> > >               -cwhere "AttributoPerSQLClip = 'valoreAttributo'"
> > >               -crop_to_cutline /path/immagine/da clippare
> > >              /path/dell'/output/
> > > il parametro -cwhere è stata una bella sopresa!
> > > grazie ancora
> > > a presto
> > > Michele
> >
> > Ciao Michele
> > senza dover pre-trattare lo shapefile, potresti farti uno script
> > Python tipo questo, veramente molto semplice (usando ogr da codice e poi
> il
> > comando gdal_translate):
> >
> > import os
> > from osgeo import ogr
> >
> > def ClipRasterWithShape(raster, shapefile):
> >     driver = ogr.GetDriverByName('ESRI Shapefile')
> >     dataset = driver.Open(shapefile, 0)
> >     layer = dataset[0]
> >     for feature in layer:
> >         geom = feature.GetGeometryRef()
> >         envelope = geom.GetEnvelope()
> >         print 'Clipping feature %s with extent: %s' %
> > (feature.GetFID(), envelope)
> >         os.system('gdal_translate -projwin %s %s %s %s  %s out%s' %
> > (envelope[0], envelope[3], envelope[1], envelope [2], raster,
> > feature.GetFID()))
> >
> > ClipRasterWithShape('myraster.tif', 'myshape.shp')
> >
> > salva tutto come myprog.py e lancia:
> > $ python myprog.py
> >
> > se hai il supporto Python per gdal dovrebbe andare senza problemi
> >
> > ciao
> > P
> >
> > --
> > Paolo Corti
> > Geospatial software developer
> > web: http://www.paolocorti.net
> > twitter: @capooti
> > _______________________________________________
> > Iscriviti all'associazione GFOSS.it:
> http://www.gfoss.it/drupal/iscrizione
> > Gfoss a lists.gfoss.it
> > http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
> > Questa e' una lista di discussione pubblica aperta a tutti.
> > Non inviate messaggi commerciali.
> > I messaggi di questa lista non rispecchiano necessariamente
> > le posizioni dell'Associazione GFOSS.it.
> > 527 iscritti al 7.7.2011
>
> --
>
>  ()   Free GIS & Flash consultant/developer
>  /\   http://strk.keybit.net/services.html
> _______________________________________________
> Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione
> Gfoss a lists.gfoss.it
> http://lists.gfoss.it/cgi-bin/mailman/listinfo/gfoss
> Questa e' una lista di discussione pubblica aperta a tutti.
> Non inviate messaggi commerciali.
> I messaggi di questa lista non rispecchiano necessariamente
> le posizioni dell'Associazione GFOSS.it.
> 527 iscritti al 7.7.2011
>
-------------- parte successiva --------------
Un allegato HTML ? stato rimosso...
URL: <http://lists.gfoss.it/pipermail/gfoss/attachments/20111020/174a360f/attachment.html>