Использование методов и средств аэрокосмического мониторинга, ГИС-технологий, дискретных ортогональных преобразований становится все более популярным в различных отраслях науки и техники [3–6, 10]. Одна из задач аэрокосмического мониторинга – детектирование областей поверхности земли с требуемыми свойствами.
ЭВМ современного поколения не достаточно развиты для того, чтобы сравниться с человеком-оператором. Это обуславливает наличие ошибок детектирования любых алгоритмов с любыми параметрами детектирования. Таким образом, для текущего поколения компьютерных технологий необходимо введение человека как основного звена космической обработки, в частности, в задаче детектирования областей. В этом случае любая обработка изображений не может быть полностью автоматической, т.е. имеет характер автоматизированной.
Постановка задачи
Число алгоритмов детектирования растет со временем. Большинство, если не все, из них имеют ошибки детектирования E, которые складываются из ошибок необнаружения истинных объектов (E1) и ошибок обнаружения ложных объектов (E2). Различным параметрам детектирования p соответствуют различные значения этих ошибок, т.е. регулирование значений этих параметров можно оптимизировать E. В этом случае область детектирования получается максимально полной (минимум E1) и четкой, т.е. с минимальным зашумлением (минимум E2). В исключительных случаях может быть E = 0. В связи с этим требуется разрабатывать специфические алгоритмы детектирования, применяемые для каждых конкретных типов объектов. При этом оптимальные параметры детектирования определяются ручным способом – варьированием p. Более того, имеет место многоэтапность обработки изображений, т.е. на каждом этапе алгоритм обрабатывает «полуфабрикат» изображения по своим информационным признакам. Рост числа этапов обработки ведет к усложнению алгоритма детектирования и к росту ошибок детектирования.
Описание метода
Предлагается метод, который условно назовем «Динамическая обработка изображений». В соответствии с этим методом человек вводится на все этапы обработки, в частности:
1) отслеживает и приостанавливает процесс обработки в любой момент времени;
2) регулированием значений параметров обработки просматривает и контролирует результат обработки;
3) обследует различные варианты результатов обработки в зависимости от тех или иных значений параметров обработки.
Для удобства и полноты управления «штурвалом самолета» детектирования предлагается идея использования принципа «кино» в обработке изображений различных видов. Данный принцип состоит в автоматическом перелистывании изображений обработки различных видов, в том числе космических снимков, каждому из которых соответствуют свои значения параметров обработки, которые изменяются во времени по некоторой временной функции. Кадры этого фильма сменяются один на другой при изменении параметров обработки за некоторый малый интервал времени.
Данная идея может быть программно реализована, например, в виде приложения мультимедиа с приведением различных внутренних механизмов обработки к «законам» мультимедиа. Например, вводится терминология, которая включает такие понятия, как видеообработка, фильм, видеоэкран, трек, кадр, проигрывание, остановка, приостановка, режим, жанр, скорость, настройка воспроизведения.
Среди возможностей метода можно выделить:
1) демонстрацию процесса динамической обработки;
2) составление карт результатов обработки, в частности детектирования [7, 9];
3) прослеживание изменений результатов обработки в зависимости от значений тех или иных параметров [1, 2, 8];
4) теоретические обследования территорий и оптимизацию параметров обработки.
Общее представление математического аппарата методики
Математическую модель фильма можно представить в виде следующей формулы:
F = U[M, p]; p = f(t), (*)
где F – фильм, который состоит из последовательности изображений обработки (видеоизображение); t – текущий момент времени; M – исходное изображение; U[M, p] – текущее изображение обработки; p – параметр (параметры) обработки над исходным изображением; U – оператор, обрабатывающий исходное изображение M по заданному алгоритму, воздействуя на M параметром обработки p; f(t) – функция параметра обработки p от времени t (параметрическая функция).
Математическая модель кадра фильма:
F0 = U[M, p0]; p0 = f(t0),
где F0 – кадр фильма, представляющий собой значение функции U0 = [M, p0] в определенной точке p = p0 для постоянного исходного изображения M; t0 – момент времени, в который отображается данный кадр (например, на экране дисплея).
Среди множества частных случаев (*) можно выделить видео пороговой фильтрации, кластеризации, дилатации, наполнения областей, выделения объектов и т.д.
Например, модель видео фильтрации полихроматического изображения:
.
где x и y – координаты абсциссы и ординаты пикселя изображений M и U; U(x, y) – КСЯ данного пикселя на изображении U; M(x, y) – КСЯ данного пикселя на исходном изображении M; U – бинарное изображение с выделением областей детектирования, которой соответствует интервал КСЯ [m1, m2].
Рис. 1. Общая блок-схема методики
Программная реализация методики
Общая блок-схема программы представлена на рис. 1. Одним из вариантов программной реализации методики является программа-мультимедиа с использованием соответствующих мультимедийных настроек, таких как Play, Stop, Pause и т.д. В программе предусматривается возможность подключения различных программных блоков, отрабатывающих те или иные алгоритмы обработки относительно регулируемых параметров обработки. Например, для программного блока «пороговая фильтрация» параметрами обработки являются средние значения коэффициентов спектральной яркости (КСЯ) и их отклонения на различных каналах исходных мультиспектральных изображений.
Программная реализация может быть в различных программных формах и на различных уровнях сложности и завершенности:
1) программного блока;
2) многооконного виртуального приложения;
3) завершенного программного комплекса.
а 
б
Рис. 2. Область наблюдения (а), программа Google Планета Земля; исходное изображение, Landsat 5 TM (июнь, 2010 г.)
Программная «философия» реализации методики, приуроченной к принципу мультимедиа, включает множество различных компонентов, в частности панель управления, строку состояния, программное тело, панель, экран, пульт навигации, пульт рисования и мн. др.
Результаты работы программы на примере космического мониторинга объектов захоронения отходов
На примере объектов захоронения отходов (ОЗО) с использованием программных блоков «пороговая фильтрация» и «наполнение областей» покажем работу методики.
Как известно [1, 8, 9], ОЗО являются структурно сложными объектами, и для их выделения недостаточно использовать только пороговую фильтрацию. В связи с этим, сперва для исходного изображения (рис. 3) просматривался фильм пороговой фильтрации на его различных каналах с регулированием пороговых значений КСЯ на других каналах. Задача состояла в приведении к выделениям эталонных ОЗО в область детектирования. При проведении эксперимента было обнаружено, что оптимизировать параметры детектирования удобно поканально, с каждым последующим каналом приближаясь к эталону. В результате была «на глаз» получена одна из многочисленных оптимальных точек параметров детектирования:
Соответственно получено одно из оптимальных изображений детектирования – рис. 4. Данное изображение приближено к эталонной области с учетом минимального значения ошибок детектирования E, но очевидно, что она сильно отличается от эталонной области (см. I–IV на рис. 2): сами ОЗО обнаруживаются на изображениях, но не выделяются. В качестве областей I-IV представлены основные полигоны ТБО окрестностей: I – Кучино, II – Савино, III – Торбеево, IV – Некрасовка.
Это предел возможностей пороговой фильтрации для ОЗО, поэтому целесообразно программный блок «пороговая фильтрация» усилить другим программным блоком, подключаемым к программе-мультимедиа. Одним из вариантов такого усиления является программный блок «наполнение областей» [7]. Кадры данного фильма – области детектирования, полученные пороговой фильтрацией с более широкой полосой пропускания для выделения самих ОЗО. При этом выделяются только те связные компоненты, которые включают области на рис. 3. Если при изменении параметров обработки в фильме площадь связных компонент меняется скачком, то предыдущее значение параметра дает уже выделение соответствующей включенной области (рис. 3).
Рис. 3. Обнаружение ОЗО в программе-мультимедиа с помощью программного блока пороговой фильтрации
Рис. 4. Выделение ОЗО в программе-мультимедиа с помощью программного блока наполнения областей для ОЗО: а – I; б – II; в – III; г – IV
На рис. 4 показаны примеры предельного наполнения областей I–IV (см. рис. 3). A – обнаружение ОЗО, Б – выделение ОЗО.
Выводы
В соответствии с предлагаемой идеей была разработана полноценная программа-мультимедиа обработки космических изображений с возможностью подключения различных программных блоков для видеообработки. Данная программа была протестирована в различных задачах космического мониторинга, в частности ОЗО, деградации почвы, обнаружения и выделения антропогенных объектов, обследования импактных областей.
Рецензенты:
Заалишвили В.Б., д.ф.-м.н., профессор, директор ФГБУН «Центр геофизических исследований», г. Владикавказ;
Марчук В.И., д.т.н., профессор, заведующий кафедрой «Радиоэлектронные системы», Институт сферы обслуживания и предпринимательства (филиал), ДГТУ, г. Шахты.
Библиографическая ссылка
Рихтер А.А., Казарян М.Л., Шахраманьян М.А. ОБРАБОТКА КОСМИЧЕСКИХ ИЗОБРАЖЕНИЙ С ИСПОЛЬЗОВАНИЕМ СРЕДСТВ МУЛЬТИМЕДИА И ЕЕ ПРИМЕНЕНИЕ В КОСМИЧЕСКОМ МОНИТОРИНГЕ ОБЪЕКТОВ ЗАХОРОНЕНИЯ ОТХОДОВ // Фундаментальные исследования. – 2015. – № 8-2. – С. 328-332;URL: https://fundamental-research.ru/ru/article/view?id=38895 (дата обращения: 19.04.2024).