[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