Molte volte e spesso ho passato ore a lottare con bizzarri bug derivanti da geometrie non valide . I sintomi variano notevolmente, tra cui:

• Rallentamenti del database
• Risultati della query \ Null errati
• File di registro sovraccarichi di errori e avvisi
• Risultato incoerente con operatori geometrici (ad es. Intersezione)

Vorrei classificare i tipi di geometrie non valide, al fine di aiutare tutti qui a disinfettare i propri dati.

Rispondi con il tuo tipo preferito di geometria non valida (un tipo per risposta). Schermate, descrizioni e soluzioni sono benvenute, ma cerca di mantenere la risposta breve.

Una miniera d'oro di 8 anni con geometrie non valide è fornita dai test di regressione isopatica di PostGIS . Molti di questi test di regressione sono stati spostati su GEOS, ad esempio qui e qui , così come altre tasche nella struttura dei test . (Siamo spiacenti, troppi da evidenziare. Scegline alcuni e incolla il WKT nel builder JTS Test per visualizzarli .)

Segmenti di lunghezza zero

Esempio:

SELECT geomFromEWKT('LINESTRING(1 1, 1 1)');

                                    geomfromewkt                                    
------------------------------------------------------------------------------------
 0102000000020000000000000000000840000000000000104000000000000008400000000000001040

Problemi: la memorizzazione del segmento di lunghezza zero può causare errori nel calcolo dei rapporti di lunghezza (divisione per zero) o nel tentativo di calcolare azimut e funzioni trigonometriche.

Rilevamento: in PostGIS è possibile rilevare segmenti di lunghezza zero utilizzando isValid().

Soluzioni: prova a mantenere segmenti di lunghezza zero come punti.

Poligoni autointersecanti

Esempio:

SELECT isValid(geomFromEWKT('POLYGON((0 0, 1 1, 1 0, 0 1, 0 0))'));
NOTICE:  Self-intersection at or near point 0.5 0.5
 isvalid 
---------
 f
(1 row)

Problemi: calcoli errati, ad esempio area():

SELECT area(geomFromEWKT('POLYGON((0 0, 1 1, 1 0, 0 1, 0 0))'));
 area 
------
    0
(1 row)

Ciò può causare errori di divisione zero con funzioni di chiamata.

Rilevamento: in PostGIS è possibile rilevare segmenti autointersecanti utilizzando isValid().

Soluzione: conversione in MULTIPOLYGON(vedi commento).