Muitas vezes, passei horas lutando com bugs estranhos resultantes de geometrias inválidas . Os sintomas variam muito, incluindo:

• Lentidão do banco de dados
• Resultados de consulta \ Null incorretos
• Arquivos de log sobrecarregados com erros e avisos
• Resultado inconsistente com operadores geométricos (por exemplo, interseção)

Gostaria de categorizar os tipos de geometrias inválidas, para ajudar todos aqui a higienizar seus dados.

Responda com seu tipo favorito de geometria inválida (um tipo por resposta). Capturas de tela, descrições e soluções são bem-vindas - mas tente manter sua resposta curta.

Uma mina de ouro de 8 anos de geometria inválida é fornecida pelos testes de regressão isométricos do PostGIS . Muitos desses testes de regressão foram movidos para o GEOS, por exemplo, aqui e aqui , além de outros bolsos na árvore de testes . (Desculpe, muitos para destacar. Escolha alguns e cole o WKT no construtor JTS Test para visualizá-los .)

Segmentos de comprimento zero

Exemplo:

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

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

Problemas: O armazenamento do segmento de comprimento zero pode causar erros ao calcular as proporções de comprimento (divisão por zero) ou ao tentar calcular azimutes e funções trigonométricas.

Detecção: no PostGIS, segmentos de comprimento zero podem ser detectados usando isValid().

Soluções: Tente manter segmentos de comprimento zero como pontos.

Polígonos com auto-interseção

Exemplo:

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)

Problemas: cálculos incorretos, por exemplo area():

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

Isso pode causar erros de divisão zero com funções de chamada.

Detecção: No PostGIS, segmentos com auto-interseção podem ser detectados usando isValid().

Solução: Conversão para MULTIPOLYGON(ver comentário).