Проверка геометрии ругается на корректный полигон

[Denis Rykov]

Сабж. Закинул его в PostGIS и проверил встроенными средствами - нет ошибок.
Подгружаю как слой в QGIS - ftools при проверке сообщает о самопересечении.

[Максим Дубинин]

похоже есть за что ругать

в полигоне 2 узла существуют в двойном экземпляре

[Denis Rykov]

Точно, сделал извлечение узлов, получилось 7 вместо 5.
Однако PostGIS на такую ситуацию не ругается, я почему то думал, что у него с кугисом общие библиотеки проверки топологической корректности.

[Максим Дубинин]

Значит надо продолжать разбираться, либо постгис тихо удаляет дубли, либо кугис излишне параноидален в такой ситуации.

Однако, ArcGIS тоже ругается. Мне все-таки кажется постгис самовольничает, так как дублированный узел это все-таки ошибка.

[Denis Rykov]

PostGIS не удаляет дубли:

Код: Выделить всё

POLYGON((317371.62563763 5731575.49932608,
317028.110300846 5731043.43103426,
317028.110300846 5731043.43103426,
316389.822682815 5731138.79331604,
316478.209591618 5731705.00321977,
316478.209591619 5731705.00321977,
317350.223258488 5731584.46661525,
317371.62563763 5731575.49932608))

Однако на запрос select ST_IsValidReason(the_geom) получаем:

Код: Выделить всё

Valid Geometry

[Denis Rykov]

А полигон то оказывается валидным согласно OGC, вот что мне ответили в рассылке:

This polygon is valid according to the OGC SFS, and thus validates correctly in PostGIS (also GEOS/JTS).

You'd have to look at the validation semantics of the other systems to determine why they don't validate it. The OGC SFS model is fairly tolerant in terms of trivial structural aspects - for instance, it allows duplicate points. I believe this is done in order to allow as a wide a range of geometries as possible to be represented, consistent with topological integrity. Some other systems are less tolerant.

Perhaps they don't allow duplicate points, or perhaps they are reducing the precision and thus causing a self-intersection (although I tried this in JTS, and reducing precision doesn't actually cause self-intersections)

[Максим Дубинин]

Продолжение, теоретическая дискуссия и результаты запроса PostGIS'овцев
ST_IsValid и повторяющиеся узлы