6876

Створення геоприв’язаної карти

Лабораторная работа

География, геология и геодезия

Створення геоприв'язаної карти Однією з базових задач в роботі із супутниковими знімками, є прив'язка їх до існуючої карти, чи створення нової карти на основі супутникових зображень. В даній роботі ми розглянемо процес створення геоп...

Украинкский

2013-01-08

3.87 MB

2 чел.

Створення геоприв’язаної карти

 Однією з базових задач в роботі із супутниковими знімками, є прив’язка їх до існуючої карти, чи створення нової карти на основі супутникових зображень.

В даній роботі ми розглянемо процес створення геоприв’язаної карти при наявності лише сканованої топографічної карти та супутникового зображення тієї ж місцевості.

  1.  Прив’язка до координат растрового зображення карти

В якості базової топографічної карти візьмемо «генштабівську» карту України в масштабі 1:100 000. Її можна скачати з [1]. Для приклада візьмемо квадрат М-36-050. (рис. 1.1)

Рис. 1.1. Карта України, квадрат М-36-050

Збережену карту відкриємо в ENVI:

  •  File/Open Image File, далі вказуємо шлях до файлу;
  •  В Available Band List (Список доступних спектральних знімків) завантажуємо зображення на екран – Load RGB (Рис.1.2)

Рис 1.2 Діалогове вікно при відкритті файлу.Available Band List

Рис. 1.3 Вікно Image to Map Registration

Створєюмо pts-файл, що містить координати чотирьох кутів.

  •  Map/Registration/select GCPs: Image to Map
  •  В діалоговому вікні Image to Map Registration зі списку Select Registration Projection (Вибір проекції) обираємо проекцію вихідного зображення GK Zone 6 (Pulkova 1942) та Datum (Рівнева поверхня)WGS-84. Натискаємо ОК. (Рис. 1.3)
  •  В діалоговому вікні Ground Control Points Selection (Вибір наземних контрольних точок) переключаємось в режим Lat/Lon кнопкою зі стрілочками (Рис. 1.4).
  •  Послідовно додаємо опорні точки в кутах карти, починаючи із верхнього лівого. Для вибору опорних точок необхідно обрати потрібних піксель в окні Zoom, ввести координати широти(Lat) та довготи (Lon) Add Point (Рис. 1.4)

Після того, як ці операції завершено, потрібно пересвідчитись, що помилка (RMS Error) дорівнює нулю.

Збережемо pts-файл (File/Save GCPs /Lat /Lon у тіє ж директорії та з тим же імям, що і файл карти (М-36-050.pts)

Рис.1.4 Вікно Ground Control Points Selection

Тепер виконаємо географічну прив’язку карти.

  •  Map/Registration/Warp from GCPs:Image to Map (Деформація зображення в карту за опорними точками).
  •  Далі обираємо pts-файл (M-36-050.pts).
  •  У вікні Registration Parameters  обираємо тип проекції, що задавався при визначенні опорних точок (GK Zone 6 (Pulkova 1942), Datum WGS-84) та вводимо розміри пікселя за осями X та Y(10м).
  •  Обираємо карту, яку треба деформувати.
  •  Вводимо наступні параметри деформації: Warp Method (Метод деформації)Polynomial (Поліноміальний), Degree 1, Resampling Cubic Convolution, Background 0.
  •  Вводимо ім’я та шлях, де будемо зберігати файл. В деяких версіях ENVI потрібно обирати вже існуючий файл, інакше програма вилетить із помилкою. Натиснемо ОК (Рис 1.5).

Рис.1.5 Вікно Registration Parameters

Тепер обріжемо поля карти за допомогою ROI (Області інтересу), використовуючи pts-файл. Обираємо – Basic Tools/Region of Interest/ROI Tool. В ROI_Type/ Input Points from ASCII, потрібно обрати збережений pts-файл  (M-36-050.pts), та проекцію, у який він був створений (Geographic Lat/Lon), рис.1.6.

Рис.1.6 Вікно Input ASCII File

Рис.1.7 Вікно Subset Data via ROI

Для того, щоб вивести на екран тільки зображення, що знаходиться в області інтересу, виберемо File/Subset Data via ROI у вікні ROI Tool.

Далі обираємо зображення, що потрібно обрізати, та за якою областю. Параметр Mask pixels outside of ROI? встановлюємо в Yes. Тепер вказуємо куди зберегти отриманий файл та натискаємо ОК (Рис.1.7).

Після перерахунку завантажимо отриманий файл в Available Band List (Рис. 1.8)

Рис.1.8 Геоприв’язанная карта з обрізаними полями

  1.  Створення карти із супутникового зображення

Закриємо всі файли, крім останнього та відкриємо заздалегідь підготовлений супутниковий знімок у форматі jpeg (файл Kyiv_LS7.jpg) тієї ж ділянки. Створення композитних знімків з окремих каналів буде детально розглянуто у наступних лабораторних роботах, а заразвикористаємо готовий знімок.

Зображення супутника Landsat 7 доступні для вільного завантаження у [2-6]

В Available Band List завантажуємо зображення на два екрани одночасно (карта – Display#1, знімок – Display#2 ) (Рис. 1.9)

Рис.1.9 Одночасно відкриті в ENVI знімок та карта

Створюємо pts-файл, що містить 10-20 опорних точок, рівномірно розподілених по карті:

  •  Map/Registration/select GCPs: Image to Image.У вікні вказати Base Image (Базове зображення) – Display#1 та Warp Image (Зображення для деформації) – Display#2.
  •  У вікні Ground Control Points Selection визначаємо опорні точки, що відповідають одне одному на обох зображеннях. Рекомендується встановлювати точки в місцях перехрестя доріг,на мостах, кутах будівель і т.п. 
  •  Після цього, натиснувши кнопку Show List, Пересвідчуємось, що помилка мінімальна (не перевищує одиниці). У разі необхідності, найбільш «грубі» точки треба замінити. (Рис. 1.10).
  •  Збережемо pts-файл File/Save GCPs to ASCII, з імям, що співпадає з імям вихідного непривязаного знімка.

Рис. 1.10 Вікно Ground Control Points Selection

Тепер виконаємо географічну прив’язку карти за допомогою інструмента Warp:

  •  Map/Registration/Warp from GCPs:Image to Image. Далі обираємо створений  pts-файл, зображення, що буде деформуватись (знімок) та базове зображення (карта). Потрібно ввести наступні параметри деформації: Warp Method – Polynomial, Degree – 1, Resampling – Cubic Convolution , Background – 0. Вводимо назву файлу та натискаємо ОК.

Звяжемо обидва зображення (Tools/Link/Link Displays) та перевіримо точність привязки.

Отримані карти можна зберегти у форматі GeoTIFF [7], що підтримується всіма ГІС та программами для роботи із супутниковими знімками. File/Save file as/TiFF/Geo TiFF.

Наприклад, комерційна версія Google Earth Pro дозволяє добавити власну карту на глобус автоматично, користуючись інформацією прив’язки (Рис. 1.11). В безкоштовній версії Google Earth доведеться власноруч підбирати масштаб карти та «висоту камери».

Рис. 1.11 Карта, додана до Google Earth Pro

Література та джерела

  1.  Сайт «Топографические карты Украины 1:100000».
  2.  Global Land Cover Facility. Landsat Imagery.


 

А также другие работы, которые могут Вас заинтересовать

47834. Ионнообменная очистка 408.5 KB
  Применяется при извлечении из сточных вод металлов цинка меди хрома свинца ртути кадмия а тагже соединений мышьяка фтора цианистых соединений и радиоактивных веществ. Метод позволяет рекуперировать ценные вещества при высокой степени очистки воды. Ионный обмен распространен при обменивании в процессе водоподготовки. Они не растворимы в воде.
47835. Информатика. Автоматизированные системы управления 2.41 MB
  Она опирается на математическую логику и включает такие разделы как теория алгоритмов и автоматов теория информации и теория кодирования теория формальных языков и грамматик исследование операций и другие. Этот раздел информатики использует математические методы для общего изучения процессов обработки информации. Примеры принципиальных ставших классическими решений в этой области неймановская архитектура компьютеров первых поколений шинная архитектура ЭВМ старших поколений архитектура параллельной многопроцессорной обработки...
47836. Правосознание и основные правовые категории 109.5 KB
  Властные отношения связанны с: Навязыванием воли властвующего субъекта подчинённому подданному; Jus – право Nomos – закон Автономия воли – делать должным для себя Право – возможность требовать должного отдавая каждому своё не чиня и не терпя вреда Должное – это возложенное на себя обязанность как проявляется автономия воли. можно дать следующую характеристику РФ: РФ является демократическим федеративным социальным правовым государством с республиканской формой правления. РФ – её цель служить человеку поскольку человек является...
47837. ФИНАНСОВО-ХОЗЯЙСТВЕННАЯ ДЕЯТЕЛЬНОСТЬ ПРЕДПРИЯТИЯ 1.01 MB
  Анализ как практика вид управленческой деятельности предшествующий принятию управленческих решений. С помощью анализа проводится обоснование управленческих решений на базе имеющейся информации. Анализ имеет свой предмет объекты и методы исследования экономические показатели и взаимосвязи между ними.
47838. АНАТОМИЯ, ФИЗИОЛОГИЯ, ПАТОЛОГИЯ ДЕТЕЙ С ОСНОВАМИ ГЕНЕТИКИ 857 KB
  Рост –- это увеличение размеров развивающегося организма. В понятие онтогенез включают все стадии развития организма от момента оплодотворения яйцеклетки до смерти человека. отрезки времени онтогенеза каждый из которых характеризуется своими специфическими особенностями организма функциональными биохимическими морфологическими и психологическими.
47839. Система бухгалтерского учета 101 KB
  В соответствии с Федеральным законом О бухгалтерском учете N 129ФЗ от 21 ноября 1996 года бухгалтерский учет представляет собой упорядоченную систему сбора регистрации и обобщения информации в денежном выражении об имуществе обязательствах организации и их движении путем сплошного непрерывного и документального учета всех хозяйственных операций. В соответствии с действующим Законом О бухгалтерском учете бухгалтерский учет ведут все организации находящиеся на территории РФ а также филиалы и представительства иностранных организаций....
47840. Економічна сутність і методи вимірювання продуктивності 110.5 KB
  Економічна сутність і методи вимірювання продуктивності Сутність і значення продуктивності праці Розрізняють поняття продуктивність і продуктивність праці. Продуктивність –- це ефективність використання ресурсів праці капіталу землі матеріалів енергії інформації – під час виробництва різних товарів і надання послуг. Тому для визначення ефективності виробництва найчастіше використовують показник продуктивності праці хоча це не означає що тільки праця є джерелом продуктивності. Продуктивність праці відбиває ступінь ефективності...
47842. Лекции по электробезопасности 247.5 KB
  Длительность протекания тока (ожоги тканей тела, нагрев внутренних органов, изменение состава крови, нарушение функций центральной нервной системы, вероятность совпадения времени протекания электрического тока с фазой Т кардиоцикла)