орторектификация растра

[hsargsyan]

Уважаемые профессионалы!
Возникла задача орторектификации снимка, необходимо создать файл в котором будет содержаться информация RPC. В наличии есть DEM.
Помогите пожалуйста разобраться с этим файлом, что входит в его состав, как его создать, какие данные нужны для этого, и как это все реализовать в ArcGIS 10.
Надеюсь на Вашу помощь...

P.S. Никогда прежде этого не делал, для меня это темная сторона ДДЗ и ГИС

[novia]

Для начала вопрос несколько вопросов:
какой снимок?
какая ЦМР (источник, размер ячейки растра)?

Если интересно почитайте информацию по ссылкам, и ответьте на вопросы, будем вместе думать.

RPC (Rational Polynomial Coefficients или Rapid Positioning Capability), представляющие собой коэффициенты полиномов рациональных функций, связывающих координаты точек местности с координатами их изображений на растре. Значения этих коэффициентов рассчитываются поставщиком данных строгим методом, с использованием геометрической модели сенсора, бортовых и телеметрических измерений элементов внешнего ориентирования, и, в некоторых случаях, опорных точек. Поскольку RPC являются аппроксимацией строгой модели, последняя теоретически способна обеспечить получение выходных продуктов большей точности. Тем не менее, практика показывает, что обычно RPC обеспечивают достижение точности уравнивания, адекватной разрешению снимков (конечно, при условии использования качественных опорных точек в процессе вычисления RPC поставщиком данных или, что более распространено, в процессе расчета поправок в RPC при уравнивании, проводимом заказчиком снимков).

Источник: Статья на сайте Ракурс

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

Файл с коэффициентами обычно поставляют вместе со снимком.

Статья из справки к ArcGIS 10.0 докладывает, что ArcGIS понимает, когда в папке со снимком лежит файл RPC-коэффициентов и пользуется им при отображении и геообработке.

ArcGIS automatically reads RPC information in the associated text file (if it can be found next to the image file) and can apply the RPC transformation (image correction) on the fly when displaying it in ArcMap. This oriented image can also be orthorectified on the fly in ArcMap by applying elevation information. You can use a DEM to perform orthorectification or use an average elevation if no DEM is available for the same area.

Источник: Веб-справка Esri

[hsargsyan]

Добрый день!
Спасибо за инфу
Снимок IKONOS-2, есть файл в котором содержатся координаты, высота и многое другое, о котором я к сожалению не имею представления. Является ли этот файл RPC коэффициентами или нет не имею представления. Вот кусок из файла..
-------------------------------------
Product Component Metadata

Number of Components: 2

Component ID: 0000000
Product Image ID: 000
Component File Name: po_2250053_pan_0000000.tif po_2250053_red_0000000.tif po_2250053_grn_0000000.tif po_2250053_blu_0000000.tif po_2250053_nir_0000000.tif
Thumbnail File Name: po_2250053_rgb_0000000_ovr.jpg
Country Code: AM
Component Geographic Corner Coordinates
Number of Coordinates: 4
Coordinate: 1
Latitude: 40.xxxxxxxxx degrees
Longitude: 44.xxxxxxxxx degrees
Coordinate: 2
Latitude: 40.xxxxxxxxx degrees
Longitude: 44.xxxxxxxxx degrees
Coordinate: 3
Latitude: 40.xxxxxxxxx degrees
Longitude: 44.xxxxxxxxx degrees
Coordinate: 4
Latitude: 40.0xxxxxxxxx degrees
Longitude: 44.xxxxxxxxx degrees
Component Map Coordinates (in Map Units)
UL Map X (Easting): 44xxxx.xx meters
UL Map Y (Northing): 445xxxx.xx meters
Pixel Size X: 1.0000000000 meters
Pixel Size Y: 1.0000000000 meters
Product Order Map Units: meters
Columns: 13808 pixels
Rows: 21072 pixels
Percent Component Cloud Cover: 0
---------------------------------------------

P.S. некоторые цифры заменены на x.Есть также какой-то XML файл. ЦМР векторный файл изолиний с шагом в 2 метра. Можно на этой основе построить GRID - ну скажем с 2 метровыми ячейками.

[novia]

Если это иконос, то нужен файл .rpb. Если его в папке нет, запросите у поставщика снимка. Должны предоставить.
Если будет этот файл то стройте цмр(в проекции как у снимка, то есть скорее всего utm) и запускайте процесс ортотрансформирования. В тулбоксе управление данными(data management tools) - растр (raster) - обработка (processing) - орто...

Фрагмент файла это метаданные. Время съемки, номер сцены, высота солнца и тп. Это тоже полезная информация, пригодится при дешифрировании.

[Сообщение с мобильного устройства]

[hsargsyan]

Многое стало ясно
Большое спасибо!

[novia]

Пишите, если что. Был рад помочь.

[hsargsyan]

Добрый день Novia
нашел файл RPB, сделал орторектификацию, применил 30м дем, но есть рахождения до 3 метров, при сравнении с другим снимком. Эти расхождения были и до орто. Может небходимо применить дем с меньшим размером пикселя?
Или проблема совсем в другом, как думаете?
Заранее спасибо

P.S. Орто делал в ENVI 4.8

[novia]

Другой снимок это какой (тип сенсора)?

Смотрите, у вас снимок ikonos, его разрешение (0.8)1 метр. Даже если бы у вас была супер точная цмр и опорные точки, их все равно пришлось бы вам указывать на снимке.
Точки должны быть с хорошей точностью измерены на местности.
Вариантов два - это результаты геодезической съемки теодолитом/тахеометром в плановом масштабе, либо это точки измеренные "дорогим" gps приемником. Скорее всего их нет.

На снимке вы глазами видите объекты от 1-2м в поперечнике и больше, то есть они как минимум занимают область 3х3 пиксела. То есть точность привязки снимка по этим этим опорным точкам все равно заведомо ниже 1 метра разрешения. и это неприятно (хочется то 0,5 пиксела) но так уж вышло и это нормально.

Вариант действия:
Если есть опорные точки - использовать их для привязки двух снимков, затем оба ортотрансформировать.

Если точек нет, то в зависимости от площади перекрытия либо (при большой площади) привязать один к другому (менее детальный к более детальному) либо (да простят меня фотограмметристы) подтянуть площадь перекрытия "резиновым листом" чтоб "красиво было".

[hsargsyan]

Спасибо за ответ
для сравнения использую снимок IKONOS (более ранний снимок той же местности). Этот снимок отлично орторектифицирован, так как многократные полевые наблюдения точно совпадают со снимком. Если использовать этот снимок для получения опорных точек, затем привязать новый снимок с этими точками и ортотрансформирвать, получится или нет? Мультиспектральный снимок имеет размер ячеек 2.4 на 2.4 м, а панхром 0.6 на 0.6 м

[novia]

0,6 и 2,4 это quickbird. Но не суть.
В целом, если все исследование делалось по первому снимку и точность пространственного положения устраивала, то да, привязывайте, потом ортотрансформируйте.
Будет не сильно лучше в пространственном плане, но зато картинка и результаты исследований будут согласоваными.

[Сообщение с мобильного устройства]

[hsargsyan]

Спасибо за помощь!