Моделирование динамики яркостной температуры земли методом инвариантного погружения и нейронных сетей

Моделирование динамики яркостной температуры земли методом инвариантного погружения и нейронных сетей

КУРСОВАЯ РАБОТА

"МОДЕЛИРОВАНИЕ ДИНАМИКИ ЯРКОСТНОЙ ТЕМПЕРАТУРЫ ЗЕМЛИ МЕТОДОМ ИНВАРИАНТНОГО ПОГРУЖЕНИЯ И НЕЙРОННЫХ СЕТЕЙ"

Введение

Необходимость усиления контроля за глобальными процессами: изменения климата, ростом негативного антропогенного воздействия на биосферу, а также потребности в прогнозировании природных и техногенных катастроф, в ведении хозяйственной деятельности выдвигают в качестве приоритетной задачи дистанционное зондирование природной среды.

Дистанционное зондирование (ДЗ) можно представить как процесс, посредством которого собирается информация об объекте, территории или явлении без непосредственного контакта с ним. Методы ДЗ основаны на регистрации в аналоговой или цифровой форме отраженного или собственного электромагнитного излучения участков поверхности в широком спектральном диапазоне (от оптического и инфракрасного до коротковолнового).

Космическое зондирование, интенсивно развивающиеся в последние десятилетия, предоставило наукам о Земле новые возможности для исследования земной поверхности. За этот период существенно возросли объем, разнообразие и качество материалов ДЗ. К настоящему времени накоплен огромный фонд (более 100 миллионов) аэрокосмических снимков, полностью покрывающих всю поверхность Земли, а для значительной части районов с многократным перекрытием. [5]

Ученые и специалисты многих стран активно разрабатывают методики изучения Земли с использованием дистанционных измерений, совершенствуются и техническая база мониторинга, и методы интерпретации данных. Радиометрические методы исследования земной поверхности являются перспективными, поскольку не зависят от естественной освещенности местности и состояния атмосферы. Наблюдения в радиоволновом диапазоне чувствительны к влажностным, геометрическим и диэлектрическим характеристикам объектов. Это делает возможным ведение наблюдения и определение характеристик объектов, которые недоступны для оптических систем наблюдения.

Один из последних выведенных на орбиту спутников – ALOS помимо сенсоров видимого и ИК-диапазона оборудован РСА PALSAR с разрешением от 10 до 100 м (L-диапазон). В 2007 году готовится к запуску спутник SMOS, имеющий на борту усовершенствованный радиометр L-диапазона для изучения влажности почв и солености океанов. [13]

Однако обработка результатов радиолокационных и радиометрических исследований является более сложной по сравнению с оптическими системами наблюдения и требует особенного подхода, адаптированных к физическому и техническому алгоритму формирования данных. Также, важной задачей является повышение оперативности получения физических данных наблюдаемого объекта. Необходимы новые методы экспресс-анализа и быстрой обработки данных ДЗ в режиме реального времени.

Целью работы является применение моделирования динамики яркостной температуры методом инвариантного погружения и нейронных сетей; решение обратной задачи радиометрии – получение физических данных исследуемого объекта (почв); обзор различных моделей нейронных сетей; оценка погрешности и выявление причин ее появления.

моделирование динамика инвариантный нейронный

1. Задачи дистанционного зондирования земли. Современные методы обработки данных

1.1 Задачи, решаемые методами радиометрии

Задачи, которые ставит перед собой радиометрическое зондирование, имеют широкий спектр в области исследований различных параметров окружающей среды (геосферы, гидросферы и атмосферы).

Исследование земных покровов – одна из важных задач нынешнего времени. Данные этих исследований очень значимы для многих областей науки, экономики и быта человечества.

Существует огромное количество различных методов, моделей, алгоритмов анализа для расчета величин и параметров исследуемой поверхности. Каждая из этих моделей допускает свои упрощения и неточности.

Задачи, решаемые радиометрическим зондированием в геосфере Земли.

Радиометрическое зондирование используют для определения влажности почвы (требует частоты в 1–6 GHz), оценки уровня грунтовых вод. Вода обладает самой большой диэлектрической проницаемостью. Таким образом, если вода попадает в почву, то она будет вносить вклад в диэлектрическую проницаемость почвы, которая состоит из воздушных полостей, плёнок воды, которые покрывают частицы грунта. Поэтому с увлажнением почвы её диэлектрическая проницаемость начинает возрастать, следовательно, увеличивается коэффициент отражения.

Динамика радиояркостной температуры в процессе испарения почвенной влажности зависит от типа почвы, что позволяет оценивать физические и механические характеристики почв, находящихся в одинаковых метеоусловиях, дистанционным радиометрическим методом [1,12].

В радиометрии легко обнаружить загрязнение почвы нефтью. Нефть является гидрофобной жидкостью, то есть она препятствует проникновению воды в почву, в результате чего образуется плотная корка или стоит лужа в низине. В виду этого, сразу после дождя загрязнённая почва кажется более сухой и более яркой. Применяя радиометр на различных длинах волн можно определить толщину загрязнённого слоя. [6,16]

Дистанционное радиометрическое зондирование дает такую полезную информацию как концентрацию в ней солей и минеральных веществ, глубину промерзания и оттаивания почв.

Фиксация границ снежного покрова лучше происходит с помощью дистанционного зондирования в видимом диапазоне, однако при дешифрировании снежного покрова лесных территорий, возникает ряд затруднений, вследствие чего приходится ориентироваться на заснеженность открытых площадок: болот, вырубок, озер, или применять радиометрические методы в ИК или микроволновом диапазоне.

Оперативное картографирование снежного покрова и скорость отступания его границ в весенний период используются для гидрологических прогнозов. Некоторые параметры можно оценить косвенно. Например, зоны, охваченные снеготаянием, выявляются в ближнем инфракрасном диапазоне спектра, а мощность снежного покрова рассчитывается по ряду последовательных снимков, скорости продвижения границ снегонакопления и температуре воздуха. [9,14,16,18]

С помощью радиометрического ДЗ ведутся наблюдения за крупными чрезвычайными происшествиями и экологическими катастрофами (наводнения, засухи, лесные пожары, загрязнения морских вод нефтепродуктами, заморозки, ураганы, туманы, гололед, пыльные бури). Например, лесные пожары хорошо обнаруживаются даже на снимках малого разрешения, за счёт того, что излучение горящих углей очень большое, то есть большая яркостная температура. Однако существуют другие проблемы: большая задымлённость приводит к маскировке теплового излучения.

Торфяные пожары страшны тем, что находятся под поверхностью. Радиометрия может определить эти очаги пожара за счёт повышения температуры поверхности. [6,8]

Задачи, решаемые радиометрическим зондированием в гидросфере Земли.

Температура поверхности земли и океана является главным, хотя и не единственным фактором, воздействующим на яркостную температуру.

Длинноволновые радиометры (длина волны 20 см и более) в принципе могут оценивать соленость морской воды, используя зависимость поглощательной способности от солености. Точность таких измерений невелика – около одного промилле (1% ₀), что сравнимо по величине с естественными вариациями солености в Мировом Океане. Однако длинноволновые радиометры могут найти применение в арктических районах, где вариации солености достигаются 10–20%.

Для слежения за изменением ледовой обстановки в морях составляют ледовые карты. Важные преимущества космической съемки – повторяемость поступления информации и оперативность обработки – дают возможность фиксировать состояние быстро изменяющихся природных явлений на различные моменты времени. Автоматизированные технологии позволяют отличать льды от облаков и разделять лед по сплоченности.

В результате, по спутниковым данным, создаются динамические карты ледовой обстановки в период навигации, а также в осенне-зимний и весенний периоды (наступление ледостава, очищение ото льда).

В результате интерференции наблюдаемая яркостная температура участка водной поверхности с разливом нефти изменяется периодически с изменением толщины пленки нефти. Сравнение наблюдаемой яркостной температуры участка разлива нефти с яркостной температурой участка чистой водной поверхности позволяет провести измерения толщины пленки нефти в пятне загрязнения. [8,13,16]

Особенностью мероприятий по метеозащите крупных городов – использование наряду с радиолокационной информацией, отражающей пространственно-временную эволюцию жидко-капельных облаков и осадков СВЧ-радиометрической информации о содержании как парообразной, так и жидко-капельной влаги в атмосфере. В ходе работ выполняется анализ временной / пространственной изменчивости характеристик влагосодержания атмосферы (водозапас облаков, влагозапас атмосферы). [6,10,13,16]

Задачи, решаемые радиометрическим зондированием для других планет Солнечной Системы.

Дистанционное зондирование радиометрическими методами проводится не только для Земли, но и для других близлежащих планет Солнечной системы. Эти методы особенно оправданы для планет с «густой» атмосферой.

Радиолокационное картирование северного полушария планеты Венера космическими аппаратами Венера-15 и Венера-16, выполненное в 1983–1984 гг. советскими учеными, заслуженно является достижением мирового уровня. Впервые в мире с борта космических аппаратов была выполнена детальная радиолокационная съемка поверхности планеты, закрытой плотной атмосферой, непрозрачной для наблюдений в оптическом диапазоне. Площадь отснятой территории, расположенной севернее 300 С.Ш., равна 115 млн. км², что составляет четверть всей поверхности Венеры и лишь на треть меньше территории всей земной суши. Идея проведения эксперимента и его научно-методическая основа разработана в ИРЭ РАН. [2,15,18]

В рамках программы по исследованию планет солнечной системы и их спутников по проекту «Марс-96» совместно с учеными ИРЭ был разработан и установлен на космическом аппарате длинноволновый радар для глубинного зондирования грунта и ионосферы Марса в диапазоне рабочих частот от 0,2 до 5 МГц. Основным препятствием к широкому использованию радиолокаторов в геологической разведке земных недр является сильное поглощение радиоволн в почве из-за наличия в ней воды. Однако на других планетах и космических телах вода практически отсутствует, поэтому глубина проникновения может быть весьма большой. Основной целью эксперимента «Марс-96» являлось исследование высотного распределения электронной концентрации ионосферы Марса, измерение диэлектрических свойств грунта на разных глубинах вдоль трассы дрейфа космического аппарата, выявление глубинной структуры полярных областей. [2,15,18]

1.2 Недостатки современных методов обработки данных дистанционного зондирования

На данный момент в мире существует достаточно большое количество теорий, которые сами по себе верны, но при некоторых условиях не выполняются, дают сбой. Наверное, ещё не существует универсального способа или метода для определения тех или иных характеристик. Также и в радиометрии. Есть множество моделей, некоторые разные, некоторые почти одинаковые. Для примера можно рассмотреть характеристики почвы, которые надо учитывать, а которыми можно пренебрегать.

В идеальном случае мы считаем, что градиент температуры неизвестен, это есть функция, которую можно только аппроксимировать более простой, но точно задать невозможно. Или разбивать почвы на слои, желательно бесконечно малые и определять температуру каждого слоя. Также следует поступать с комплексной диэлектрической проницаемостью (КДП) и волновым числом, которое выражается через КДП. Зависит от КДП и коэффициент отражения каждого слоя, но эту зависимость можно определить с помощью простого математического аппарата. Влажность каждого слоя также нужно учитывать и определять.

Не надо забывать, что почва – это неоднородная среда. В почве встречаются камни, песчинки, живые организмы, мёртвые организмы, вода, трава, огромное количество природных ресурсов и многое другое. Каждая составляющая почвы имеет свою КДП, волновое число, влажность, температуру, коэффициент отражения. То есть для задания алгоритма по расчёту яркостной температуры и отражающей способности надо каждый слой разбивать на много подслоёв. Такое задание параметров почвы очень сильно осложняет работу и сводит её к нереализуемой, так как всё учесть невозможно, причём наука далеко ещё не всё знает об окружающей нас природе (среде). Также на поверхности почвы существуют шероховатости, неровности, растительность, органические соединения (живые существа). В зависимости от географического расположения изменяются характеристики среды, а также угол под которым растут растения на поверхности, вид растений (насколько сильно они поглощают тепло, как глубоко корни приникают в почву и прочие факторы). Даже если учесть все известные зависимости и закономерности максимально точно, то всегда будет какое-то допущение, упрощение. Для этого и создано это множество моделей, которые учитываю одно, досконально просчитывая все возможные варианты, и не учитывают другое.

Модели, используемые при тематической обработке радарных или радиометрических изображений, можно условно поделить на два класса: физические и статистические. Физические модели строятся на основе знания закономерностей собственного излучения или рассеяния волн. Они содержат функциональные соотношения, связывающие совокупность геофизических параметров исследуемого природного объекта с измеряемыми характеристиками принимаемого микроволнового излучения. Это позволяет построить количественные алгоритмы восстановления с использованием математических методов решения обратных задач, соответствующих форме найденных функциональных соотношений. Статистические модели не содержат априорных функциональных соотношений. Они рассчитаны на получение статистических оценок геофизических параметров с помощью выборочных значений для конкретной совокупности характеристик электромагнитных полей, которые получаются в процессе зондирования, и геофизической информации, собираемой с тестовых участков. Этот подход широко использует обучение по выборкам и нейронные сети. Построение статистических моделей трудоемко, и они обычно справедливы для конкретных природных объектов. Однако ввиду сложности процессов собственного излучения и рассеяния электромагнитных волн при построении физических моделей для большинства реальных природных объектов часто встречаются непреодолимые трудности. Поэтому в настоящее время используются оба вида моделей в зависимости от сложности зондируемых природных комплексов. В ряде случаев применяются комбинированные модели, в которых используются как элементы статистических оценок, так и физические закономерности взаимодействия волн с природной средой.

2. Необходимость использования нейронных сетей в моделировании динамики яркостной температуры земли при дистанционном зондировании земли

Одним из направлений дальнейшего совершенствования алгоритмов обработки данных ДЗ является использование искусственных нейронных сетей, обладающих способностью реализовать широкий класс функций без априорных допущений о законах их распределения. На основе их обучения может быть обеспечено более точное и оперативное получение искомых параметров в реальном масштабе времени. Разработка эффективных алгоритмов, основанных на искусственных нейронных сетях, для выполнения задач подобного класса представляет собой актуальную проблему данного исследования.

2.1 Нейронные сети и их применение

Нервная система и мозг человека состоят из нейронов, соединенных между собой нервными волокнами. Нервные волокна способны передавать электрические импульсы между нейронами. Все процессы передачи раздражений от нашей кожи, ушей и глаз к мозгу, процессы мышления и управления действиями – все это реализовано в живом организме как передача электрических импульсов между нейронами. Рассмотрим строение биологического нейрона. Каждый нейрон имеет отростки нервных волокон двух типов – дендриты, по которым принимаются импульсы, и единственный аксон, по которому нейрон может передавать импульс. Аксон контактирует с дендритами других нейронов через специальные образования – синапсы, которые влияют на силу импульса.

Рис. 1. Биологическая нейронная сеть

Можно считать, что при прохождении синапса сила импульса меняется в определенное число раз, которое мы будем называть весом синапса. Импульсы, поступившие к нейрону одновременно по нескольким дендритам, суммируются. Если суммарный импульс превышает некоторый порог, нейрон возбуждается, формирует собственный импульс и передает его далее по аксону. Важно отметить, что веса синапсов могут изменяться со временем, а значит, меняется и поведение соответствующего нейрона.

Нетрудно построить математическую модель описанного процесса.

Рис. 2. Схематичное изображение обработки сигнала нейроном

На рисунке изображена модель нейрона с тремя входами (дендритами), причем синапсы этих дендритов имеют веса w₁, w₂, w₃. Пусть к синапсам поступают импульсы силы x₁, x₂, x₃ соответственно, тогда после прохождения синапсов и дендритов к нейрону поступают импульсы w₁x₁, w₂x₂, w₃x₃. Нейрон преобразует полученный суммарный импульс x=w₁x₁+ w₂x₂+ w₃x₃ в соответствии с некоторой передаточной функцией f(x). Сила выходного импульса равна y=f(x)=f(w₁x₁+ w₂x₂+ w₃x₃). Таким образом, нейрон полностью описывается своими весами w_k и передаточной функцией f(x). Получив набор чисел (вектор) x_k в качестве входов, нейрон выдает некоторое число y на выходе.

Как работает нейронная сеть

Искусственная нейронная сеть (ИНС, нейронная сеть) – это набор нейронов, соединенных между собой. Как правило, передаточные функции всех нейронов в нейронной сети фиксированы, а веса являются параметрами нейронной сети и могут изменяться. Некоторые входы нейронов помечены как внешние входы нейронной сети, а некоторые выходы – как внешние выходы нейронной сети. Подавая любые числа на входы нейронной сети, мы получаем какой-то набор чисел на выходах нейронной сети. Таким образом, работа нейронной сети состоит в преобразовании входного вектора в выходной вектор, причем это преобразование задается весами нейронной сети.

Биологическая эволюция, которая привела к столь эффективным решениям, шла по пути от образов к логике. Так и человек после рождения сначала учится распознавать образы, а только потом приобретает умение рассуждать логически и строить алгоритмы. Компьютеры же, напротив, начав с логики, лишь спустя несколько десятилетий осваивают распознавание образов за счет создания специальных программ для компьютеров традиционной архитектуры или благодаря созданию специализированных аппаратных нейропроцессоров.

Искусственные нейронные сети, подобно биологическим, являются вычислительной системой с огромным числом параллельно функционирующих простых процессоров с множеством связей. Несмотря на то что при построении таких сетей обычно делается ряд допущений и значительных упрощений, отличающих их от биологических аналогов, искусственные нейронные сети демонстрируют удивительное число свойств, присущих мозгу, – это обучение на основе опыта, обобщение, извлечение существенных данных из избыточной информации.

Нейронные сети могут менять свое поведение в зависимости от состояния окружающей их среды. После анализа входных сигналов (возможно, вместе с требуемыми выходными сигналами) они самонастраиваются и обучаются, чтобы обеспечить правильную реакцию. Обученная сеть может быть устойчивой к некоторым отклонениям входных данных, что позволяет ей правильно «видеть» образ, содержащий различные помехи и искажения.

В 50-х годах прошлого века группа исследователей объединила биологические и физиологические подходы и создала первые искусственные нейронные сети. Тогда казалось, что ключ к искусственному интеллекту найден. Но, хотя эти сети эффективно решали некоторые задачи из области искусственного зрения – предсказания погоды и анализа данных, иллюзии вскоре рассеялись. Сети были не в состоянии решать другие задачи, внешне похожие на те, с которыми они успешно справлялись. С этого времени начался период интенсивного анализа. Были построены теории, доказан ряд теорем. Но уже тогда стало понятно, что без привлечения серьезной математики рассчитывать на значительные успехи не следует.

С 70-х годов в научных журналах стали появляться публикации, касающиеся искусственных нейронных сетей. Постепенно был сформирован хороший теоретический фундамент, на основе которого сегодня создается большинство сетей. В последние два десятилетия разработанная теория стала активно применяться для решения прикладных задач. Появились и фирмы, занимающиеся разработкой прикладного программного обеспечения для конструирования искусственных нейронных сетей. К тому же 90-е годы ознаменовались приходом искусственных нейронных сетей в бизнес, где они показали свою реальную эффективность при решении многих задач. В числе задач, решение которых доверяют искусственным нейронным сетям, можно назвать следующие – это распознавание текста и речи, семантический поиск, экспертные системы и системы поддержки принятия решений, предсказание курсов акций, системы безопасности, анализ текстов. Рассмотрим несколько особенно ярких и интересных примеров использования нейронных сетей в разных областях.

Техника и телекоммуникации

В 1996 году фирмой Accurate Automation Corp (#"479434.files/image003.jpg">

Рис. 3. Схема передвижного штатива с радиометрами

Рис. 4. Передвижение рамы с радиометрами с одного участка на другой

Рис. 5. Экспериментальные участки

Конструкция и размеры передвижного штатива (рамы) были выбраны исходя из удовлетворения следующим критериям [1]:

1. Эталоны и исследуемые участки должны находиться в дальней зоне антенн радиометров

2. Ширина главного лепестка по уровню 0,5 дб должна быть меньше геометрических размеров, эталонов и участков.

3. Угол наклона антенн должен составлять 10⁰

Участки №1, №2, №3, №4 являются экспериментальными участками, на которых расположены следующие виды почв: земля, загрязненная нефтью, песок, глина, чернозем, соответственно. На участке №5 расположен лист металла, а на участке №6 – емкость с водой. Участки №5 и №6 используются для калибровки оборудования, а измерения снятые с этих участков являются опорной точкой (калибровкой) для расчета яркостных температур почвы участков №1 – №4.

Штатив (рама) передвигается с одного на другой участок, измерения проводятся последовательно радиометрами с частотой 5 Gh, затем 11 Gh, затем 3,6 Gh.

При калибровке радиометров используются следующие эталоны излучения: излучения неба, отраженное металлическим листом, излучение гладкой водной поверхности, излучение поглощающего покрытия. Размеры эталонов должны превышать размеры пятна, излучающего в главный лепесток, чтобы дифракционными явлениями на краях образцов можно пренебречь.

В ходе эксперимента замерялись: время, температура окружающей среды, температура слоя почвы на глубине 0,5 см и 2 см, а также измерялась яркостная температура почв.

Дважды в день брались пробы почв на влажность: поверхностный слой 0–1 см, 1–2 см, и 3–4 см.

Исследуемыми объектами являлись участки №2 и №3, песчаная и глинистая почва соответственно (рис. 4; рис. 5).

Рис. 6. Песчаная почва. Участок №2

Рис. 7. Средний суглинок. Участок №3

Из полученных данных видно, что почти одинаковые в оптическом диапазоне почвы, кардинально отличаются по физическим свойствам.

Таблица 1. Гранулометрический состав почв (% от массы сухой почвы)

Был проведен цикл измерений в течении четырех дней яркостной температуры почвы «сухая-влажная-сухая». Цикл измерений состоит из 39 точек, для трех радиометров разной частоты (и соответственно разной глубины проникновения в почву). Измерения на участках проводятся между двумя опорными точками: калибровкой металлом и водой.

3.2 Постановка задачи. Поиск алгоритма решения

Задачей построенной нейронной сети должно быть восстановление параметров почв по экспериментальным данным значения яркостной температуры с радиометров различной частоты в момент экспериментального получения проб влажности почв.

Для моделирования методом Нейронных сетей используются несколько программ, среди них прикладной пакет для MathLab 6.1, Statistica Neural NetWorks, NNMath 3.1 и др. Будем использовать программы MathLab 6.1, Statistica Neural NetWorks, выбранные за простоту работы и возможности настройки нужной модели.

В общем случае, для получения алгоритмов и моделей выполняется следующее:

1. Определяется изучаемый класс объектов (представленный входными и выходными данными).

2. Для этого класса выбирается настраиваемая модель (модель, параметры которой можно изменить), удовлетворяющая определенным критериям и требованиям.

3. Выбирается оценка качества идентификации (потери, характеризующие различие между выходными величинами объекта и модели).

4. Формируется алгоритм идентификации, который, изменяя параметры модели, минимизирует потери.

Выбор и разработка моделей и алгоритмов требуют серьезных усилий для экспериментального исследования и сравнения с уже ранее предложенными. В то же время это предоставляет большую свободу в выборе направления в науке, знания которой будут привлекаться для создания новых моделей и алгоритмов.

Будем пользоваться двумя подходами в решении поставленной задачи:

1. Используем модель двухслойного персептрона, и поэтапно увеличиваем количество нейронов на каждом слое с 3 до 5, параметры которых оптимально подобраны в исследовании[15].

2. Функция автоматического подбора параметров, которая присутствует в программе Statistica Neural NetWorks и работает на основе анализа количества и вида входных и выходных данных.

Рис. 8. Двухслойный персептрон с одним скрытым слоем и 3 нейронами на каждом слое

Рис. 9. Двухслойный персептрон с одним скрытым слоем и 5 нейронами на каждом слое

Рис. 10. Нейронная сеть построенная методом автоматического подбора параметров

Радиометрическое исследование природной среды и восстановление требуемых параметров выполняется в работах. С учетом предлагаемой в рефракционной модели КДП почв связь между радиоизлучательной способностью почвогрунта и его КДП может быть разрешена в явном виде относительно объемной влажности. При этом решение содержит ряд априорно задаваемых параметров. Входными данными будут являться яркостные температуры, полученные при теоретическом расчете методом инвариантного погружения. Выходными данными являться параметры почв: объемная влажность, уровень слоя, и КДП почвы (при анализе которой можно определить класс почвы).

В решении поставленной задачи будем считать отсутствие как шума, создаваемого растительностью (т.к. участки в эксперименте были подготовлены и очищены от растительности), так и техногенного шума. При наличии же таковых, представляется возможным создание комплекса нейронных сетей, одни из которых будут отвечать за отчистку радиометрического сигнала от шумов, другая – решать основную задачу, в противном случае целесообразно использование вейвлет фильтров для отчистки сигнала.

3.3 Метод инвариантного погружения. Теоретический расчет. Режим обучения

Метод инвариантного погружения берет за основу слоистость почвы и различное влагосодержание в различных слоях почвы. Слоистая модель по Башаринову А.Е.

Данная модель содержит следующие допущения:

1. Излучение считается некогерентным.

2. Нет ослабления излучения между поверхностью и антенной.

3. Яркость неба считается изотропной и имеет значение 3К.

4. Влажность и температура считаются функциями только глубины.

5. Диэлектрические и тепловые свойства почв постоянны в слоях определённой толщины.

6. Поверхность почвы считается гладкой. (растительность отсутствует)

Также будем считать, что диэлектрическая проницаемость имеет комплексный вид, тем самым рассмотрим наиболее общий случай.

С увеличением глубины, диэлектрическая проницаемость, а также влажность увеличивается.

Эта модель рассматривает тепловое излучение слоисто-неоднородных неизотермических детерминированных сред с произвольной величиной поглощения. Диэлектрическая проводимость имеет общий комплексный вид:

Пусть температура слоёв имеет следующую зависимость:

,

где z – глубина почвенных слоёв

Для расчета КДП почвы (e) в работе используется рефракционная модель диэлектрических свойств смесей. Она представляет собой описание влажностной зависимости диэлектрической проницаемости почв. Модель является наиболее успешной при работе в диапазонах СВЧ, учитывает двойственность диэлектрических свойств почв, определяемых содержанием связанной и свободной воды [15].

Рис. 11. Слоистая модель по Башаринову А.Е.

Данная модель учитывает связь свободных и связанных молекул воды в почвах.

Для обучения нейронной сети были использованы данные теоретических расчетов по программе на базе Microsoft Excel, рассчитывающей яркостную температуру по данным КДП полученным, при моделировании методом инвариантного погружения в лаборатории Радиометрического зондирования Земли Омского Государственного Педагогического Университета на основе исследования [15].

Коэффициент диэлектрической проницаемости почв также используется в режиме обучения и восстанавливается по яркостным температурам радиометров различной частоты в нейронной сети с целью классификации типа почв.

Таблица 2. Пример обучающей выборки для соотношений параметров песчаной почвы и яркостных температур при различных частотах радиометра

Количество обучающих данных было 120 измерений, что в 5 крат больше количества связей в двухслойном персептроне с 3 нейронами на каждом слое. Количество эпох обучения равнялось 100 в один период установления весов.

Рис. 12. Расчет весов в зависимости от количества эпох при обучение нейронной сети, созданной методом автоматического подбора параметров

Рис. 13. Расчет весов в зависимости от количества эпох обучения при обучение нейронной сети двухслойного персептрона методом обратного распространения ошибки

Рис. 14. Расчет весов в зависимости от количества эпох обучения при обучение нейронной сети двухслойного персептрона методом сопряженного градиента

Из рисунков видно, что наиболее быстро «обучаемой» нейронной сетью, является двухслойный персептрон, устанавливающий параметры весов методом сопряженного градиента. При увеличении количества нейронов на слоях нейронной сети период становления увеличивается во всех вышеприведенных методах обучения нейронной сети.

3.4 Решение поставленной задачи на экспериментальных данных методом нейронных сетей

Из экспериментальных данных полученных на производственной практике делается тестовая выборка яркостных температур трех радиометров в момент получения проб влажности почв. Таким образом для тестирования нейронной сети у нас будут два эталона: теоретический расчет яркостной температуры по методу инвариантного погружения на основе экспериментальных данных влажности почв и экспериментальные данные влажности почв.

Рис. 15. Критические или противоречивые данные

Ошибки и фиксации противоречий возникали при восстановлении параметров почв по яркостным температурам, на тех данных, которые были получены в момент эксперимента, когда один из радиометров отключался или выдавал нечеткий результат.

В обученную нейронную сеть, созданную методом автоматического подбора модели и параметров при решении прямой задачи (получение яркостной температуры почв с заданными параметрами) вводим входные данные: экспериментально полученные влажность, глубину слоя.

Получаем на выходе из нейронной сети яркостную температуру для трех радиометров определенные, в среднем, с точностью ±5К.

Таблица 3

В обученную нейронную сеть, созданную методом автоматического подбора модели и параметров при решении обратной задачи (получение параметров почв при заданной яркостной температуре с трех радиометров) вводим входные данные: экспериментально полученные яркостную температуру с трех радиометров разной частоты.

Получаем на выходе из нейронной сети данные об объемной влажности почв определенные, в среднем, с точностью ±0,04.

Таблица 4

В обученную нейронную сеть двухслойный персептрон при решении прямой задачи (получение яркостной температуры почв с заданными параметрами) вводим входные данные: экспериментально полученные влажность, глубину слоя.

Получаем на выходе из нейронной сети яркостную температуру для трех радиометров определенные с точностью ± 3К для двухслойного персептрона с 3 нейронами на каждом слое и ± 2К для двухслойного персептрона с 5 нейронами на каждом слое.

Таблица 5

В обученную нейронную сеть двухслойный персептрон при решении обратной задачи (получение параметров почв при заданной яркостной температуре с трех радиометров) вводим входные данные: экспериментально полученные яркостную температуру с трех радиометров разной частоты.

Получаем на выходе из нейронной сети данные об объемной влажности почв определенные с точностью ± 0,03 для двухслойного персептрона с 3 нейронами на каждом слое и ± 0,02 для двухслойного персептрона с 5 нейронами на каждом слое.

Таблица 6

3.5 Полученные результаты

В процессе обучения ИНС формирует некоторое скрытое, неявное знание о способе решения поставленной задачи, которое заложено в виде значений весов ее адаптивных параметров. Если требуется лишь получение результата, удовлетворяющего требованиям поставленной задачи, сама ИНС и является конечным решением задачи.

Таблица 7. Погрешность определения данных методом нейронных сетей различных моделей

Таблица 8. Скорость обучения нейронных сетей различных моделей

Из приведенных выше таблиц видно что наиболее оптимальной моделью точной моделью для решения данной задачи является двухслойный персептрон с 5 нейронами на каждом слое, а оптимальным методом обучения является метод сопряженного градиента.

Из полученных данных наиболее достоверными считаются результаты, определенные с точностью ±2К для T_я и ±0.02 для W и Wsl. Обученная ИНС, показывающая наилучшие результаты на тестовой выборке и имеющая необходимый набор служебных процедур для работы с измеряемыми данными и выводом искомых параметров, является в итоге рабочим нейросетевым решателем[15].

Полученная погрешность объясняется следующими допущениями:

·   Недостаточный объем данных в режиме обучения

·   Трехканальность входных и выходных данных

·   Неучет шероховатости поверхности

·   Неучет динамики температуры внешней среды

·   Пренебрежение шумом растительности

·   Пренебрежение техногенным шумом

Выводы

Построены нейронные сети для определения параметров почв с погрешностью в 14–19% и классификации почв на основе разработанного нейросетевого способа определения неоднородности распределения влаги в приповерхностном слое почвы по данным многоканальных измерений в СВЧ-диапазоне [15].

Наиболее удачной в решении как прямой так и обратной задачей можно считать двухслойный персептрон с 5 нейронами на каждом слое (дальнейшее увеличение количества нейронов приводит к ухудшению результатов из-за большого количества связей и малого объёма входных данных в режиме обучения).

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

В настоящее время перспективы дальнейшего развития данных методик обработки данных радиометрического зондирования не вызывает сомнения. Сейчас уже работают комплексы нейронных сетей по определению: влажности, температуры почв.

Серия работ Л.Е. Назарова посвящена вопросам нейросетевой классификации земных объектов (лесов, водоемов).

Разработана методика нейросетевого определения участков лесного пожара по данным RADARSAT-1.

Список литературы

1. Под ред. акад. РАН В.А. Садовничьего. Космическое землеведение. – М.: Изд-во МГУ, 1992. – Ч. 1. – 269 с.; 1998. – 4.2. – 571 с.

2. Башаринов А.Е., Гурвич А.С., Егоров СТ. Радиоизлучение Земли как планеты. – М.: Наука, 1974. – 207 с.

3. Шутпко A.M. СВЧ-радиометрия водной поверхности и почвогрунтов. – М.: Наука, 1986. – 190 с.

4. Арманд Н.А., Крапивин В.Ф., Мкртчян Ф.А. Методы обработки данных радиофизического исследования окружающей среды. – М.: Наука, 1987. – 270 с.

5. Шанда Э. Физические основы дистанционного зондирования: Пер. с англ. – М.: Недра, 1990. – 208 с.

6. Кондратьев К.Я. Ключевые проблемы глобальной экологии // Теоретические и общие вопросы географии. – М.: ВИНИТИ, 1990. – 454 с. – (Итоги науки и техники; Т. 9.)

7. Аэрокосмические методы в почвоведении и их использование в сельском хозяйстве. – М.: Наука, 1990. – 247 с.

8. Сост. М. Назиров, А.П. Пичугин, Ю.Г. Спиридонов. Под ред. Л.М. Митника, СВ. Викторова. Радиолокация поверхности Земли из космоса. – Л.: Гидрометеоиздат, 1990. – 200 с.

9. Баранов Д.В., Бобров П.П. Моделирование и экспериментальное исследование собственного радиотеплового излучения влажных почв. // Дипломная работа на соискание степени бакалавра радиофизики – Омск – 2006 – 30 с.

10. Караваев Д.М., Щукин Г.Г. СВЧ-радиометрические исследования влагозапаса атмосферы и водозапаса облаков. Тезисы докладов региональной XXIII конференции по распространению радиоволн. С-Петербург, 1997, с. 76.

11. «Потенциальные возможности бистатического радиометра для наблюдения поверхности Земли с высоким разрешением» // А.П. Верещак, В.В. Пискорж. – Журнал Радиоэлектроники – 2003 – №3

12. Баррет Э., Куртис Л. «Введение в космическое землеведение.» – пер. с англ. – М – Прогресс – 1979 г.

13. «Наблюдение океана из космоса при помощи микроволновых радиометров» Ю.А. Кравцов // Соросовский Образовательный Журнал – 1999 – (44)№7.

14. Медведев В.С., Потемкин В.Г. «Нейронные сети MatLab 6» М – ДиалогМИФИ – 2002 г.

15. Мансуров А.В. дисс. канд. ф.-м.н. «Алгоритмы обработки данных радиоволнового дистанционного зондирования поверхности Земли на основе искусственных нейронных сетей», Алтайский Государственный Университет, Барнаул – 2006