Математическое и программное обеспечение системы анализа динамики земной поверхности по данным дистанционного зондирования Земли

Замятин, Александр Владимирович

Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей

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

кандидата технических наук: Замятин, Александр Владимирович
город: Томск
год: 2005
специальность ВАК РФ: 05.13.11

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

Автореферат диссертации по теме "Математическое и программное обеспечение системы анализа динамики земной поверхности по данным дистанционного зондирования Земли"

На правах рукописи

Замятин Александр Владимирович

Математическое и программное обеспечение системы анализа динамики земной поверхности по данным дистанционного зондирования Земли

05.13.11— Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей

Автореферат

диссертации на соискание ученой степени кандидата технических наук

Томск — 2005

Работа выполнена в Томском политехническом университете

Научный руководитель:

доктор технических наук, профессор Н.Г. Марков

Официальные оппоненты:

доктор технических наук, профессор А.М. Кори ков доктор технических наук, доцент А.В. Скворцов

Ведущая организация:

Институт вычислительной математики и математической геофизики СО РАН, г. Новосибирск.

Защита состоится « 16 » марта 2005 г. в 15 ч. в ауд. 214 на заседании диссертационного совета Д 212.269.06 при Томском политехническом университете по адресу: 634034, г. Томск, ул. Советская, 84, Институт «Кибернетический центр» ТПУ.

С диссертацией можно ознакомиться в библиотеке Томского политехнического университета по адресу: 634034, г. Томск, ул. Белинского, 53.

Автореферат разослан февраля 2005 г.

Ученый секретарь

диссертационного Совета /)

к.т.н., доцент {/^/¿(¿^л/ МА Сонькин

Общая характеристика работы

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

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

За последние десятилетия XX века отечественными учеными (Берлянт A.M., Виноградов Б. В., Тикунов B.C. и другие) был выполнен большой объем исследований, посвященных, в основном, теоретическим и в меньшей степени практическим аспектам решения проблем анализа динамики земной поверхности, в том числе и с использованием данных ДЗЗ. Первые практически интересные результаты при решении этих проблем получены зарубежными учеными (Агравал Ц., Бэйкер В., Белл Э., Вербург П., Хаггет П. и многими другими). Однако в их работах не уделяется должного внимания этапу автоматизированной интерпретации данных ДЗЗ, являющемуся ключевым этапом при анализе динамики земной поверхности.

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

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

Исследования и разработки по теме диссертационной работы проводились в соответствии с утвержденными планами НИР Института «Кибернетический центр» ТПУ в 2000-2005 гг., а также в рамках проекта Минпромнауки РФ по гранту №00-15-98478 поддержки ведущих научных школ России, проектов РФФИ по грантам №00-07-90124, №03-07-90124, № 03-07-06024 и гранту Федерального агентства по образованию №432.

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

Для достижения поставленной цели необходимо последовательное решение следующих задач:

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

2. Разработка математического обеспечения (методов и алгоритмов) системы анализа динамики земной поверхности. Решение данной задачи предполагает также исследование эффективности предлагаемых методов и алгоритмов.

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

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

Научную новизну полученных в работе результатов определяют:

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

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

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

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

Практическая ценность и реализация результатов работы. Практически значимыми являются созданные методы, алгоритмы и программное обеспечение (ПО) системы анализа динамики земной поверхности, включая ПО для интерпретации АИ, имеющее самостоятельную практическую ценность. Программные средства этой системы функционируют на компьютерах типа IBM PC под управлением ОС Windows 98/NT/2000. Объем исходного кода системы составляет более 9000 строк на языках C++ и Object Pascal.

Созданные программные средства были использованы при выполнении х/д №8-21/03 между Институтом «Кибернетический центр» Томского политехнического университета и ОАО «Востокгазпром» и внедрены в ОАО «Востокгаз-пром» при решении задач автоматизированного получения тематических карт территории Северо-Васюганского газоконденсатного месторождения Томской области по панхроматическим и многозональным космическим снимкам системы ДЗЗ 1RS. Результаты внедрения подтверждены соответствующим актом.

Основные положения, выносимые на защиту:

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

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

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

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

Апробация работы. Основные результаты работы докладывались и обсуждались на следующих конгрессах, симпозиумах, конференциях и семинарах: Международные конференции Association of Geoinformation Laboratories of Europe AGILE V, VI, VII (г. Пальма, Испания, 2002 г.; г. Лион, Франция, 2003 г.; г. Ираклион, Греция, 2004 г.), Русско-американский семинар «Studies of Socio-Natural Co-Evolution from Différent Parts of the World» (г. Новосибирск, 2002 г.), IX Международный симпозиум "Remote Sensing 2002" SPIE, (г. Ираклион, Греция, 2002 г.), V Международный семинар «Computer Science and Information Technologies» (г. Уфа, 2003 г.), V и VI Всероссийские научно-

технические конференции «Нейроинформатика— 2003» (г. Москва) и «Нейро-информатика — 2004» (г. Москва), XX Международный конгресс ISPRS (г. Стамбул, Турция, 2004 г.), X Юбилейная международная научно-практическая конференция студентов и молодых ученых «Современные техника и технологии» (г. Томск, 2004 г.), I Международная конференция «Земля из космоса — наиболее эффективные решения» (г. Москва, 2004 г.).

По результатам диссертационных исследований опубликовано 14 работ, в том числе 11 статей.

Личный вклад:

1. Постановка задач исследования и разработка концепции построения системы анализа динамики земной поверхности с использованием аэрокосмических изображений выполнены автором совместно с Н.Г. Марковым.

2. Математическое обеспечение системы анализа динамики земной поверхности разработано автором.

3. Постановки задач исследования эффективности предложенных методов и алгоритмов выполнены автором совместно с Н.Г. Марковым. Результаты исследования этих алгоритмов получены лично автором.

4. Разработка программного обеспечения системы анализа динамки земной поверхности выполнена автором, за исключением подсистемы формирования текстурных признаков, разработанной совместно с Напрюшкиным А.А. Программная реализация алгоритмов параметрической оценки плотности и оценки плотности распределения Розенблатта-Парзена выполнена Напрюшкиным А.А.

Объем и структура работы. Диссертация состоит из введения, четырех глав, заключения, списка использованных источников из 144 наименований и приложения. Объем основного текста диссертации составляет 148 страниц машинописного текста, иллюстрированного 46 рисунками и 11 таблицами.

Содержание работы

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

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

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

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

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

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

Этап усовершенствованной, более точной тематической интерпретации АИ нацелен на преодоление недостатков традиционной интерпретации и базируется на ряде сформулированных в работе принципов:

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

2. Раздельное использование спектральных и текстурных признаков при классификации, направленное на улучшении детальности интерпретации.

3. Использование для проведения классификации всей имеющейся информации без какого-либо предварительного редуцирования признакового пространства.

4. Совместное использование статистических (параметрических и непараметрических) и нейросетевых подходов к классификации АИ.

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

При реализации предложенных методов классификации формирование ВК может осуществляться двумя способами. В статистическом методе классификации текстурные признаки предлагается учитывать с помощью марковских случайных полей первого порядка и наиболее распространенных текстурных характеристик Харалика— энтропия, энергия и вариация, определяемых по гистограммам частот первого и второго порядка для анализируемого окна размером 5x5 и 9x9. Такой способ формирования ВК позволяет учитывать пространственные особенности АИ в различных масштабах. Указанные текстурные характеристики рассчитываются для каждого из каналов исходного (панхроматического или многозонального) АИ, поэтому ВК представляет собой признаковое пространство большой размерности, требующее дополнительного преобразования.

Формирование ВК в нейросетевом методе классификации предлагается осуществлять более просто — путем формирования признакового вектора контекстно-спектральным способом. Суть этого способа заключается в том, что каждый компонент вектора представляет собой значения яркостей пикселей в некоторой выбранной окрестности для всех каналов АИ. Такой способ позволяет учитывать текстурные особенности типов земной поверхности без их специального расчета и оптимизации признакового пространства.

Моделирование изменений земной поверхности в рамках многоэтапного подхода предложено выполнять с использованием статистического аппарата марковских цепей. В этом случае вероятностная информация об изменениях земной поверхности содержится в матрице вероятностей переходов Р = \р>], каждый элемент которой является вероятностью перехода типа земной поверхности оо, в тип со,, где у - 1,...М, М— число типов земной поверхности на изучаемой территории, выявленное на этапе интерпретации АИ. Количественная информация об изменениях содержится в матрице фактических переходов территорий Мф>кт = [т,^] и матрице ожидаемых переходов территорий Можид = [и,/*]. Матрица Мф^ содержит количество элементов изображения (пикселей), перешедших из типа земной поверхности о, в тип а Можид — ко-

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

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

Первый принцип декларирует использование априорной информации о вероятности развития (угнетения) каждого из типов земной поверхности на исследуемой территории. Такую информацию удобно представить в виде А/ вероятностных карт, называемых картами соответствия, где М— число типов земной поверхности на исследуемой территории. В каждой точке ьй карты исследуемой территории содержится вероятность замещения типа поверхности со, на тип (Яр обозначаемая р,/"". Формироваться такие карты могут с использованием функций пространственного анализа геоинформационных систем (ГИС), позволяющих строить вероятностные растровые изображения.

Второй принцип декларирует использование при формировании правил функционирования КА вероятностных и пространственных характеристик каждого из типов земной поверхности исследуемой территории, рассчитываемых по исходной тематической карте, созданной на этапе интерпретации АИ. Вероятностные характеристики вычисляются на основе подхода, суть которого в том, что в некоторой анализируемой окрестности точки (пикселя) тематической карты, вероятность перехода типа поверхности ш, в тип со, зависит не только от вероятности р9 из матрицы Р, но и от количества элементов типа со, в этой окрестности. Таким образом, для каждого типа земной поверхности в анализируемой окрестности вычисляется вероятность где т — общее количество типов земной поверхности в анализируемой окрестности, а п — количество элементов типа <а; в этой окрестности.

Использование особенностей пространственных характеристик типов земной поверхности основано на применении критерия насыщенности F, значение которого определяется как отношение частоты встречаемости конкретного типа земной поверхности в анализируемой окрестности точки h к встречаемости этого типа в целом на изображении (тематической карте)

р _ пк,*М / "<и,

где характеризует насыщенность анализируемой окрестности с расстоянием от центрального элемента скользящего окна d (например, при d =1 размер окна 3x3) в точке /г типа поверхности о*; ,/,/>— число элементов земной поверхности типа в анализируемой окрестности, —общее число элементов в анализируемой окрестности, N. — число элементов типа на всем изобра-

жении, a N— общее число элементов изображения. Расчет усредненных характеристик критерия насыщенности Fhkd для каждого типа Ш/ проводится согласно выражению:

. ( 1 Ï

(2)

где L — множество элементов изображения, относящихся к типу <й/, a N/— общее число элементов типа to/ на изображении.

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

Учитывая проблемы, затрудняющие формирование вероятностных карт соответствия, и, следовательно, реализацию первого принципа, результирующая вероятность перехода p/f" может учитывать только указанные вероятностную и пространственную составляющие и будет иметь вид Если удается выявить дополнительные факторы, влияющие на развитие процессов на земной поверхности, и, следовательно, построить вероятностные карты соответствия, то результирующая вероятность в этом случае будет иметь вид рif3 Рц""). Именно такой набор вероятностных составляющих

предложено использовать для определения правил функционирования КА.

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

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

Предложен алгоритм, реализующий метод статистической классификации. На первом шаге этого алгоритма проводится формирование исходного множества текстурных признаков согласно описанному выше способу формирования ВК. На втором шаге с помощью критерия «хи-квадрат» осуществляется проверка каждого элемента ВК на согласованность с нормальным распределением. Это позволяет на третьем шаге преобразовать исходную ВК в ряд подмножеств, согласованных и не согласованных с нормальным законом распределения. В свою очередь, это дает возможность обрабатывать соответствующие подмножества параметрическим или непараметрическим байесовским классификатором и тем самым повышать вычислительную эффективность обработки ВК. Апостериорные вероятности pjifù^Xj"9) (/'= \ количество подмно-

жеств), полеченные с помощью обработки каждого подмножества «своим» классификатором, агрегируются в общую апостериорную вероятность jо(о); I Л*1р), которая используется в дальнейшем. Такое агрегирование может проводиться путем усреднения всех или выбором максимального

1 На четвертом шаге полученные значения р{о, | Л"1*) трансформи-

руются в вероятностные карты для каждого типа поверхности На

заключительных пятом и шестом шагах алгоритма, представляющих собой реализацию второго этапа двухэтапной классификации, признаковое пространство, представленное в ПК, проверяется на согласованность с нормальным распределением. Затем проводится байесовская классификация с признаковым вектором Хт<'> в которой полученная на четвертом шаге каждая /-я вероятностная карта несет информацию об априорной вероятности принадлежности к типу (классу) поверхности со,. Отметим, что в байесовском решающем правиле, так же как и на четвертом шаге алгоритма, используется либо традиционная параметрическая оценка условной плотности распределения (УПР) признаков, либо непараметрическая оценка УПР в зависимости от результатов проверки гипотезы о согласованности с нормальным распределением.

В настоящее время наиболее распространенными являются два способа непараметрической оценки УПР. Один способ основан на оценке УПР типа Ро-зенблатта-Парзена (УПР РП)

х,. -х,,

С,.

1 ,...№,

(3)

где п, — количество наблюдений в обучающей выборке V, типа ©,; Р — количество признаковых каналов изображения; с/— сглаживающие параметры обучающей выборки V класса ю,; Ф(и) - ядерная (гауссова) функция оценки плотности распределения.

Другой способ основан на УПР по методу к-ю ближайшего соседа (для краткости будем приводить англоязычную аббревиатуру названия метода — к-КЫ), которая определяется исходя из выражения

где к„ — параметр близости соседа, N— величина выборки, \(к„, М,Х) — объем множества всех точек обучающей выборки, расстояние которых до точки Хв п-мерном пространстве меньше или равно значению Я". В случае использования евклидова расстояния функция' V(к,,, ДО, X) вычисляется из выражения

(5)

где Г — гамма-функция, А — единичная матрица.

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

Повышение вычислительной эффективности оценки УПР по методу £-ЫЫ может быть получено за счет более быстрого вычисления расстояний Я в выражении (5), что достигается путем использования методов пространственного индексирования.

Предложены способ индексирования многомерного пространства, учитывающий несколько вариантов обхода этого пространства с помощью 2-кривой

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

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

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

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

• определение рациональных параметров модифицированного алгоритма оценки УПР по методу включая вид индексирования многомерного пространства и количество обходов признакового пространства;

• поиск пределов применимости различных подходов к непараметрической оценке УПР;

• определение рационального способа преобразования ВК большой размерности;

• определение эффективности статистического и нейросетевого метода классификации в целом.

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

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

нии со статистическим методом. В то же время показано, что более простой способ формирования ВК в нейросетевом методе ведет к снижению точности на 3%-5%. Доверительные интервалы средней ошибки распознавания для статистического и нейросетевого метода (уровень доверия 95%, размер каждой выборки 30 наблюдений) оказались более широкими для нейросетевого метода, что также подтверждает большую надежность и точность результатов, доставляемых статистическим методом На основании результатов исследований предложен комбинированный алгоритм непараметрической оценки УПР, основанный на модифицированном алгоритме оценки УПР по методу к-ММ и известном алгоритме оценки УПР РП Он обеспечивает более эффективную классификацию в пространстве признаков большой и малой размерности.

0 ЯО 400 ОФО >00 1000 1200 Ш 1100 • МО 400 000 000 1ЙОО 1100 1400 1000

раммр выборки, кикс разиар выбора* пик

Рис 1 Эффективность алгоритмов классификации, составляющих основу метода статистической классификации, в зависимости ог размера обучающей выборки 1 — байесовская классификация с параметрической оценкой УПР, 2 байесовская классификация с непараметрической оценкой УПР по методу без индексирования 3 — байесовская классификация с модифицированным алгоритмом оценки УПР РП, 4 — байесовская классификация с модифицированным алгоритмом оценки УПР по методу А-ЫЫ

' > • . 9 <

Количество Шилов изображения Квпичкпо нищим мзобравония

Рис 2 Эффективность распознавания различных классификаторов в зависимости от числа каналов изображения 1 — байесовская классификация с использованием параметрической оценки УПР, 2 — байесовская классификация с использованием УПР РП, 3 — байесовская классификация с использованием УПР по методу к-1Ж , 4 — классификация нейросетевым методом, 5 — классификация статистическим методом

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

Моделирование изменений земной поверхности основано на использовании КА с описанным выше вероятностным подходом к формированию правил их функционирования Каждый КА представлен в виде квадратной матрицы

Мкд = [С/у] порядка d, в которой значение центрального элемента с^ (к = h = nil, здесь «/» — целочисленное деление) зависит определенным образом от значений всех элементов этой матрицы. Матрица МКд в этом случае формируется путем сканирования растровой матрицы изображения I скользящим окном размера (2d+l) х (2d+l) пикселей и сохранения, результата обработки Ckh = fi.C\\,C\2,...,Chh,-.;C,j<j) в текущей позиции центрального пикселя скользящего окна в качестве значения пикселя нового (результирующего) изображения I.

Исследование

точности

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

лг л _ ДОП

рования КА, вероятность pv не учитывалась. На рис. 3 в качестве примера представлена часть результатов этих исследований. На оси абсцисс рис. 3 показаны различные виды комбинаций четырех модельных изображений, имитирующих разновременные АИ. Так, запись 1,2->3 означает, что расчет параметров и построение прогнозной карты проводилось на основе модельных изображений №1 и №2, а точность полученной прогнозной карты оценивалась по изображению №3. Результаты показали значимость каждой вероятностной характеристики предлагаемых правил функционирования КА и более высокую точность полученных прогнозных карт при учете вероятностных характеристик чем точность карт, созданных с помощью модуля

CAMarkov, реализованного в ГИС Idiisi32.

Четвертая глава посвящена разработке ПО системы и ее апробации при решении практических задач анализа динамики земной поверхности.

Разработана структура ПО системы, укрупненная схема которой приведена на рис. 4. В состав ПО системы включено также ПО других производителей:

• программные средства предварительной обработки АИ и библиотека базовых функций работы с изображениями системы ER Mapper 5.5;

• СУБД для работы с БД1 и БД2 которые хранят обученные нейросети, АИ, результаты их обработки и интерпретации, а также дополнительную атрибутивную информацию;

• программные средства растровой ГИС Idrisi32 для проведения пространственного анализа.

Оригинальной программной частью созданной системы являются подсистемы, объединенные на рис. 4 в рамках ПО для решения прикладных задач. Человеко-машинный интерфейс системы (интерфейс пользователя) разработан в

среде Borland Delphi 6.0, а ПО оригинальных подсистем, учитывая требования к их высокой вычислительной эффективности,— в среде Microsoft Visual C++ 6.0. Взаимодействие между программными модулями, созданными с использованием разных сред и языков программирования, осуществляется посредством вызова соответствующих функций из динамически связываемых библиотек. Серверная компонента, необходимая для хранения различных данных, реализована с помощью языка Transact-SQL на основе СУБД MS SQL Server 2000, а в качестве механизма доступа к БД используется технология Microsoft ActiveX Data Objects.

Рис. 4. Укрупненная структурная схема ПО системы При реализации функций нейросетевой имитации, необходимых для реализации нейросетевого метода классификации, в качестве основы была использована известная библиотека классов с базовыми классами Neuron, Layer и Net, связанными иерархически.

Приведены результаты апробации разработанного ПО системы при решении задачи анализа динамики земной поверхности фрагмента территории Ханты-Мансийского автономного округа. Для ее решения были использованы космические снимки системы ДЭЗ Landsat ETM+ на 10 августа 1999 г., 6 сентября 2000 г. и 3 июля 2002 гг. Для каждого из них были получены соответствующие тематические карты (рис. 5а для снимка 2002 г.). На основании карт 1999 и 2000 гг. были получены прогнозные карты на 2001 и 2002 гг. (рис. 5в, карта на 2002 г.), что позволило с использованием программных средств пространственного анализа оценить изменение границ и площадей и выявить динамику изменений тематических классов земной поверхности.

Так, например, на 2000 г. пойменные участки исследуемого фрагмента земной поверхности занимали общую площадь около 257.6 км2, а к 2002 г. —

около 205.9 км2, то есть общая площадь пойменных участков уменьшилась приблизительно на 51.6 км2. В то же время мшистые болота на 2000 г. занимали общую площадь того же фрагмента земной поверхности около 1198.5 км2, а к 2002 г. — около 1258.1 км2, то есть прослеживается увеличение площади данного типа земной поверхности на 59.6 км2. Водные объекты исследуемой земной поверхности, а также границы кустарников практически не претерпели никаких изменений.

Для численной оценки точности прогнозной карты на 2002 г., полученной с помощью ПО системы (рис. 5в), было проведено ее сравнение с реальной тематической картой на 2002 г. (рис. 5а). Значение КИС в этом случае составило 65% совпадения с реальной тематической картой. Для сравнения на основе аналогичных разновременных тематических карт была получена прогнозная карта этой же территории на 2002 г. с помощью программного модуля СА Markov (рис. 56). Точность этой карты оказалась значительно более низкой и КИС составила 55%. Это подтверждает результаты исследований на модельных разновременных н^бражешж,^^ ^^

Рис. 5. Фрагменты тематических карт на 2002 г. а) — оригинальная, б) — прогнозная, полученная с помощью модуля СА Markov, в) — прогнозная, полученная с помощью системы

Описаны результаты апробации ПО системы при решении задачи тематического картирования территории Северо-Васюганского газоконденсатного месторождения ОАО «Востокгазпром» в Томской области по данным космической съемки системой ДЗЗ IRS. Из них следует, что созданное ПО системы эффективно не только при решении задач анализа динамики земной поверхности, но и при решении задачи тематического картирования. Это позволяет также сделать вывод о самостоятельной практической ценности программных модулей, обеспечивающих автоматизированную интерпретацию АИ.

Все приведенные выше результаты апробации системы с использованием данных систем ДЗЗ Landsat ETM+ и IRS подтвердили обоснованность и достоверность результатов выполненных исследований, а также показали высокую эффективность созданного математическою и программного обеспечения системы при решении прикладных задач.

В заключении приведены основные выводы и результаты диссертационной работы.

В приложение вынесен акт о внедрении полученных результатов.

Основные результаты и выводы

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

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

2. Предложена концепция создания системы анализа динамики земной поверхности. В ее основу положен многоэтапный подход, суть которого в последовательном применении более точных методов интерпретации АИ и моделирования изменений земной поверхности, основанного на использовании КА с вероятностным определением правил их функционирования. Сформулирован ряд принципов улучшенной автоматизированной интерпретации АИ и функционирования КА.

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

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

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

6. Разработана структура ПО системы анализа динамики земной поверхности и в соответствии с ней на базе предложенных методов и алгоритмов разработано ПО системы. Объем исходного кода системы составляет более 9000 строк на языках C++ и Object Pascal.

7. Проведена апробация разработанного математического и программного обеспечения системы при решении задачи анализа динамики земной поверхности фрагмента территории Ханты-Мансийского автономного округа по разновременным космическим снимкам системы ДЗЗ Landsat ETM+. Проведена также апробация системы при решении задачи тематического картирования территории Северо-Васюганского газоконденсатного месторождения ОАО «Восток-газпром» (Томская область) по данным системы ДЗЗ IRS. Результаты апробации подтвердили работоспособность и эффективность созданного математического и программного обеспечения системы.

Основные публикации по теме диссертации

1. Замятин А. В. Использование искусственных нейронных сетей при классификации многозональных аэрокосмических изображений // Материалы VI Всероссийской научно-технической конференции «Нейроинформатика — 2004», 4.2. М.: Изд-во МИФИ. —2004. — С. 239-246.

2. Замятин А.В. Классификация разновременных снимков в задачах прогнозирования изменений ландшафта земной поверхности // Сборник тезисов I Международной конференции «Земля из космоса— наиболее эффективные решения», М.: ООО «Бином-пресс». — 2004. — С. 116-117.

3. Замятин А.В., Марков Н.Г. Интеллектуальная информационная система для решения задач интерпретации объектов земной поверхности на аэрокосмических изображениях // Труды X Юбилейной Международной научно-практической конференции студентов и молодых ученых «Современные техника и технологии». Томск: Изд-во ТПУ. — 2004. — С. 140-142.

4. Замятин А.В., Марков Н.Г., Напрюшкин А.А. Хранение и поиск искусственных нейронных сетей в информационной системе обработки и интерпретации данных дистанционного зондирования Земли // Материалы V Всероссийской научно-технической конференции «Нейроинформатика— 2003», 4.2. М.: Изд-во МИФИ. — 2003. — С. 178-185.

5. Замятин А.В., Марков Н.Г., Напрюшкин А.А, Адаптивный алгоритм классификации с использованием текстурного анализа для автоматизированной интерпретации аэрокосмических изображений // Исследование Земли из космоса, 2004.— №2.— С. 32-40.

6. Марков Н.Г., Напрюшкин А.А., Замятин А.В. Применение нейронных сетей при решении задач экологического мониторинга// Межвуз. научно-технический сборник «Кибернетика и вуз», вып. 30.- Томск: Изд-во ТПУ. — 2003.—С. 76-81.

7. Markov N.G., Napryushkin A.A., Zamyatin A.V. Adaptive system for extracting information about spatial objects from aerospace imageries // Abstracts of Russian-American Workshop «Studies of Socio-Naturai Co-Evolution from Different Parts of the World», Novosibirsk, 2002, http://ivvep.ab.ru/~workshop/ (01.12.2004).

Марков Н.Г., Напрюшкин А.А., Замятин А.В. Адаптивная система для извлечения информации о пространственных объектах на аэрокосмических изображениях // Тезисы Российско-Американского семинара «Изучение природно-социального развития в различных частях света».

8. Markov N.G., Napryushkin A.A., Zamyatin A.V. Application of Neural Network Methods in RS-based Thematic Mapping // Proceedings of the 5-th AGILE conference on Geographic science, Palma de Mallorca, Spain, 2002, pp. 485-490.

Марков Н.Г., Напрюшкин А.А., Замятин А.В. Применение нейросетевых методов для тематического картирования на основе дистанционного зондирования // Труды 5-ой конференции по геоинформатике Ассоциации ГИС-лабораторий Европы, г. Пальма, Испания.

9. Markov N.G., Napryushkin A.A., Zamyatin A.V., Vertinskaya E. V. Adaptive Procedure of RS Images Classification with Use of Extended Feature Space // Proceedings of SPIE, 2003, Vol. 4885, pp. 489-500.

Марков Н.Г., Напрюшкин А.А., Замятин А.В., Вертинская Е.В. Адаптивная процедура классификации аэрокосмических изображений с использованием расширенного признакового пространства // Труды Международного общества оптического приборостроения.

10. N.G. Markov, A.A. Napryushkin, A.V. Zamyatin, E. V. Vertinskaya, Adaptive Procedure of RS Images Classification with Use of Extended Feature Space //Abstracts of 9th International Symposium on Remote Sensing 2002, Aghia Pelagia, Crete, Greece, 2002, p. 129.

Марков Н.Г., Напрюшкин А.А., Замятин А.В., Вертинская Е.В. Адаптивная процедура классификации аэрокосмических изображений с использованием расширенного признакового пространства // Тезисы 9-го международного симпозиума по дистанционному зондированию, г. Ираклион, Греция.

11. Markov N.G., Zamyatin A.V., Napryushkin A.A., Vertinskaya E.V., Advanced Thematic Mapping Approach for Forecasting Landscape Change Using GIS // Proceedings of the 6th AGILE Conference on Geographic Information Science, Lyon, France. AGILE, 2003, pp. 687-693.

Марков Н.Г., Напрюшкин А.А., Замятин А.В. Вертинская Е.В. Подход к улучшенному тематическому картированию для прогнозирования изменений ландшафта с использованием ГИС // Труды 6-ой конференции по геоинформатике Ассоциации ГИС-лабораторий Европы, г. Лион, Франция.

12. Zamyatin A.V., Markov N.G. Advanced GIS Tool for Assessment of Land Use Change // Proceedings of the 5th International Workshop on Computer Science and Information Technologies, Ufa, Russia, 2003, pp. 115-118.

Замятин А.В., Марков Н.Г. Улучшенный ГИС-инструмент для оценки изменений территорий землепользования //Труды 5-го международного семинара по компьютерным наукам и информационным технологиям, г.Уфа, Россия.

13. Zamyatin A.V., Markov N.G. Designing Forecast Thematic Maps Using Time Series Remotely Sensed Images // XX ISPRS Congress. Proceedings of The International Archives of The Photogrammetry, Remote Sensing and Spatial Information Sciences, vol. XXXV, part B, Istanbul, Turkey, 2004, pp. 492-497.

Замятин А.В., Марков Н.Г. Создание прогнозных карт с использованием разновременных аэрокосмических изображений // XX Международный конгресс. Труды международного архива фотограмметрии, дистанционного зондирования и наук о пространственно распределенной информации, г. Стамбул, Турция.

14. Zamyatin Alexandr V., Markov Nikolay G. Multi-stage Processing of Time Series Aerospace Images for Obtaining Enhanced Forecast Land Cover Maps // Proceedings of the 7th AGILE Conference on Geographic Information Science, Herak-lion, Greece, AGILE, 2004, pp. 457-463.

Замятин Александр В., Марков Николай Г. Многоэтапная обработка разновременных аэрокосмических изображений для формирования улучшенных прогнозных карт земной поверхности // Труды 7-ой конференции по геоинформатике Ассоциации ГИС-лабораторий Европы, г. Ираклион. Греция.

Подписано к печати 02 02 2005 г Тираж 100 экз Заказ №8 Бумага офсетная Услпечл 1,10 Уч-издл 1,00 Печать Ш80 Отпечатано в типографии ООО «РауШ йбХ»

Лицензия Серия ПД № 12-0092 от 03 05 2001 г , 634034, г Томск, ул Усова 7, ком 052 тел (3822) 56-44-54

951

Оглавление автор диссертации — кандидата технических наук Замятин, Александр Владимирович

ВВЕДЕНИЕ.

ГЛАВА 1. АНАЛИЗ ДИНАМИКИ ЗЕМНОЙ ПОВЕРХНОСТИ С ИСПОЛЬЗОВАНИЕМ ДАННЫХ ДИСТАНЦИОННОГО ЗОНДИРОВАНИЯ ЗЕМЛИ.

1.1. Математическое и программное обеспечение для анализа динамики земной поверхности.

1.1.1. Классификация моделей изменения земной поверхности.

1.1.2. Программное обеспечение систем для решения задач анализа динамики земной поверхности.

1.2. Проблемы анализа динамики земной поверхности с использованием данных дистанционного зондирования Земли.

1.2.1. Общие положения.

1.2.2. Анализ динамики земной поверхности.

1.2.3. Особенности использования аэрокосмических изображений.

1.3. Проблема автоматизированной интерпретации аэрокосмических изображений.

1.4. Цель и задачи исследования.

Введение 2005 год, диссертация по информатике, вычислительной технике и управлению, Замятин, Александр Владимирович

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

Все это, как правило, ведет к появлению серьезных локальных и глобальных экологических изменений, которые необходимо анализировать и оценивать, а также прогнозировать их развитие с целью комплексного исследования изучаемой территории земной поверхности и повышения эффективности использования ее природных ресурсов. Существенная пространственная неоднородность исследуемых природно-территориальных комплексов, их динамичность и большая площадь зачастую не позволяют применять такие традиционные методы слежения за их состоянием как стационарные исследования, тематическое картирование на местности и т.п. [11,51,72]. Поэтому как альтернатива таким методам за последние десятилетия развивалась концепция экосистем-ного мониторинга как системы наблюдения, регистрации, контроля динамики, прогнозирования изменений и, наконец, управления и оптимизации экосистем [11]. Для наблюдения за состоянием экосистем и их отдельных компонент с 30-х годов XX века используются аэрометоды, а с 70-х годов XX века и космические методы. В настоящее время аэрокосмические методы изучения экосистем называют методами дистанционного зондирования Земли (ДЗЗ) [21,31].

За последние десятилетия XX века отечественными учеными, такими как Берлянт A.M., Виноградов Б. В., Тикунов B.C. и другими был выполнен большой объем исследований, посвященных, в основном, теоретическим [4,11,12,14,16,47,58,72] и в меньшей степени практическим [44,46,48] аспектам решения проблем анализа динамики земной поверхности, в том числе и с использованием данных ДЗЗ.

При этом неоправданно мало внимания отечественными исследователями уделялось и уделяется автоматизации решения подобных задач, использованию постоянно возрастающих возможностей компьютерной техники. Первые практически интересные результаты при решении проблем анализа динамики земной поверхности получены зарубежными учеными, такими как Агравал Ц., Бэйкер В., Белл Э., Вербург П., ХаггетП. и многими другими [65,7375,80,81,104,116,121,130,133,137]. Однако в их работах не уделяется должного внимания этапу автоматизированной интерпретации данных ДЗЗ, являющемуся ключевым этапом при анализе динамики земной поверхности.

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

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

Для достижения поставленной цели необходимо последовательное решение следующих задач:

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

Методы исследований. В работе использованы методы многомерного статистического анализа, теории распознавания образов, нейросетевого анализа, пространственной индексации данных, объектно-ориентированного программирования. Проверка эффективности разработанных методов и алгоритмов, предназначенных для интерпретации АИ и моделирования изменений земной поверхности, а также оценка достоверности полученных при этом результатов осуществлялись путем проведения численных экспериментов с использованием модельных и реальных АИ.

2002 г.), IX Международный симпозиум "Remote Sensing 2002" SPIE, (г. Ираклион, Греция, 2002 г.), V Международный семинар «Computer Science and Information Technologies» (г. Уфа, 2003 г.), V и VI Всероссийские научно-технические конференции «Нейроинформатика — 2003» (г. Москва) и «Нейро-информатика — 2004» (г. Москва), XX Международный конгресс ISPRS (г. Стамбул, Турция, 2004 г.), X Юбилейная международная научно-практическая конференция студентов и молодых ученых «Современные техника и технологии» (г. Томск, 2004 г.), I Международная конференция «Земля из космоса — наиболее эффективные решения» (г. Москва, 2004 г.).

По результатам диссертационных исследований опубликовано 14 работ, в том числе 11 статей.

Кратко изложим основное содержание работы.

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

На основе приведенного анализа проблем формулируются цель и задачи исследования.

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

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

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

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

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

В четвертой главе рассматривается созданное ПО системы анализа динамики по данным ДЗЗ и приводятся результаты его апробации при решении практических задач анализа динамики земной поверхности и интерпретации по данным с различных спутниковых систем ДЗЗ.

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

Приводятся результаты апробации разработанного ПО для решения прикладной задачи анализа динамики земной поверхности. Для этого по разновременным космическим снимкам системы ДЗЗ Landsat (ЕТМ+) на 1999 и 2000 гг. были получены прогнозные тематические карты на 2001 и 2002 гг., отражающие динамику изменений земной поверхности фрагмента территории Ханты-Мансийского автономного округа. Описываются результаты применения разработанного ПО для решения задачи тематического картирования территории Северо-Васюганского газоконденсатного месторождения Томской области по данным системы ДЗЗ IRS. По результатам апробации сделан вывод о работоспособности и эффективности математического и программного обеспечения системы. Указывается, что разработанное ПО системы и полученные с его помощью результаты решения этой задачи внедрены в ОАО «Востокгазпром».

Научную новизну полученных в работе результатов определяют:

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

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

Практическая ценность и реализация результатов работы. Практически значимыми являются созданные методы, алгоритмы и ПО системы анализа динамики земной поверхности, включая ПО для интерпретации АИ, имеющее самостоятельную практическую ценность. Программные средства этой системы функционируют на компьютерах типа IBM PC под управлением ОС Windows 98/NT/2000. Объем исходного кода системы составляет более 9000 строк на языках С++ и Object Pascal.

Созданные программные средства были использованы при выполнении х/д №8-21/03 между Институтом «Кибернетический центр» Томского политехнического университета и ОАО «Востокгазпром» и внедрены в ОАО «Востокгаз-пром» при решении задач автоматизированного получения тематических карт территории Северо-Васюганского газоконденсатного месторождения Томской области по панхроматическим и многозональным космическим снимкам системы ДЗЗIRS. Результаты внедрения подтверждены соответствующим актом.

Личный вклад:

1. Постановка задач исследования и разработка концепции автоматизированной интерпретации и анализа динамики с использованием аэрокосмических изображений выполнены автором совместно с Н.Г. Марковым.

2. Математическое обеспечение системы автоматизированной интерпретации и анализа динамики земной поверхности с использованием аэрокосмических изображений разработано автором.

3. Постановки задач исследования эффективности предложенных алгоритмов классификации и моделирования выполнены автором совместно с Н.Г. Марковым. Результаты исследования этих алгоритмов получены лично автором.

4. Разработка программного обеспечения системы интерпретации и анализа динамки земной поверхности выполнена автором, за исключением модуля построения текстурных характеристик, разработанного совместно с Напрюш-киным A.A. Программная реализация алгоритмов параметрической оценки плотности и оценки плотности Розенблатта-Парзена выполнена Напрюшки-ным A.A.

Основные положения, выносимые на защиту:

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

Автор выражает глубокую благодарность научному руководителю доктору технических наук, профессору Н.Г. Маркову за большую помощь в подготовке диссертационной работы, ценные замечания и советы. Автор также благодарит за плодотворные дискуссии кандидата технических наук П.М. Острасть, доцентов Томского политехнического университета, кандидатов технических наук A.B. Кудинова, Р.В. Ковина, Е. А. Мирошниченко, A.A. Напрюшкина и A.B. Сарайкина. Автор выражает отдельную благодарность Югорскому НИИ информационных технологий и лично руководителю центра ДЗЗ, кандидату технических наук Копылову В.Н. за предоставленные разновременные космические снимки Landsat ЕТМ+.

Заключение диссертация на тему "Математическое и программное обеспечение системы анализа динамики земной поверхности по данным дистанционного зондирования Земли"

4.5. Основные результаты и выводы по главе

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

2. Осуществлен выбор базового ПО в виде программных модулей других производителей — это система ER Mapper 5.5, ГИС Idrisi32 и СУБД MS SQL

Server 2000. Проанализированы несколько вариантов реализации программных средств нейросетевой имитации. В качестве базового выбран вариант на основе существующей библиотеки базовых классов С. Короткого.

3. Разработан удобный интерфейс пользователя, позволяющий упростить человеко-машинное взаимодействие. Приведены особенности разработанной подсистемы хранения данных, которая обеспечивает автоматизированный поиск и хранение ранее обученных нейросетей. Созданы с использованием объектно-ориентированного подхода программные средства нейросетевой имитации. Объем исходного кода программ составляет более 9000 строк на языках С++ и Object Pascal.

4. Проведена апробация разработанного ПО при решении задачи анализа динамики земной поверхности— фрагмента территории Ханты-Мансийского автономного округа по разновременным космическим снимкам системы ДЗЗ Landsat (ЕТМ+). Получены прогнозные карты и выявлены тенденции изменений типов поверхностей исследуемой территории. В частности, по результатам пространственного анализа полученных тематических и прогнозных карт на 1999-2002 гг. сделан вывод о преимущественной замене пойменных участков на мшистые болота на всей исследуемой территории.

5. Проведена апробация разработанного ПО при решении задачи автоматизированного тематического картирования территории Северо-Васюганского га-зоконденсатного месторождения по космическим снимкам системы ДЗЗ IRS.

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

7. Результаты, полученные при решении задачи тематического картирования, а также разработанное математическое и программное обеспечение улучшенной, более точной интерпретации внедрены в ОАО «Востокгазпром», о чем получен соответствующий акт.

ЗАКЛЮЧЕНИЕ

Диссертационная работа посвящена созданию математического и программного обеспечения системы анализа динамики земной поверхности по данным дистанционного зондирования Земли (ДЗЗ).

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

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

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

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

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

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

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

7. Разработана структура ПО системы анализа динамики земной поверхности и в соответствии с ней на базе предложенных методов и алгоритмов разработано ПО системы, создан удобный пользовательский интерфейс, позволяющий упростить человеко-машинное взаимодействие. Объем исходного кода разработанных программ составляет более 9000 строк на языках С++ и Object Pascal.

8. Проведена апробация разработанного математического и программного обеспечения системы при решении задачи анализа динамики земной поверхности фрагмента территории Ханты-Мансийского автономного округа по разновременным космическим снимкам системы ДЗЗ LandsatETM+. Проведена апробация системы при решении задачи тематического картирования территории Северо-Васюганского газоконденсатного месторождения ОАО «Востокгаз-пром» по данным системы ДЗЗ 1RS. Результаты апробации подтвердили работоспособность и эффективность созданного математического и программного обеспечения системы, получен также соответствующий акт о внедрении.

Библиография Замятин, Александр Владимирович, диссертация по теме Математическое и программное обеспечение вычислительных машин, комплексов и компьютерных сетей

1. Айвазян С А. Классификация многомерных наблюдений. — М.: Изд-во Статистика, 1978. — 232 с.

2. Баруча-Рид А.Т. Элементы теории Марковских процессов и их приложения. — М.: Изд-во «Наука», 1969. — 512 с.

3. Белолипецкий В.М., Шокин Ю.И. Математическое моделирование в задачах охраны окружающей среды. — Новосибирск: Изд-во «ИНФОЛИО-пресс», 1997. — 240 с.

4. Берлянт A.M. Геоинформационное картографирование // Картография и геоинформатика. Итоги науки и техники. Сер. Картография. М.: ВИНИТИ АН СССР, 1991. — Т.14.— С. 80-117.

5. Берлянт A.M. Геоиконика. — М.: Астрея, 1996. — 208 с.

6. Бутусов О.Б. Алгоритмы текстурной классификации типов лесов на основе анализа космических снимков с ИСЗ «Landsat-7» // Исследование Земли из космоса. — 2002. — №5. — С. 87-96.

7. Буч Г. Объектно-ориентированный анализ и проектирование с примерами приложений на С++. Пер. с англ. под ред. И. Романовского и Ф. Андреева, М.: Rational, 1998. — 420 с.

8. Вапник В.Н. Восстановление зависимостей по эмпирическим данным. — М.: Изд-во «Наука», 1979. — 488 с.

9. Вапник В.Н., Стефанюк А.Р. Непараметрические методы восстановления плотности вероятности // Автоматика и телемеханика. — 1978.— № 8. — С. 38-52.

10. Ватолин Д., Ратушняк А., Смирнов М., Юкин В. Методы сжатия данных. Устройство архиваторов, сжатие изображений и видео. — М.: Изд-во «ДИАЛОГ-МИФИ», 2002. — 384 с.

11. Виноградов Б. В. Аэрокосмический мониторинг экосистем. — М.: Изд-во Наука, 1984. — 320 с.

12. Виноградов Б. В. Основы ландшафтной экологии. — М.: Изд-во «Геос», 1998. —418 с.

13. Виноградов Б. В., Кондратьев К.Я. Космические методы землеведения. — JL: Гидрометеорологическое изд-во, 1971. — 190 е.: ил.

14. Виртуальные методы дешифрирования.— М.: Изд-во Недра, 1990.—341 с.

15. Гарбук C.B., Гершензон В.Е. Космические системы дистанционного зондирования Земли. — М.: Издательство А и Б, 1997. — 296 е.: ил.

16. География из космоса / В.Т. Савиных, В.А. Малинников, С.А. Сладкопевцев, Э.Н. Цыпина. — Изд-во Московского государственного университета геодезии и картографии, 2000. — 224 с.

17. Гмурман В.Е. Теория вероятностей и математическая статистика.— М.: Изд-во Высш. шк., 2000. — 479 с.

18. Горелик А.Л., Скрипкин В.А. Методы распознавания. — 3-е изд. перераб. и доп. — М.: Изд-во Высш. шк., 1989. — 232 с.

19. Дантеман Д., Мишел Д., Тейлор Д. Программирование в среде Delphi: Пер. с англ./ — К.: НИПФ «ДиаСофт Лтд.», 1995. — 608 с.

20. Джеффри Рихтер Windows для профессионалов: Программирование для Windows 95 и Windows NT 4 на базе Win32 API/Пер. с англ. М.: Издательский отдел «Русская редакция» ТОО «Channel Trading Ltd.», 1997.— 712 с.

21. Дистанционное зондирование: количественный подход / Ш.М. Дейвис, ДА. Ландгребе, Т.Л. Филлипс и др. Пер с англ. М., Недра, 1983.— с. 415.

22. Дуда Р., Харт П. Распознавание образов. Пер. с англ. — М.: Изд-во Наука, 1981. —450 с.

23. Епанечников В.А. Непараметрическая оценка многомерной плотности вероятности // Теория вероятностей и ее применения.— 1969.— Т. 14, вып. 1. — С. 156-161.

24. Забавин А.Б. Использование искусственных нейронных сетей в задачах изучения Земли из космоса // Исследование Земли из космоса. — 2000. — №6. — С. 79-93.

25. Загоруйко Н.Г. Методы распознавания и их применение.— М.: Изд-во Наука, 1984. —256 с.

26. Замятин А.В. Использование искусственных нейронных сетей при классификации многозональных аэрокосмических изображений // Материалы VI Всероссийской научно-технической конференции «Нейроинформатика — 2004», 4.2. М.: МИФИ. — 2004. — С. 239-246.

27. Замятин А.В., Марков Н.Г., Напрюшкин А.А, Адаптивный алгоритм классификации с использованием текстурного анализа для автоматизированной интерпретации аэрокосмических изображений// Исследование Земли из космоса. — 2004.— №2.— С. 32-40.

28. Книжников Ю.Ф. Аэрокосмическое зондирование.— М.: Изд-во МГУ, 1997. — 129 с.

29. Компания Advanced Information Systems и др. Oracle8. Энциклопедия пользователя.: Пер. с англ./ К.: Издательство «ДиаСофт Лтд.», 1998. — 864 с.

30. Круглински Д., Уингоу С., Шефер Дж. Программирование на Microsoft Visual С-н- 6.0 для профессионалов / Пер. с англ. — СПб: Питер; М.: Изда-тельско-торговый дом «Русская редакция», 2002. — 864 с.

31. Лапко A.B., Ченцов C.B. Непараметрические системы обработки информации: Учебное пособие. — М.: Изд-во Наука, 2000. — 350 с.

32. Мания Г.М. Статистическое оценивание распределения вероятностей.— Тбилиси: Изд-во Тбил. ун-та, 1974. — 238 с.

33. Марков Н.Г., Напрюшкин A.A., Замятин A.B. Применение нейронных сетей при решении задач экологического мониторинга // Межвузовский на-учно-техн. сб. "Кибернетика и вуз", вып. 30- Томск: Изд-во ТПУ. — 2003.— С. 76-81.

34. Минаси М. Графический интерфейс пользователя: секреты проектирования: Пер. с англ. —М.: Изд-во Мир, 1996. —160 с.

35. Назаров JI.E. Применение многослойных нейронных сетей для классификации земных объектов на основе анализа многозональных сканерных изображений // Исследование Земли из космоса. — 2000. — №6. — С.41-50.

36. Назимова Д.И., Поликарпов Н.П. Возможен ли прогноз лесного покрова Сибири на XXI век? // Природа. — 2001. — №4. — С. 55-62.

37. Напрюшкин A.A. Алгоритмическое и программное обеспечение системы интерпретации аэрокосмических изображений для решения задач картирования ландшафтных объектов: диссертация на соискание ученой степени к.т.н. Томск. — 2002. — 168 с.

38. Наумов JI.A., Шалыто A.A. Клеточные автоматы. Реализация и эксперименты // Мир ПК. — 2003. — №8. С. 64-71.

39. Нейронные сети. Statistica Neural Networks: Пер. с англ. — М.: Горячая линия-Телеком. 2000.—182 с.

40. Нейроинформатика / А.Н. Горбань, B.JI. Дунин-Барковский, А.Н. Кир дин и др.— Новосибирск: Изд-во Наука. Сибирское предприятие РАН. — 1998. —296 с.

41. Новый метод картографического представления информации о динамике экосистем // Аэрокосмические методы исследований при мелиоративном и водохозяйственном строительстве. М.: Союзводпроект. — 1990.— С. 9299.

42. Орлов А. И. Какие гипотезы можно проверять с помощью двухвыборочно-го критерия Вилкоксона? // Заводская лаборатория. — 1999.— №1.— С. 51-55.

43. Принципы и методика составления прогнозных металлогенических карт. Сб. статей. Отв. Ред. Е.Т. Шаталов. М.: Изд-во Недра, 1966.

44. Принципы и методы геосистемного мониторинга / А.М. Грин, Н.Н. Клюев, В.Д. Утехин и др. — М.: Изд-во Наука, 989. — 168 с.

45. Проблемы создания региональных геоинформационных комплексов и опыт решения прикладных задач на основе картографической информации / И.Е. Бруни, В.Ю. Вознесенский, А.Ю. Воробьев и др. Под ред. В.В. Лебедева.— М.: Изд-во Наука, 2002. — 239 с.

46. Прэтт У. Цифровая обработка изображений: Пер. с англ.— М.: Мир, 1982. — Кн.1— 456 е.: ил.

47. Прэтт У. Цифровая обработка изображений: Пер. с англ.— М.: Мир, 1982. —Кн.2 —480 е.: ил.

48. Ревзон А.А., Камышев А.П. Природа и сооружения в критических ситуациях // Дистанционный анализ.— М.: Изд-во «Триада Лтд», 2001.—2008 с.

49. Региональный географический прогноз — управление природопользованием. Сб. статей. — М.: Изд-во Наука, 1989. — 80 с.

50. Розанов Ю.А. Стационарные случайные процессы.—2-е изд., доп.—М.: Изд-во Наука. Гл.ред.физ.-мат.лит., 1990.—272 с.

51. Саймон P. Microsoft Windows 2000 API: Энциклопедия программиста (пер. с англ.). К.: ДиаСофт, 2001. — 1088 с.

52. Скворцов А.В. Триангуляция Делоне и ее применение.— Томск: Изд-во Том. ун-та, 2002.—128 с.

53. Смирнов Н.В., Дунин-Барковский И.В. Краткий курс теории вероятностей и математической статистики для технических приложений.— М.: Изд-во Наука. Физматлит, 1969. — 436 с.

54. Тексейра С., Пачеко К. Borland Delphi 6. Руководство разработчика.: Пер.с англ.—М.: Издательский дом «Вильяме», 2002.—1120 с.

55. Теория и методы географического прогнозирования: возможности и пути: Сб. науч. тр. / Научный совет по проблемам биосферы.— М.: Наука, 1992.—151 с.

56. Тикунов B.C. Моделирование в картографии. — М.: Изд-во МГУ, 1997. — 405 с.

57. Ту Д., Гонсалес Р. Принципы распознавания образов. Пер. с англ.— М.: Изд-во Мир, 1978. — 412 с.

58. Уоссерман Ф. Нейрокомпьютерная техника.— М.: Изд-во Мир, 1992.

59. Фу К. Последовательные методы в распознавании образов и обучении машин. Пер. с англ. — М.: Изд-во Наука, 1971. — 256 с.

60. Фу К. Структурные методы в распознавании образов. Пер. с англ. — М.: Изд-во Мир, 1977. — 320 с.

61. Фукунага К. Введение в статистическую теорию распознавания образов / Пер. с англ. М.: Изд-во Наука, 1979. — 368 с.

62. Хаггет П. Пространственный анализ в экономической географии. — М.: Изд-во Прогресс, 1968.

63. Харалик P.M. Статистические и структурные подходы к описанию текстур.—ТИИЭР.—Т. 67.—№56.— 1979.

64. Харин Ю.С., Степанова М.Д. Практикум на ЭВМ по математической статистике. — Мн.: Изд-во «Университетское», 1987. — 304 с.

65. Харламов Б.П. Непрерывные полумарковские процессы.— Спб.: Изд-во Наука, 2001. —432 с.

66. Ховард Р.А. Динамическое программирование и марковские процессы. — М.: 1964.—302 с.

67. Шамис В.A. Borland С++ Builder. Программирование на С++ без проблем. М.: «Нолидж», 1997. 266 с.

68. Шапиро Е.И. Непараметрические оценки плотности вероятности в задачах обработки результатов наблюдений // Зарубежная радиоэлектроника.— 2000.—№2.—С. 3-22.

69. Экоинформатика. Теория. Практика. Методы и системы. Под ред. РАН В.Е. Соколова, С.-Пб., Изд-во Гидрометеоиздат, 1992.— 520 с.

70. Baker W.L. A review of models of landscape change // Landscape Ecology, no. 2, 1989, pp. 111-133.

71. Bell E.J. Markov analysis of Land Use Change — An Application of Stochastic Processes to Remotely Sensed Data. J. of Socieoeconomic Planning Sciences №8. 1974, pp. 311-316.

72. Berry M., Comiskey J., Minser K. Parallel Analysis of Clusters in Landscape Ecology // IEEE Computational Science and Engineering, 1994, vol. 1, no. 2, pp. 24-38.

73. Bischof H., Scheider W., Pinz A.J. Multispectral Classification of LandsatImages Using Neural Networks // IEEE Transactions on Geoscience and Remote Sensing, 1992, pp. 482-493.

74. Bishop С. M. Neural Networks for Pattern Recognition. Oxford Univ. Press, 1995.

75. Briassoulis, H. Analysis of Land Use Change: Theoretical and Modeling Approaches. The Web Book of Regional Science, Regional Research Institute, West Virginia University, 1999 — http://www.rri.wvu.edu /WebBook/ Briassoulis/ contents.htm (01.12.2004).

76. Brown D.G., Goovaerts P., Burnicki A., Li M.Y. Stochastic simulation of landcover change using geostatistics and generalized additive models // Photogram metric Engineering and Remote Sensing, 68(10), 2002, pp. 1051-1061.

77. Brown D.G., Pijanowski B.C., Duh J.D. Modeling the relationships between land use an land cover on private lands in the Upper Midwest, USA. J. of Environmental Management, 2002, vol. 59, pp. 1-17.

78. Butz A.R. Convergence with Hilbert's Space Filling Curve // J. Computer and System Sciences, 1969, vol. 3, pp. 128-146.

79. Clark Labs. Геоинформационная система Idrisi 32. — http:// www.clarklabs.org/ IdrisiSoftware/ (01.12.2004).

80. Clarke K.C., Hoppen S., Gaydos L. A self-modifying cellular automaton model of historical urbanization in the San Francisco Bay area // Environment and Planning B: Planning & Design: 24, 1997, pp. 247-261.

81. Earth Resource Mapping, ER Mapper 5.5 —http://www.ermapper.com (01.12.2004).87. eCognition. Система для объектно-ориентированной классификация аэрокосмических изображений. —http://www.definiens-imaging.com/(01.12.2004)

82. ERDAS Corp. Система обработки данных ДЗЗ ERDAS Imagine. — http:// www.erdas.com (01.12.2004).

83. ESRI Software. — http://www.esri.com (01.12.2004).

84. Faloutsos С. Gray Codes for Partial Match and Range Queries. IEEE Transactions on Software Engineering (TSE) 14(10), 1988. — pp. 1381-1393.

85. Fitz H.C., DeBellevue E.B., Costanza R., Boumans R., Maxwell Т., Wainger L., Sklar F.H. Development of a General Ecosystem Model for a Range of Scales and Ecosystems // Ecological Modeling,vol. 88, 1996, pp. 263-295.

86. Geman S., Geman D. Stochastic relaxation, Gibbs distributions and the Bayesian restoration of images // IEEE Transaction on Pattern Analysis and Machine Intelligence, vol. 6, 1984, pp.721-741.

87. Giacinto G., Roli F. Adaptive Selection of Image Classifiers // ICIAP '97, 9th International Conference on Image Analysis and Processing, Florence, Italy, Springer Verlag Ed. 1997, pp.38-45.

88. Giacinto G., Roli F. Dynamic Classifier Selection Based on Multiple Classifier Behaviour // Pattern Recognition, 34(9), 2001, pp. 179-181.

89. Giacinto G., Roli F., Bruzzone L., Combination of Neural and Statistical Algorithms for Supervised Classification of Remote-Sensing Images // Pattern Recognition Letters, 2000, vol. 21, no. 5, pp. 385-397.

90. Giacinto G., Roli F., Fumera G., Selection of Image Classifiers // Electronics Letters, 36(5), 2000, pp. 420-422.

91. GotAI.NET — Материалы — Нейронные сети http://www.gotai.net/ documents/ doc-nn-004.aspx (01.12.2004).

92. Haralick R.M., Shapiro L.G. Computer and Robot Vision, Addison-Wesley, Massachusetts, 1992.

93. Hervas J., Barredo J.I., Rosin P., Passuto A., Mantovani F., S. Silvano. Monitoring landslides from optical remotely sensed imagery: the case history of Tes-sina landslide, Italy // Geomorphology, 2003, vol. 54, pp. 63-75.

94. Hecht-Nielsen R. Kolmogorov's Mapping Neural Network Existence Theorem // IEEE First Annual Int. Conf.on Neural Networks, San Diego, 1987, vol. 3, pp. 11-13.

95. Jain A., Zongker D. Feature Selection: Evaluation, Application, and Small Sample Performance // IEEE Transactions on Pattern Analysis and Machine Intelligence, , 1997, vol. 19, no. 2pp. 153-158.

96. Jenerette G.D. and J. Wu. Analysis and Simulation of Land Use Change in the Central Arizona — Phoenix region // Landscape Ecology, 2001, vol. 16,m pp. 611-626.

97. Kaneko I.S., Igarashi S. Combining Multiple k-Nearest Neighbour Classifiers Using Feature Combinations // J. IECI, 2000, vol. 2, no. 3, pp. 23-31.

98. Katkovnik V., Shmulevich I. Nonparametric density estimation with adaptive varying window size // Conference on Image and Signal Processing for Remote Sensing VI, European Symposium on Remote Sensing, Barcelona, Spain,—• 2000.

99. Kumar A. A Study of Spatial Clustering techniques // DEXA 1994. pp. 57-71.

100. Maplnfo Corporation. GIS Maplnfo — http://www.mapinfo.com (01.12.2004).

101. Markov N.G., Napryushkin A.A., Zamyatin A.V. Adaptive system for extracting information about spatial objects from aerospace imageries // Abstracts of Russian-American Seminar «Studies of Socio-Natural Co-Evolution from Differentft

102. Parts of the World», Novosibirsk, 2002, http://iwep.ab.ru/~workshop/ (01.12.2004).

103. Markov N.G., Napryushkin A.A., Zamyatin A.V. Application of Neural Network Methods in RS-based Thematic Mapping// Proceedings of the 5-th AGILE conference on Geographic science, Palma de Mallorca, Spain, 2002, pp. 485-490.

104. Markov N.G., Napryushkin A.A., Zamyatin A.V., Vertinskaya E. V., Adaptive

105. Procedure of RS Images Classification with Use of Extended Feature Space // Proceedings of SPIE, 2003, Vol. 4885, pp. 489-500.

106. Markov N.G., Zamyatin A.V., Napryushkin A.A., Vertinskaya E.V., Advanced

107. Thematic Mapping Approach for Forecasting Landscape Change Using GIS //

108. Proceedings of the 6th AGILE Conference on Geographic Information Science, Lyon, France, AGILE, 2003, pp. 687-693.

109. Microimages Inc. Система обработки данных ДЗЗ TNT MIPS. — http://www.microimages.com (01.12.2004).

110. Microsoft Windows 2000 Server. Учебный курс MSCE: Пер. с англ. — 2-е изд., перераб.—М.: Издательско-торговый дом «Русская редакция»,2001.—912 с.

111. Muller М. R., Middleton J. A Markov model of land-use change dynamics in the Niagara Region, Ontario, Canada// Landscape Ecology, 1994, vol. 9 no.2, pp. 151—157.

112. Official GRASS GIS homepage. ГИС GRASS —http://grass.itc.it/.(01.12.2004)

113. Omachi S., H. Aso. A Fast Algorithm for k-NN Classifier Based on Branch and Bound Method and Computational Quantity Estimation // Transaction IEICE DII, J82-D-II, 4, 1999, pp. 641-649.

114. O'Sullivan D., P. M. Torrens. Cellular models of urban systems// Theoretical and Practical Issues on Cellular Automata. London: Springer-Verlag, 2000, pp. 108-117.

115. Park S., Wagner D. F. Incorporating Cellular Automata simulators as analytical engines in GIS // Transaction in GIS, 1997, vol. 2, no. 3, pp.213-231.

116. Pijanowski, B.C., D. Brown, B. Shellito and G. Manik. Using neural networks and GIS to forecast land use changes: A Land Transformation Model //Computers, Environment and Urban Systems, 2004, vol. 26, pp. 553-576.

117. Podest E., Saatchi S. Application of Multiscale Texture in Classifying JERS-1 Radar Data Over Tropical Vegetation. International Journal of Remote Sensing,2002, vol. 23, no 7, pp. 1487-1506.

118. Purdue/LARS MultiSpec—http://dynamo.ecn.purdue.edu/~biehl/MultiSpec/ (01.12.2004).

119. Raymer M.L., Punch W.F., Goodman E.D., Kuhn L.A., and Jain L.C. Dimensionality reduction using genetic algorithms. IEEE Trans, on Evolutionary Computation, 4(2), 2000, pp. 164-171.

120. Richards J.A., Xiuping Jia. Remote Sensing Digital Image Analysis: An Introduction. Berlin: Springer, 1999.

121. Richardson D., van Oosterom P. Urban Simulation Using Neural Networks and Cellular Automata for Land Use Planning // Advances in Spatial Data Handling. Berlin: Springer-Verlag, 2002. pp. 451-464.

122. Scanex NeRIS. Программа нейросетевой интерпретации данных ДЗЗ. — www.scanex.ru/rus/tematica/Program/neris.htm (01.12.2004).

123. SNNS — Stuttgart Neural Network Simulator. — http://www-ra.informatik.uni-tuebingen.de/SNNS/ (01.12.2004).

124. Steele B.M., J.C. Winne, R.L. Redmond. Estimation and mapping of misclassi-fication probabilities for thematic land cover maps // Remote Sensing of the Environment, vol. 66, 1998. pp. 192-202.

125. Turner M.G. Spatial simulation of landscape changes in Georgia: a comparison of 3 transition models //Landscape Ecology. 1987, vol 1, pp. 29-36.

126. Tutorial on ER Mapper SDK (Earth Resource Mapping). ER Mapper 5.5 user's manual. —1997.

127. Veldkamp A., Verbürg P.H., et al., Modeling the Spatial Dynamics of Regional Land Use: The CLUE-S Model // Environmental Management, 2002, vol. 30, no. 3,pp. 391-405.

128. Verbürg P., Schot P., Dijst M, Veldkamp A. Land use change modelling: current practice and research priorities // Geo Journal: 2002, в печати in press.

129. Verbürg P.H. et al., A method to analyse neighborhood characteristics of land use patterns // Computers, Environment and Urban Systems. 2003, pp.354-369.

130. Verbürg P.H., Soepboer W., Limpiada R., Espaldón V., Mastura S., Veldkamp S.A. Modelling the spatial dynamics or regional land use: the CLUE-S model // Environmental management, Springer-Verlag New-York Inc. 30(3), 2002, pp. 391-405.

131. Voinov A., Costanza R., Wainger L., Boumans R., Villa F., Maxwell T., Voinov H. Patuxent Landscape Model: Integrated Ecological Economic Modeling of a Watershed// Environmental Modeling and Software, 1999, vol. 14, pp. 473-491.

132. Weng Q. Land Use Change Aanalysis in the Zhujiang Delta of China Using Satellite Remote Sensing, GIS, and Stochastic Modeling // J. of Environmental Management, 2002, vol. 64, pp. 273-284.

133. Westin L. K. Department of computer Science Umee University. Receiver operating characteristic (ROC) analysis. — http://www.cs.umu.se/ research/ reports/ 2001/018/partl.pdf (01.12.2004)

134. Widrow B., Lehr M.A. 30 Years of Adaptive Neural Networks: Perceptron, Madaline, and Backpropagation // Proceedings of the IEEE, 1990, vol. 78, no. 9, p. 1415-1442.

135. Zadeh L.A. Fuzzy Sets // Information control, 1965, pp. 338-353.

136. Zamyatin A.V., Markov N.G. Advanced GIS Tool for Assessment of Land Use Change // Proceedings of the 5th International Workshop on Computer Science and Information Technologies, Ufa, Russia, 2003, pp. 115-118.

Похожие работы

Информатика, вычислительная техника и управление
05.13.00