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

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
мНет описания правки
Нет описания правки
Строка 110: Строка 110:
В продуктовой линейке нас интересует продукты Basic Imagery и Standard Imagery. По поводу возможности ортокоррекции продуктов Basic Imagery не скажу точно (вроде в поставке имеются необходимые метаданные), а вот Standard Imagery наиболее подходящий вариант.
В продуктовой линейке нас интересует продукты Basic Imagery и Standard Imagery. По поводу возможности ортокоррекции продуктов Basic Imagery не скажу точно (вроде в поставке имеются необходимые метаданные), а вот Standard Imagery наиболее подходящий вариант.
Комплект поставки продукции DigitalGlobe включает в себя (на примере [http://www.digitalglobe.com/downloads/product_samples/ortho/Ortho_Ready_Stnd_Bundle_16bit.zip QuickBird: Ortho ready standard Bundle 16bit]):
Комплект поставки продукции DigitalGlobe включает в себя (на примере [http://www.digitalglobe.com/downloads/product_samples/ortho/Ortho_Ready_Stnd_Bundle_16bit.zip 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
Как вы видите структура директорий и состав файлов отличаются от продуктов компании GeoEye.
== RPC ==
== RPC ==


Версия от 19:46, 29 января 2012

Эта страница является черновиком статьи.


Введение

Начнем с определения, что же такое ортокоррекция.
Ортотрансформирование (ортокоррекция) изображения (снимка) – математически строгое преобразование исходного изображения (снимка) в ортогональную проекцию и устранение искажений, вызванных рельефом, условиями съемки и типом камеры.[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

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

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=1741