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

Материал из GIS-Lab
Перейти к навигации Перейти к поиску
Нет описания правки
мНет описания правки
Строка 7: Строка 7:
На самом деле '''orthorectification''' — ОИ  ортотрансформирование, ортоисправление
На самом деле '''orthorectification''' — ОИ  ортотрансформирование, ортоисправление
ортокоррекция (orthocorrection) с трансформированием изображения в заданную проекцию<ref name="orthorectification">[http://www.dataplus.ru/dict/show.asp?ID=1566 Англо-русский толковый словарь по геоинформатике]. Андрианов В. Ю. - М.: ДАТА+, 2001. - 122 с. </ref>
ортокоррекция (orthocorrection) с трансформированием изображения в заданную проекцию<ref name="orthorectification">[http://www.dataplus.ru/dict/show.asp?ID=1566 Англо-русский толковый словарь по геоинформатике]. Андрианов В. Ю. - М.: ДАТА+, 2001. - 122 с. </ref>
Существует еще и такое определение. '''Ортотрансформирование, орторектификация''' (orthorectification, ortho-transformation, orthofototransformation) — устранение на изображении геометрических искажений, вызванных рельефом, для создания ортофото-снимков, ортофотокарт, ортофотопланов и др. ортотрансформированных (орторектифицированных) изображений и продуктов<ref name="орторектификация">[http://edu-knigi.ru/tikunov/geoinformatika.php?id=140 Словарь терминов по геоинформатике]. Геоинформатика. Тикунов В. С. - М.: Академия, 2005. — 480 с</ref>.


Что такое ортотрансформация? Это процесс геометрической коррекции изображений, при котором устраняются перспективные искажения, развороты, искажения вызванные дисторсией объектива и другие. Изображение при этом приводится к плановой проекции, то есть такой при которой каждая точка местности наблюдается строго вертикально, в надир. Чтобы выполнить такое преобразование необходимо устранить искажения, вызванные рельефом. Следовательно, для трансформации нужна модель рельефа, нужно знать высоту местности для каждой точки снимка.<ref name="Ортотрансформирование">[http://www.vinmap.net/?act=google_maps Google maps и другие. Использование картографических веб-сервисов в Digitals]. Бондарец А., НПП "Геосистема", Винница</ref>
Что такое ортотрансформация? Это процесс геометрической коррекции изображений, при котором устраняются перспективные искажения, развороты, искажения вызванные дисторсией объектива и другие. Изображение при этом приводится к плановой проекции, то есть такой при которой каждая точка местности наблюдается строго вертикально, в надир. Чтобы выполнить такое преобразование необходимо устранить искажения, вызванные рельефом. Следовательно, для трансформации нужна модель рельефа, нужно знать высоту местности для каждой точки снимка.<ref name="Ортотрансформирование">[http://www.vinmap.net/?act=google_maps Google maps и другие. Использование картографических веб-сервисов в Digitals]. Бондарец А., НПП "Геосистема", Винница</ref>
Строка 60: Строка 62:


Как можно заметить все файлы начинаются на '''po_344780''' - это номер заказа по которому подготовлен это рабочий набор.
Как можно заметить все файлы начинаются на '''po_344780''' - это номер заказа по которому подготовлен это рабочий набор.
Если посмотреть на содержимое каталога с использованием [http://code.google.com/p/wxgis/downloads/list wxGIS Catalog] то можно увидеть что в комплекте поставки имеются три шейп-файла (район заказа, проекции изображений - я так понимаю до и после предварительной обработки, во всяком случае, обычно они совпадают), текстовый файл метаданных, собственно изображение в формате TIFF, overview в формате JPEG с привязкой и файл с RPC данными.
Если посмотреть на содержимое каталога с использованием [[wxGIS Catalog]] (последнюю версию можно скачать по [http://code.google.com/p/wxgis/downloads/list ссылке]), то можно увидеть что в комплекте поставки имеются три шейп-файла (район заказа, проекции изображений - я так понимаю до и после предварительной обработки, во всяком случае, обычно они совпадают), текстовый файл метаданных, собственно изображение в формате TIFF, overview в формате JPEG с привязкой и файл с RPC данными.


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


=== DigitalGlobe ===
=== DigitalGlobe ===
Продуктовая линейка DigitalGlobe включает в себя<ref name="DigitalGlobe_products">[http://digitalglobe.com/products#imagery Products], © 2012 DIGITALGLOBE, INC.</ref>:
* [http://www.digitalglobe.com/downloads/BasicImagery-DS-BASIC-Web.pdf Basic Imagery]
* [http://www.digitalglobe.com/downloads/StandardImagery-DS-STAND-Web.pdf Standard Imagery]
* [http://www.digitalglobe.com/downloads/Stereo-DS-STER-Web.pdf Stereo Imagery]
* [http://www.digitalglobe.com/downloads/AdvancedOrtho-DS-AOS-Web.pdf Advanced Ortho Series]
В продуктовой линейке нас интересует продукты 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]):


== RPC ==
== RPC ==
== DEM ==
воспользоваться одной из общедоступных грубых ЦМР: [http://gis-lab.info/qa/srtm.html SRTM] (разрешение 30-90 м) и [http://gis-lab.info/qa/aster-gdem.html ASTER GDEM] (разрешение (15-90 м)


== Ортокоррекция космических снимков в wxGIS ==
== Ортокоррекция космических снимков в wxGIS ==
[[Участник:Bishop|Bishop]] Указать как ортокорректировать pan и multi
http://gis-lab.info/qa/wxgis-intro.html
http://gis-lab.info/qa/wxgis-intro.html
[[Участник:Bishop|Bishop]] Перенести статью в WiKi
[[Участник:Bishop|Bishop]] Перенести статью в WiKi
Строка 110: Строка 124:
== Литература ==
== Литература ==
<references />
<references />
== Ссылки ==
# [http://www.sovzond.net/articles/article_ortosoft.html Как провести ортотрансформирование изображений IKONOS Ortho-Kit в программном комплексе ENVI]
# [http://www.sovzond.net/articles/article_eros.html Ортотрансформирование снимков, полученных со спутника EROS-A]
# [http://www.sovzond.net/articles/article_geom.html Ортотрансформирование данных со спутника OrbView-3 в программной среде PCI Geomatica]
# [http://www.sovzond.net/articles/article_ortospot.html Использование программного комплекса ENVI для ортотрансформирования аэрофотоснимков и изображений со спутника SPOT]
# [http://www.sovzond.net/articles/article_mnogoch.html Использование коэффициентов рационального многочлена (RPCs) для ортотрансформирования спутниковых изображений IKONOS/QuickBird]
# [http://www.sovzond.net/articles/article_posib.html Возможности данных Ikonos и QuickBird для точного трехмерного позиционирования, получения ортоисправленных изображений и создания цифровых моделей поверхности (DSM)]
# [http://www.sovzond.net/articles/article_geomobr.html Геометрическая обработка данных со спутника QuickBird]
# [http://www.sovzond.net/products/ikonos.html Продукты с КА IKONOS]
# [http://www.sovzond.net/articles/article_estim.html Оценка геометрической точности данных дистанционного зондирования, полученных со спутника IKONOS]
# [http://www.sovzond.net/articles/article_stereot.html Точность стереосъемки со спутника IKONOS без наземной привязки]
# [http://www.racurs.ru/?page=300 Определение точностных характеристик снимков QuickBird]

Версия от 19:33, 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):

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