Геоинформационные технологии в системах территориального управления

Содержание

Слайд 2

МЕТОДЫ ТЕМАТИЧЕСКОГО ДЕШИФРИРОВАНИЯ

МЕТОДЫ ТЕМАТИЧЕСКОГО ДЕШИФРИРОВАНИЯ

Слайд 3

СХЕМА ИНФОРМАЦИОННЫХ ПОТОКОВ в АСУРТ

СХЕМА ИНФОРМАЦИОННЫХ ПОТОКОВ в АСУРТ

Слайд 4

ТЕХНОЛОГИЧЕСКАЯ СХЕМА КОСМИЧЕСКОГО МОНИТОРИНГА

ТЕХНОЛОГИЧЕСКАЯ СХЕМА КОСМИЧЕСКОГО МОНИТОРИНГА

Слайд 5

Данные ДЗЗ являются косвенными наблюдениями, так как регистрируются не значения свойств

Данные ДЗЗ являются косвенными наблюдениями, так как регистрируются не значения свойств

реальных объектов, а интенсивность собственного или отраженного электромагнитного излучения участков подстилающей поверхности. Чтобы получить "объективную" картину НМ космический снимок нужно правильно обработать, т.е. перейти от яркостных характеристик поверхности к свойствам объектов.
Для определения соответствия между изображением объектов на космических снимках и их наземными образами в
СКМ подсистема наблюдения разделена на две части:
космический сегмент включает КА ДЗЗ с бортовой целевой аппаратурой;
подспутниковые полигоны на Земле.
Слайд 6

ТЕХНОЛОГИЧЕСКАЯ СХЕМА ОБРАБОТКИ ДАННЫХ ДЗЗ В СИСТЕМАХ КОСМИЧЕСКОГО МОНИТОРИНГА Чтобы перейти

ТЕХНОЛОГИЧЕСКАЯ СХЕМА ОБРАБОТКИ ДАННЫХ ДЗЗ В СИСТЕМАХ КОСМИЧЕСКОГО МОНИТОРИНГА

Чтобы перейти от

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

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

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

Слайд 8

ПРЕДОБРАБОТКА Радиометрическая калибровка аппаратных погрешностей Геометрическая коррекция изображений с учетом суточного

ПРЕДОБРАБОТКА

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

и кривизны Земли, перемещения КА по орбите в процессе проведения съемки, нестабильности ориентации КА во время съемки, неравномерности расположения сенсоров, и т.п.
Слайд 9

ПРИВЕДЕНИЕ ИЗОБРАЖЕНИЙ К ОДИНАКОВЫМ УСЛОВИЯМ СЪЕМКИ Атмосферная коррекция - рассеяние и

ПРИВЕДЕНИЕ ИЗОБРАЖЕНИЙ К ОДИНАКОВЫМ УСЛОВИЯМ СЪЕМКИ Атмосферная коррекция

- рассеяние и поглощение излучения

в атмосфере
- дополнительная «подсветка»
Слайд 10

ПРИВЕДЕНИЕ ИЗОБРАЖЕНИЙ К ОДИНАКОВЫМ УСЛОВИЯМ СЪЕМКИ (нормализация яркости) Тестовые полигоны –

ПРИВЕДЕНИЕ ИЗОБРАЖЕНИЙ К ОДИНАКОВЫМ УСЛОВИЯМ СЪЕМКИ (нормализация яркости)

Тестовые полигоны – участки, для

которых известны яркостные характеристики
Подспутниковые радиометрические измерения
Слайд 11

ГЕОМЕТРИЧЕСКАЯ КАЛИБРОВКА И ПРОСТРАНСТВЕННОЕ СОВМЕЩЕНИЕ СНИМКОВ видимый диапазон Географическая привязка: -

ГЕОМЕТРИЧЕСКАЯ КАЛИБРОВКА И ПРОСТРАНСТВЕННОЕ СОВМЕЩЕНИЕ СНИМКОВ

видимый диапазон

Географическая привязка:
- уровень 1

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

ТЕМАТИЧЕСКОЕ ДЕШИФРИРОВАНИЕ Процедуру построения карты территории Т по космическому снимку называют

ТЕМАТИЧЕСКОЕ ДЕШИФРИРОВАНИЕ

Процедуру построения карты территории Т по космическому снимку называют тематическим

дешифрированием.
Следует различать два вида задач тематического дешифрирования:
- задача объектного дешифрирования - отнесение изображений наземных объектов на снимках к одному из заданных классов;
- задача параметрического дешифрирования - определение количественных значений свойств наземных объектов по космическим изображениям.
Для корректного решения задач тематического дешифрирования необходимы соответствующие наземные данные
Слайд 13

ОБЪЕКТНОЕ ДЕШИФРИРОВАНИЕ Формально задачу объектного дешифрирования можно определить как задачу перехода

ОБЪЕКТНОЕ ДЕШИФРИРОВАНИЕ

Формально задачу объектного дешифрирования можно определить как задачу перехода от

изображения территории Т, заданного в виде матрицы “элементарных” пикселей F (X,Y), где X =(x1, … , xN), Y= (y1, … ,yM), обладающих размерами (δxi, δyj) и яркостью F(xi, yj), к описанию территории Т в виде совокупности "объектов" {Тk : k=1,…, K} различных классов с размерами (Δxk, Δyk) и отношениями R (Tk1, Tk2).


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

Слайд 14

Для корректного решения задачи объектного дешифрирования необходимо предварительно построить классификатор объектов

Для корректного решения задачи объектного дешифрирования необходимо предварительно построить классификатор объектов

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

Классификатор объектов определяется целями мониторинга. Следует использовать известные классификаторы, такие как LCCS (Land Cover Classification System).

Слайд 15

В результате объектного тематического дешифрирования формируется объектная карта территории (ОКТ )

В результате объектного тематического дешифрирования формируется объектная карта территории (ОКТ )

А

- исходный космический снимок (RGB-синтез), В – объектная карта территории (ОКТ), С – легенда ОКТ

А

В

С

Слайд 16

Возможны два варианта постановки этой задачи: - вариант «объект – пиксель»,

Возможны два варианта постановки этой задачи:
- вариант «объект – пиксель», когда

размеры пикселя изображения соответствуют размерам наземных объектов, т.е. (Δxk, Δyk) = (δxi, δyj) и К=M*N. В этом случае задача дешифрирования сводится к классической задаче распознавания образов. В результате ее решения каждому пикселю изображения присваивается номер соответствующего класса k;
- вариант «объект – много пикселей», когда пространственные размеры наземных объектов существенно больше пикселя, т.е. (Δxk, Δyk) >> (δxi, δyj) и К << M*N. В этом случае задача дешифрирования объединяется с задачей районирования.
Слайд 17

Постановка задачи объектного дешифрирования предполагает, что в пределах территории Т, выбраны

Постановка задачи объектного дешифрирования предполагает, что в пределах территории Т, выбраны

типичные участки, так называемые, тестовые полигоны {ΩTl: l =1,…,L}, внутри которых проведены наземные наблюдения и выделены объекты разных классов, т.е. построен классификатор объектов и определены решающие правила их "прямого" распознавания. Используя эти данные, нужно провести распознавание объектов на территории Т по космическим снимкам.
Важно, что суммарная площадь тестовых полигонов существенно меньше площади исходной территории, т.е. {∑ΩTl << S (Т)}. На практике это соотношение обычно изменяется от 5 до 20 % в зависимости от степени изученности и неоднородности территории Т.
С учетом того, что затраты на наземные измерения на несколько порядков выше затрат на космическую съемку, процедура тематического дешифрирования позволяет заменить «прямой» способ картирования объектов («точный», но «дорогой») - на «косвенный» («дешевый», но «приближенный»).
Слайд 18

РЕШЕНИЕ ЗАДАЧИ ОБЪЕКТНОГО ДЕШИФРИРОВАНИЯ : 1. Обучение 2. Тестирование (Экзамен) Экстраполирование

РЕШЕНИЕ ЗАДАЧИ ОБЪЕКТНОГО ДЕШИФРИРОВАНИЯ :

1. Обучение
2. Тестирование (Экзамен)
Экстраполирование


Слайд 19

ОБУЧЕНИЕ Обучение предусматривает определение критериев и правил распознавания объектов по их

ОБУЧЕНИЕ

Обучение предусматривает определение критериев и правил распознавания объектов по их изображениям

с использованием наземных наблюдений на тестовых полигонах.
Результаты наземных наблюдений делятся на две части: материал «обучения» и материал «экзамена». Материал «обучения» используется для выявления соответствия между «спектральными образами» на снимках и наземными «образами» реальных объектов и построения решающих правил.
Слайд 20

Сложность решения задачи объектного дешифрирования определяется сложностью построения решающих правил, которая

Сложность решения задачи объектного дешифрирования определяется сложностью построения решающих правил, которая

в свою очередь зависит от особенностей объектов распознавания, точнее - от специфики их «образов».
Для распознавания статических образов (СО) применяются статические решающие правила (СРП), а для распознавания динамических образов (ДО) – динамические решающие правила (ДРП).
СРП могут быть как простыми (однокритериальными) так и сложными (многокритериальными).
ДО можно представить в виде упорядоченного ряда СО, каждый из которых отвечает определенному периоду жизненного цикла ДО. Соответственно ДРП представляет собой упорядоченную последовательность СРП, каждое из которых применяется в определенный период времени.
Слайд 21

ПРИМЕРЫ РЕШАЮЩИХ ПРАВИЛ ::= Пример : Red = 0.73; Пороговый критерий:

ПРИМЕРЫ РЕШАЮЩИХ ПРАВИЛ

<простой критерий > ::= <имя свойства > < арифметическое

условие > <значение >
Пример : Red = 0.73;
Пороговый критерий: 0.23<= NDVI <= 0.34
<СРП > ::= <простой критерий > [ <логическая связка > <простой критерий > ]
Пример: Red >= 0.73 & 0.23<= NDVI <= 0.34
<ДРП > ::= {[ < период - 1> ] : (< СРП-1 >) ; . . .; [ < период - n> ] : (< СРП-n>)}
<период > ::= <дата начала – дата конца >
Пример: [01.05. – 10.06 ] : (0.13<= NDVI <= 0.23); [10.06. – 10.08 ] : ( 0.23<= NDVI <= 0.74); [10.08. – 10.09 ] : (0.13<= NDVI <= 0.34)
Слайд 22

Объектное дешифрирование без обучения На практике при объектном дешифрировании космических снимков

Объектное дешифрирование без обучения

На практике при объектном дешифрировании космических снимков часто

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

Корректность применения процедур распознавания без обучения в системах космического мониторинга вызывает серьезные сомнения, хотя результаты ее решения могут быть полезны для построения классификатора объектов, формирования «образов» и уточнения мест расположения тестовых полигонов.

Слайд 23

В настоящее время большое число алгоритмов классификации и распознавания образов, включая

В настоящее время большое число алгоритмов классификации и распознавания образов, включая

и нейросетевые, встроены в стандартные пакеты обработки космических снимков, такие как Erdas Imagine, PCI, Er Mapper, ENVI, Imagine Pro и др. При этом
результаты дешифрирования представляются в форматах совместимых с форматами наиболее распространенных ГИС-оболочек (ARC GIS, MapInfo и др.).
Слайд 24

Матрица погрешностей и ошибки распознавания Построенные правила применяются для контрольного распознавания

Матрица погрешностей и ошибки распознавания

Построенные правила применяются для контрольного распознавания объектов

из материала «экзамена» и оценивается качество распознавания.
Слайд 25

ОШИБКИ РАСПОЗНАВАНИЯ Еij – объект принадлежащий классу j при распознавании отнесен

ОШИБКИ РАСПОЗНАВАНИЯ
Еij – объект принадлежащий классу j при распознавании отнесен к

классу i
Точность распознавания объектов классу j
Δj= Еjj / Σi Еij
Вероятность ошибки первого рода ( пропуски объекта класса j )
Λ1j = (Σi Еij – Еjj) / Σi Еij
Вероятность ошибки второго рода («ложные» объекты отнесены к классу j)
Λ2j = ( Σj Еkj - Еjj ) / (ΣjΣi Еij - Σi Еij )
При суммировании в числителе индекс k обозначает строку,соответствующую классу j
Риск потерь определяется двумя компонентами: вероятностью ошибок распознавания и ценой ошибки или величиной "штрафа" за ошибку.
Риск «потерь» для объектов класса j
Rj = Λ1j х S1j + Λ2j x S2j
где S1j и S2j цена ошибок первого и второго рода для объектов класса j
Суммарный риск потерь для всей матрицы
R = Σj R j
Слайд 26

Пример матрицы погрешностей и ошибки распознавания

Пример матрицы погрешностей и ошибки распознавания

Слайд 27

ОЦЕНКА КАЧЕСТВА РАСПОЗНАВАНИЯ «Обучения» продолжается до тех пор, пока на «экзамене»

ОЦЕНКА КАЧЕСТВА РАСПОЗНАВАНИЯ

«Обучения» продолжается до тех пор, пока на

«экзамене» не будут минимизированы вероятности ошибок распознавания 1 и 2-го рода (или риски потерь) для всех классов объектов.
Рекомендуется:
В процессе обучения периодически «перемешивать» материал «обучения» и материал «экзамена».
.
Слайд 28

ЭКСТРАПОЛИРОВАНИЕ. Правила, построенные для объектов тестовых полигонов, применяются для распознавания и

ЭКСТРАПОЛИРОВАНИЕ.

Правила, построенные для объектов тестовых полигонов, применяются для распознавания и картирования

объектов на остальной территории Т.
Рекомендуется:
Периодически проводить дополнительные исследования по верификации результатов дешифрирования за пределами тестовых участков.
Слайд 29

Параметрическое дешифрирование - перехода от изображения территории Т, заданного в виде

Параметрическое дешифрирование - перехода от изображения территории Т, заданного в виде

“элементарных” пикселей F (X,Y), где X =(x1, … , xN), Y= (y1, … ,yM), обладающих размерами (δxi, δyj) и яркостью F(xi, yj), к представлению территории Т в виде матрицы значений искомого параметра G (X,Y), где X =(x1, … , xN), Y= (y1, … ,yM).
Примерами параметров: влажность почв, содержание гумуса в почве, рельеф, температура поверхности Земли, объем надземной биомассы и т.п.
Важно, что результат параметрического дешифрирования представляется в сильных шкалах.
Если для решения задачи объектного дешифрирования необходимо предварительно построить классификатор объектов, то для корректного решения задачи параметрического дешифрирования нужна шкала градуировки значений параметра, построенная по контактным измерениям.
Слайд 30

РЕШЕНИЕ ЗАДАЧИ ПАРАМЕТРИЧЕСКОГО ДЕШИФРИРОВАНИЯ 1. Построение калибровочной функции R(F,G), устанавливающую соответствие

РЕШЕНИЕ ЗАДАЧИ ПАРАМЕТРИЧЕСКОГО ДЕШИФРИРОВАНИЯ

1. Построение калибровочной функции R(F,G), устанавливающую соответствие между

значениями яркостей F и численными значениями параметра G.
Нужно иметь снимки, на которых представлены все градации значений параметра, включая пиксели с минимальным и максимально возможным значением.
2. Тестирование калибровочной функции
По аналогии с объектным дешифрированием данные контактных измерений следует разделить на две части: материал «обучения» и материал «экзамена». Материал «обучения» используется для построения калибровочной функции, а материал «экзамена» для оценки точности калибровки.
3. Применение калибровочной функции для определения значений параметра
Примечание. В ряде случаев вид калибровочной функции может зависеть от периода съемки (например: зима, весна, лето, осень). Кроме того для больших и неоднородных территорий построить единую калибровочную функцию часто не удается. В этих случаях следует провести предварительное районирование территории , разделить ее на однородные «районы» и строить свои калибровочные функций для каждого «района».
Слайд 31

Карта температуры подстилающей поверхности.

Карта температуры подстилающей поверхности.

Слайд 32

Частным случаем параметрического дешифрирования является задача построения рельефа местности по данным

Частным случаем параметрического дешифрирования является задача построения рельефа местности по данным

ДЗЗ.

а) растровое представление, б) карта изолиний (вектор)

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

Слайд 33

Пример неправильного использования результатов дешифрирования (распознавание дат сева)

Пример неправильного использования результатов дешифрирования (распознавание дат сева)

Слайд 34

ТЕХНОЛОГИЧЕСКАЯ СХЕМА ОБРАБОТКИ ДАННЫХ ДЗЗ ДЛЯ БОЛЬШИХ ТЕРРИТОРИЙ Предварительная обработка отдельных

ТЕХНОЛОГИЧЕСКАЯ СХЕМА ОБРАБОТКИ ДАННЫХ ДЗЗ
ДЛЯ БОЛЬШИХ ТЕРРИТОРИЙ

Предварительная обработка отдельных сцен
Формирование

групповых сцен (мозаик)
Тематическое дешифрирование (построение карты)
Формирование временных рядов результирующих карт
Слайд 35

ОПИСАНИЕ ТЕРРИТОРИЙ В ГИС Результаты дешифрирования могут быть преобразованы в векторные

ОПИСАНИЕ ТЕРРИТОРИЙ В ГИС Результаты дешифрирования могут быть преобразованы в векторные

карты и сохранены в ГИС . При этом дешифрированные объекты представляют собой тематические слои ГИС.
Слайд 36

БД ГИС- ТЕРРИТОРИИ БД ГИС- территории представляет собой временной ряд ОКТ

БД ГИС- ТЕРРИТОРИИ

БД ГИС- территории представляет собой временной ряд ОКТ

описывающих состояния территории, начиная с некоторого стартового момента t 0.
Актуализация БД ГИС осуществляется на основе данных мониторинга с определенным временным шагом δt.
Таким образом, содержание БД ГИС- схематично можно представить в виде
{ОКТ (t 0), ОКТ (t 0 +δt), … ОКТ (t n - δt), ОКТ (t n), … }.
Слайд 37

Временной ряд значений декадных VCI в вегетационном сезоне 2010 года май

Временной ряд значений декадных VCI в вегетационном сезоне 2010 года

май

июнь

июль

август

сентябрь

1 декада

3

декада

2 декада

Слайд 38

Динамика декадных значений VCI за 2012 г.

Динамика декадных значений VCI за 2012 г.

Слайд 39

Динамика IVСI территории Казахстана за период 2000-2013гг

Динамика IVСI территории Казахстана за период 2000-2013гг

Слайд 40

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

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

обработки, включая, распознавание изменений, анализ и диагностика тенденций, прогноз развития и выработка рекомендаций.
Слайд 41

Распознавание изменений в городской застройке

Распознавание изменений в городской застройке

Слайд 42

Картирование изменений Аман-Карагайского лесного массива на территории Костанайской области за 2005-2009

Картирование изменений Аман-Карагайского лесного массива на территории Костанайской области за 2005-2009

гг. (данные MODIS – разрешение 250 м, дневные пролеты, зимний период)

18 января 2005

MODIS

Слайд 43

Изменение площади посевных земель и паров на территории Костанайской области

Изменение площади посевных земель и паров на территории Костанайской области

Слайд 44

А – система озер Тоболо-Убаганского водороздела; В – система озер Тургайской

А – система озер Тоболо-Убаганского водороздела; В – система озер Тургайской ложбины;

С – система озер Сыпсынагашской ложбины

Пример сравнительного анализа состояния озерных систем (июль 1986 – июнь 2001)

Динамика изменения площади озер

Уменьшение

Увеличение

Слайд 45

Геопортал ИТЦ Сканэкс http://www.kosmosnimki.ru/ http://osm.kosmosnimki.ru/ https://scihub.copernicus.eu/

Геопортал ИТЦ Сканэкс http://www.kosmosnimki.ru/ http://osm.kosmosnimki.ru/

https://scihub.copernicus.eu/