Как взять векторный слой с перекрывающимися полигонами и превратить его в растр, где каждая ячейка считает количество полигонов, содержащих эту ячейку?
Мы смотрим на полигоны областей горения огня в период с 1976 по 2000 год и хотим узнать, сколько пожаров будет содержаться в каждой растровой ячейке.
Ответы:
Сделайте это в три этапа: разбейте полигоны на составные части, подсчитайте перекрытия и преобразуйте в растр. Это позволяет избежать потенциально больших вычислительных затрат на отдельное преобразование каждого многоугольника в растр и объединение этих растров.
Union(в Geoprocessingменю) разбивает многоугольники на их части.
К сожалению, каждое перекрытие дублируется в выходных данных: у него есть одна идентичная копия для каждого исходного полигона, покрывающего его. Следовательно
Dissolve(снова в Geoprocessingменю) объединит перекрывающиеся части, при условии, что вы сможете найти способ их уникальной идентификации. Прочитайте диалог: ближе к концу у вас будет возможность вычислить «статистику». Выберите любое поле, которое, возможно, идентифицировало исходные полигоны, и запросите счет.
Во многих случаях сочетание площади многоугольника и периметра позволяет однозначно идентифицировать детали. Если нет, вы можете добавить дополнительные геометрические свойства в дополнительные поля, такие как координаты центроида, до тех пор, пока вы не накопите достаточно информации, чтобы различить каждый объект.
Результирующий слой имеет одну особенность для каждого перекрытия полигонов и какое-то поле «считать», считающее количество перекрытий.
Преобразуйте это в растр, используя поле «count» для атрибутов.
Например, вот некоторые перекрывающиеся полигоны и их идентификаторы с показанной таблицей атрибутов:
После второго шага у нас есть одна запись для каждой перекрывающейся области вместе со счетчиком, который уже можно использовать для обозначения величины перекрытия:
Остальное просто - и это всего лишь одна операция растеризации.
Unionзаключается в том, что при почти одинаковом рабочем процессе он может разрешить ситуацию, когда полигоны находятся в нескольких наборах данных (что обычно является плохим дизайном базы данных, но, к сожалению, распространенным): просто объединяются все входные наборы данных одновременно.
union) с последующей растеризацией? Векторные операции прекратятся, когда объекты будут оцифрованы с большим количеством деталей, чем необходимо, наделив их слишком большим количеством вершин. В этих экстремальных условиях растровый подход мог бы быть лучше (хотя первое упрощение многоугольников могло бы быть лучшим вариантом). Однако во всех других ситуациях растеризация каждого многоугольника в отдельности является огромной тратой компьютерного и человеческого времени.
В следующем посте есть несколько похожий вопрос, который вы могли бы найти связанное решение: Создание растровой поверхности перекрытий из шейп-файлов векторного полигона? ,
Для этого быстродействующего и простого в вычислительном отношении растрового подхода вам необходимо (1) разделить перекрывающиеся полигоны на отдельные слои (предположительно по годам для пожарных полигонов), используя Select By Attributes в ModelBuilder с итератором или инструментом сценария, (2) преобразовать каждый Полигон в растр с назначением ячейки MAXIMUM_COMBINED_AREA (обеспечение того же размера ячейки , растров привязки и того, что экстент остается тем же, что и весь набор полигонов) - с использованием постоянного значения поля (например, используйте поле года или строки, которые имеют все 1) для преобразования (снова используйте ModelBuilder с итератором или скриптом Python для автоматизации), а затем (3) примените следующий инструмент Spatial Analyst: Статистика ячеек - используйте статистический тип VARIETY, если каждый растр имеет уникальное значение, такое как год или SUM, если значения всех ячеек растра равны 1 с - обязательно установите флажок, чтобы игнорировать NoData.
Промежуточные растры (из предыдущего преобразования) могут быть удалены или готовы к использованию в последующих растровых анализах.