RPC-файл

[AndreyL]

Дамы и Господа!

Кто нибудь знает, как сделать RPC-файл? Понятно, что для этого нужно знать конкретные параметры съемки конкретного сенсора – т.е. модель камеры. Но есть подозрение, что простенькую модельку можно получить и без модели сенсора, зная только геометрические параметры съемки и имея несколько хорошо привязанных точек. А как это корректно сделать? Самое главное: как правильно записать все 80 коэффициентов – нигде не нашел полного описания формата: на каком месте какой коэффициент?

[AndreyL]

Интересная штука!

Есть снимок SPOT 5 Pan. Открываю его в ЭНВИ 4.3 как SPOT DIMAP. Строю три слоя геопривязанных точек в соответствии с моделью SPOT (слои 10х10, высота 500, 700 и1000 м), сшиваю слои в один PTS файл, проставляя каждой точке соответствующую высоту. Далее строю RPC по этому PTS файлу. Все параметры в соответствии с инструкцией Хелпа к ЭНВИ: Pushbroom Sensor, Focal Length 580, XY Pixel Size 0.0065, 0.0065, along track incidence angle = <INCIDENCE_ANGLE>, across track incidence angle =<VIEW_ANGLE>, для всех параметров полином второго порядка – точек хватает, Principal Points по нулям. Все нормально, точность – сотые метра. По этим RPC делаю орторектификацию по высоте с использованием ДЕМ (3-секундного, максимальная высота 1200, минимальная 450) и своих GSP (их штук 30) – получаю картинку. Кроме того, делаю орторектификацию с использованием модели SPOT напрямую – с той же ДЕМ и теми же GSP – картинка другая, причем сдвиг в отдельных местах пикселей на 20, а это порядка 100 м. Кто нибудь работал с функцией Building RPCs в ЭНВИ 4.3? Есть подозрение, что косяк именно там.

[AndreyL]

Косяк нашел, не знаю, как обойти. Оказывается при расчете RPC программа теряет информацию о датуме, т.е. устанавливает датум <none>. По идее можно бы было привести ДЕМ к такому датуму, а потом перепроектировать обратно. Но в ЭНВИ нет как такового датума <none>. В любом случае какой то датум используется, только как понять какой.