CC BY
17
Взаимное ориентирование аэроснимков с оптимальным сочетанием
угловых элементов в стереопаре
Н. Ф. Добрынин, Т.М. Пимшина Ростовский государственный университет путей сообщения
Аннотация: Рассмотрен возможный вариант взаимного ориентирования аэроснимков с оптимальным выбором угловых элементов.
Ключевые слова: Фотограмметрия, условия компланарности, оптимальное сочетание элементов взаимного ориентирования, корреляционные связи.
Отрасль геодезии и картографии переживает непростые времена. Состояние предприятий, особенно финансовое, объединенных в ОАО «Роскартография» еще в 2012 году [1], из-за допущенных промахов проработки Концепции развития отрасли до 2020 года [2] достаточно тяжелое. По признанию работников Росреестра [3] реализовано лишь 15% от всех предусмотренных Концепцией мероприятий.
Однако научные исследования в отрасли геодезии и картографии, может быть в меньшей степени из-за отсутствия дополнительного финансирования со стороны ОАО «Роскартографии», продолжаются. Остановимся лишь на совершенствовании методов взаимного ориентирования аэроснимков и построении одиночной модели.
Все настойчивее традиционному методу построения модели по паре снимков ставится в альтернативу предложенный еще в конце прошлого столетия проф. М.С. Урмаевым [4] новый принцип реализации уравнений коллинеарности и компланарности на основе аппарата алгебры кватернионов. Например, в работах В.М. Безменова [5, 6] (Казанский федеральный университет) приведены алгоритмы реализации некоторых уравнений фотограмметрии названным способом. Однако, как уже указывалось в работе авторов данной статьи [7], пока далее возможности использования алгебры кватернионов в фотограмметрии по нашим данным
сведений нет. Остается ждать дальнейших практических результатов, чтобы выявить достоинства этого подхода к решению фотограмметрических задач.
Сделана попытка «реанимировать» предложенный докт. техн. наук, проф. В.Б. Дубиновским (Московский государственный университет геодезии и картографии) метод калибровки снимков реализацией только условий компланарности, т.е. без использования опорных данных. Позднее было установлено, что такой подход полностью не решает поставленную задачу. Однако в работе [8] предпринята новая попытка отказаться от создания дорогостоящего полигона с опорными геодезическими точками, а найти более простое решение. Авторы данной публикации полностью согласны с позицией докт. техн. наук, проф. И.Т. Антипова (Сибирская государственная геодезическая академия), изложенной в его работе [9], в которой он последовательно и убедительно показал несостоятельность приведенных в статье [8] выводов и предложений.
К взаимному ориентированию аэроснимков на протяжении всей истории развития фотограмметрии уделялось и уделяется в настоящее время пристальное внимание. Во-первых, в результате взаимного ориентирования снимков стереопары создается модель сфотографированного объекта. Во-вторых, от надежности решения этой фундаментальной задачи фотограмметрии зависит точность окончательных результатов обработки стереопары.
Обращаясь к точности определения элементов взаимного ориентирования (ЭВзО), необходимо заметить, что кроме известных факторов, влияющих на точность взаимного ориентирования снимков (дисторсия объектива АФА, деформация, погрешности измерительного прибора и т.д.), замечены корреляционные связи между ошибками самих ЭВзО. Причем коэффициенты корреляции варьируют в зависимости от выбора фотограмметрической системы координат, в пространстве которой
ведется определение ЭВзО. Связано это с различным сочетанием угловых элементов, которое обусловлено как раз выбором фотограмметрической системы координат. Например, в базисной системе элементов наибольшая корреляционная связь замечена между углами а и % как левого, так и правого снимков стереопары. Предлагаемый ниже вариант взаимного ориентирования снимков практически свободен от взаимных корреляционных связей ЭВзО.
Для взаимного ориентирования аэроснимков используем следующие ЭВзО: ш1, %1 для левого и а2 и для правого снимков стереопары, а также угол наклона базиса фотографирования (проектирования) V. Остальные углы (т, а1 и ю2) должны быть нулевыми. Такая ситуация возникнет в том случае, если фотограмметрическую систему координат, в пространстве которой будет осуществляться взаимное ориентирование аэроснимков, назначить следующим образом:
- начало системы в левой точке фотографирования (51);
- плоскость 7Х совместить с главной базисной плоскостью правого снимка (проходит через базис фотографирования и главный луч правого снимка стереопары);
- плоскость направить через главный луч левого снимка перпендикулярно плоскости 7Х.
При таком выборе фотограмметрической системы координат базисные компоненты (ВХ, Ву, В£) получат следующие значения:
ВХ = Х52 - Х51 = В ^^
В = V — V = 0-
V 1 8 2 В2 = 2 - = В ^^
(1)
где В - базис фотографирования;
V - угол его наклона относительно плоскости XV.
J
Традиционное условие пересечения соответственных (одноименных) проектирующих лучей в координатной форме в виде определителя третьего порядка примет вид
1 1 0 tgv
х; у' 7; = 0,
Х 2 У2 7 2
или, раскладывая его по элементам первой строки,
(XY - y[x2 )tgv + (YZ2 - z;y;) . (2)
Здесь X ', Y ', Z ' - пространственные координаты точки на левом и правом снимках стереопары в пространственных прямоугольных системах с началами в точках фотографирования (S1 и S2) и осями, параллельными одноименным осям фотограмметрической системы.
Направляющие косинусы, являющиеся функциями углов ю1, %ь а2, х2, примут следующие значения:
a11 = cosхj; a21 = cosa2 cosx2; а12 = - sin X1; a22 =- cos a2 sin X2 ;
a13 = 0; a23 = sina2;
bu = cos sin x 1; ¿21 = sin x 2; b12 = cos®1cos x 1; ¿22 = cos X2;
b13 = - sin a>1; b23 = 0;
cu = sin®1 sin x 1; c21 = - sina2cosx2; c12 = sin®1 cos x 1; c22 = sina2sinx2; c13 = cos®1; c23 = cosa2,
а связь между пространственными и плоскими (х, у) координатами точки на стереопаре будет выглядеть так:
:
Х1 = X! С0Б X1 - у1 ^п х
У' = х1 о х! + у1 соб о соб х! + / бш о ^ 2' = хх о х! + У! о соб х! - / соб о 1; Х2 = х2 соб а2 соб х2 - У2 соб а2 х2 -1 ^п а2; У2 = х2 бШ Х2 + У2 СОБ %2;
22 = -х2 а2 соб х2 + У2 вт а2 Бт х2 - / соб а2,
(3)
где / - фокусное расстояние снимков стереопары.
Выражение (2) является уравнением взаимного ориентирования снимков стереопары в назначенной (целенаправленно выбранной) фотограмметрической системе координат. Уравнение строгое и может использоваться для определения любых значений угловых ЭВзО. В то же время оно тригонометрическое (нелинейное) и содержит пять неизвестных, для вычисления которых необходимо в пределах продольного перекрытия аэроснимков измерить плоские координаты (х1, у1, х2, у2) как минимум пяти точек стереопары. Решают такую систему уравнений, как правило, методом Ньютона через приближенные значения определяемых неизвестных. При избыточном числе уравнений (измерены плоские прямоугольные координаты более пяти точек стереопары) применяют метод наименьших квадратов. Рассмотрим параметрический способ реализации метода Ньютона.
Пусть нам известны предварительные (приближенные) величины ЭВзО
(V), (^1), (х1), (а2) и (хг) и измерены плоские координаты соответственных
(одноименных) точек стереопары. При этом п > 5 (п - число точек для
определения ЭВзО), а поправки к предварительным значениям неизвестных
обозначим Sv, Ью 1, §х1, §а2 и 8%2. Тогда традиционно уравнения поправок,
соответствующих исходному выражению (2), примут следующий вид:
а1Ьу + Ь1до1 + с1Ьх 1 + dlSa2 + е1Ьх 2 +11 = и1; а2Ьу + Ъ2д(о1 + с2Ьх 1 + d2дa2 + е2Ьх 2 +12 = и2;
(4)
апЬу + Ьпдсо1 + спЬх 1 + dnдaг + епЬх 2 + I =ип,
п п 1 п /М п 2 п г* 2 1 п"
J
где а, Ь, с, d, е - коэффициенты, являющиеся частными производными первого порядка от выражения (2) по переменным V, Ш\, %ь а2 и а свободные члены Д вычисляются подстановкой в уравнение (2) приближенных значений определяемых ЭВзО и вычисленных по ним (формулы (3)) пространственных координат соответственных точек стереопары.
Вычислим несколько коэффициентов для уравнений поправок (4), например, а и Ь.
Коэффициент а вычисляется просто:
^ = {х; Y2- Yi' Х2) sec2 V, dv
так как во втором слагаемом функции (2) нет переменной v (частная производная будет равна нулю), нет ее и в пространственных координатах соответственных точек стереопары, которые здесь играют роль лишь коэффициентов при указанном выше определяемом неизвестном.
Чтобы вычислить коэффициент b сначала найдем частные производные для первых трех выражений в формулах (3), в которые входит ЭВзО т1:
дХ' 0 ( Х, )
—1 = 0 1не участвует элемент a при вычислении Х1 I ;
Sai 1 '
dY' ... „
-— = -X' Sine ' SinX1 - У1 Sine ' cosX1 + f cose ' = -Z{;
1
dZ' = da
= x1 sin a sin a1 sin x1 + y1 cos a cos x1 + f sin a = Y'.
1
Пространственные координаты точки на правом снимке стереопары будут выступать константами в общем дифференциале функции (2), поскольку в вычислениях этих координат не используется ю1.
Таким образом, для коэффициента Ь получим следующий результат:
Ь = 0- = 1[ (х2 tgv-12)-у;У'.
да
:
Аналогично вычисляются частные производные (коэффициенты уравнений (4)) по остальным переменным. В результате получим недостающие для уравнений поправок (4) коэффициенты:
[- у; ( МП X! + У: сое X!)- Х'2 ( сое X! - У! вт х! )со8 0 1 ] ^ +
dF
с = ■
dX i
+ (x1 cos х 1 - у1 sin х 1) (z2 cos с1 - Y2' sin со1 )
d = а = у;Z2 tgv-Y1X2; да
dF
е =
дХ
= (2 cos ^2 - У2 sin ^2 ) (x1 tgv - ZÖ+ (x2 sin X2 + У2 COs X2 ) ( COs «2 + Y1 sin «2 )
2
Составленную систему уравнений поправок (4) при n > 5 решают под условием
[[[] = min,
т.е. с использованием метода наименьших квадратов.
Поскольку в результате взаимного ориентирования аэроснимков, вообще говоря, создается модель сфотографированного объекта, приведем формулы для вычисления пространственных фотограмметрических координат точек этой модели:
X = N X' ; Y = N Y''; Z = N Z[,
где множитель N для любой (правой) фотограмметрической системы
координат можно вычислить по следующей универсальной формуле [10]:
N = Bx(Y2-Z2)+By (z2-x2)+bz (X2-y;) X' (y; - z2)+Y1 (z2 - x2)+Z1 (X2 - Y;) '
Процесс взаимного ориентирования с новыми коэффициентами в уравнениях поправок (4) проверен по макетным снимкам. Максимальный коэффициент корреляции, как и ожидалось, выявлен между ошибками в
углах а2 и хх и составил +0,23, т.е. по нашим данным ниже, чем при других сочетаниях ЭВзО.
Литература:
1. Указ Президента Российской Федерации «Об открытом акционерном обществе «Роскартография»» от 12.03.2012 г. № 296 // Собрание актов Президента и Правительства Российской Федерации.
2. Распоряжение Правительства Российской Федерации «Об утверждении Плана мероприятий по реализации Концепции развития отрасли геодезии и картографии до 2020 г.» от 07.07.2011 г. № 1177-р // Собрание актов Президента и Правительства Российской Федерации.
3. Васильев И.В., Коробов А.В., Побединский Г.Г., Приданкин А.Б. Топографо-геодезическое и картографическое обеспечение Российской Федерации. Состояние и перспективы развития отрасли геодезии и картографии // Геодезия и картография. 2014. №12. С. 2-11.
4. Урмаев М.С. Применение алгебры кватернионов в фотограмметрии // Изв. вузов «Геодезия и аэрофотосъемка». 1986. №2. С. 81-90.
5. Безменов В.М. Элементы кватернионов в фотограмметрии // Изв. вузов «Геодезия и аэрофотосъемка». 2009. №4. С. 48-56.
6. Безменов В.М. Применение кватернионов в фотограмметрии // Изв. вузов «Геодезия и аэрофотосъемка». 2014. №5. С. 22-27.
7. Добрынин Н.Ф., Пимшина Т.М. Взаимное ориентирование снимков с новым сочетанием угловых элементов в стереопаре // «Инженерный вестник Дона», 2014, № 2 URL: ivdon.ru/magazine/archive/n2y2014/2347.
8. Семенцов А.В. Калибровка камер без использования твердых опорных данных // Геодезия и картография. 2014. №4. С. 26-30.
9. Антипов И.Т. Калибровка камер без использования твердых опорных данных (рецензия) // Геодезия и картография. 2014. №7. С. 59-60.
10. Добрынин Н.Ф., Пимшина Т.М. Использование космических средств позиционирования при обработке аэро- и космической информации // «Инженерный вестник Дона», 2013, № 3 URL: ivdon.ru/magazine/archive/n3y2013/1835.
11. Lu Dengsheng, Moran Emilio, Hetrick Scott. Detection of impervious surface change with multitemporal Landsat images in an urban- rural frontier // ISRRS J. Photogramm. and Retome Sens, 2011. 66, № 3. pp. 298-306.
12. Balci Murat, Foroosh Hassan. Subpixel estimation of shifts directly in the Fourier domain //IEEE Trans. Image Process, 2006. 15, - № 7. pp. 1965-1972.
References
1. Ukaz Prezidenta Rossijskoj Federacii ot 12.03.2012 g. № 296. Sobranie aktov Prezidenta i Pravitel'stva Rossijskoj Federacii.
2. Rasporjazhenie Pravitel'stva Rossijskoj Federacii ot 07.07.2011 g. № 1177-r. Sobranie aktov Prezidenta i Pravitel'stva Rossijskoj Federacii.
3. Vasil'ev I.V., Korobov A.V., Pobedinskij G.G., Pridankin A.B. Geodezija i kartografija. 2014. №12. pp. 2-11.
4. Urmaev M.S. Izv. vuzov «Geodezija i ajerofotos'emka». 1986. №2. pp. 8190.
5. Bezmenov V.M. Izv. vuzov «Geodezija i ajerofotos'emka». 2009. №4. pp. 48-56.
6. Bezmenov V.M. Izv. vuzov «Geodezija i ajerofotos'emka». 2014. №5. pp. 22-27.
7. Dobrynin N.F., Pimshina T.M. Inzenernyj vestnik Dona (Rus), 2014, № 2. URL: ivdon.ru/magazine/archive/n2y2014/2347.
8. Semencov A.V. Geodezija i kartografija. 2014. №4. pp. 26-30.
9. Antipov I.T. Geodezija i kartografija. 2014. №7. pp. 59-60.
10. Dobrynin N.F., Pimshina T.M. Inzenernyj vestnik Dona (Rus), 2013, №3 URL: ivdon.ru/magazine/archive/n3y2013/1835.
11. Lu Dengsheng, Moran Emilio, Hetrick Scott. Detection of impervious surface change with multitemporal Landsat images in an urban- rural frontier. ISRRS J. Photogramm. and Retome Sens, 2011. 66, № 3. pp. 298-306.
12. Balci Murat, Foroosh Hassan. Subpixel estimation of shifts directly in the Fourier domain. ЯЕЕЕ Trans. Image Process, 2006. 15, № 7. pp. 1965-1972.
