<<На Главную         <В меню Советы по ArcGIS

   Довелось недавно решать следующую задачу - нужно было определить площадь по интервалам агрохимического показателя (pH, если уж совсем конкретно) на каждом из контуров сельхозугодий. Исходные данные данные были такие: имелась интерполированная поверхность в виде растра, построенная по показателю кислотности почв pH (слой Idw_shp2) и собственно слой участков в виде полигонов (parc_region, как построить интерполированную поверхность можно узнать на нашем сайте здесь и здесь).

 

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

 

 Поначалу карта выглядела следующим образом:

Zonal Statistics ArcGIS

Решение этой задачи было найдено с использованием инструментов Зональной статистики (Zonal) из набора инструментов Пространственного анализа (Spatial Analyst).

Для начала нужно было переклассифицировать исходную растровую поверхность с помощью инструмента Reclassify:

Интерфейс утилиты Reclassify выглядит сл. образом:

Здесь на вход утилиты подаём растровую поверхность и заменяем интервалы параметра каким-то условным значением. В данном случае дискретным интервалам параметра pH присваиваем значения от 1 до 6 (6 интервалов).

После отработки указанной утилиты получаем новую, переклассифицированную поверхность - Reclass_Idw_4:

Использование Zonal Statistics

Так, поверхность переклассифицировали - заменили интервалы параметра pH значениями классов. Теперь можно подсчитать искомую статистику по площади по каждому контуру. Сделать это можно с помощью утилиты Tabulate Area:

Дважды жмём на указанный пункт и видим сл. форму:

Здесь в качестве входных параметров используется переклассифицированная поверхность (Reclass_Idw_4) и слой наших учебных участков (parc_region).

После отработки утилита создаёт таблицу с нужной нам статистикой по каждому участку:

Зональная статистика ArcGIS

В ячейках таблицы видим значение площади в размерностях карты по каждому контуру по всем интервалам ( в нашем случае от 1 до 6).

Далее с помощью инструмента Zonal Statistics можно узнать полную статистику по каждому участку:

Интерфейс указанной утилиты имеет сл. вид:

Здесь в качестве входных параметров используется переклассифицированная поверхность Reclass_Idw_4 и слой наших учебных участков parc_region.

По опции ALL параметра Statistics type в результирующую таблицу будут помещены все статистики: разброс, минимум, максимум и т.д.

Жмём кнопку OK и получаем приблизительно сл. картинку:

На рисунке выше красные стрелки соответствия записи в таблице контуру на карте  были нанесены для наглядности.

©Simashkov.2011

10.07.2011 г.

Назад