Мы должны помнить, что эти данные являются образцами дискретных литологических доменов. Часто граница между двумя такими доменами не может быть идентифицирована в поле, и поэтому не стоит ожидать, что многие местоположения выборки будут лежать точно вдоль границ. Правильным решением будет раздел исследуемой области, и каждый многоугольник в этом разделе может (и часто будет) выходить за пределы местоположений образцов, которые его определяют. За исключением грубых приближений, это исключает любой подход, который использует местоположения выборки в качестве вершин результирующих многоугольников .
Для высококачественной работы наилучшим методом является подбор обобщенной линейной пространственной модели для многочленного процесса. Это процедура, которая требует значительных знаний и усилий. В качестве замены вы можете рассмотреть возможность расширения каждой точки выборки в ее полигон влияния (он же полигон Тиссена, полигон Вороного или ячейка Дирихле). Ограничение расширения на земельные участки является хорошей идеей; это можно сделать с помощью сетки маски.
Чтобы проиллюстрировать это, рассмотрим этот намного меньший набор данных (из 14 136 точек), представляющий 12 литологических классов, различаемых по цвету:
Вот деталь из центра восточной доли, показывающая неправильные положения точек и относительно быстрые изменения литологии там. Отслеживание этого вручную будет сложной и произвольной процедурой:
Я выполнил расширение, преобразовав эти точки в сетку (около 800 строк и 1000 столбцов) и вычислив их евклидово распределение , используя маску, которая ограничивала вычисление неглазированной землей. (Цветовая схема на следующих двух рисунках отличается от предыдущей.)
Для сравнения приведена подробная литологическая карта того же района, нарисованная в том же масштабе с той же символикой:
С действительно большим набором данных или извилистой областью исследования может оказаться целесообразным разбить область на части и выполнить эту процедуру отдельно для каждой плитки, объединяя результаты в один выходной растр, если это необходимо. Чтобы это работало, плитки должны слегка перекрывать друг друга, чтобы избежать краевых эффектов (а затем должны быть равномерно обрезаны перед нанесением мозаики).
Основными причинами перехода к растровому представлению являются (1) его быстрое и простое вычисление и (2) трудно будет найти точные векторные решения. Если вы попробуете буферы, выпуклые оболочки, вогнутые оболочки или что-то еще, вы обнаружите, что они все взаимно пересекаются и все еще оставляют пробелы: другими словами, они не будут создавать топологически согласованное разбиение пространства на отдельные литологические области.
Один векторный метод, который будет работать, состоит в том, чтобы вычислить ограниченную тесселяцию Вороного точек ( хорошие методы требуют времени O (n * log (n)) для n точек), пространственно объединить ячейки Вороного в соответствии с литологическими атрибутами их связанных точки, а затем разделить получающиеся мульти-полигоны на их связанные компоненты (если вы хотите). Однако, если все, что вам нужно, это векторный вывод , проще сгруппировать растровый результат и преобразовать его в векторный формат.