Форумы GIS-Lab.info

Получил от коллеги векторные данные. Он уверял, что данные кристально чистые.

Изначально ошибка топологии возникла при работе функции rgeos::gBuffer в R. Там полно полигонов, которые имеют дыры, замыкающиеся на границу полигонов, т.е. граница сама на себя перехлёстывается.
См. рисунок: жёлтый цвет - полигон с ошибкой топологии, прозрачные круги - составные точки суррогатного полигона, который нужно удалить. Я в курсе, что ArcGIS for Desktop Advanced умеет фиксить ошибки топологии, но только для своего формата данных - GDB. Мне это не совсем подходит, лень конвертить шейп.

Перегонку через растр не предлагать, т.к. это очень затратно как в плане времени, так и ресурсов.

Пробовал v.clean в GRASS: не помогло. Странно, но в GRASS именно такая ошибка топологии почему-то не учтена, поправьте если ошибаюсь.

Пробовал плагин repair for Processing для QGIS. По-моему, он вообще не работает. Оригинальный repair, который требует доустановки софта от Hugo Ledoux, чем отличается?

PostGIS (ST_IsValid()) назвал это Ring self-intersection. Функция ST_MakeValid() не дала результатов.

Можно ли исправить именно эту ошибку топологии автоматом?

А ST_MakeValid вообще ничего не сделал с геометрией или превратил её в корректную геометрию с дыркой, касающейся внешней границы? Ещё ST_Buffer(geom, 0) можно попробовать.

ST_MakeValid() странно сработал - некоторые полигоны растянулись на несколько км в разные стороны.

Правильно ли я понял, что нужно удалить верхний треугольник серого цвета, вклинивающийся в желтый? А получившееся пустое место присоединить к желтому?

Если так, то по-моему в GRASS все-таки можно такое сделать, хотя и не полностью автоматом (потребуется играть параметрами). Нужно:
1. Создать узлы в пересечениях линий (есть инструмент для этого), если этих узлов еще нет.
2. Удалить области, меньшие по площади заданного порога.

Это если говорить о готовых инструментах. Возможно, что получится написать скрипт, который делает то, что нужно, если сформулировать по пунктам, как должен вести себя этот скрипт для разных входных данных.

Что-то наподобие этого?

KolesovDmitry писал(а):Правильно ли я понял, что нужно удалить верхний треугольник серого цвета, вклинивающийся в желтый? А получившееся пустое место присоединить к желтому?

Проблема двояка. Есть треугольники, которые просто суррогаты, т.е. возникли из-за плохоруких операторов, а есть треугольники с осмысленными данными, которые трогать нельзя.
Некоторые полигоны вместо серого треугольника не имеют ничего, т.е. буквально дыра. Значит, инструмент с удалением маленьких полигонов не решит проблему полностью.

Да, что-то похожее. Только я бы использовал break, чтобы создать узел на пересечении линий: Но на самом деле я не проверял код, нужно смотреть, что там есть среди инструментов

Делал двойное преобразование в Mapinfo. Сначала все полигоны преобразовал в полилинии; все линии объединил в одну; разъединить линии; создал полигоны из линий; отсортировал полигоны по площади; удалил самые мелкие. Помогло ~ на 90%.

По-хорошему, чтобы оценить возможность обработки таких полигонов в GRASS или ещё в каком-либо ПО, нужны тестовые данные. Было бы интересно поэкспериментировать.

Согласен, только данных на 1,5 Гб.
Временно, решил проблему через растр.

rhot писал(а):Согласен, только данных на 1,5 Гб.
Временно, решил проблему через растр.

Круто, конечно, только вот почему бы сразу не написать здесь это решение?

Потому что это не решение, а костыль.

rhot писал(а):Проблема двояка. Есть треугольники, которые просто суррогаты, т.е. возникли из-за плохоруких операторов, а есть треугольники с осмысленными данными, которые трогать нельзя.

Всё-таки нужен критерий, позволяющий отличить осмысленные треугольники от бессмысленных, и тестовый пример с хотя бы несколькими "хорошими" и "плохими" полигонами. Иначе задача не поставлена корректно.