Многокритериальный анализ ГИС: взвешивание, нечёткая логика и анализ чувствительности

Карта пригодности может выглядеть безупречно — аккуратная компоновка, соблюдённая семиология Бертена, понятная легенда — и при этом целиком опираться на решения, которые не подкреплены ни одной цифрой: веса, присвоенные критериям, и пороги, выбранные для перехода от благоприятного к неблагоприятному. Я убеждался в этом снова и снова — сначала на примере поиска благоприятного участка в Веркоре в QGIS, затем в реальных проектах. Поэтому интересен не вопрос как построить многокритериальную карту — освоить инструмент несложно, — а вопрос, как выбрать правильный метод и доказать, что результат не рассыпается при первом же смещении веса.

Я руковожу CODRUM, бюро геоматики и веб-разработки в Кэлэрэшь, Республика Молдова. Эта статья — методическое руководство: стандартизация, взвешивание и, главное, анализ чувствительности — этап, который превращает красивую карту в результат, который можно отстоять. Она опирается на две статьи, которые я считаю отправными: взвешивание по методу AHP Саати в ArcGIS Pro для «как именно взвешивать» и кейс Веркора для булевой логики через ограничения.

Настоящая проблема не в инструменте, а в допущениях

Weighted overlay или растровый калькулятор выполняются за несколько минут. Трудность в другом — в решениях, принятых до того, как нажать кнопку. Первое и самое определяющее — разделить два семейства критериев, которые часто путают. Ограничения — это запреты, булева маска «да или нет», без полутонов: зона Natura 2000, аэронавигационный сервитут, граница парка. Факторы — это градуированные предпочтения, более или менее благоприятные, которые компенсируют друг друга: уклон, инсоляция, расстояние до дороги. Если трактовать ограничение как фактор или наоборот, искажается вся остальная цепочка.

Ограничение или фактор: решение, которое сдвигает гектары

Возьмём порог уклона. Свыше 30 % — это абсолютный запрет (ограничение) или всего лишь штраф, который могут перекрыть другие достоинства (фактор)? Ответ не технический, а профессиональный — и от него зависят сотни гектаров. В кейсе Веркора я трактовал уклон свыше 15° как жёсткое ограничение; превращение его в постепенный фактор вновь открыло бы часть территории для оценки. Ни тот, ни другой вариант не «верен» в абсолютном смысле: важно осознанно сделать выбор и иметь возможность его переиграть.

Ловушка компенсации

Взвешенная сумма складывает. Участок, катастрофический по жизненно важному критерию — скажем, по риску затопления, — может поэтому «подняться» в рейтинге за счёт отличных оценок в другом месте. Математически верно и операционно абсурдно. Решение сводится к простому правилу: то, что является отсекающим, остаётся ограничением (умножение на ноль, исключение на входе), а во взвешенную сумму входят только действительно компенсируемые критерии. Сначала фильтруешь по ограничениям, затем ранжируешь по факторам: эти два этапа дополняют друг друга, а не конкурируют.

Стандартизируйте перед взвешиванием

Нельзя взвешивать несопоставимые единицы. Уклон в процентах, расстояние в метрах, тип землепользования в категориях — всё это нужно сначала привести к общей шкале пригодности, как правило 0–1 или 0–255. Сосуществуют три семейства функций стандартизации, и сам выбор функции — это допущение, а не формальность.

Булева, по классам или нечёткая: компромисс, а не иерархия

Булев подход бинарен — благоприятно или нет — и непревзойдён по наглядности, но груб на порогах: участок при 14,9° проходит, его сосед при 15,1° отсекается, хотя на местности их ничто не различает. Переклассификация по классам (например, по Дженксу) сглаживает ступень дискретизацией на уровни ценой внутренних пороговых эффектов. Нечёткая логика заменяет резкую границу непрерывным переходом: она точно отвечает на произвольность жёстких порогов. Она не «лучше» в абсолютном смысле — это компромисс между обоснованностью и наглядностью. Когда я представляю результат депутату, булев вариант комментируется одной фразой; нечёткий честнее, но требует схемы, чтобы быть понятым.

Функции принадлежности нечёткой логики на практике

Конкретно: функцию принадлежности выбирают по форме предпочтения. Линейная возрастающая — для «чем больше, тем лучше»; убывающая — для обратного; трапециевидная или гауссова — когда есть ограниченный идеальный интервал. Числовой пример для виноградника: благоприятная высота от 800 до 1 800 м, плато максимальной пригодности между 1 100 и 1 400 м, линейный спад по обе стороны. Из инструментов это Fuzzy Membership в ArcGIS Pro (Spatial Analyst) и плагин FuzzyAttributes или растровый калькулятор в QGIS. Ничего эзотерического: одна задокументированная функция на критерий.

Назначайте веса, не обманывая себя

Три варианта в порядке возрастающей сложности. Равные веса: честное значение по умолчанию, когда нет ни одной веской причины выстраивать иерархию — лучше признать это, чем выдумывать ложную точность. Прямое оценивание по шкале важности: быстро, прозрачно, достаточно для трёх-четырёх очевидных критериев. Наконец, AHP — попарное сравнение по шкале Саати от 1 до 9, которое я подробно разбираю в статье, посвящённой ArcGIS Pro, чтобы не дублировать здесь.

Когда AHP действительно что-то даёт

AHP блистает, когда нужно проследить и обосновать разногласие между экспертами: матрица сравнений делает рассуждение проверяемым. На трёх критериях с очевидной иерархией она лишь создаёт ложное впечатление точности — три знака после запятой, чтобы прикрыть приблизительное суждение. Моё мнение из практики: AHP — как только группе нужно прийти к согласию или результат будет оспорен; прямое оценивание — когда взвешивание единодушно.

Проверьте согласованность: отношение согласованности (RC)

AHP даёт защиту, которую две существующие статьи лишь упоминают: отношение согласованности. Целевое значение — RC < 0,10. Выше — значит, эксперт противоречит сам себе: он рассудил A > B, B > C, но C > A. На практике, когда RC поднимается до 0,15 или 0,20, я возвращаюсь к матрице, нахожу несогласованную тройку (часто достаточно одного ошибочного суждения, чтобы сбить весь набор) и корректирую её вместе с экспертом. RC ниже 0,10 не гарантирует, что веса верны, — лишь что они согласованы между собой. Это различие важно.

Этап, который почти никто не показывает: анализ чувствительности

Вот сердцевина темы и точка, в которой большинство анализов останавливаются слишком рано. Карту пригодности без проверки устойчивости невозможно отстоять, потому что она зависит от весов, выбранных с долей произвола. Анализ чувствительности состоит в том, чтобы изменять каждый вес на ±10–20 % и смотреть, сохраняется ли иерархия лучших участков или она перестраивается. Одна карта не решает вопрос о собственной надёжности; анализ чувствительности — решает.

Воспроизводимый протокол

Я отталкиваюсь от базового сценария — номинальных весов. Затем строю альтернативные сценарии: минимальный и максимальный вес каждого критерия в диапазоне его неопределённости и сценарий с поочерёдным удалением каждого критерия (вес ноль), чтобы измерить его реальное влияние. Отслеживать нужно не абсолютное значение оценки — оно всегда немного меняется, — а стабильность вершины рейтинга: остаётся ли первая десятка той же первой десяткой, когда я встряхиваю веса? Именно такую цепочку ModelBuilder в ArcGIS Pro или графический модельер QGIS позволяют переигрывать, не переделывая всё вручную.

Чтение и представление результата

В итоге я различаю два профиля. Устойчивый участок остаётся на вершине рейтинга при любом сценарии — его я рекомендую в первую очередь. Неустойчивый участок появляется только при определённом допущении: он заслуживает упоминания, но никогда — твёрдой рекомендации. В формулировке для депутата или для инвестора это даёт чёткую фразу: «эти четыре участка держатся при любом раскладе, а эти два выходят вперёд, только если переоценить доступность». Именно этот нюанс позволяет отстоять исследование в комиссии или на публичных слушаниях.

Сквозная нить: полная цепочка на реальном кейсе

Чтобы всё стало конкретным, я возвращаюсь к структуре поиска благоприятного участка типа Веркора — кейсу, которым владею. Всё начинается с перепроецирования в EPSG:2154 (RGF93 / Lambert-93), без которого километровые буферы или пороги в один гектар не имеют смысла. На 83 коммунах исследования я выстраиваю цепочку: стандартизация уклона (gdal:slope) и экспозиции (gdal:aspect), булевы исключения (парк, застройка, существующая инфраструктура), затем взвешивание оставшихся факторов. Базовый сценарий давал 406 благоприятных участков площадью не менее гектара, во главе с Лан-ан-Веркором — 871 гектар. Добавленная ценность приходит потом: переигрывая анализ с переоценённой, а затем недооценённой доступностью, я проверяю, что первые коммуны не меняются местами, — иначе рейтинг не более чем артефакт весов.

QGIS или ArcGIS Pro для этой цепочки

Оба справляются нативно. В QGIS я выстраиваю gdal:slope, gdal:aspect, растровый калькулятор, FuzzyAttributes и заключаю всё в графический модельер. В ArcGIS Pro это Spatial Analyst (Reclassify, Fuzzy Membership, Weighted Overlay), ModelBuilder и ArcPy. Картографический результат строго эквивалентен — меняется лишь грамматика инструментов. Мой критерий выбора остаётся прежним — оснащение клиента: QGIS, когда нет лицензии ESRI, ArcGIS Pro, когда его специалисты будут поддерживать анализ.

Сохранение анализа воспроизводимым и проверяемым

Автоматизация — не удобство для разработчика: именно она делает анализ чувствительности осуществимым. Переигрывать десять сценариев вручную — рутина, которую никогда не делаешь; через ArcPy или модельер это вопрос нескольких минут. Поэтому результат включает карту и таблицу допущений: журналируемые веса, версионируемые пороги, явно описанные функции стандартизации. Прослеживаемость — не деталь в конце отчёта, это сам результат.

Ошибки, которые дискредитируют многокритериальный анализ

Несколько ловушек повторяются, и я совершил почти все, прежде чем научился их избегать. Коррелированные критерии — например, высота и температура — дважды взвешивают одну и ту же реальность и молча переоценивают фактор. Веса, подогнанные постфактум, до тех пор пока не получится ожидаемая карта: это и есть отрицание метода. Произвольные оценки, принятые за истину, хотя они лишь кодировка предпочтений. И отсутствие полевой проверки, которая пропускает «идеальные» на экране, но непригодные в реальности участки. Ни один из этих изъянов не виден на итоговой карте — отсюда важность документирования цепочки.

Перенос метода и рамки услуги

Рассуждение переносится как есть: размещение ветровых или фотоэлектрических установок, участок ZAC, зонирование PLUi, поиск земли, локализация чувствительных объектов. То, что выигрывает муниципалитет или проектное бюро, передавая эту работу на аутсорс, — не карта (её сделает кто угодно), а анализ, который можно отстоять: прослеженный метод, обоснованное взвешивание, проверенная чувствительность, карта, имеющая силу в споре. Это требование вписывается в мою работу по пространственному анализу для муниципалитетов и в прикладное многокритериальное исследование из моего портфолио.

О тарифах скажу прямо: картографический аудит технической возможности начинается от 750 € (≈ 14 600 леев) (первый проход по двум-трём критериям); полное многокритериальное ГИС-исследование — стандартизация, взвешивание AHP, анализ чувствительности, сводная карта 300 dpi и поставляемый ГИС-слой — от 2 300 € (≈ 45 000 леев). Итоговая цена зависит от размера зоны и числа критериев. Базируясь в Кэлэрэшь (Республика Молдова), я работаю по всей Республике Молдова и удалённо — в ЕС и по всему миру. Обсудим ваш проект — первичная консультация бесплатно, смета в течение 48 часов.

Полный перечень услуг и диапазоны цен — на странице «Услуги и тарифы».