УДК 681.786.4
К. Г. Араканцев, Д. В. Жуков, И. А. Коняхин
МЕТОДЫ ОБРАБОТКИ ИЗМЕРИТЕЛЬНОЙ ИНФОРМАЦИИ И ОЦЕНКИ ПОГРЕШНОСТЕЙ В СТЕРЕОСКОПИЧЕСКОЙ СИСТЕМЕ КОНТРОЛЯ ПРОСТРАНСТВЕННОГО ПОЛОЖЕНИЯ ОБЪЕКТОВ
Рассмотрены методы обработки измерительной информации в стереоскопической системе контроля пространственного положения объектов. Обоснована возможность применения теории возмущений для оценки систематических и случайных погрешностей измерений.
Ключевые слова: перспективное преобразование, стереоскопическая система, сингулярное разложение, теория возмущений, погрешность измерения.
Введение. Для дистанционного контроля пространственного положения объектов широко используются стереоскопические системы технического зрения. В частности, оптико-электронная система стереоскопического типа может быть использована для контроля смещений железнодорожного полотна относительно проектного положения [1]. В этом случае задача измерительной системы сводится к определению расстояния до рельсов в профиле и плане относительно реперных меток, расположенных вдоль пути на столбах геодезической сети. Диапазон смещений реперной метки относительно начала приборной системы координат должен составлять 2000—10 000 мм в плане и ±150 мм — в профиле, при условии, что среднеквадратическое отклонение (СКО) измерений по каждой координате не превысит 1 мм.
В настоящей работе описана модель стереоскопической системы, основанной на перспективной проекции, представлены два метода численного расчета положения точечного объекта относительно приборной системы координат, также показана возможность использования теории возмущений для анализа систематических и случайных погрешностей в стереоскопических оптико-электронных системах.
Модель стереоскопической системы. Стереоскопическая измерительная система, построенная на основе двух видео- или фотокамер с матричными приемниками изображения и перекрывающимися полями зрения объективов, позволяет получить четыре координаты изображений объекта, которые совместно с центрами проекции камер определяют в пространстве положение двух прямых, на пересечении которых находится точечный объект.
Координаты изображений точки Р в плоскостях анализа стереоскопической системы, полученные в результате прямого перспективного преобразования [2], описываются переопределенной системой нелинейных уравнений, где в качестве неизвестных выступают координаты объекта хо, Уо и хо.
Хр^ — Оц
' А1 Хо + А2 Уо + Аз хо - Ъ1
А4 ХО + А5 Уо + 4 хО -— А7 хо + А8 Уо + 4 хо -
А10 хо + А11 Уо + А12 хо - tZ11
' А13 хо + А14 Уо + А15 хо -УР — а\-,
I А16 хо + А17 Уо + А18 хо - tZ 1
, А19 хо + А20Уо + А21хо - \ Ур — а11-—
II А22 хо + А23 Уо + А24 хо - tZll
(1)
Методы обработки информации и оценки погрешностей пространственного положения объектов 11
Здесь Л1—А 24 — коэффициенты, зависящие от внешних и внутренних параметров камер; а'Т, а'ц — задние отрезки объективов; х^ , Хртт , ур , ур — координаты изображений объекта в
плоскостях анализа камер рТ и рц; ¿у, — линейные координаты центра проекции камеры в глобальной (О) системе координат. Для нахождения координат объекта, удовлетворяющих системе уравнений (1), предлагается использовать один из предложенных методов.
Метод скрещивающихся прямых основан на оценке вектора координат Р точки Р, сумма квадратов расстояний от которой до прямых, задаваемых точками изображений и центрами проекции камер стереоскопической системы, минимальна. Точка Р лежит на середине отрезка р Рц, который является общим перпендикуляром для двух указанных прямых [3].
Показанные на рис. 1 направляющие векторы г и гп проецирующих лучей задаются выражениями
ГРТ
О 7
'II
рТТ
О 7
где Ррт, Рртт — векторы, состоящие из координат изображений объекта в глобальной системе координат; О7, О^ — векторы, определяющие положение начал отсчета локальных (РТ и РТТ) систем координат камер в глобальной системе координат.
О
7
^Ор,
р>
рт \ ..-• \
р
Т.;
\
рТТ
/
.........
/
/
О
\
/
\
\
2
V \
/
/
/
/
\
/
\
о ь1
I I
\
/
/
\
О,
■/о
\
Р
о
Х(
К
/р \
/
о
(2)
/
Рис. 1
Уравнения проецирующих лучей запишутся в виде
Ррт = О7Т - ЧтГТ , Рртт = О7ТТ - ЧПГП , где чт, чтт — свободные параметры.
Значения свободных параметров Чт и Чп , соответствующие минимальному расстоянию
между скрещивающимися прямыми, находятся из выражений
\Т Г/ \Г
(О7тт - О7т ) ГТ - (О7ТТ - О7т ) ГТТ (ГТТГТТ )
Чт =
1-
( ГТТ)
1 -(■/ * )
После подстановки найденных коэффициентов 1 и 1П в уравнения прямых (2) получим координаты концов искомого отрезка в глобальной системе координат
гст
О,
¿¿т гт, Роп = О Ьтт
■ 1тт гтт ,
полусумма которых даст оценку Р вектора координат точечного объекта Р
Рот + Р
Отт
го
2
Метод сингулярного разложения основан на решении переопределенной системы линейных уравнений с помощью сингулярного разложения [4]. На рис. 2 приведена геометрическая схема метода, в котором для нахождения координат точечного объекта при статических измерениях для каждой из камер выбираются две плоскости (I' и I", П' и П"), задающие своим пересечением две прямые, каждая из которых содержит точку изображения объекта и центр проекции камеры.
1а
Ха
О
От
Рис. 2
Уравнения выбранных плоскостей образуют следующую систему:
АотРох + В'оРт + СотРог + АОт = 0, ' АО т Рох + во т Рот + СОт Рог + Щ, = 0, АОтт Рох + Воп Рот + Со,, Рог + ^ = 0, Аоп Рох + воп Рот + Соп Рог + ^ = 0,
(3)
п
Методы обработки информации и оценки погрешностей пространственного положения объектов 13 где А,—AGn, AG— AGn, B'Gl— BGa, B^-BGa, CGI — Cbn, C'^-CGa, D'Gi-DGa, Dg — DGn — коэффициенты, зависящие от внешних и внутренних параметров камер, а также измеренных координат изображений объекта.
Для каждой камеры в плоскости анализа задаются две прямые, пересекающиеся в точке изображения
yf = kf X' + /,', y'i = kf, X', + /¿,
~ff 1 rr~rr . Iff ~ff 1 ft ~ff . Iff
y, = kf Xf + /f, yn = k" xn + ^ где kf, kf, kfi, kfj — коэффициенты наклона прямых в плоскостях анализа Ff и Fn; /f, /' , /', /Ц — коэффициенты линейного смещения прямых в плоскостях анализа Fi и Fn. Затем для камер выбираются четыре плоскости, каждая из которых содержит одну из заданных прямых и центр проекции камеры. Для каждой камеры получается по два линейных уравнения плоскостей, коэффициенты в которых равны
А^ = /'a', В'к = -af, C'li = /', D'k = 0;
A"li = /faf, BLj =-af, CLj = /f, DLj = 0;
A'T = /IIan, B'r = -ah, C'r = HI, DL = 0;
'nun
.|ff _ Ttt ff
Alп = 'И aii
73/f ff /"fff Iff ТЛ" Г\
BL" =-aII, CLu = /II, DLn =
В матричном виде уравнения плоскостей запишутся следующим образом:
H l Tg R g' PPg = 0, H L' tgl1 R G PPg = 0,
H LfI Tg11 R G11 PPg = 0, H LfI Tl11 R£ Pg = 0
или
HGi Pg = 0, HGi Pg = 0, H Gn Pg = 0, H G,, Pg = 0,
где Н, Н^ , Н^ , Н— столбцы коэффициентов плоскостей для локальных систем координат камер; , Н GI, Н^ , Н^ — столбцы коэффициентов плоскостей для глобальной
"" Т
системы координат; I>G = ((, ^, % ,1) — искомый вектор оценок однородных координат
объекта; RG, Т^^ — матрицы поворота и переноса локальных систем координат камер относительно глобальной системы координат.
Окончательно система линейных уравнений запишется в виде
МР^ = N,
где
M
Система уравнений (3) является переопределенной. Для ее решения предлагается использовать алгоритм сингулярного разложения БУО [4], в соответствии с которым любая (ЖхК)-матрица вещественных чисел М может быть представлена в виде
М=иЕУТ,
H ' T Л HGj Ah, BGi f DGI '
HGjT AGi BS, C&i , N = - DGi
HGT Gii Ah', B4i CG'i DGI,
HG T Gii J Ah', Bh'i C" IDGii J
где и — ортогональная (Жх Ж)-матрица, V — ортогональная (КхК)-матрица, £ — (ЖхЖ)-матрица, на главной диагонали которой находятся сингулярные числа матрицы М, расположенные в порядке убывания, а все внедиагональные элементы равны нулю. Для нахождения подматриц V, и существуют итерационные алгоритмы, реализованные в системах МаШСаё и Ма1ЬаЬ.
Сингулярное разложение матрицы М позволяет найти псевдообратную ей матрицу М*, умножение которой слева на вектор N дает вектор искомых координат объекта, соответствующий минимальному значению СКО:
М* = .
Тогда искомая оценка Ре глобальных координат объекта может быть найдена в виде вектора
Ро = М * N .
Анализ погрешностей. Предложенные методы позволяют численно рассчитывать погрешности измерений в стереоскопической системе. Рассмотрим алгоритм скрещивающихся прямых для случая наиболее простой структуры стереоскопической системы, когда оптические оси камер параллельны горизонту, находятся в вертикальной плоскости, а глобальная система координат совмещена с системой координат нижней камеры.
На рис. 3 представлен результат стохастического моделирования измерений, когда номинальные значения координат объекта хо=уо=150 мм, а расстояние хо изменялось в диапазоне 2000—10 000 мм. На графике отражена зависимость СКО ахо измерения координаты хо объекта от дистанции при СКО измеренных координат изображения в плоскостях анализа камер <5Хр=<5ур=0,1р [5], где р=2,775 мкм — размер пиксела. При этом были приняты следующие значения параметров системы: база (см. рис. 1 и 2) Ъ=300 мм, а'1=а'ц=35 мм.
3
2
1
2000 4000 6000 8000
Рис. 3
Для стохастического моделирования погрешностей стереоскопической системы с использованием метода сингулярного разложения воспользуемся работами [6, 7], в которых доказывается возможность применения теории возмущений при анализе погрешности решения переопределенной системы линейных уравнений. В соответствии с этими работами погрешность оценки искомых координат объекта можно представить в виде
АРе = Ем ЛМ + Е N AN,
Методы обработки информации и оценки погрешностей пространственного положения объектов 15
где АРо = Ро - Ро, ДМ — вектор погрешностей коэффициентов матрицы М системы линейных уравнений (3), ДN — вектор погрешностей свободных членов той же системы уравнений.
Матрицы Ем и описывающие влияние ошибок коэффициентов системы уравнений (3) на решение в виде вектора координат объекта, можно найти из выражений
Ем = V (||Р0||) Wг, Е N = УЕ-1ит,
где ||Ре || — норма вектора Ре .
Матрица W состоит из столбцов, элементы которых рассчитываются по формуле
(Ра ) • ,
^ = НР 11 , / = 1, P,
11Ро II
где р — ранг матрицы М, а и^ — элементы столбцов матрицы и.
Рассмотренный метод позволяет оценить погрешности первого порядка, т.е. погрешности измерения координат объекта, которые описываются первыми производными коэффициентов матрицы М по исследуемым параметрам.
На рис. 4 показана зависимость СКО измерения расстояния до объекта в зависимости от СКО измерения координат изображений на матрицах. Стохастическое моделирование было проведено для двух методов обработки измерительной информации (кривая 1 — для алгоритма сингулярного разложения, 2 — для алгоритма скрещивающихся прямых). СКО измерения координат изображений схр и сур изменялись в диапазоне 0—100р, номинальное положение объекта было задано координатами хо=уо=150 мм, га=10 000 мм. Из рис. 4 видно, что зависимость для алгоритма, основанного на теории возмущений (кривая 1), носит линейный характер. Разница между зависимостями 1 и 2 становится существенной при очень больших значениях охр и оур, которые никогда не имеют места на практике. На практике алгоритм поиска энергетического центра изображения обеспечивает СКО порядка 0,1 пиксела [5]. При этом разница СКО, вычисленных с использованием двух методов, не превышает 2-10-6 мм для ха, 4-10-6 — дляуо и 2-10-4 — для го, что для измерений в диапазоне 2000—10 000 мм несущественно.
Ыо
3000
2000
1000
0 0,05 0,1 0,15 0,2 0,25 схо, схо
Рис. 4
Результаты исследования. Предложен метод сингулярного разложения матрицы коэффициентов переопределенной системы линейных уравнений для расчета координат объекта в стереоскопических и гиперстереоскопических оптико-электронных системах.
Обосновано применение теории возмущений для анализа систематических и случайных погрешностей измерений в стереоскопической системе. Такой подход удобен вследствие
1 V*"" * ■
Л—/" /V 2
использования матрицы влияния ошибок, которая зависит только от конфигурации стереоскопической системы и описывает, как ошибка оценивания того или иного параметра влияет на погрешность результата измерений. Данный алгоритм может найти применение для анализа влияния числа камер на погрешность измерения положения объекта в системе с п числом измерительных преобразователей. Метод исследования влияния случайной ошибки в стереоскопической системе [8] позволяет обойтись без стохастического моделирования при анализе случайных погрешностей измерений.
Работа осуществлялась при финансовой поддержке Федерального агентства по науке и инновациям РФ и была проведена в рамках аналитической ведомственной целевой программы „Развитие научного потенциала высшей школы (2009—2010 гг.)" и Федеральной целевой программы „Научные и научно-педагогические кадры инновационной России" на 2009—2013 гг.
список литературы
1. Arakantsev K. G., Konyakhin I. А., Timofeev А. N. Precision system for motion path parameters measurement of wheel and rail transport // Intern. Symp. on Instrum. Sci. and Techn. Institute of Physics Publishing. J. of Phys. 2006. Vol. 48. P. 998—1002.
2. Дуда Р., Харт П. Распознавание образов и анализ сцен. М.: Мир, 1976. 521 с.
3. Форсайт Д., Понс Ж. Компьютерное зрение. Современный подход. М.: Вильямс, 2004. 928 с.
4. Тыртышников Е. Е. Матричный анализ и линейная алгебра. М.: ФИЗМАТЛИТ, 2007. 480 с.
5. ЯкушенковЮ. Г. Проектирование оптико-электронных приборов. М.: Машиностроение, 1981. 263 с.
6. Wedin P., Wikstrom G. First Order Error Analysis of a Linear System of Equations by use of Error Propagation Matrices connected to the Pseudo Inverse Solution // Proc. of the Householder Symp. XV. 2005. 31 р.
7. Wedin P. Perturbation Theory for Pseudoinverses // BIT. 1973. Vol. 13. P. 217—232.
8. Haralick R. M. Propagating Covariance in Computer Vision [Электронный ресурс]: <http://www.citr.auckland. ac.nz/techreports/1998/CITR-TR-29.pdf>.
Сведения об авторах
Константин Геннадьевич Араканцев — аспирант; Санкт-Петербургский государственный университет
информационных технологий, механики и оптики, кафедра оптико-электронных приборов и систем; младший научный сотрудник; E-mail: kostya3312@mail.ru, kostya3312@yandex.ru Дмитрий Валерьевич Жуков — студент; Санкт-Петербургский государственный университет
информационных технологий, механики и оптики, кафедра оптико-электронных приборов и систем; E-mail: dzhukov@gmail.com Игорь Алексеевич Коняхин — д-р техн. наук, профессор; Санкт-Петербургский государствен-
ный университет информационных технологий, механики и оптики, кафедра оптико-электронных приборов и систем; E-mail: igor@grv.ifmo.ru
Рекомендована кафедрой Поступила в редакцию
оптико-электронных приборов и систем 26.04.10 г.