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

   В п. 6  меню Советы по MapInfo нашего сайта уже было рассказано о том, как можно построить интерполированную поверхность (surface) используя ГИС MapInfo - эту информацию можно посмотреть здесь. Однако, если вы располагаете  ArcGIS в комплектации ArcMap, то можно решить такую же задачу используя и эту ГИС. Сразу необходимо оговориться, что лицензии ArcView на выполнение этой задачи не хватит (некоторые инструменты будут недоступны). Но, с другой стороны, возможностей ArcMap вполне достаточно, чтобы решить подобный вопрос   от начала и до конца. Итак, начнём ...

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

 

 

 

Исходная  картинка может выглядеть так (как пример):

Задача - построить интерполированную поверхность (растр) по одному из показателей. Ну, допустим, по кислотности почв - pH.

Для выполнения этого задания необходимо, чтобы наши данные были представлены в виде слоя точек (points). Имея исходный слой в виде линий,  такое требование проще всего выполнить путём 'развала' линий на узлы (vertices). Для этого воспользуемся соответствующей утилитой Feature Vertices To Points из набора инструментов Feature :

Загружаем в утилиту нужный нам набор полилиний и смело жмём кнопку OK:

Утилита Feature Vertices To Points

В результате получаем слой точек, с прикреплённой к нему атрибутивной информацией :

Таак ... Предварительную работу провели. Теперь можно приступать к построению поверхности. Для этого выбираем инструмент создания интерполированной поверхностей методом IDW (метод обратных взвешенных расстояний) из набора инструментов Spatial Analyst (Пространственный анализ):

Метод IDW выбран исключительно из учебных соображений. Никто не мешает, естественно, поэкспериментировать с другими методами интерполяции - методом Кригинга(Kriging),  методом ближайшего соседства (Natural Neighbor) или сплайн-интерполяции(Spline). Вообще, выбор наиболее подходящего метода интерполяции в зависимости от характера исходных данных - это целое искусство. О сравнении методов интерполяции, реализованных в ArcGIS, можно посмотреть  вот здесь

 

Интерфейс у этой утилиты следующий:

Утилита IDW

 Здесь в качестве параметра Z выбираем значения pH. Размер сетки для интерполяции 40 - в измерениях карты ( в данном случае - метры). В качестве входного слоя - слой полученных после "развала" узлов.

 После отработки указанной утилиты наша картинка может выглядеть так :

Прим. Существует ещё один способ построения интерполированной поверхности. Если интересно - вам сюда

Теперь для придания нашей карте более профессионального вида построим изолинии. Для этого выбираем инструмент Contour List из набора инструментов Surface:

Вид у этой утилиты такой:

Утилита Contour List Здесь с помощью кнопки "+" добавляем значения кривых  равного уровня.

В результате получаем такую картинку:

Но результат построения изолиний в некоторых локализациях (обведено красным на рисунке выше) нас не совсем устраивает. Получились не очень гладкие кривые. Чтобы выправить этот недостаток, можно воспользоваться инструментом Smooth Line (Сглаживание линий):

Интерфейс этого инструмента таков:

Утилита Smooth Line в ArcGIS

В качестве параметра выглаживания выбираем 400 единиц ( в данном случае, это метры-координаты карты).

Окончательный итог всей нашей работы с нанесёнными значениями изолиний может выглядеть так - сравните результат до и после выглаживания:

Создание изолиний и поверхностей в ArcGIS

Ссылки к статье

1. Как  ещё можно построить интерполированный растр - смотрите здесь

2. Как решить аналогичную задачу по построению тематической поверхности используя ГИС MapInfo, можно посмотреть на нашем сайте здесь

3. Как сделать тематическую поверхность в виде "отмытого"  рельефа - смотрим  тут

4. Информацию о построении тематической карты в программе Surfer смотрите на нашем сайте по этой ссылке

Назад

 

©Simashkov.2010

31.12.2010 г.