Описание работы со скриптом получения данных по полигональным объектам в Arcview GIS для сравнения приводится также описание как выполнить ту же операцию в ERDAS IMAGINE.
Обсудить в форуме Комментариев 0
Для некоторых задач, связанных со статистическим анализом тренировочных данных, таких как классификации при помощи деревьев регрессии (decision - regression trees) или метода максимального подобия (maximum likelihood) требуется извлечение данных о тренировочных (эталонных) объектах не по полигону в целом как в зональной статистике, а по каждому пикселю изображения попадающему в этот полигон. Как правило, для этих задач используются многоканальные растровые данные, и, соответственно, каждый пиксель имеет N+2 значений, N = число каналов + координата X и координата Y.
В качестве профиля используется полигон или полигоны в формате shape-файла.
Статистика подсчитывается для растрового файла, загруженного как Image Theme в содержание вида. Для того чтобы получить данные из растров в формате ERDAS IMAGINE нужно загрузить их как Image Theme предварительно включив расширение IMAGINE Image Support.
Заметим, если растровые уже находятся в формате GRID и профили представляют из себя набор точек, то проще всего воспользоваться другой схемой работы и другим расширением.
В процессе работы, для каждого полигонального объекта в shape файле, для каждого канала растра, будет произведена операция получения данных из растрового слоя. Растр будет разбит поканально, каждый канал будет представлен в формате GRID, далее для каждого полигона будут определены попадающие в него ячейки растра (поканально) и полученные данные будут экспортированы в новый файл. При получении данных, если одна из точек попадает на пиксел со значением NODATA значение получено не будет. Сгенерированные в процессе работа наборы данных в формате GRID будут удалены в случае успешного выполнения операции.
Если в полигональной теме присутствует поле Id и в нем есть значения, они будут использованы как названия классов и добавлены дополнительной, последней, колонкой в результирующий файл. Рассмотрим пример: для некого изображения оцифрованы области (полигоны) представляющие из себя участки растительности типа A и участки растительности типа B. Тип растительности указан в атрибутивной таблице кодом 1 или 2. С помощью данного расширения можно для каждого полигона в наборе эталонов получить значения всех пикселей попадающих в этот полигон поканально + добавить колонку с кодом из атрибутивной таблицы темы эталонов. Такие данные потом очень просто импортировать и использовать в любом статистическом пакете, например R.
Скрипт может работать как со спроектированным, так и неспроектированным видом.
Таким образом, для получения пиксельных данных в Arcview GIS необходимо:
Примечания:
ERDAS IMAGINE
Операцию экстракции попиксельных значений во всех каналов для набора полигонов-эталонов можно также осуществить в ERDAS IMAGINE. Utilities\Convert Pixels to ASCII...
По сравнению с предложенным нами способом этот инструмент обладает следующими недостатками:
Обсудить в форуме Комментариев 0
Последнее обновление: September 09 2021
Дата создания: 06.11.2007
Автор(ы): Максим Дубинин
© GIS-Lab и авторы, 2002-2021. При использовании материалов сайта, ссылка на GIS-Lab и авторов обязательна. Содержание материалов - ответственность авторов. (подробнее).