[Gfoss] problema postgis

Iacopo Zetti iacopozetti a gmail.com
Ven 9 Maggio 2008 15:10:58 CEST


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



Maggiori informazioni sulla lista Gfoss