Коллектив авторов - Океанография и морской лед
Свои недостатки и преимущества в зависимости от времени года и конкретной ледовой обстановки имеют и остальные из перечисленных алгоритмов.
Важно отметить, что микроволновые радиометры в настоящее время позволяют автоматически ежесуточно различать границу льда и общую сплоченность морского льда в течение сезонных периодов вне летнего таяния и осеннего замерзания. Однако автоматическая оценка сплоченности многолетнего льда все еще нестабильна и менее точна, чем оценки, основанные на экспертном анализе SAR или видимых/ИК изображений. В периоды летнего таяния и осеннего замерзания величина ошибок для автоматически получаемых продуктов увеличивается в два раза и более (Andersen, Tonboe et al.,2007). Для снижения погрешности необходима дополнительная экспертная проверка качества последовательных серий информационных продуктов.
Возраст льда. Спутниковые наблюдения в различных диапазонах спектра позволяют определить возрастной состав льда – важнейшую характеристику ледяного покрова. При этом по изображениям видимого диапазона опытный ледовый эксперт может при визуальном анализе определить до 8 градаций развития (возраста) льдов; причем ниласовые и молодые льды разделяются на несколько возрастных градаций, а более старые льды объединяются в одну категорию, не разделяемую по градациям. По ИК-изображениям могут быть определены толщины ниласовых, молодых и однолетних тонких и средних льдов. По радиолокационным изображениям дешифрируются старые льды. В условиях Арктики в зимний период диагностировать возрастной состав льдов можно путем совместного использования ИК– и СВЧ-диапазонов.
Автоматизированные методы картирования морских льдов по возрастным градациям с использованием радиолокационных данных. Технологии использования данных микроволнового диапазона для картирования льда по возрастным градациям развиваются с 1980-х гг. Одна из первых методик оперативного картирования ледовой обстановки была разработана для российского ИСЗ «Океан», единственного на тот период спутника, оснащенного радиолокатором бокового обзора. В процедурах картирования применялся как метод визуальной экспертной интерпретации изображений, так и алгоритмы автоматизированной количественной оценки льдов (Александров, Лощилов, 1985; Александров и др., 1989). Данные спутника «Океан» широко использовались в ААНИИ для ежедневного анализа оперативной ледовой обстановки по изображениям на отдельных орбитах и для составления еженедельных композитных ледовых карт для всей трассы СМП (Асмус и др., 2002; Бушуев и др., 2004).
Одним из методов объективного оценивания возрастных градаций ледяного покрова является метод нейронных сетей (НС) (Bogdanov, 2005; Hara et al., 1994), широко применяемый при обработке данных дистанционного зондирования благодаря ряду преимуществ перед традиционно используемыми статистическими методами классификации. Нейросетевые алгоритмы показывают лучшие результаты при наличии спекл-шума в спутниковых радиолокационных данных. Обученные НС имеют более высокое быстродействие, что позволяет использовать их для обработки больших массивов информации. Это особенно важно при оперативной работе со спутниковыми изображениями с высоким пространственным разрешением. Наиболее сложной задачей является настройка (тренировка) НС с использованием обучающей выборки. При распознавании образов, где обучающая выборка может не охватывать всех возможных состояний данных, важнейшим свойством НС становится способность классифицировать вектора данных, которые не использовались при тренировке алгоритма, то есть способность к обобщению. Это свойство позволяет использовать нейронные сети как универсальный классификатор в задачах дистанционного зондирования.
Метод НС был применен в ААНИИ для разработки метода картирования ледяного покрова по спутниковой радиолокационной информации SAR. На начальном этапе были определены текстурные характеристики изображения. Текстура изображения определяется статистическими взаимосвязями значений соседних пикселей радиолокационного изображения, что выражается в виде «узора» или «рисунка», воспринимаемого глазом при визуальном анализе изображения. Особенности формирования ледяного покрова, стадии его развития, состояние поверхности и другие факторы определяют внешний вид текстуры изображения. Характеристики текстур различных типов поверхностей вычисляются на основе матрицы совместной встречаемости уровней яркости.
Статистический анализ полученных текстурных характеристик показал, что корреляция, инерция и выпуклость кластера – это наиболее значимые признаки. Энергия, гомогенность и энтропия менее значимы, однако являются источниками дополнительной информации. Целесообразность использования этих характеристик подтверждена корреляционным анализом. Яркости SAR-изображения, как известно, подвержены значительному краевому эффекту: уменьшение удельной эффективной площади рассеяния (УЭПР) электромагнитных импульсов морским льдом с увеличением угла зондирования приводит к значительному уменьшению сигнала в дальней части полосы обзора относительно ближней. Для получения равноконтрастного изображения по всей полосе обзора была разработана методология приведения УЭПР морских льдов к одному углу зондирования (Александров, Пиотровская, 2008а; 2008б). Алгоритм приведения изображения к фиксированному углу зондирования включает в себя пересчет яркости SAR-изображения в значение УЭПР для данного угла зондирования, пересчет полученного значения для угла зондирования 25° и пересчет скорректированной УЭПР в новое значение яркости по ранее определенным коэффициентам, уникальным для каждого класса льда (однолетнего, однолетнего деформированного и старого льда).
Верификация алгоритма оценки возрастных градаций льда по методу НС проводилась в период Международного полярного года, на акватории российской Арктики, в частности, по району дрейфа СП-35. Были использованы SAR-изображения с европейского спутника Envisat с разрешением 150 м и полосой обзора 400 км. Благодаря высокоширотному положению траектории дрейфа станции, ее изображение (точнее морских льдов в области станции) можно было получать практически ежедневно. Согласно экспертной оценке, ледяное поле, на котором был построен лагерь дрейфующей станции, состояло из смерзшихся обломков полей старого льда различных размеров, толщины и конфигурации. Монолитные обломки полей старого льда разделены полосами однолетнего льда с включениями битого старого льда. Поэтому сеть настраивалась на выделение на спутниковой сцене трех классов: старого льда, однолетнего ровного льда и однолетнего деформированного льда. Для обучения этой сети использовался стандартный алгоритм обратного распространения ошибки. Сеть была настроена на классификацию SAR-изображений центрального района Арктики в зимний период года. Полученные результаты классификации соответствовали результатам экспертного дешифрирования. Рассчитанные ошибки классификации составили: для однолетнего ровного льда 15 %, для однолетнего деформированного льда 17 % и для старого льда 20 %. Обученная сеть также была использована для автоматического картирования этих же возрастных градаций льда по данным Envisat в Арктическом бассейне (рис. 1).
Рис. 1. Картирование возрастных стадий льда в Арктике с помощью метода нейронных сетей по данным Envisat. 16.01.2008 (слева – исходный снимок, в центре – снимок с угловой коррекцией, справа – классифицированное изображение) 1 – старый лед; 2 – однолетний ровный лед; 3 – однолетний деформированный лед
Для использования метода НС в автоматическом режиме, с целью оперативного картирования возрастных стадий льда, необходима настройка сети в направлении выделения большего количества классов ледовых образований.
Определять возрастной состав морских льдов по спутниковым данным в автоматизированном режиме можно также с помощью методов теории вероятности, в частности, применяя формулу Байеса для расчета апостериорной вероятности p(ωj/xi)
p(ωj/xi) = p(xi/ωj) × p(ωj)/p(xi), p(xi) = p(xi/ωj) × p(ωj),
где p(ωj) – априорная вероятность, p(xi/ωj) – условная плотность распределения величины xi в состоянии ωj.
Метод байесовской классификации обеспечивает оптимальное решение с точки зрения минимума вероятности ошибки. При классификации ледяного покрова на спутниковом снимке принимается решение в пользу того вида льда, для которого величина апостериорной вероятности p(ωj/xi) максимальна.
