Scusate se mi auto-rispondo, ma magari l'esperienza può essere utile ad altri.
Ho risolto semplificando i poligoni: CREATE TABLE nomenuovolayer AS SELECT SnapToGrid(vecchiolayer.the_geom, 0.001) AS geom, altricampi FROM vecchiolayer e poi rifacendo l'intersezione come prima, ma sulle nuove tabelle. CREATE TABLE intersezione AS SELECT intersection(f.the_geom, c.the_geom) as the_geom, f.altricampi, c.altricampi FROM layer1 f, layer2 c WHERE f.the_geom && c.the_geom AND Intersects(f.the_geom, c.the_geom); Se qualcuno che ne sa più di me ha voglia di completare l'informazione mi rimangono 2 domande: 1) si può fare lo SnapToGrid senza creare nuovi layer, ma direttamente dentro il comando Intersects? 2) se i poligoni non sono validi si può correggere l'errore della topologia con postgis o tocca passare per grass o altro? Buon fine settimana a tutti. Iacopo _______________________________________________ Iscriviti all'associazione GFOSS.it: http://www.gfoss.it/drupal/iscrizione [email protected] http://www.faunalia.com/cgi-bin/mailman/listinfo/gfoss Questa e' una lista di discussione pubblica aperta a tutti. I messaggi di questa lista non rispecchiano necessariamente le posizioni dell'Associazione GFOSS.it.
