Выделение класса(ов) из серии одинаково индексированных растров

Задача: Выделить один или несколько классов из нескольких одинаково индексированных растров. Например, класс лесов из одинаково индексированных (цветоделенных по одной легенде) топокарт.

Главным условием успешного выполнения является одинаковая индексация растров (т.е. универсальная легенда), помимо возможности использования таких расширений, как TRLegend, такая индексация позволяет производить в автоматическом режиме операции над классами, одновременно для неограниченного количества растров.

Сама по себе процедура достаточно проста, для Arcinfo она выглядит так (выделенным показаны команды, которые необходимо набрать в Arcinfo Workstation):

1. Конвертация растра в GRID
   imagegrid L53-10 L53-10
2. Переход в модуль GRID
   grid
3. Переклассификация полученного GRID'а
   L53-10r = reclass(L53-10, remap, nodata)
   remap - файл таблицы переклассификации созданный по определенному шаблону, например
   /*Присвоить индексу 0 и 1 значение 0
   0 1 : 0
   /*Присвоить индексу 2 значение 2
   2 : 2
   3 4 : 1
   /*Присвоить индексам с 5 по 16 значение 0
   5 16 : 0
4. Установка значения NODATA (по умолчанию отображается как прозрачный класс в Arcview 3.x, необязательный этап)
   L53-10z = setnull(L53-10r <= 0, L53-10r)
5. Склейка полученных GRID'ов в единый (необязательный этап)
   result = merge(L53-10z,L53-11z)

Небольшая программа Extract-Class позволяет автоматизировать этот процесс, ей достаточно показать директорию где лежат файлы TIF или GIS, из которых вы хотите выделить класс и она создаст sml-файл (файл набора команд для Arcinfo Workstation), который нужно будет запустить в Arcinfo. Так же вам понадобится создать свою таблицу переклассификации так как показано выше (или воспользоваться помощью программы, чтобы она создала образец).

У программы есть следующие настройки:

Convert zeroes to NODATA - Нули из таблицы обычно показываться как прозрачный класс при загрузке растров в Arcview 3.x. Усли вам нужно сохранить нулевые значения, снимите флажок.

Remove intermediate grids - Удалять промежуточные GRID'ы. В процессе работы Arcinfo Workstation создается некоторое количество промежуточных файлов, которые потом ею же и удаляются. Например результирующие GRID'ы склеиваются в один, а сами удаляются. Если вам необходимо их оставить - снимите флажок

Create example remap table - Создать образец таблицы переклассификации, который можно использовать для начала. Файл таблицы переклассификации называется remap, расположен в той же папке что и входные растры. Если он уже существует, то программа не будет его заменять и выдаст предупреждение.

Запустить полученный sml-файл на исполнение можно так:

1. Задать рабочую директорию, имя директории где лежат все исходные файлы и куда программа помещает созданный файл extractclass.sml, далее, после запуска Arcinfo Workstation:

2. Запустить пакетный файл, файл должен лежать там же где и HGT файлы

Результатом работы данного пакета является одно покрытие в формате GRID (result).

Обсудить в форуме (Комментариев - 0)

Последнее обновление: August 01 2007
(Наверх)