Научная статья на тему 'Взаимное ориентирование снимков с новым сочетанием угловых элементов в стереопаре'

Взаимное ориентирование снимков с новым сочетанием угловых элементов в стереопаре Текст научной статьи по специальности «Математика»

CC BY
996
118
i Надоели баннеры? Вы всегда можете отключить рекламу.
Область наук
Ключевые слова
ФОТОГРАММЕТРИЯ / ЭЛЕМЕНТЫ ВЗАИМНОГО ОРИЕНТИРОВАНИЯ / КОМПЛАНАРНОСТЬ / АЛЬТЕРНАТИВНАЯ ГРУППА УГЛОВЫХ ЭЛЕМЕНТОВ / ПРЯМАЯ ФОТОГРАММЕТРИЧЕСКАЯ ЗАСЕЧКА / ФОТОГРАММЕТРИЧЕСКАЯ СИСТЕМА КООРДИНАТ / PHOTOGRAMMETRY / ELEMENTS OF THE RELATIVE ORIENTATION / COMPLANARITY / ALTERNATIVE GROUP OF THE ELEMENTS / STRAIGHT PHOTOGRAMMETRY MARCH / PHOTOGRAMMETRY CO-ORDINATES

Аннотация научной статьи по математике, автор научной работы — Добрынин Николай Федорович, Пимшина Татьяна Михайловна

Взаимное ориентирование аэроснимков фундаментальная задача фотограмметрии. От точности ее решения зависит достоверность и надежность окончательного результата, характеризующего пространственное положение сфотографированного объекта (местности). Производственное применение в нашей стране нашли две группы элементов взаимного ориентирования (далее ЭВзО). Однако они не исчерпывают всех возможных сочетаний угловых элементов как по их сущности, так и по распределению между снимками стереопары. К тому же известно, что в обеих системах взаимного ориентирования аэроснимков (базисной и линейно-угловой) обнаружены заметные корреляционные связи между ошибками ЭВзО, приводящие к дополнительному прогибу маршрутной модели в координатных плоскостях ХУ и ХZ. В данной статье рассмотрен один из возможных альтернативных вариантов использования в группе ЭВзО иных угловых элементов. Пять элементов взаимного ориентирования аэроснимков представлены на левом снимке стереопары двумя углами: поперечным углом наклона ω1 и продольным углом наклона α1, двумя углами (α2 и χ2) на правом и азимутом направления базиса проектирования τ. Такое распределение ЭВзО заметно уменьшило корреляционную связь между ошибками угловых элементов. Однако возможны и другие варианты (например, τ, ν, χ1, α2 и ω2). Для вычисления фотограмметрических координат определяемых точек предложена общая формула прямой фотограмметрической засечки, которая может быть использована как для снимков воздушной съемки по любому направлению, так и при фототопографической съемке с точек земной поверхности.

i Надоели баннеры? Вы всегда можете отключить рекламу.
iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.
i Надоели баннеры? Вы всегда можете отключить рекламу.

Relative orientation of the aerial views with the new combination of the corner piece within the stereo mate

Relative orientation of the aerial views is the fundamental task of the photogrammetry. The final results depends on the accuracy and reliability of its solving, characterizing the spatial position of the photographed object (area). Two groups of the elements of the relative orientation found their industrial application in our country (then ERO). However they don’t use up all the possible combinations of the corner pieces res ipsa loquitur and also between the photogrammes. Moreover it is known that in both systems of the relative orientation of the aerial views (base and linear-angular) visible correlation of the mistakes of the ERO are found, which leads to the additional of the route trough model within the coordinate space XY and XZ. One of the possible alternative variants of the application within the group ERO of another corner piece is examined in this article. Five elements of the relative orientation of the aerial views are presented on the left view of the stereo mate as two angles: transverse inclination angle ω1 and longitudinal slope α1, and two angles (α1 and χ2) on the right and the azimuth of the direction of the basis of the design τ. Such a diviosn of the ERO has reduced much the correlation between the mistakes of the angle elements. However other variants are possible (e.g., τ, υ, χ1, α2 and ω2). To calculate the photogrammetry co-ordinates of the general formula of the straight photogrammetry mark is proposed, which can be used both for the aerial shots in any direction and for the shooting from the points of the earth surface. Relative orientation of the aerial views is the fundamental task of the photogrammetry. The final results depends on the accuracy and reliability of its solving, characterizing the spatial position of the photographed object (area). Two groups of the elements of the relative orientation found their industrial application in our country (then ERO). However they don’t use up all the possible combinations of the corner pieces res ipsa loquitur and also between the photogrammes. Moreover it is known that in both systems of the relative orientation of the aerial views (base and linear-angular) visible correlation of the mistakes of the ERO are found, which leads to the additional of the route trough model within the coordinate space XY and XZ. One of the possible alternative variants of the application within the group ERO of another corner piece is examined in this article. Five elements of the relative orientation of the aerial views are presented on the left view of the stereo mate as two angles: transverse inclination angle ω1 and longitudinal slope α1, and two angles (α1 and χ2) on the right and the azimuth of the direction of the basis of the design τ. Such a diviosn of the ERO has reduced much the correlation between the mistakes of the angle elements. However other variants are possible (e.g., τ, ν, χ1, α2 and ω2). To calculate the photogrammetry co-ordinates of the general formula of the straight photogrammetry mark is proposed, which can be used both for the aerial shots in any direction and for the shooting from the points of the earth surface.

Текст научной работы на тему «Взаимное ориентирование снимков с новым сочетанием угловых элементов в стереопаре»

Взаимное ориентирование аэроснимков с новым сочетанием угловых элементов в стереопаре

Н.Ф. Добрынин, Т.М. Пимшина

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

На практике широкое распространение получили две системы взаимного ориентирования стереопары [1]. Их названия базировались на аналоговых принципах обработки снимков. В первой системе взаимное ориентирование осуществлялось на аналоговых приборах при неподвижном базисе проектирования (отсюда и название «базисная»). Вторая система (линейно-угловая) предполагала выполнять взаимное ориентирование пары снимков при неподвижном левом снимке. Такое название она получила из-за того, что на аналоговых приборах отечественного производства взаимное ориентирование велось не только угловыми движениями правой камеры фотограмметрического прибора, но и линейными (движения bz и Ьу). В настоящее время первую систему по-прежнему называют базисной, а вторую - системой левого снимка [2].

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

элементов т и и оказались равными нулю, так как за ось абсцисс фотограмметрической системы координат принят базис проектирования (ву = В'2 = 0), а ю = 0, поскольку плоскость Х2 этой системы

совмещена с главной базисной плоскостью левого снимка, в которой лежит главный луч левой связки, что и привело к отсутствию угла ю в данной ситуации. Однако невозможно «придумать» такую систему, чтобы в ней для стереопары осталось меньше пяти угловых элементов, ибо взаимное положение двух плоскостей в пространстве с точностью до параллельного переноса определяют пять параметров, которые в фотограмметрии называют элементами взаимного ориентирования пары перекрывающихся снимков. Не удастся, разумеется, избавиться одновременно и от любых двух одноименных элементов левого и правого снимков стереопары. Тем не менее, можно предложить целый ряд новых сочетаний элементов взаимного ориентирования (ЭВзО) снимков посредством варьирования направлениями координатных осей фотограмметрической системы координат, в которой решается данная задача. К этому побуждают и следующие обстоятельства.

Погрешности ЭВзО, вызванные известными причинами, приводят к деформациям в плане и по высоте созданной с их использованием пространственной фототриангуляционной сети (построенной, например, способом независимых или частично зависимых моделей [1]). И хотя эти погрешности вполне подчиняются нормальному закону их распределения (обладают свойством компенсации), тем не менее, кроме уже известных факторов ситуация усугубляется наличием корреляционных зависимостей между ошибками ЭВзО смежных моделей. На это обратил внимание еще в середине прошлого века доктор технических наук И.Т. Антипов [3]. Согласно его исследованиям эти зависимости (корреляционные связи) могут приводить к накоплению ошибок в угловых элементах внешнего ориентирования аэроснимков, что вызовет дополнительный прогиб маршрутной модели в координатных плоскостях Х2 и ХУ.

Поскольку используемые на практике базисная и линейно-угловая группы ЭВзО получены без учета корреляционных связей [3], можно предположить, что могут существовать сочетания элементов и с меньшими коэффициентами корреляции между их ошибками. Современные тенденции развития фотограмметрических принципов измерений, о которых шла речь в работе [4], также могут повлиять на взаимозависимости ошибок ЭВзО в цифровой фотограмметрии. Например, способ идентификации соответственных точек на цифровых снимках существенно отличается от бинокулярного принципа создания искусственного стереоэффекта на аналоговых фотограмметрических приборах. Может быть и сам «пиксельный» подход внесет определенные коррективы в решение задачи взаимного ориентирования аэроснимков (см. например, [5, 6]). Появились публикации (например, [7]), которые возвращают нас к использованию во взаимном ориентировании аэроснимков алгебры кватернионов. Этот метод предложен профессором М.С. Урмаевым [8] в восьмидесятых годах прошлого столетия. Однако, судя по публикациям, преимущества его пока еще не доказаны.

Рассмотрим один их возможных вариантов нового сочетания угловых элементов в стереопаре. Он не исследован на корреляционные связи, а призван обратить внимание ученых и специалистов на возможно новое направление в исследовании технологии взаимного ориентирования аэроснимков.

Обращаясь к теоретическим основам взаимного ориентирования аэроснимков с новой группой угловых элементов, вместо традиционной угловой ориентации снимков [9] а, ю, % воспользуемся новой трактовкой их последовательности, то есть ю, а, %. Тогда связь направляющих косинусов с угловыми элементами внешнего ориентирования аэроснимка будет представлена так:

а; = cosa cosx; а2 = - cosa sin x;

a3 = sin a;

b; = sin® sin a cos x + cos® sin x;

b2 =-sin® sina sinx + cos® cosx;\ (1)

b3 = - sin® cosa;

c; = - cos ® sin a cos x + sin ® sin x; c2 = cos® sina sinx + sin® cosx; c3 = cos® cosa.

Пусть взаимное ориентирование аэроснимков ведется в фотограмметрической системе координат (Х, У, Z), начало которой совмещено с центром проекции (S1) левого снимка. При этом плоскость ХУ проходит через базис проектирования (В') параллельно оси ординат (у2) правого снимка, а координатная плоскость 2У параллельна оси ординат (у1) левого снимка.

При таком выборе фотограмметрической системы координат, в которой будет осуществляться взаимное ориентирование аэроснимков, шесть элементов внешнего ориентирования пары снимков примут нулевые значения. Так как начало системы совпадает с S1, то XS1 = yS1 = ZS1 = 0. Базисный компонент B'Z также примет нулевое значение, поскольку базис проектирования оказался в плоскости ХУ. Равны нулю и угловые элементы %ь первый вследствие взаимной параллельности осей У и у1, второй - из-за параллельности оси у2 плоскости ZУ. Оставшиеся не нулевые угловые элементы и составят группу ЭВзО:

т, mi, ai, a2 и Х2.

Известное условие пересечения соответственных лучей стереопары [1] для назначенной выше системы координат примет следующий вид:

В'х ВУ 0

х; у z'

х 2 у 2 z 2

или

ВХ

1 Щт 0

х' у1 7 Х2 у2 7'

= 0.

где

Х1 , У1 , 71; Х2, У2, 7 2

пространственные координаты

соответственных точек на левом и правом аэроснимках. Поскольку ВХ * 0, следовательно,

1 tgт 0

х; у; 7;

х2 у2 72

=0

или, раскладывая данный определитель по элементам первой строки,

(у; 72 - 7; У2) + (7; х2 - х; 72) tgт = 0.

(2)

Выражение (2) представляет собой уравнение взаимного ориентирования снимков в назначенной выше фотограмметрической системе координат. Заметим, что в данном случае пространственная координатная система левого снимка совпадает с фотограмметрической, а правого снимка стереопары (с началом в 52) ей параллельна. Однако фотограмметрические координаты (Х, У, 7) относятся к точке модели (точке пересечения соответственных лучей после взаимного ориентирования снимков), а пространственные х; , у;, 7;' ; Х2, У2, 72 характеризуют положение этой же

точки на левом и правом аэроснимках.

Поскольку при взаимном ориентировании непосредственно измеряемыми величинами являются плоские координаты соответственных точек на снимках (х;, у1; х2, у2), то для перехода к пространственным воспользуемся формулами связи [9], имеющими следующий вид:

Х' = а; х + а2 у - а3 /; У ' = Ь х + Ь2 у - Ьз /; 7' = с; х + с2 у - с3 /.

Уравнения пространственных преобразований (3) совместно с направляющими косинусами (1) позволяют включить в уравнение взаимного ориентирования стереопары (2) определяемые неизвестные т, a;, a2 и Х2. Учитывая, что угловые элементы и m2 равны нулю, направляющие косинусы (1) примут новые значения. Тогда уравнения (3) в окончательном виде будут выглядеть так:

Х; = х; cos a; - f sin a;;

У' = х; sin ® ; sin a; + y; cos ® ; + f sin ® ; cos a;; Z; = -x; cos ® ; sin a; + y; sin ® ; - f cos ® ; cos a;; Х2 = х2 cos a2 cos x2 - y 2 cos a2 sin x2 - f sin a2;

у2 = x2sin x2 + y 2cos x2;

Z2 = -x2 sin a2 cos x2 + y2 sin a2 sin x2 - f cos a2,

(4)

где / - фокусное расстояние аэроснимка.

После подстановки пространственных координат (4) конкретной точки с плоскими ее координатами х1, у1; х2, у2 в уравнение взаимного ориентирования (2) получим тригонометрическое (нелинейное) его выражение с пятью определяемыми неизвестными ЭВзО. Следовательно, для их определения необходимо измерить плоские координаты хотя бы пяти точек. Такую систему уравнений обычно решают методом Ньютона через приближенные значения неизвестных. При избыточном числе точек (п > 5) используют метод наименьших квадратов (МНК). Рассмотрим параметрический способ его реализации.

Пусть известны приближенные значения определяемых ЭВзО (т), (юд (а1), (а2) и (х2), и измерены плоские координаты соответственных точек стереопары (п > 5). Тогда уравнения поправок, соответствующих функции (2), примут следующий вид:

а1 8т + Ъ1 8о) + с1 а + 8аг + е1 8х2 +11 = и1;

а2 8т + Ъ2 8а 1 + с2 8ах + ё2 да2 + е2 д%2 +12 = и2; (5)

8т + bn 8® + cn 8a; + dn 8a2 + en 8x2 +1 = un,

n n i n 1 n 2 n Л 2 n n5

где a¡, bu ..., et - коэффициенты, которыми являются частные производные от функции (2) по переменным т, а;, а2, Х2, а свободные члены (l) подсчитываются посредством подстановки в нее приближенных значений ЭВзО, используя выражение (4) и плоские координаты соответственных точек стереопары.

Вычислим несколько коэффициентов, например, а и b.

Коэффициент а вычисляется просто:

а = — = (z; X2 - X[ Z2) sec2 т. дт

Для определения коэффициента b сначала вычислим частные производные от пространственных координат левого снимка, в которые входит угловой элемент

dx; 0 ( )

—L = 0 1не используется элемент а при вычислении Х J ;

да i 1 '

дУ^ i у r~z t

-- = xcose sma; - y1smc + f cose cosa; =-Z;;

да 1 1 1 1 1 1 1 1

i

dZ'

-- = x;sinc sina; + y;cose + f sine cosa; = У'.

да 1 11

1

Остальные пространственные координаты соответственных точек для переменной будут считаться константами. Подставляя теперь в общий дифференциал функции (2) значения полученных частных производных, для коэффициента b получим следующее значение:

b = с = у; (x2 tgT-у2)-z; z2.

да

i

Аналогично вычисляются и остальные коэффициенты уравнений поправок (5):

X; (у 2 cos с + Z 2 sine ')+ [z 2 ( sin«' + f cosa;)- X; X 2 cos с ' ] tgT;

dF

с =-= Xi У2 cose + Z2 sine J+ |

da;

е =

d = — = (x; x2 + z; Z2) tgT - у; X2;

oa2

dF í \

-= У' У2 sin a2 - Z; (x2 cos x2 - y2 sin x2) - (x; У2 sin a2 + Z; У2 cos a2) tgT.

Полученную систему уравнений поправок (5) при n > 5 решают под условием

[[[] = min,

то есть используя МНК. Тогда от n уравнений поправок переходят к пяти нормальным уравнениям, из решения которых получают поправки к приближенным ЭВзО. Поскольку они, как правило, далеки от их истинных значений, решение задачи ведется методом последовательных приближений. При этом в каждом приближении коэффициенты уравнений поправок (5) уточняются новыми значениями ЭВзО, вычисляемые по следующим формулам:

т = (т )+ 5т' + 5т;2 +...;

a = (с )+ Sa 1 + 5a 2 +...;

1 1 1 1

а1 = (а )+ 5а11 + 5а;2 +...;

а2 = (а2)+ 5а21 + 5а22 +...;

X = (Х ) + S%2 +5Х22 + ... .

Обычно достаточно 4-6 приближений до полной сходимости итерационного процесса.

Известно, что в результате взаимного ориентирования стереопары создается модель сфотографированного объекта. В компьютерных технологиях вычисляются фотограмметрические координаты (Х, У, Z) на основе прямой фотограмметрической засечки [1]:

Х = N X'; У = N У'; Z = N Z'.

Множитель N вычисляют по разным формулам в зависимости от вида съемки (фототопографическая, воздушная, космическая). Предлагается общая формула для всех упомянутых видов съемки:

N = B'x (у2 - z2)+ву (z2 - x2)+bz (x2 - У2) х; (у2 - z2)+у; (у2 - x2)+z; (x2 - У2).

iНе можете найти то, что вам нужно? Попробуйте сервис подбора литературы.

Для воздушной съемки данная формула повышает точность определения N за счет использования большего числа измеренных величин при его вычислении.

При полной фотограмметрической обработке отдельной стереопары (одиночной модели) для перехода от фотограмметрической системы координат к геодезической применяются известные [10] уравнения пространственных преобразований, для чего используются опорные геодезические точки.

Литература:

1. Лобанов, А.Н. Фотограмметрия [Текст]: Учебник для вузов / А.Н. Лобанов. - М.: Недра, 1984. - 552 с.

2. Обиралов, А.И. Фотограмметрия [Текст]: Учебник для вузов / А.И. Обиралов, А.Н. Лимонов, Л. А. Гаврилова. - М.: КолосС, 2004. - 240 с.

3. Антипов И.Т. О корреляционных связях ошибок элементов взаимного ориентирования в рядах аналитической пространственной фототриангуляции [Текст] // Труды НИИГАиК, 1968. - т. XXI. - С. 39-46.

4. Тюфлин Ю.С. Фотограмметрия - вчера, сегодня и завтра [Текст] // Изв. Вузов. «Геодезия и аэрофотосъемка», 2011. - №2. - С. 3-8.

5. 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. - P. 298-306. Англ.

6. Balci Murat, Foroosh Hassan. Subpixel estimation of shifts directly in the Fourier domain //IEEE Trans. Image Process, 2006. 15, - № 7. - P. 1965-1972. Англ.

7. Безменов В.М. Элементы кватернионов в фотограмметрии [Текст] // Изв. Вузов. «Геодезия и аэрофотосъемка». 2009. - №4. - С. 48-56.

8. Урмаев М.С. Применение алгебры кватернионов в фотограмметрии [Текст] // Изв. Вузов. «Геодезия и аэрофотосъемка». 1986. - №2. - С. 81-90.

9. Добрынин Н.Ф. , Пимшина Т.М. Использование космических средств позиционирования при обработке аэро- и космической информации [Электронный ресурс] // «Инженерный вестник Дона», 2013, № 3. - Режим доступа: http://www.ivdon.ru/magazine/archive/n3y2013/1835 (доступ свободный) - Загл. с экрана. - Яз. рус.

10. Куштин, В.И. Преобразование координат, используемых в геодезии [Электронный ресурс] // «Инженерный вестник Дона», 2012, №4-2. - Режим доступа: http://www.ivdon.ru/magazine/archive/n4p2y2012/1240 (доступ свободный) - Загл. с экрана. - Яз. рус.

i Надоели баннеры? Вы всегда можете отключить рекламу.