Ортокоррекция космических снимков с использованием RPC

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Эта страница является черновиком статьи.


Введение

Начнем с определения, что же такое ортокоррекция.
Ортотрансформирование (ортокоррекция) изображения (снимка) – математически строгое преобразование исходного изображения (снимка) в ортогональную проекцию и устранение искажений, вызванных рельефом, условиями съемки и типом камеры.[1]

При этом, иногда употребляют не совсем корректный термин орторектификация, который по сути является англоизмом термина orthorectification.
На самом деле orthorectification — ОИ ортотрансформирование, ортоисправление ортокоррекция (orthocorrection) с трансформированием изображения в заданную проекцию[2]

Существует еще и такое определение. Ортотрансформирование, орторектификация (orthorectification, ortho-transformation, orthofototransformation) — устранение на изображении геометрических искажений, вызванных рельефом, для создания ортофото-снимков, ортофотокарт, ортофотопланов и др. ортотрансформированных (орторектифицированных) изображений и продуктов[3].

Что такое ортотрансформация? Это процесс геометрической коррекции изображений, при котором устраняются перспективные искажения, развороты, искажения вызванные дисторсией объектива и другие. Изображение при этом приводится к плановой проекции, то есть такой при которой каждая точка местности наблюдается строго вертикально, в надир. Чтобы выполнить такое преобразование необходимо устранить искажения, вызванные рельефом. Следовательно, для трансформации нужна модель рельефа, нужно знать высоту местности для каждой точки снимка.[4]

Перераспределение пикселов на изображении в результате ортокоррекции

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

Зачем проводить ортокоррекцию космических снимков, если все и так на изображении можно различить, а в результате дополнительных операций качество ухудшится? Если целью стоит только выявление фактов по изображению и никаких требований к точности позиционирования, измерению длин и площадей не предъявляется, то проводить ортокоррекцию вовсе не обязательно. Но если нужны измерительные и позиционные свойства изображения, а также если необходимо точное совмещение разновременные изображения (или даже стыковка перекрывающихся включений), то ортокоррекция крайне рекомендуется.

Необходимые данные

Для проведения ортокоррекции космического снимка необходимо:

  • Космический снимок в одном из форматов поставки (обычно TIFF)
  • Коэффициенты RPC к нему
  • Информация о рельефе в виде DEM (Digital Elevation Model)

Рассмотрим подробнее комплекты поставки основных компаний рынка детальной космической съемки: GeoEye и DigitalGlobe (образы продукции собраны на отдельной странице).

GeoEye

Продуктовая линейка GeoEye включает в себя[5]:

  • Geo
  • GeoProfessional
  • GeoStereo

Подробнее о продуктах компании можно почитать скачав брошюру Product Guide (требует регистрацию, но потом дает прямую ссылку).

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

Комплект поставки продукции GeoEye включает в себя (на примере Transportation GeoEye-1 Sample):

Состав поставки данных с КА GeoEye-1
GeoEye_logo.tif
po_344780_aoi.dbf
po_344780_aoi.prj
po_344780_aoi.shp
po_344780_aoi.shx
po_344780_component.dbf
po_344780_component.prj
po_344780_component.shp
po_344780_component.shx
po_344780_image.dbf
po_344780_image.prj
po_344780_image.shp
po_344780_image.shx
po_344780_metadata.txt
po_344780_rgb_0000000.hdr
po_344780_rgb_0000000.tfw
po_344780_rgb_0000000.tif
po_344780_rgb_0000000_ovr.jgw
po_344780_rgb_0000000_ovr.jpg
po_344780_rgb_0000000_rpc.txt
SingleOrganization_license.txt

Как можно заметить все файлы начинаются на po_344780 - это номер заказа по которому подготовлен это рабочий набор. Если посмотреть на содержимое каталога с использованием wxGIS Catalog (последнюю версию можно скачать по ссылке), то можно увидеть что в комплекте поставки имеются три шейп-файла (район заказа, проекции изображений - я так понимаю до и после предварительной обработки, во всяком случае, обычно они совпадают), текстовый файл метаданных, собственно изображение в формате TIFF, overview в формате JPEG с привязкой и файл с RPC данными.

wxGIS Catalog может работать с ZIP архивами, как с обычными папками. Если посмотреть метаданные изображения (щелкнуть правой клавишей мыши на po_344780_rgb_0000000.tif и выбрать "Свойства"), то можно убедиться, что изображение географически привязано (имеется описание системы координат, значения экстента представляют собой географически координаты, а не пиксельные) и корректно определяется наличие RPC данных.

Диалог вывода метаданных изображения

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

Следует отметить, что программное обеспечение воспринимает изображение не как один файл в формате TIFF, а как набор из файлов, в который входят, для нашего случая, еще и файлы с расширением hdr, tfw, _rpc.txt. wxGIS Catalog корректно работает с таким набором и при переименовании будет переименовывать все связанные файлы, а при удалении, копировании, перемещении - соответственно, удалять, копировать и перемещать.

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

Аналогичный состав имеет комплект поставки материалов съемки с КА Ikonos.

А вот комплект поставки с третьего КА компании - OrbView-3, который стал свободно доступен (подробнее можно почитать на странице описания каталога данных OrbView-3), имеет ряд нюансов.

Комплект поставки продукции GeoEye с КА OrbView-3 включает в себя (на примере территории в Белоруссии): 

3v050909p0000897861a520004700712m_001631680.att
3v050909p0000897861a520004700712m_001631680.dbf
3v050909p0000897861a520004700712m_001631680.eph
3v050909p0000897861a520004700712m_001631680.jgw
3v050909p0000897861a520004700712m_001631680.jpg
3v050909p0000897861a520004700712m_001631680.prj
3v050909p0000897861a520004700712m_001631680.pvl
3v050909p0000897861a520004700712m_001631680.shp
3v050909p0000897861a520004700712m_001631680.shx
3v050909p0000897861a520004700712m_001631680.tif
3v050909p0000897861a520004700712m_001631680_aoi.dbf
3v050909p0000897861a520004700712m_001631680_aoi.prj
3v050909p0000897861a520004700712m_001631680_aoi.shp
3v050909p0000897861a520004700712m_001631680_aoi.shx
3v050909p0000897861a520004700712m_001631680_rpc.txt
3v050909p0000897861a520004700712m_001631680_src.dbf
3v050909p0000897861a520004700712m_001631680_src.prj
3v050909p0000897861a520004700712m_001631680_src.shp
3v050909p0000897861a520004700712m_001631680_src.shx
unrestricted_license.txt

Как видим названия файлов также имеют общую часть 3v050909p0000897861a520004700712m_001631680 - здесь мы видим дату съемки 050909 (2005-09-09T09:32:38.745600Z). номер в каталоге и уровень обработки 897861a (89786 и 1А), номер заказа 001631680 (OPS001631680). А вот состав файлов отличается. Но главное отличие, что поставляемые TIFF файлы не имеют географической привязки и описания СК, внедренные в файл.

wxGIS на лету подставляет необходимые данные в заголовок TIFF файла при его открытии, но это накладывает ограничение на работу с ZIP архивами - данные не будут записаны в сам файл, поэтому для ортокоррекции материалов съемки с КА OrbView-3 их необходимо предварительно распаковать. Для ортокоррекции в других программных продуктах может понадобится предварительная подготовка таких данных. Особенности ортокоррекции с использованием библиотеки GDAL изложены в отдельной статье.

DigitalGlobe

Продуктовая линейка DigitalGlobe включает в себя[6]:

В продуктовой линейке нас интересует продукты Basic Imagery и Standard Imagery. По поводу возможности ортокоррекции продуктов Basic Imagery не скажу точно (вроде в поставке имеются необходимые метаданные), а вот Standard Imagery наиболее подходящий вариант. Комплект поставки продукции DigitalGlobe включает в себя (на примере QuickBird: Ortho ready standard Bundle 16bit): 

|-005554608030_01
| |-005554608030_01_P001_MUL
| | |-05JUL04180116-M2AS-005554608030_01_P001-BROWSE.JPG
| | |-05JUL04180116-M2AS-005554608030_01_P001.IMD
| | |-05JUL04180116-M2AS-005554608030_01_P001.RPB
| | |-05JUL04180116-M2AS-005554608030_01_P001.TIF
| | |-05JUL04180116-M2AS-005554608030_01_P001.TIL
| | |-05JUL04180116-M2AS-005554608030_01_P001.XML
| | |-05JUL04180116-M2AS-005554608030_01_P001_README.TXT
| | |-DEMONSTRATION.TXT
| |-005554608030_01_P001_PAN
| | |-05JUL04180116-P2AS-005554608030_01_P001-BROWSE.JPG
| | |-05JUL04180116-P2AS-005554608030_01_P001.IMD
| | |-05JUL04180116-P2AS-005554608030_01_P001.RPB
| | |-05JUL04180116-P2AS-005554608030_01_P001.TIF
| | |-05JUL04180116-P2AS-005554608030_01_P001.TIL
| | |-05JUL04180116-P2AS-005554608030_01_P001.XML
| | |-05JUL04180116-P2AS-005554608030_01_P001_README.TXT
| | |-DEMONSTRATION.TXT
| |-GIS_FILES
| | |-005554608030_01_ORDER_SHAPE.dbf
| | |-005554608030_01_ORDER_SHAPE.shp
| | |-005554608030_01_ORDER_SHAPE.shx
| | |-005554608030_01_PRODUCT_SHAPE.dbf
| | |-005554608030_01_PRODUCT_SHAPE.shp
| | |-005554608030_01_PRODUCT_SHAPE.shx
| | |-005554608030_01_STRIP_SHAPE.dbf
| | |-005554608030_01_STRIP_SHAPE.shp
| | |-005554608030_01_STRIP_SHAPE.shx
| | |-005554608030_01_TILE_SHAPE.dbf
| | |-005554608030_01_TILE_SHAPE.shp
| | |-005554608030_01_TILE_SHAPE.shx
| |-005554608030_01_LAYOUT.JPG
| |-005554608030_01_README.TXT
| |-005554608030_01_README.XML
|-005554608030_01.MAN
|-005554608030_01_EOT.TXT

Как вы видите структура директорий и состав файлов в поставке с КА QuickBird отличаются от продуктов компании GeoEye.

В данному продукте имеются как панхроматический, так и мультиспектральный снимки. Панхроматический имеет размер пиксела 0,6 м, а мультиспектральный - 2,4 м. Каждый из снимков лежит в соответствующей папке 005554608030_01_P001_PAN и 005554608030_01_P001_MUL. Посмотрим какие данные составляют снимок: файл, который оканчивается на -BROWSE.JPG, представляет собой overview, файлы с расширениями IMD (image metadata), XML, _README.TXT содержат метаданные изображения , файл с расширением RPB содержит RPC данные, файл с расширением TIF - это само изображение, файл с расширением TIL содержит в себе схему нарезки на тайлы.

Следует обратить внимание, что если изображение разделено на тайлы, то комплект RPC данных будет предоставляться для всех файлов. Что бы корректно выполнить ортокоррекцию такого набора, необходимо выполнить мозаику тайлов в единое изображение с тем же именем как и RPB файл или работать с TIL файлом (не все ПО корректно работает с TIL файлами).

Для проверки, что комплект поставки материалов съемки содержит необходимую информацию для выполнения ортокоррекции, посмотрим метаданные изображения в wxGIS Catalog. Как видим? изображение географически привязано (имеется описание системы координат, значения экстента представляют собой географически координаты, а не пиксельные) и корректно определяется наличие RPC данных.

RPC

DEM

воспользоваться одной из общедоступных грубых ЦМР: SRTM (разрешение 30-90 м) и ASTER GDEM (разрешение (15-90 м)

Ортокоррекция космических снимков в wxGIS

Bishop Указать как ортокорректировать pan и multi http://gis-lab.info/qa/wxgis-intro.html Bishop Перенести статью в WiKi Остальные операции по ортокоррекции космических снимков в wxGIS являются аналогичными. Подробнее о

Литература

  1. ↑ Словарь терминов. Журнал Геоматика, 2009, №2. С. 119
  2. ↑ Англо-русский толковый словарь по геоинформатике. Андрианов В. Ю. - М.: ДАТА+, 2001. - 122 с.
  3. ↑ Словарь терминов по геоинформатике. Геоинформатика. Тикунов В. С. - М.: Академия, 2005. — 480 с
  4. ↑ Google maps и другие. Использование картографических веб-сервисов в Digitals. Бондарец А., НПП "Геосистема", Винница
  5. ↑ Satellite Imagery Products, Copyright © 2012 GeoEye, Inc.
  6. ↑ Products, © 2012 DIGITALGLOBE, INC.

Ссылки

  1. Как провести ортотрансформирование изображений IKONOS Ortho-Kit в программном комплексе ENVI
  2. Ортотрансформирование снимков, полученных со спутника EROS-A
  3. Ортотрансформирование данных со спутника OrbView-3 в программной среде PCI Geomatica
  4. Использование программного комплекса ENVI для ортотрансформирования аэрофотоснимков и изображений со спутника SPOT
  5. Использование коэффициентов рационального многочлена (RPCs) для ортотрансформирования спутниковых изображений IKONOS/QuickBird
  6. Возможности данных Ikonos и QuickBird для точного трехмерного позиционирования, получения ортоисправленных изображений и создания цифровых моделей поверхности (DSM)
  7. Геометрическая обработка данных со спутника QuickBird
  8. Продукты с КА IKONOS
  9. Оценка геометрической точности данных дистанционного зондирования, полученных со спутника IKONOS
  10. Точность стереосъемки со спутника IKONOS без наземной привязки
  11. Определение точностных характеристик снимков QuickBird
Источник — https://wiki.gis-lab.info/index.php?title=Ортокоррекция_космических_снимков_с_использованием_RPC&oldid=1742