Grazie,<div>aspettando postgis 2.0 proverņ anche questa soluzione......</div><div><br></div><div>vi farņ sapere.</div><div><br></div><div>Michele<br><br><div class="gmail_quote">2011/10/20 Sandro Santilli <span dir="ltr"><<a href="mailto:strk@keybit.net">strk@keybit.net</a>></span><br>

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