Обработка многозональных космоснимков в MultiSpec: различия между версиями
(Новая страница: «{{Статья|Черновик}} {{Аннотация|Использование выражений на языке Python для операций с полям…») |
Нет описания правки |
||
Строка 1: | Строка 1: | ||
{{Статья|Черновик}} | {{Статья|Черновик}} | ||
{{Аннотация| | {{Аннотация|В статье описываются основные функции бесплатной программы MultiSpec, используемые для обработки, интерпретации и анализа многозональных космоснимков.}} | ||
Программа MultiSpec была разработана и поддерживается Дэвидом Лэндгрибом (David Landgrebe) и Ларри Бьелем (Larry Biehl) из университета Пердью (США). MultiSpec является результатом продолжающихся многолетних исследований, направленных на получение мощного и надёжного инструмента для анализа многозональных и гиперспектральных данных. MultiSpec вместе с документацией доступна для свободного скачивания с официального сайта программы . | |||
К основным функциям программы относятся [1]: | |||
* Импорт данных, как из двоичных, так и текстовых данных. | |||
* Визуализация многозональных изображений с помощью множества черно-белых и цветных режимов; показ создаваемых типизированных растров. Наложение векторных данных поверх изображений. | |||
* Расчёт гистограмм по многозональным изображениям. | |||
* Переформатирование данных, включающее изменение каналов, объединение файлов, создание мозаик, изменение геометрии и др. | |||
* Создание новых каналов из существующих с помощью извлечения информации из многозонального изображения или в результате попиксельных алгебраических вычислений. | |||
* Кластерный анализ (классификация без обучения) в отношении многозональных изображений и создание в результате типизированных растров. | |||
* Классификация с обучением после подготовки набора эталонов для основных классов наземных объектов, видимых на многозональном изображении. | |||
* Построение графика спектральных значений для выделенных участков многозонального изображения. | |||
* Экспорт исходных данных и их трансформант в текстовый файл. | |||
* Экспорт изображений в графический файл и буфер обмена. | |||
== Объединение файлов отдельных каналов в единое многозональное изображение == | |||
Многозональные спутниковые снимки поступают из хранилища ГС США в виде упакованных в архив монохромных файлов формата GeoTIFF (*.tif). Причём каждый файл представляет отдельный канал многозонального изображения. Такой подход к предоставлению материалов имеет как плюс – пользователь может ограничить объём обрабатываемых данных только нужными ему каналами, так и минус – для дальнейшей работы каналы должны быть предварительно объединены в многоканальный файл. | |||
Создание многозонального изображения выполняется в программе MultiSpec следующим образом [2]: | |||
1. Выполнить команду из главного меню File Open Image… или щёлкнуть по кнопке Open на панели инструментов. Появится диалоговое окно Open Image File. | |||
2. Выделить в окне те файлы каналов, которые должны быть включены в многозональное изображение (рис. 1). | |||
Рис. 1. Выделение нескольких одноканальных файлов перед объединением | |||
Можно обратить внимание на то, что при выделении более чем одного файла опция в нижнем левом углу переключается в Link selected files to new image window (Привязать выделенные файлы к новому окну изображения). | |||
3. Щёлкнуть по кнопке Открыть. Появится диалоговое окно Set Display Specifications for (Параметры визуализации) (рис. 2). Это окно позволяет устанавливать большое число параметров, используемые при показе собираемого изображения. | |||
Рис. 2. Окно установки параметров визуализации изображения | |||
По умолчанию для показа представляется полный размер сцены (параметры Area to Display) и тип изображения выбран 3-Channel Color (3-канальный цветной). | |||
4. Щёлкнуть по кнопке OK. В появившемся окне Histogram (Гистограмма) тоже щёлкнуть по кнопке OK. В результате в отдельном графическом окне появится многозональное изображение, оформленное в соответствии с установленными параметрами. | |||
== Управление визуализацией многозонального изображения == | |||
Открытое в программе MultiSpec с параметрами по умолчанию многозональное изображение можно визуализировать как-нибудь иначе с помощью группы опций, доступных с помощью команды Processor Display Image… главного меню. К этим опциям относятся Bits of color (глубина цвета), Stretch (кривая яркости), Min-max (пределы), Treat ‘0’ (представление нуля) и Display levels per channel (число уровней яркости) в диалоговом окне Set Display Specification (рис. 2.2). | |||
Назначение этих параметров следующее [3]: | |||
* Параметры Bits of color и Display levels per channel по умолчанию имеют максимально возможные значения (24 и 256 соответственно). При желании их можно уменьшить, что может быть полезно для выполнения простейшей кластеризации изображения (разделения на объекты, принципиально различающиеся по своим свойствам: пашня, леса, вода, строения и т.п.) (рис. 3). | |||
Рис. 3. Визуализация изображения с 256 уровнями яркости (миллионы цветов) и 2 уровнями яркости (8 цветов) | |||
* Параметр Treat ‘0’ позволяет управлять интерпретацией нулевых значений в исходных файлах каналов материалов Landsat. Изначально нулевые значения обозначают фон, который не используется при дальнейших вычислениях и может быть показан чёрным (black) или белым (white) цветом. Однако при желании можно нулевые значения считать в качестве данных с нулевой яркостью (as data) (рис. 4). | |||
Рис. 2.4. Изображение фона (нулевых данных) чёрным или белым цветом | |||
* Параметры Stretch и Min-max обычно используются для улучшения визуализации отдельных частей изображения. Они контролируют соотношением между значениями яркости в исходных файлах каналов и уровнями яркости на экране. Есть три варианта параметра Stretch: Linear (линейная), Equal Area (равноплощадная) и Gaussian (гауссова) (рис. 5). | |||
Рис. 5. Три варианта кривой яркости (линейная, равноплощадная и гауссова) | |||
** В случае линейной кривой яркости интервалы серого цвета на экране равномерно распределены между значениями яркости в исходных данных. | |||
** В то же время равноплощадная кривая автоматически подбирает такое соотношение между ними, чтобы каждый уровень яркости занимал на экране одинаковое число пикселей. | |||
** Гауссова кривая яркости использует гауссово распределение уровней серого цвета на экране, что позволяет улучшить изображения слабосветящихся объектов. | |||
* Опция Min-max позволяет установить предельные значения (и диапазон между ними) для исходных данных, которые будут использоваться при показе изображения на экране. Уровни яркости и кривая яркости будут ограничены этими пределами. | |||
** Вариант Entire Range (весь диапазон) для параметра Min-maxes ставит в соответствии минимальное значение исходных данных (0) чёрному цвету, а максимальное (255) – белому. Однако иногда требуется не показывать значения, например, меньше 50 и больше 150. | |||
** Вариант 2 Percent Tails Clipped (обрезать двухпроцентные края) позволяет сделать так, чтобы визуализировались значения исходных данных только в диапазоне 2-98% от их предельных значений. Это позволяет увеличить контрастность изображения за счёт исключения из показа редко встречающихся самых тёмных и самых светлых точек. | |||
** Если выбрать вариант User Specified..., то появится диалоговое окно Set Specifications for Display Enhancement (Установить параметры кривой яркости). Это окно позволит задать свой собственный процент отсечения краёв (рис. 6), который может отличаться от 2%. | |||
Рис. 6. Окно установки пределов кривой яркости | |||
1.3. ВЫДЕЛЕНИЕ УЧАСТКОВ НА МНОГОЗОНАЛЬНОМ ИЗОБРАЖЕНИИ | |||
Если в программе MultiSpec открыто изображение, то часть его можно выделить для того чтобы проводить дальнейшие преобразования только с этой частью, а не со всей картой. Существует 3 способа создания выделения [4]: | |||
1. Выполнить команду Edit Selection Rectangle (Рамка выделения) главного меню программы. Появится соответствующее диалоговое окно. В этом окне надо ввести или отредактировать значения, задающие минимальное и максимальное значение координат в выбранной системе координат. Имеется 3 варианта систем координат: | |||
• Lines-Columns (строки и столбцы) – локальная система координат данного изображения, основанная на нумерации его пикселей от верхнего левого угла вниз (строки) и направо (столбцы). | |||
• Transverse Mercator-Meters (глобальная прямоугольная поперечная проекция Меркатора с метрами в качестве единиц измерения). | |||
• Latitude-Longitude (Широта-Долгота) – глобальная сферическая система координат с единицами измерения – градусами (Рис. 2.7). | |||
Рис. 2.7. Окно задания параметров рамки выделения | |||
Кнопка Preview (просмотр) позволяет увидеть на карте сделанные изменения выделения. | |||
2. С меньшей точностью, но быстрее рамку выделения можно нарисовать по карте прямо с помощью указателя мыши, ориентируясь на значения координат, которые можно видеть на панели Coordinate View (координаты) (Рис. 2.8). | |||
Рис. 2.8. Панель координат и выделение с помощью указателя мыши | |||
1.4. СОХРАНЕНИЕ В ЕДИНОМ ФАЙЛЕ СОБРАННОГО МНОГОКАНАЛЬНОГО ИЗОБРАЖЕНИЯ | |||
Постоянно собирать многоканальное изображение из разрозненных файлов каналов при каждой необходимости использования – нерационально. Кроме того часто после сборки требуется продолжить работу не со всей 185-км сценой, а только с областью интереса, которая обычно значительно меньше. Существенно упростить работу может сохранение собранного многозонального изображения или его части в отдельный файл. Для этого потребуется следующее [1]: | |||
1. Собрать многозональное изображение из тех каналов, которые планируется использовать в дальнейшем. | |||
2. Выделить часть изображения, представляющую область интереса или удалить выделение, если должна быть сохранена сцена целиком. | |||
3. Выполнить команду Processor Reformat Rectify Image… (Сохранить прямоугольный фрагмент). Появится диалоговое окно Set Image Rectification Parameters (Установить параметры прямоугольного фрагмента). Все параметры этого окна имеют по умолчанию значения, изменять которые обычно нет необходимости. | |||
4. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. Надо ввести имя и расширение “.tif” для этого файла. По неизвестной причине MultiSpec сама в данном случае расширение к имени файла не добавляет, что впоследствии мешает ей же открывать такие файлы. | |||
1.5. ПОЛУЧЕНИЕ СВОДНОЙ СТАТИСТИЧЕСКОЙ ИНФОРМАЦИИ ПО МНОГОЗОНАЛЬНОМУ ИЗОБРАЖЕНИЮ | |||
Пункт меню Histogram image (гистограмма изображения) позволяет вычислить и вывести в текстовом окне сводную статистическую информацию по всем каналам многозонального изображения или его части. Причём это может быть как пространственная часть изображения, так и выборка из спектральных каналов [1]. | |||
1. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение. | |||
2. Выполнить команду Processor Histogram image. Появится диалоговое окно Set Histogram Specifications (Установить параметры гистограммы) (Рис. 2.9). | |||
Рис. 2.9. Окно установки параметров вычисления гистограммы изображения | |||
3. Выбрать метод (Method) вычисления гистограммы – Compute new histogram (Вычислить новую гистограмму). | |||
4. С помощью списка Channels (Каналы) можно выбрать каналы, используемые для вычислений. По умолчанию установлен вариант All (Все). | |||
5. Для получения сводной информации надо включить опцию List histogram summary (вывести сводную гистограмму) и выключить List histogram (вывести всю гистограмму). | |||
6. При нажатии на кнопку OK в текстовом окне появится информация похожая на приведённый ниже пример (Табл. 2.1). | |||
Табл. 2.1. Пример таблицы сводной статистической информации | |||
Channel Channel | |||
Description Data Range Mean Median Standard | |||
Deviation | |||
1 0.45-0.52 um 39 to 50 43.1 43 1.6 | |||
2 0.53-0.61 um 15 to 21 17.8 18 1.2 | |||
3 0.63-0.69 um 12 to 21 15.6 15 2.0 | |||
4 0.78-0.90 um 24 to 59 41.7 41 6.4 | |||
5 1.55-1.75 um 22 to 54 38.2 39 6.0 | |||
6 2.09-2.35 um 10 to 23 15.4 15 2.4 | |||
7 10.4-12.5 um 103 to 111 106.7 107 1.8 | |||
1.6. СОЗДАНИЕ СХЕМ РАСПРЕДЕЛЕНИЯ NDVI | |||
В программе MultiSpec можно создавать трансформанты многозональных снимков – изображения, представляющие собой алгебраические комбинации значений исходных файлов каналов материалов Landsat. При этом чаще всего создают схемы распределения нормализованного дифференциального вегетационного индекса (NDVI). Назначение этого индекса – улучшить читаемость и представительность изображения путём подчёркивания особенностей распределения растительности и почвы. | |||
Создаётся такая схема следующим образом [5]: | |||
1. В меню Processor (процессор) выбрать команду Reformat Change Image File Format… (изменить формат файла изображения). В результате появится соответствующее диалоговое окно (Рис. 2.10). | |||
Рис. 2.10. Окно установки параметров изменения формата файла | |||
2. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение. | |||
3. Отметить опцию Transform Data… (преобразовать данные). В результате появится соответствующее диалоговое окно. | |||
4. Выбрать опцию New Channel from General Algebraic Transformation (создать новый канал с помощью алгебраического преобразования). Окно изменит свой интерфейс (Рис. 2.11). | |||
Рис. 2.11. Ввод формулы алгебраического преобразования каналов многозонального снимка | |||
Ячейки для ввода значений в этом окне позволяют задать формулу алгебраической комбинации исходных каналов. С их помощью можно создать новый канал, например NDVI, если ввести после символа “C” номера каналов, соответствующие красному и ближнему ИК каналу. | |||
Если многозональный снимок содержит все каналы сцены Landsat, то это будут: канал 1 – синий, канал 2 – зелёный, канал 3 – красный, канал 4 – ближний инфракрасный и т.д. Тогда в ячейки окна надо вводить формулу следующим образом: | |||
5. Далее щёлкнуть по кнопке OK, что применит введённую формулу преобразования, закроет окно Set Reformat Transform Parameters и вернёт фокус в окно Image File Format Change Options. | |||
6. В окне Image File Format Change Options для параметра Data value type (тип значений данных) задать вариант 32-bit Real (вещественный). Программа может и сама установить этот тип, если определит, что результат вычислений будет не целочисленным. | |||
7. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. Надо ввести имя и расширение “.tif” для этого файла. | |||
8. После этого можно просмотреть получившейся результат, если при открытии файла задать тип изображения 1-channel thematic (1-канальный тематический) (Рис. 2.12). | |||
Рис. 2.12. Схема распределения NDVI | |||
На данном изображении сине-зелёные оттенки соответствуют воде, красно-оранжевые – растительности разной пышности. Остальное – почво-грунты. | |||
1.7. ЭКСПОРТ ИЗОБРАЖЕНИЯ В ТЕКСТОВЫЙ ФАЙЛ | |||
Нередко изобразительных свойств программы MultiSpec бывает недостаточно. Возникает желание передать информацию в другую, более подходящую, программу, например, Golden Software Surfer. Для этого надо экспортировать открытое в MultiSpec изображение в текстовый файл, который впоследствии может быть открыт в Surfer. | |||
Экспорт в текстовый файл может быть выполнен так [1]: | |||
1. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение. | |||
2. Выполнить команду Processor List Data… Появится диалоговое окно Set List Data Specifications (Установить параметры экспорта данных) (Рис. 2.13). | |||
Рис. 2.13. Параметры экспорта данных | |||
3. Включить переключатель Include latitude and longitude values (включить значения широты и долготы). | |||
4. В группе Output results to (Вывести результаты в) убрать отметку опции Text window (Текстовое окно) и включить отметку опции Disk file (Дисковый файл). | |||
5. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. | |||
Этот файл имеет простой текстовый формат. В нём сформирована таблица, содержащая как координаты каждой точки изображения, так и значения всех каналов в ней. | |||
Рис. 2.14. Пример построения распределения NDVI средствами Surfer | |||
На данном изображении благодаря богатым возможностям визуализации двумерных данных у Surfer применены более естественные цвета (зелёный – пышная растительность, жёлто-зелёный – редкая растительность, коричневый – грунт, синий – вода) с точно подобранными рубежами перехода между цветами. | |||
Рассмотрение технологии использования Surfer для подобных построений находится за пределами данного пособия. | |||
1.8. НАЛОЖЕНИЕ ВЕКТОРНЫХ ДАННЫХ В ОКНЕ ИЗОБРАЖЕНИЯ | |||
Векторные данные – это подготовленные в других программах наборы графических примитивов (полигонов, линий, точек), в совокупности представляющие собой чаще всего элементы топографической карты (реки, дороги, населённые пункты и т.п.). Векторные данные, наложенные на многозональное или классифицированное изображение, помогают легче ориентироваться по нему и быстро находить области интереса. | |||
MultiSpec имеет возможность открывать векторные данные из shape-файлов (*.shp) [6]. После открытия хотя бы одного shape-файла (при условии, что до этого был открыт хотя бы один растровый файл) на панели инструментов станет активной кнопка наложений (оверлеев) , предо¬ста¬вля¬ю¬щая доступ к списку открытых shape-файлов (Рис. 2.15). Ставя или убирая галочку в соответствующей строек можно управлять видимостью накладываемых слоёв. | |||
Рис. 2.15. Пример наложения векторных данных на многозональный снимок | |||
На этом рисунке показан пример, где маленькими чёрными квадратиками отмечены пробные площадки в лесном массиве, окружающем апатитово-нефелиновую обогатительную фабрику № 2 в г. Апатиты. | |||
Для открытия shape-файла требуется в диалоговом окне Открыть в списке фильтров Тип файлов выбрать вариант Shape (*.shp) (Рис. 2.16). | |||
Рис. 2.16. Установка фильтра для открытия shape-файла | |||
Имеется возможность произвольно менять цвет и толщину линий, с помощью которых прорисовываются векторные данные. Для этого надо щёлкнуть по кнопке оверлеев при нажатой клавише Shift. В выпавшем списке открытых shape-файлов выбрать строку с нужным слоем. Появится диалоговое окно Set Overlay Parameters (Установка параметров оверлея) (Рис. 2.17). | |||
Рис. 2.17. Установка параметров линий shape-файла | |||
Для удаления наложенных векторных данных можно использовать команду Edit Clear Overlays (Убрать оверлеи). | |||
1.9. РАЗБИЕНИЕ ИЗОБРАЖЕНИЯ НА КЛАССЫ ПО ОБРАЗЦАМ | |||
Программа MultiSpec позволяет выполнять автоматическое выделение на многозональном изображении участков, сходных с эталонными по спектральным характеристикам [7]. Такая операция называется классификация с обучением. Когда алгоритм классификации разобьёт всё изображение на ограниченное число классов (соответствующих типам наземных объектов), то это будет большим подспорьем при дешифрировании и интерпретации космоснимка. Предварительно требуется произвести обучение алгоритма на эталонных образцах – участках изображения с характерными (и известными интерпретатору) спектральными портретами. | |||
1.9.1. Выбор эталонных участков | |||
Подготовка набора эталонных участков производится следующим образом. | |||
1. Открыть ранее собранный файл многоканального изображения. Если просто открыть несколько файлов отдельных каналов одной сцены, то необходимая команда меню будет недоступной. | |||
2. В меню Processor выбрать команду Statistics. В появившемся диалоговом окне Set Project Options (Рис. 2.18) щёлкнуть по кнопке OK. Значения по умолчанию, применённые в таком случае будут вполне уместны. | |||
Рис. 2.18. Окно установки параметров проекта | |||
3. Новое окно проекта (с надписью Project в заголовке) появится в правой части окна. | |||
4. Для выбора эталонных участков по каждому классу потребуется просто нарисовать прямоугольную рамку выделения поверх космоснимка. Или (при выборе опции Polygon Enter) нарисовать многоугольник, щёлкая мышкой в вершинах ломаной линии. | |||
Желательно так выделить эталонный участок, чтобы пространство внутри него было максимально однородно по спектральным характеристикам. Также важно отступать на пару пикселей внутрь от границы участка. Это позволит избежать краевых эффектов. | |||
5. В случае ошибки просто повторить процесс рисования границы эталонного участка сначала. | |||
6. Щёлкнуть по кнопке "Add To List…" (Рис. 2.19). | |||
Рис. 2.19. Пример выделения эталонного участка озера | |||
7. Появится диалоговое окно Define Class and/or Field Description (Определить описание эталонного класса и/или участка) (Рис. 2.20). В этом окне можно определить название класса и дать наименование участка эталона. | |||
Рис. 2.20. Диалоговое окно описания эталонного участка | |||
8. Завершить с помощью кнопки OK. | |||
Таким образом надо создать набор обучающих эталонов для всех существенных классов объектов, встречающихся на изображении (Рис. 2.21). Для одного класса можно создавать несколько эталонов. | |||
Рис. 2.21. Пример результата создания набора обучающих эталонов | |||
Окно проекта (Project) может находиться в четырёх различных режимах: | |||
• Выделение обучающих эталонов (кнопка >Select). | |||
• Список классов (кнопка >Classes). | |||
• Список эталонов (кнопка >Fields). | |||
• Список координат (кнопка >Coord.). | |||
Пользователь может, переключившись в соответствующий режим, удалить ненужный класс. Для этого надо выделить его в списке классов и выполнить команду Edit Cut Class (Вырезать класс). Аналогичным образом можно удалить и эталон. | |||
Имеется возможность задавать границы эталонных участков не с помощью прямоугольника, как по умолчанию, а в виде полигона. Для этого надо включить опцию Polygon Enter, когда окно проекта находится в режиме выделения. Затем надо щёлкать мышью по изображению, задавая таким образом вершины полигона. Двойной щелчок завершает ввод полигональной границы. | |||
1.9.2. Выполнение классификации | |||
Когда набор эталонных участков разных классов подготовлен, можно приступать к выполнению самой классификации. Для этого надо выполнить следующие действия: | |||
1. Выполнить команду Processor Classify… (Классификация). Появится диалоговое окно Set Classification Specifications (Установка параметров классификации) (Рис. 2.22). | |||
Рис. 2.22. Диалоговое окно установки параметров классификации | |||
2. На первом этапе отключить параметр Image Selection (Выделение изображения). Это необходимо для того, чтобы проводить классификацию только по обучающим эталонам и получить начальную оценку качества определения классов. | |||
3. Оставить все остальные параметры со значениями по умолчанию. | |||
4. Щёлкнуть по кнопке OK в диалоговом окне параметров классификации, а затем по кнопке OK в окне Update Project Statistics (Обновить статистику проекта) (Рис. 2.23) для запуска процесса классификации. | |||
Рис. 2.23. Предупреждение об обновлении статистики проекта | |||
5. После окончания вычислений можно посмотреть результаты классификации (Табл. 2.2). Для этого выполнить команду Window Text Output (Текстовое окно). | |||
Табл. 2.2. Пример таблицы результатов классификации | |||
Project | |||
Class | |||
Name Reference | |||
Class | |||
Number Accuracy+ | |||
(%) Number of Samples in Class | |||
Number | |||
Samples 1 2 3 4 5 6 7 8 | |||
Озеро Хвостохранилище Облако Тень Тайга Луг Грунт Скалы | |||
Озеро 1 100.0 1064 1064 0 0 0 0 0 0 0 | |||
Хвостохранилище 2 100.0 720 0 720 0 0 0 0 0 0 | |||
Облако 3 100.0 60 0 0 60 0 0 0 0 0 | |||
Тень 4 100.0 330 0 0 0 330 0 0 0 0 | |||
Тайга 5 100.0 156 0 0 0 0 156 0 0 0 | |||
Луг 6 100.0 182 0 0 0 0 0 182 0 0 | |||
Грунт 7 100.0 36 0 0 0 0 0 0 36 0 | |||
Скалы 8 100.0 132 0 0 0 0 0 0 0 1321 | |||
TOTAL 2680 1064 720 60 330 156 182 36 132 | |||
Reliability Accuracy (%)* 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0 | |||
Таблица, напечатанная в текстовом окне, содержит информацию о том, насколько качественно были классифицированы пиксели каждого класса. Параметр Accuracy (Точность) должен быть около 100% для всех обучающих эталонов. Если же этот параметр имеет низкое значение, меньше 50% у какого-то класса, то это значит, что для него надо переназначить эталонные участки. | |||
Если результаты классификации являются удовлетворительными, то можно приступать к классификации всего изображения. | |||
6. Снова выполнить команду Processor Classify… Появится диалоговое окно Set Classification Specifications. | |||
7. Отключить параметр Training (resubstitution) (Обучение (повторная подстановка)). | |||
8. Включить параметр Image Selection. | |||
9. Убедиться, что будет классифицироваться всё изображение, для чего нажать на квадратную кнопку в группе Area to Classify (Область классификации), если она активна. | |||
10. Включить параметр Disk file в группе Write classification results to:, что позволит сохранить результаты для дальнейшего использования. По умолчанию сохранение будет производиться в формате ERDAS.GIS. | |||
11. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла с результатами классификации. | |||
12. Сохранить проект с помощью команды File Save Project (Сохранить проект). Это приведёт к сохранению набора обучающих эталонных участков и соответствующих им классов для последующего анализа. | |||
1.9.3. Просмотр классифицированного растра | |||
Результат классификации в виде типизированного растра был сохранён в файле «*.gis». Если открыть его, то результат будет подобен приведённому ниже. | |||
Рис. 2.24. Предупреждение об обновлении статистики проекта | |||
Выполнить команду Project Add as Associated Image (Добавить как связанное изображение). Это приведёт к наложению границ обучающих эталонов с подписями поверх изображения. | |||
Есть возможность изменить цвет этих линий и надписей на чёрный, если выполнить команду Processor Statistics… и выбрать пункт Black (Чёрный) в списке Color (Цвет) группы Outline selected areas: (Линии выделенных областей). | |||
В помощь для более удобного анализа графических результатов классификации пользователю предоставляется две полезные возможности: | |||
• Если при нажатой клавише Shift щёлкнуть левой кнопкой мыши по цветному прямоугольнику интересуемого класса в легенде, то это приведёт к окрашиванию в белый цвет всех пикселей этого класса на изображении. | |||
• Если одновременно нажать клавиши Shift и Ctrl и щёлкнуть левой кнопкой мыши по цветному прямоугольнику класса, то все классы, кроме этого, поменяют свой цвет на белый. | |||
• Двойной щелчок мыши по цветному прямоугольнику класса позволит изменить цвет, ассоциированный с ним по умолчанию на более подходящий по смыслу. | |||
1.10. ЭКСПОРТ ГРАФИЧЕСКОГО ИЗОБРАЖЕНИЯ В РАСТРОВЫЙ ФАЙЛ | |||
Многозональные изображения, создаваемые в MultiSpec, особенно, если они содержат более 3 каналов, не могут быть открыты для просмотра многими графическими программами или импортированы текстовыми редакторами. Поэтому для передачи полученных изображений в другие программы можно, помимо копирования через буфер обмена, использовать возможность сохранения в простой растровый файл [1]. | |||
При экспорте информация сохраняется «как есть», т.е. просто делается мгновенный снимок изображения при текущих значениях параметров визуализации космоснимка. Дальнейшее изменение внешнего вида, такие как смена комбинации каналов, использование разных кривых яркости и т.п. становится недоступным. | |||
Одновременно с созданием растрового файла (*.tif), если это многозональный, а не тематический растр, создаётся одноимённый kml-файл, позволяющий сохранить также пространственную привязку изображения. | |||
Для экспорта в растровый файл необходимо выполнить команду File Save Image To GeoTIFF As… (Сохранить изображение в формате GeoTIFF) и в появившемся стандартном диалоговом окне задать имя экспортируемого фала. | |||
Если предварительно выделить часть изображения, то команда для экспорта преобразится в Save Selection Image To GeoTIFF As… (Сохранить выделение в формате GeoTIFF). | |||
ЛИТЕРАТУРА | |||
1. An introduction and reference for MultiSpec. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/MultiSpec_Intro_9_11.pdf. – 2011. – 189 pp. (дата обращения: 31.08.2015). | |||
2. Tutorial: Combining separate image files into a single multispectral image file. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Tutorial_5.pdf. – 2013. – 10 pp. (дата обращения: 31.08.2015). | |||
3. MultiSpec exercise: Image enhancement. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_2.pdf. – 2009. – 5 pp. (дата обращения: 31.08.2015). | |||
4. MultiSpec exercise: Selecting areas and the coordinate view. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_7.pdf. – 2009. – 2 pp. (дата обращения: 31.08.2015). | |||
5. MultiSpec exercise: Creating vegetation indices images. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_8.pdf. – 2010. – 7 pp. (дата обращения: 31.08.2015). | |||
6. MultiSpec exercise: Overlay shape files on image window. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_6.pdf. – 2009. – 2 pp. (дата обращения: 31.08.2015). | |||
7. MultiSpec exercise: Supervised classification. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_4.pdf. – 2009. – 7 pp. (дата обращения: 31.08.2015). |
Версия от 19:05, 31 августа 2015
В статье описываются основные функции бесплатной программы MultiSpec, используемые для обработки, интерпретации и анализа многозональных космоснимков.
Программа MultiSpec была разработана и поддерживается Дэвидом Лэндгрибом (David Landgrebe) и Ларри Бьелем (Larry Biehl) из университета Пердью (США). MultiSpec является результатом продолжающихся многолетних исследований, направленных на получение мощного и надёжного инструмента для анализа многозональных и гиперспектральных данных. MultiSpec вместе с документацией доступна для свободного скачивания с официального сайта программы . К основным функциям программы относятся [1]:
- Импорт данных, как из двоичных, так и текстовых данных.
- Визуализация многозональных изображений с помощью множества черно-белых и цветных режимов; показ создаваемых типизированных растров. Наложение векторных данных поверх изображений.
- Расчёт гистограмм по многозональным изображениям.
- Переформатирование данных, включающее изменение каналов, объединение файлов, создание мозаик, изменение геометрии и др.
- Создание новых каналов из существующих с помощью извлечения информации из многозонального изображения или в результате попиксельных алгебраических вычислений.
- Кластерный анализ (классификация без обучения) в отношении многозональных изображений и создание в результате типизированных растров.
- Классификация с обучением после подготовки набора эталонов для основных классов наземных объектов, видимых на многозональном изображении.
- Построение графика спектральных значений для выделенных участков многозонального изображения.
- Экспорт исходных данных и их трансформант в текстовый файл.
- Экспорт изображений в графический файл и буфер обмена.
Объединение файлов отдельных каналов в единое многозональное изображение
Многозональные спутниковые снимки поступают из хранилища ГС США в виде упакованных в архив монохромных файлов формата GeoTIFF (*.tif). Причём каждый файл представляет отдельный канал многозонального изображения. Такой подход к предоставлению материалов имеет как плюс – пользователь может ограничить объём обрабатываемых данных только нужными ему каналами, так и минус – для дальнейшей работы каналы должны быть предварительно объединены в многоканальный файл. Создание многозонального изображения выполняется в программе MultiSpec следующим образом [2]: 1. Выполнить команду из главного меню File Open Image… или щёлкнуть по кнопке Open на панели инструментов. Появится диалоговое окно Open Image File. 2. Выделить в окне те файлы каналов, которые должны быть включены в многозональное изображение (рис. 1).
Рис. 1. Выделение нескольких одноканальных файлов перед объединением
Можно обратить внимание на то, что при выделении более чем одного файла опция в нижнем левом углу переключается в Link selected files to new image window (Привязать выделенные файлы к новому окну изображения). 3. Щёлкнуть по кнопке Открыть. Появится диалоговое окно Set Display Specifications for (Параметры визуализации) (рис. 2). Это окно позволяет устанавливать большое число параметров, используемые при показе собираемого изображения.
Рис. 2. Окно установки параметров визуализации изображения
По умолчанию для показа представляется полный размер сцены (параметры Area to Display) и тип изображения выбран 3-Channel Color (3-канальный цветной). 4. Щёлкнуть по кнопке OK. В появившемся окне Histogram (Гистограмма) тоже щёлкнуть по кнопке OK. В результате в отдельном графическом окне появится многозональное изображение, оформленное в соответствии с установленными параметрами.
Управление визуализацией многозонального изображения
Открытое в программе MultiSpec с параметрами по умолчанию многозональное изображение можно визуализировать как-нибудь иначе с помощью группы опций, доступных с помощью команды Processor Display Image… главного меню. К этим опциям относятся Bits of color (глубина цвета), Stretch (кривая яркости), Min-max (пределы), Treat ‘0’ (представление нуля) и Display levels per channel (число уровней яркости) в диалоговом окне Set Display Specification (рис. 2.2). Назначение этих параметров следующее [3]:
- Параметры Bits of color и Display levels per channel по умолчанию имеют максимально возможные значения (24 и 256 соответственно). При желании их можно уменьшить, что может быть полезно для выполнения простейшей кластеризации изображения (разделения на объекты, принципиально различающиеся по своим свойствам: пашня, леса, вода, строения и т.п.) (рис. 3).
Рис. 3. Визуализация изображения с 256 уровнями яркости (миллионы цветов) и 2 уровнями яркости (8 цветов)
- Параметр Treat ‘0’ позволяет управлять интерпретацией нулевых значений в исходных файлах каналов материалов Landsat. Изначально нулевые значения обозначают фон, который не используется при дальнейших вычислениях и может быть показан чёрным (black) или белым (white) цветом. Однако при желании можно нулевые значения считать в качестве данных с нулевой яркостью (as data) (рис. 4).
Рис. 2.4. Изображение фона (нулевых данных) чёрным или белым цветом
- Параметры Stretch и Min-max обычно используются для улучшения визуализации отдельных частей изображения. Они контролируют соотношением между значениями яркости в исходных файлах каналов и уровнями яркости на экране. Есть три варианта параметра Stretch: Linear (линейная), Equal Area (равноплощадная) и Gaussian (гауссова) (рис. 5).
Рис. 5. Три варианта кривой яркости (линейная, равноплощадная и гауссова)
- В случае линейной кривой яркости интервалы серого цвета на экране равномерно распределены между значениями яркости в исходных данных.
- В то же время равноплощадная кривая автоматически подбирает такое соотношение между ними, чтобы каждый уровень яркости занимал на экране одинаковое число пикселей.
- Гауссова кривая яркости использует гауссово распределение уровней серого цвета на экране, что позволяет улучшить изображения слабосветящихся объектов.
- Опция Min-max позволяет установить предельные значения (и диапазон между ними) для исходных данных, которые будут использоваться при показе изображения на экране. Уровни яркости и кривая яркости будут ограничены этими пределами.
- Вариант Entire Range (весь диапазон) для параметра Min-maxes ставит в соответствии минимальное значение исходных данных (0) чёрному цвету, а максимальное (255) – белому. Однако иногда требуется не показывать значения, например, меньше 50 и больше 150.
- Вариант 2 Percent Tails Clipped (обрезать двухпроцентные края) позволяет сделать так, чтобы визуализировались значения исходных данных только в диапазоне 2-98% от их предельных значений. Это позволяет увеличить контрастность изображения за счёт исключения из показа редко встречающихся самых тёмных и самых светлых точек.
- Если выбрать вариант User Specified..., то появится диалоговое окно Set Specifications for Display Enhancement (Установить параметры кривой яркости). Это окно позволит задать свой собственный процент отсечения краёв (рис. 6), который может отличаться от 2%.
Рис. 6. Окно установки пределов кривой яркости
1.3. ВЫДЕЛЕНИЕ УЧАСТКОВ НА МНОГОЗОНАЛЬНОМ ИЗОБРАЖЕНИИ
Если в программе MultiSpec открыто изображение, то часть его можно выделить для того чтобы проводить дальнейшие преобразования только с этой частью, а не со всей картой. Существует 3 способа создания выделения [4]:
1. Выполнить команду Edit Selection Rectangle (Рамка выделения) главного меню программы. Появится соответствующее диалоговое окно. В этом окне надо ввести или отредактировать значения, задающие минимальное и максимальное значение координат в выбранной системе координат. Имеется 3 варианта систем координат:
• Lines-Columns (строки и столбцы) – локальная система координат данного изображения, основанная на нумерации его пикселей от верхнего левого угла вниз (строки) и направо (столбцы).
• Transverse Mercator-Meters (глобальная прямоугольная поперечная проекция Меркатора с метрами в качестве единиц измерения).
• Latitude-Longitude (Широта-Долгота) – глобальная сферическая система координат с единицами измерения – градусами (Рис. 2.7).
Рис. 2.7. Окно задания параметров рамки выделения
Кнопка Preview (просмотр) позволяет увидеть на карте сделанные изменения выделения. 2. С меньшей точностью, но быстрее рамку выделения можно нарисовать по карте прямо с помощью указателя мыши, ориентируясь на значения координат, которые можно видеть на панели Coordinate View (координаты) (Рис. 2.8).
Рис. 2.8. Панель координат и выделение с помощью указателя мыши
1.4. СОХРАНЕНИЕ В ЕДИНОМ ФАЙЛЕ СОБРАННОГО МНОГОКАНАЛЬНОГО ИЗОБРАЖЕНИЯ
Постоянно собирать многоканальное изображение из разрозненных файлов каналов при каждой необходимости использования – нерационально. Кроме того часто после сборки требуется продолжить работу не со всей 185-км сценой, а только с областью интереса, которая обычно значительно меньше. Существенно упростить работу может сохранение собранного многозонального изображения или его части в отдельный файл. Для этого потребуется следующее [1]:
1. Собрать многозональное изображение из тех каналов, которые планируется использовать в дальнейшем.
2. Выделить часть изображения, представляющую область интереса или удалить выделение, если должна быть сохранена сцена целиком.
3. Выполнить команду Processor Reformat Rectify Image… (Сохранить прямоугольный фрагмент). Появится диалоговое окно Set Image Rectification Parameters (Установить параметры прямоугольного фрагмента). Все параметры этого окна имеют по умолчанию значения, изменять которые обычно нет необходимости.
4. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. Надо ввести имя и расширение “.tif” для этого файла. По неизвестной причине MultiSpec сама в данном случае расширение к имени файла не добавляет, что впоследствии мешает ей же открывать такие файлы.
1.5. ПОЛУЧЕНИЕ СВОДНОЙ СТАТИСТИЧЕСКОЙ ИНФОРМАЦИИ ПО МНОГОЗОНАЛЬНОМУ ИЗОБРАЖЕНИЮ
Пункт меню Histogram image (гистограмма изображения) позволяет вычислить и вывести в текстовом окне сводную статистическую информацию по всем каналам многозонального изображения или его части. Причём это может быть как пространственная часть изображения, так и выборка из спектральных каналов [1].
1. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение.
2. Выполнить команду Processor Histogram image. Появится диалоговое окно Set Histogram Specifications (Установить параметры гистограммы) (Рис. 2.9).
Рис. 2.9. Окно установки параметров вычисления гистограммы изображения
3. Выбрать метод (Method) вычисления гистограммы – Compute new histogram (Вычислить новую гистограмму). 4. С помощью списка Channels (Каналы) можно выбрать каналы, используемые для вычислений. По умолчанию установлен вариант All (Все). 5. Для получения сводной информации надо включить опцию List histogram summary (вывести сводную гистограмму) и выключить List histogram (вывести всю гистограмму). 6. При нажатии на кнопку OK в текстовом окне появится информация похожая на приведённый ниже пример (Табл. 2.1). Табл. 2.1. Пример таблицы сводной статистической информации Channel Channel Description Data Range Mean Median Standard Deviation 1 0.45-0.52 um 39 to 50 43.1 43 1.6 2 0.53-0.61 um 15 to 21 17.8 18 1.2 3 0.63-0.69 um 12 to 21 15.6 15 2.0 4 0.78-0.90 um 24 to 59 41.7 41 6.4 5 1.55-1.75 um 22 to 54 38.2 39 6.0 6 2.09-2.35 um 10 to 23 15.4 15 2.4 7 10.4-12.5 um 103 to 111 106.7 107 1.8 1.6. СОЗДАНИЕ СХЕМ РАСПРЕДЕЛЕНИЯ NDVI В программе MultiSpec можно создавать трансформанты многозональных снимков – изображения, представляющие собой алгебраические комбинации значений исходных файлов каналов материалов Landsat. При этом чаще всего создают схемы распределения нормализованного дифференциального вегетационного индекса (NDVI). Назначение этого индекса – улучшить читаемость и представительность изображения путём подчёркивания особенностей распределения растительности и почвы. Создаётся такая схема следующим образом [5]: 1. В меню Processor (процессор) выбрать команду Reformat Change Image File Format… (изменить формат файла изображения). В результате появится соответствующее диалоговое окно (Рис. 2.10).
Рис. 2.10. Окно установки параметров изменения формата файла
2. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение. 3. Отметить опцию Transform Data… (преобразовать данные). В результате появится соответствующее диалоговое окно. 4. Выбрать опцию New Channel from General Algebraic Transformation (создать новый канал с помощью алгебраического преобразования). Окно изменит свой интерфейс (Рис. 2.11).
Рис. 2.11. Ввод формулы алгебраического преобразования каналов многозонального снимка
Ячейки для ввода значений в этом окне позволяют задать формулу алгебраической комбинации исходных каналов. С их помощью можно создать новый канал, например NDVI, если ввести после символа “C” номера каналов, соответствующие красному и ближнему ИК каналу. Если многозональный снимок содержит все каналы сцены Landsat, то это будут: канал 1 – синий, канал 2 – зелёный, канал 3 – красный, канал 4 – ближний инфракрасный и т.д. Тогда в ячейки окна надо вводить формулу следующим образом:
5. Далее щёлкнуть по кнопке OK, что применит введённую формулу преобразования, закроет окно Set Reformat Transform Parameters и вернёт фокус в окно Image File Format Change Options. 6. В окне Image File Format Change Options для параметра Data value type (тип значений данных) задать вариант 32-bit Real (вещественный). Программа может и сама установить этот тип, если определит, что результат вычислений будет не целочисленным. 7. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. Надо ввести имя и расширение “.tif” для этого файла. 8. После этого можно просмотреть получившейся результат, если при открытии файла задать тип изображения 1-channel thematic (1-канальный тематический) (Рис. 2.12).
Рис. 2.12. Схема распределения NDVI
На данном изображении сине-зелёные оттенки соответствуют воде, красно-оранжевые – растительности разной пышности. Остальное – почво-грунты. 1.7. ЭКСПОРТ ИЗОБРАЖЕНИЯ В ТЕКСТОВЫЙ ФАЙЛ Нередко изобразительных свойств программы MultiSpec бывает недостаточно. Возникает желание передать информацию в другую, более подходящую, программу, например, Golden Software Surfer. Для этого надо экспортировать открытое в MultiSpec изображение в текстовый файл, который впоследствии может быть открыт в Surfer. Экспорт в текстовый файл может быть выполнен так [1]: 1. Выделить часть изображения, по которому будут производиться вычисления. Или удалить выделение с помощью клавиши Delete, чтобы обработать всё изображение. 2. Выполнить команду Processor List Data… Появится диалоговое окно Set List Data Specifications (Установить параметры экспорта данных) (Рис. 2.13).
Рис. 2.13. Параметры экспорта данных
3. Включить переключатель Include latitude and longitude values (включить значения широты и долготы). 4. В группе Output results to (Вывести результаты в) убрать отметку опции Text window (Текстовое окно) и включить отметку опции Disk file (Дисковый файл). 5. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла. Этот файл имеет простой текстовый формат. В нём сформирована таблица, содержащая как координаты каждой точки изображения, так и значения всех каналов в ней.
Рис. 2.14. Пример построения распределения NDVI средствами Surfer
На данном изображении благодаря богатым возможностям визуализации двумерных данных у Surfer применены более естественные цвета (зелёный – пышная растительность, жёлто-зелёный – редкая растительность, коричневый – грунт, синий – вода) с точно подобранными рубежами перехода между цветами. Рассмотрение технологии использования Surfer для подобных построений находится за пределами данного пособия. 1.8. НАЛОЖЕНИЕ ВЕКТОРНЫХ ДАННЫХ В ОКНЕ ИЗОБРАЖЕНИЯ Векторные данные – это подготовленные в других программах наборы графических примитивов (полигонов, линий, точек), в совокупности представляющие собой чаще всего элементы топографической карты (реки, дороги, населённые пункты и т.п.). Векторные данные, наложенные на многозональное или классифицированное изображение, помогают легче ориентироваться по нему и быстро находить области интереса. MultiSpec имеет возможность открывать векторные данные из shape-файлов (*.shp) [6]. После открытия хотя бы одного shape-файла (при условии, что до этого был открыт хотя бы один растровый файл) на панели инструментов станет активной кнопка наложений (оверлеев) , предо¬ста¬вля¬ю¬щая доступ к списку открытых shape-файлов (Рис. 2.15). Ставя или убирая галочку в соответствующей строек можно управлять видимостью накладываемых слоёв.
Рис. 2.15. Пример наложения векторных данных на многозональный снимок
На этом рисунке показан пример, где маленькими чёрными квадратиками отмечены пробные площадки в лесном массиве, окружающем апатитово-нефелиновую обогатительную фабрику № 2 в г. Апатиты. Для открытия shape-файла требуется в диалоговом окне Открыть в списке фильтров Тип файлов выбрать вариант Shape (*.shp) (Рис. 2.16).
Рис. 2.16. Установка фильтра для открытия shape-файла
Имеется возможность произвольно менять цвет и толщину линий, с помощью которых прорисовываются векторные данные. Для этого надо щёлкнуть по кнопке оверлеев при нажатой клавише Shift. В выпавшем списке открытых shape-файлов выбрать строку с нужным слоем. Появится диалоговое окно Set Overlay Parameters (Установка параметров оверлея) (Рис. 2.17).
Рис. 2.17. Установка параметров линий shape-файла
Для удаления наложенных векторных данных можно использовать команду Edit Clear Overlays (Убрать оверлеи). 1.9. РАЗБИЕНИЕ ИЗОБРАЖЕНИЯ НА КЛАССЫ ПО ОБРАЗЦАМ Программа MultiSpec позволяет выполнять автоматическое выделение на многозональном изображении участков, сходных с эталонными по спектральным характеристикам [7]. Такая операция называется классификация с обучением. Когда алгоритм классификации разобьёт всё изображение на ограниченное число классов (соответствующих типам наземных объектов), то это будет большим подспорьем при дешифрировании и интерпретации космоснимка. Предварительно требуется произвести обучение алгоритма на эталонных образцах – участках изображения с характерными (и известными интерпретатору) спектральными портретами. 1.9.1. Выбор эталонных участков Подготовка набора эталонных участков производится следующим образом. 1. Открыть ранее собранный файл многоканального изображения. Если просто открыть несколько файлов отдельных каналов одной сцены, то необходимая команда меню будет недоступной. 2. В меню Processor выбрать команду Statistics. В появившемся диалоговом окне Set Project Options (Рис. 2.18) щёлкнуть по кнопке OK. Значения по умолчанию, применённые в таком случае будут вполне уместны.
Рис. 2.18. Окно установки параметров проекта
3. Новое окно проекта (с надписью Project в заголовке) появится в правой части окна. 4. Для выбора эталонных участков по каждому классу потребуется просто нарисовать прямоугольную рамку выделения поверх космоснимка. Или (при выборе опции Polygon Enter) нарисовать многоугольник, щёлкая мышкой в вершинах ломаной линии. Желательно так выделить эталонный участок, чтобы пространство внутри него было максимально однородно по спектральным характеристикам. Также важно отступать на пару пикселей внутрь от границы участка. Это позволит избежать краевых эффектов. 5. В случае ошибки просто повторить процесс рисования границы эталонного участка сначала. 6. Щёлкнуть по кнопке "Add To List…" (Рис. 2.19).
Рис. 2.19. Пример выделения эталонного участка озера
7. Появится диалоговое окно Define Class and/or Field Description (Определить описание эталонного класса и/или участка) (Рис. 2.20). В этом окне можно определить название класса и дать наименование участка эталона.
Рис. 2.20. Диалоговое окно описания эталонного участка
8. Завершить с помощью кнопки OK. Таким образом надо создать набор обучающих эталонов для всех существенных классов объектов, встречающихся на изображении (Рис. 2.21). Для одного класса можно создавать несколько эталонов.
Рис. 2.21. Пример результата создания набора обучающих эталонов
Окно проекта (Project) может находиться в четырёх различных режимах: • Выделение обучающих эталонов (кнопка >Select). • Список классов (кнопка >Classes). • Список эталонов (кнопка >Fields). • Список координат (кнопка >Coord.). Пользователь может, переключившись в соответствующий режим, удалить ненужный класс. Для этого надо выделить его в списке классов и выполнить команду Edit Cut Class (Вырезать класс). Аналогичным образом можно удалить и эталон. Имеется возможность задавать границы эталонных участков не с помощью прямоугольника, как по умолчанию, а в виде полигона. Для этого надо включить опцию Polygon Enter, когда окно проекта находится в режиме выделения. Затем надо щёлкать мышью по изображению, задавая таким образом вершины полигона. Двойной щелчок завершает ввод полигональной границы. 1.9.2. Выполнение классификации Когда набор эталонных участков разных классов подготовлен, можно приступать к выполнению самой классификации. Для этого надо выполнить следующие действия: 1. Выполнить команду Processor Classify… (Классификация). Появится диалоговое окно Set Classification Specifications (Установка параметров классификации) (Рис. 2.22).
Рис. 2.22. Диалоговое окно установки параметров классификации
2. На первом этапе отключить параметр Image Selection (Выделение изображения). Это необходимо для того, чтобы проводить классификацию только по обучающим эталонам и получить начальную оценку качества определения классов. 3. Оставить все остальные параметры со значениями по умолчанию. 4. Щёлкнуть по кнопке OK в диалоговом окне параметров классификации, а затем по кнопке OK в окне Update Project Statistics (Обновить статистику проекта) (Рис. 2.23) для запуска процесса классификации.
Рис. 2.23. Предупреждение об обновлении статистики проекта
5. После окончания вычислений можно посмотреть результаты классификации (Табл. 2.2). Для этого выполнить команду Window Text Output (Текстовое окно). Табл. 2.2. Пример таблицы результатов классификации Project Class Name Reference Class Number Accuracy+ (%) Number of Samples in Class Number Samples 1 2 3 4 5 6 7 8 Озеро Хвостохранилище Облако Тень Тайга Луг Грунт Скалы Озеро 1 100.0 1064 1064 0 0 0 0 0 0 0 Хвостохранилище 2 100.0 720 0 720 0 0 0 0 0 0 Облако 3 100.0 60 0 0 60 0 0 0 0 0 Тень 4 100.0 330 0 0 0 330 0 0 0 0 Тайга 5 100.0 156 0 0 0 0 156 0 0 0 Луг 6 100.0 182 0 0 0 0 0 182 0 0 Грунт 7 100.0 36 0 0 0 0 0 0 36 0 Скалы 8 100.0 132 0 0 0 0 0 0 0 1321
TOTAL 2680 1064 720 60 330 156 182 36 132 Reliability Accuracy (%)* 100.0 100.0 100.0 100.0 100.0 100.0 100.0 100.0
Таблица, напечатанная в текстовом окне, содержит информацию о том, насколько качественно были классифицированы пиксели каждого класса. Параметр Accuracy (Точность) должен быть около 100% для всех обучающих эталонов. Если же этот параметр имеет низкое значение, меньше 50% у какого-то класса, то это значит, что для него надо переназначить эталонные участки. Если результаты классификации являются удовлетворительными, то можно приступать к классификации всего изображения. 6. Снова выполнить команду Processor Classify… Появится диалоговое окно Set Classification Specifications. 7. Отключить параметр Training (resubstitution) (Обучение (повторная подстановка)). 8. Включить параметр Image Selection. 9. Убедиться, что будет классифицироваться всё изображение, для чего нажать на квадратную кнопку в группе Area to Classify (Область классификации), если она активна. 10. Включить параметр Disk file в группе Write classification results to:, что позволит сохранить результаты для дальнейшего использования. По умолчанию сохранение будет производиться в формате ERDAS.GIS. 11. Щёлкнуть по кнопке OK. Это закроет окно и вызовет стандартный диалог для сохранения создаваемого файла с результатами классификации. 12. Сохранить проект с помощью команды File Save Project (Сохранить проект). Это приведёт к сохранению набора обучающих эталонных участков и соответствующих им классов для последующего анализа. 1.9.3. Просмотр классифицированного растра Результат классификации в виде типизированного растра был сохранён в файле «*.gis». Если открыть его, то результат будет подобен приведённому ниже.
Рис. 2.24. Предупреждение об обновлении статистики проекта
Выполнить команду Project Add as Associated Image (Добавить как связанное изображение). Это приведёт к наложению границ обучающих эталонов с подписями поверх изображения. Есть возможность изменить цвет этих линий и надписей на чёрный, если выполнить команду Processor Statistics… и выбрать пункт Black (Чёрный) в списке Color (Цвет) группы Outline selected areas: (Линии выделенных областей). В помощь для более удобного анализа графических результатов классификации пользователю предоставляется две полезные возможности: • Если при нажатой клавише Shift щёлкнуть левой кнопкой мыши по цветному прямоугольнику интересуемого класса в легенде, то это приведёт к окрашиванию в белый цвет всех пикселей этого класса на изображении. • Если одновременно нажать клавиши Shift и Ctrl и щёлкнуть левой кнопкой мыши по цветному прямоугольнику класса, то все классы, кроме этого, поменяют свой цвет на белый. • Двойной щелчок мыши по цветному прямоугольнику класса позволит изменить цвет, ассоциированный с ним по умолчанию на более подходящий по смыслу. 1.10. ЭКСПОРТ ГРАФИЧЕСКОГО ИЗОБРАЖЕНИЯ В РАСТРОВЫЙ ФАЙЛ Многозональные изображения, создаваемые в MultiSpec, особенно, если они содержат более 3 каналов, не могут быть открыты для просмотра многими графическими программами или импортированы текстовыми редакторами. Поэтому для передачи полученных изображений в другие программы можно, помимо копирования через буфер обмена, использовать возможность сохранения в простой растровый файл [1]. При экспорте информация сохраняется «как есть», т.е. просто делается мгновенный снимок изображения при текущих значениях параметров визуализации космоснимка. Дальнейшее изменение внешнего вида, такие как смена комбинации каналов, использование разных кривых яркости и т.п. становится недоступным. Одновременно с созданием растрового файла (*.tif), если это многозональный, а не тематический растр, создаётся одноимённый kml-файл, позволяющий сохранить также пространственную привязку изображения. Для экспорта в растровый файл необходимо выполнить команду File Save Image To GeoTIFF As… (Сохранить изображение в формате GeoTIFF) и в появившемся стандартном диалоговом окне задать имя экспортируемого фала. Если предварительно выделить часть изображения, то команда для экспорта преобразится в Save Selection Image To GeoTIFF As… (Сохранить выделение в формате GeoTIFF). ЛИТЕРАТУРА 1. An introduction and reference for MultiSpec. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/MultiSpec_Intro_9_11.pdf. – 2011. – 189 pp. (дата обращения: 31.08.2015). 2. Tutorial: Combining separate image files into a single multispectral image file. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Tutorial_5.pdf. – 2013. – 10 pp. (дата обращения: 31.08.2015). 3. MultiSpec exercise: Image enhancement. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_2.pdf. – 2009. – 5 pp. (дата обращения: 31.08.2015). 4. MultiSpec exercise: Selecting areas and the coordinate view. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_7.pdf. – 2009. – 2 pp. (дата обращения: 31.08.2015). 5. MultiSpec exercise: Creating vegetation indices images. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_8.pdf. – 2010. – 7 pp. (дата обращения: 31.08.2015). 6. MultiSpec exercise: Overlay shape files on image window. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_6.pdf. – 2009. – 2 pp. (дата обращения: 31.08.2015). 7. MultiSpec exercise: Supervised classification. [Электронный ресурс] // PURDUE.EDU: сервер университета Пердью. URL: https://engineering.purdue.edu/~biehl/MultiSpec/tutorials/MultiSpec_Exercise_4.pdf. – 2009. – 7 pp. (дата обращения: 31.08.2015).