Как найти пересечение объектов

При анализе картографических данных часто бывает нужно найти объекты из одного слоя, относящиеся к определенным объектам в другом слое. Например, у нас есть слой c административными районами и слой с дорожной сетью, и мы хотим выделить все дороги, относящиеся к одному из районов.

Для этого в MapInfo Pro можно написать небольшой универсальный скрипт на MapBasic, который можно будет сохранить и использовать в дальнейшем во всех таких случаях. Скрипт будет сравнивать выбранные в данный момент объекты (текущую выборку) с объектами из указанной таблицы, и будет отбирать только те записи, которые пересекаются.

Создание скрипта

Откроем окно «SQL», добавленное в MapInfo v2019, из меню на вкладке «Объекты» и создадим в нем следующий скрипт.

Давайте рассмотрим запрос немного подробнее.

Таблица, которую мы указали для выбора — sTabToQuery, является не конкретной таблицей, а переменной. Создадим эту переменную с помощью «Редактирования переменных» из раскрывающего списка «Переменные» в окне SQL.

Тип переменной — Таблица. Это означает, что при запуске скрипта MapInfo Pro предложит выбрать таблицу из списка всех открытых таблиц. В поле «Значения» укажем следующее выражение: TableInfo(ScriptContext(SCRIPT_CONTEXT_TABLEALIAS), TAB_INFO_MAPPABLE), чтобы ограничить список таблиц только таблицами с возможностью отображения. В выражение также можно добавить и другие дополнительные условия на ваше усмотрение, например, чтобы в списке были только TABы или нельзя было выбрать растры.

Сохраним новую переменную и вернемся назад в окно SQL.

Обратите внимание на две вещи в «Where» части запроса. Первый заключается в том, что мы указали Intersect в качестве пространственного оператора. Если вы знаете, что будете сравнивать точки с полигонами, вы также можете использовать оператор Within/Contains. Почему не стоит использовать Within/Contains для полигонов и полилиний? Потому, что эти операторы сравнивают центроиды объектов, а центроиды полигона и полилинии зачастую не дают хорошего представления об объекте, т.к. находятся вне его пределов.

Второе, на что следует обратить внимание, это то, что мы используем подвыборку, в которой объединяем объекты из текущей выборки в один объект с помощью пространственного агрегата AggregateCombine. Объединение выбранных объектов обычно выполняется очень быстро. Особенно, если вы используете небольшие полигоны, такие как зоны транспортной доступности или административные границы. Однако, если вы используете полигоны большего размера — например, огромные полигоны покрытия телекоммуникационных компаний, можно заметить, что запрос выполняется медленее обычного.

Не забудьте указать имя для скрипта, а также сохранить его.

Использование скрипта

Выберем интересующий нас площадной объект, например, парковую зону, как в данном примере, состоящую из нескольких полигонов.

Далее из раскрывающегося списка «SQL» на вкладке «Объекты» выберем и запустим созданный нами скрипт.

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

В результате выполнения скрипта создастся новый запрос с объектами из указанного слоя, пересекающими изначальный выбранный объект.

Заключение

Несмотря на то что запрос составлен правильно, он может не дать вам именно то, что вы ищете. Данный запрос возвращает любую запись, которая каким-либо образом пересекает выборку. Это означает, что возвращаемые записи могут лишь частично пересекаться с выборкой. Может быть и такое, что они вообще имеют только один общий узел, формально это тоже считается пересечением.

Поэтому обратите внимание: если вы используете результат этого запроса, то будете использовать объект целиком, а не только ту его часть, которая пересекает выборку. То же самое относится и к значениям атрибутов. Возможно, вы захотите воспользоваться функцией ProportionalOverlap для вычисления пропорциональных значений. Если вы знаете, какие столбцы есть в вашей таблице, вы можете также включить эти пропорциональные вычисления непосредственно в запрос. То есть вместо использования простого оператора Intersect можно добавить более сложные методы для принятия решения о включении или невключении записи в набор результатов.

Другие полезные статьи:

Инструмент Пересечение (Intersect) вычисляет геометрическое пересечение любого количества классов пространственных объектов и векторных слоёв. Пространственные объекты или части пространственных объектов, общие для всех входных данных (т.е. они пересекаются), будут записаны в выходной класс объектов.

Пересечение делает следующее:

• Определяет пространственную привязку для обработки. Она также будет пространственной привязкой выходного класса объектов. Более подробно о том, как это происходит, см. в разделе Пространственная привязка (Spatial Reference). Все входные объекты проецируются в эту пространственную привязку для обработки.
• Выполняет разбиение и кластеризацию пространственных объектов. Разбиение добавляет вершины в местах пересечении ребер пространственных объектов; кластеризация состыковывает вершины в пределах допуска xy.
• Обнаруживает геометрические отношения (пересечения) между пространственными объектами из всех классов пространственных объектов или слоёв.
• Записывает эти пересечения как пространственные объекты (точки, линии или полигоны) в выходные данные.

Если в списке входных объектов задано несколько классов пространственных объектов или слоёв, порядок объектов в списке не влияет на тип выходных объектов, но пространственная привязка верхнего объекта в диалоговом окне инструмента (первый объект в скриптах) в списке будет использоваться в процессе обработки и будет установлена на выходные данные.

Чтобы контролировать выходную пространственную привязку (систему координат и домены), установите правильно параметры среды: Выходные Z значения (Output Z Aware) и Выходные M значения (Output M Aware).

Входными данными может быть любая комбинация типов геометрии (точки, мультиточки, линии или полигоны). Геометрия выходного класса объектов должна быть того же измерения, что и входной класс объектов, или на порядок ниже (точка = 0-е измерение, линия = 1-е измерение, полигон = 2-е измерение). Определение разных типов геометрии для выходных данных дает разные типы пересечений входных классов объектов. Это не разное представление одних и тех же пересечений; это пересечения, которые могут быть представлены соответствующим типом геометрии (точки, линии или полигоны).