CC BY
217
МЕТОД ОПРЕДЕЛЕНИЯ СУММАРНОЙ ДИСТОРСИИ ЦИФРОВЫХ ИЗОБРАЖЕНИЙ
Жимбуева Л.Д.
Восточно-Сибирский государственный технологический университет
Аннотация
Описывается метод определения суммарной дисторсии цифровых изображений. Разработан алгоритм решения задачи, который состоит из следующих этапов: создание идеального изображения тестового объекта по его реальному изображению, построение векторной диаграммы дисторсии изображения, нелинейный регрессионный анализ аппроксимируемых функций, решение задачи интерполяции, корректировка изображения детали и восстановление плоского контура детали. В качестве тестового объекта использована квадратная сетка. Модуль «Нелинейное оценивание системы STATISTICA» позволяет разработать математическую модель дисторсии изображения. При использовании полинома третьего порядка величина достоверности аппроксимации исследуемой функции составляет 0,98 - 0,99. Разработанный метод позволяет повысить точность восстановления контуров плоских деталей, что актуально при использовании систем технического зрения в интеллектуальных ро-бототехнических системах, при распознавании объектов и т.д.
Ключевые слова: аберрация изображений, дисторсия, нелинейный регрессионный анализ, корректировка изображения, точность восстановления контура.
Введение
Оснащение робототехнических комплексов системами технического зрения обеспечивает существенное расширение их технологических и функциональных возможностей, наделяет их адаптивными возможностями и элементами искусственного интеллекта [1 - 2]. Промышленные системы технического зрения широко используются для бесконтактного контроля геометрических размеров, поиска и обнаружения самых разнообразных дефектов продукции бесконтактным способом, позволяют заменить человека на операциях сортировки и отбраковки продукции, полностью автоматизировать контроль сборки деталей в процессе производства, контроль качества упаковки готовой продукции. Системы технического зрения и видеодатчики являются одним из наиболее эффективных инструментов в задачах идентификации и автоматизации учёта продукции в производственном процессе.
В системах технического зрения можно использовать относительно недорогие цифровые камеры. Конечно, такие камеры по точности результатов не в состоянии конкурировать со специальной аппаратурой, используемой в фотограмметрии, но имеется большой круг задач, где они могут успешно применяться.
При реализации любой методики с использованием реальной камеры требуется знание внутренних и внешних параметров калибровки камеры. Паспортные данные таких сведений не содержат.
Одним из значимых внутренних параметров калибровки камеры является дисторсия цифрового изображения, которая в камерах, не предназначенных для измерительных целей, может достигать заметных величин.
Измерению дисторсии изображения посвящено множество работ [3 - 31]. Условно их можно разделить на две группы. Первая группа работ посвящена аттестации оптических систем на предмет дистор-сии при расчёте оптических систем и оценке их качества. Вторая группа работ посвящена задаче ка-
либровки камеры с целью определения её элементов внутреннего ориентирования (ЭВО), в том числе дисторсии изображения.
В первой группе нормативным документом является стандарт [3], устанавливающий два метода измерения дисторсии съёмочных объективов, один -для измерения дисторсии на оптической скамье для предмета, расположенного в бесконечности, и второй - для измерения дисторсии объективов, рассчитанных для работы с конечного расстояния.
При аттестации и юстировке объективов на предмет дисторсии важно разделение дисторсии на центрированную и нецентрированную и аппроксимация дисторсии полиномами различных видов. В большинстве случаев рассматривается математический аппарат обработки измерений координат изображений тестовых марок, и расчёт коэффициентов глобальной аппроксимации дисторсии по полиномам Цернике, известным в оптике как базис разложения волновых аберраций. Например, в [4] предусмотрено получение до 98 коэффициентов глобального разложения волновой аберрации, вызванной дисторсией. При разложении волновой аберрации по полиномам Цернике коэффициенты этого разложения, как известно [5], показывают непосредственный вклад аберрации каждого вида в среднеквадратичную волновую аберрацию. Для определения коэффициентов дисторсии по Цернике полученные после измерений составляющие дисторсии трактуются в качестве поперечных аберраций как частные производные волновой аберрации [6].
Известны методы аттестации и юстировки объективов на предмет дисторсии, основанные на измерении искажений квадратной сетки [7] по двум ортогональным осям. Формируемые при этом элементы двух массивов чисел ДА'и А У суммируют в себе расчётную дисторсию и дисторсию, вызванную разными типами погрешностей изготовления и сборки деталей объектива, погрешностями установки объ-
ектива при контроле дисторсии, погрешностями измерений. Это обстоятельство не позволяет эффективно использовать массивы АХ, А У непосредственно для аттестации и юстировки объективов.
В задаче калибровки камеры это обстоятельство наоборот имеет преимущество, так как учитывает суммарную погрешность измерений, независимо от причин их возникновения. Суммарные погрешности искажений изображения многими авторами рассматриваются как систематическая погрешность объектива камеры.
При калибровке камеры в основном используются тестовые объекты, как плоские, так и пространственные. Отмечается, что для определения ЭВО камеры достаточны плоские тестовые объекты. По искажениям изображений тестовых объектов формируются массивы для измерения дисторсии камеры. Сейчас имеется опыт успешных аэросъёмок цифровой камерой и попытки определения параметров цифровых камер [8 - 11].
Известны способы калибровки аэрофотографической системы с использованием снимков звёзд [12]. Преимущества этого метода - в простоте реализации, в отсутствии необходимости дополнительного оборудования, кроме точных часов. Основным недостатком метода является то, что условия калибровки отличаются от реальной съёмки. Этот факт ограничивает применение метода калибровки по снимкам звёзд.
Более эффективным по сравнению с предыдущими методами выглядит метод калибровки камеры по снимкам специального стенда или испытательного полигона. Данный метод калибровки является наиболее распространённым и легко реализуемым на практике [13 - 19]. Сущность метода - получение калибровочных поправок, учитывающих влияние всех источников систематических искажений, существующих при реальной съёмке. При этом дистор-сия может описываться полиномами Эбнера, Якобсона, Грюна, Брауна и т.д. [19].
Использование данного метода калибровки даёт заметное повышение точности измерений. Сложность практической реализации данного метода связана с необходимостью создания дорогостоящих полигонов в различных областях страны. И до настоящего времени не сформулированы оптимальные параметры калибровки камеры по тестовому полигону.
В работах других авторов [19 - 27] решение задачи калибровки с необходимой точностью достигалось путём вывода новых математических зависимостей и алгоритмов. Так, в [21] поле дисторсии определяется по изображению плоского тест-объекта и описывается полиномом третьего порядка. Полиномы Цернике использованы в работе [27] для разделения направления децентрировок. Затем с учётом порядка разложения сформированы наборы косинусных и синусных полиномов, что позволяет оценить присутствие в оптической системе ошибок сборки и изготовления.
Также известны работы, в которых описывается аналитический метод самокалибровки [12, 28]. Сущ-
ность метода состоит в том, что одновременно определяются аналитически поправки к приближённым элементам внешнего ориентирования, координатам точек местности, элементам внутреннего ориентирования и суммарные поправки к координатам точек за влияние всех источников систематических ошибок.
Преимущество данного метода перед методом испытательного полигона состоит в том, что он не требует специального полигона, дополнительных измерений. Поправки в координаты определяются для снимков, по которым строится сеть фототриангуляции, а не по специальным снимкам полигона.
Очевидно, что метод самокалибровки - эффективный способ учёта систематических ошибок снимка, но нет единодушных практических рекомендаций. В работах представлены самые разнообразные полиномы с количеством членов до 20 и со степенью от 1 до 12. В работе [12] отмечается, что нахождение оптимального полинома для описания систематических искажений снимка невозможно, что полином должен формироваться внутри программы при решении конкретной задачи, а не даваться неизменным извне.
Несмотря на многообразие полиномов [12], все они составлены по двум основным идеям. Первая -интегральная систематическая ошибка ф, описываемая полиномом, представлена как сумма членов отдельных систематических искажений (дисторсии, деформации и т. д.).
Вторая идея состоит в том, что представление полиномом интегральной ошибки не связывается с отдельными видами искажений. Интегральная ошибка описывается степенным многочленом.
Этот способ представления дополнительных параметров, вероятно, более правильный, как отмечается в работе [12], так как не связан с определёнными видами искажений. Полиномы же первого вида могут не учесть искажений, которые ещё не изучены или просто не описаны математически.
Другая проблема, с которой сталкиваются исследователи при выборе дополнительных параметров, -это корреляция коэффициентов полинома между собой и другими неизвестными, элементами внешнего ориентирования, координатами точек местности. Это снижает точность фототриангуляции, а иногда сводит на нет попытки увеличения точности за счёт использования самокалибровки. Эта трудность может быть преодолена путём отбраковки сильно коррелированных коэффициентов в процессе уравнивания сети фототриангуляции.
Не сложилось ещё единого мнения по выбору оптимальной длины полинома. Здесь много противоречий. С одной стороны, увеличение количества членов полинома детализирует описание систематических искажений, с другой стороны, ухудшается обусловленность матрицы нормальных уравнений, загружается память ЭВМ. Очевидно, что существуют определённые пределы количества членов полинома, обусловленные исходными материалами, а также памятью ЭВМ. Поэтому длина и вид полинома могут изменяться для различных исходных данных.
Изложенные в работах [6 - 28] способы определения оптических аберраций камеры в основном реализованы в виде программных комплексов.
В настоящее время в программных комплексах реализованы два основных способа: полиномиальный, например, работы [19, 21], и способ зон [29]. Для исправления изображения с помощью зон осуществляется поиск формулы преобразования в виде нелинейной функции Е от наблюдаемых координат. Предлагается разбивка поля изображения на отдельные клеточные области с предположением линейности закона внутри каждой клетки.
Кроме того, на рынке присутствует больше десятка цифровых фотограмметрических станций (ЦФС) («Дельта» - Н1111 «Геосистема»; Photomod - компания «Ракурс»; «ЦФС Талка», «Талка ГИС», «Талка КПК» - компания «Талка» и ряд других). Каждая из данных программ обладает своими достоинствами и недостатками.
Область применения ЦФС - фотограмметрия, где используются очень дорогие и точные камеры для формирования изображений, в частности, это космические изображения, предназначены они в основном для картографических работ и обеспечивают возможности вР8-навигатора [30].
Наиболее близким по назначению для нашей работы является программное обеспечение РТЪепБ [31]. Оно предназначено для исправления искажений оптики цифровых камер. Чтобы исправить искажение, программа считывает ЕХШ-информацию из ДРЕв-файла, находит ближайшие соответствия в базе данных калибровок камер и автоматически исправляет искажения. Результаты использования программы при корректировке контура цифрового изображения квадратной сетки показали, что программа довольно точно исправляет искажения отдельных линий сетки, но неточно определяет местоположение линии на изображении. Поэтому программа РТЬеш, предназначенная в основном для корректировки изображений архитектурных сооружений, не приемлема для калибровки камеры.
Таким образом, нужно отметить, что в задачах аттестации и юстировки объективов важно разделение дисторсии на центрированную и нецентриро-ванную, определение причин их возникновения. В задаче же калибровки камеры при определении её внутренних параметров необходимо учитывать суммарную дисторсию, вызванную разными типами погрешностей изготовления и сборки деталей объектива, погрешностями установки объектива при контроле дисторсии, погрешностями измерений, погрешностями установки ПЗС-матрицы камеры. При этом суммарные погрешности искажений изображения надо рассматривать как систематическую погрешность объектива камеры.
Дисторсия изображений
Целью данной работы является разработка метода определения суммарной дисторсии цифровых изображений независимо от причин их возникновения, с дальнейшим исследованием зависимостей
дисторсии от фокусного расстояния камеры, расстояния от камеры до исследуемого объекта.
Дисторсия характерна тем, что её величина нелинейно зависит от величины предмета, и приводит к искажению прямых линий, не проходящих через ось [32 - 35]. При этом квадратный предмет изображается в виде подушки при положительной дистор-сии и в виде бочки при отрицательной дисторсии.
Для вычисления дисторсии необходимо определить точку идеального изображения, в которой должно находиться изображение по законам гауссовой оптики. Если кроме дисторсии нет других аберраций, то точка А трёхмерного пространства изображается в виде точки А (а1, У), не совпадающей с точкой идеального изображения А'0 (Х0, у0) (рис. 1).
Рис. 1. Дисторсия
Тогда отклонения координат точки А - точки пересечения реального луча с плоскостью изображения - от координат точки А'0 идеального изображения и есть дисторсия ДУ, Ду соответственно по осям X и у [32]:
Дх = х 0 — х, (1)
Ду = У о— У'-
Величина дисторсии различна для разных точек поля изображения и является функцией от реальных координат точек изображения хр, ур:
Дх' = Дх'(хр, Ур X
ДУ' = ДУ'(хр, Ур).
(2)
Допустимая относительная дисторсия (т.е. дис-торсия, при которой нет ощущения, что изображение искажено) около 5-10%. Величина максимально допустимой дисторсии в оптических системах зависит от области их применения. Например, в оптической литографии допуск на абсолютную дисторсию не должен превышать 20 нм [32]. В фотограмметрии погрешность измерений не должна превышать 1 мм. Недорогие цифровые камеры со значительной дис-торсией не предназначены для измерений, потому для них нет требований по погрешностям измерений. Тем не менее, в системах технического зрения довольно широко стали использовать неметрические камеры, и известны методы калибровки камер, позволяющие значительно уменьшить дисторсию до 0,3 - 1 пикселя.
Общее описание метода
В приосевой, так называемой параксиальной области оптическая система близка к идеальной. Поэтому в центральной области реального изображения искажения близки к нулю. В этой области точка изображается точкой, прямая линия - прямой. На основе анализа структуры центральной области реального изображения возможно формирование идеального изображения объекта.
В нашей работе определение суммарной дистор-сии цифровых изображений основано на измерении искажений квадратной сетки (далее в тексте просто сетки) по двум ортогональным осям.
Сформированные при этом массивы чисел Да1, ДУ позволяют разработать математическую модель нелинейной зависимости дисторсии от реальных координат узловых точек изображения сетки. Алгоритм решения задачи следующий:
- восстановление контура реального изображения сетки;
- создание идеального изображения сетки;
- определение отклонений узловых точек сетки на реальном изображении от соответствующих точек на её идеальном изображении и построение векторной диаграммы дисторсии изображения;
- представление функций Да1, Ay дисторсии в матричной форме;
- нелинейный регрессионный анализ аппроксимируемых функций методом наименьших квадратов, подбор вида функции, достоверно описывающей характер изменения дисторсии;
- вычисление значений функций Да1, Ay в точках, отличных от узловых, т. е. решение задачи интерполяции;
- корректировка реального изображения детали с учётом дисторсии и восстановление плоского контура детали.
Создание идеального изображения сетки
Реальный контур сетки представляет собой линии опредёленной толщины, причём горизонтальные линии сетки имеют верхнюю и нижнюю, вертикальные - левую и правую границы. Восстановление граничных линий сетки выполняется сравнением яркостей пикселей изображения с выбранным значением порога бинаризации.
Далее определяются осевые линии между соответствующими граничными линиями и на пересечении вертикальных и горизонтальных осевых линий определяются узловые точки сетки.
Для восстановления контура реального изображения сетки использованы алгоритмы обработки изображений [34], на основе которых разработана программа в среде Delphi.
Контур изображения сетки в его центральной области зависит от взаимного расположения плоскостей изображения и предмета. Если плоскости параллельны, то перспективные искажения на изображении отсутствуют. Изображение будет представлять собой в общем случае прямоугольную сетку (если коэффици-
енты масштабирования разные по осям X и У). Тогда по размерам прямоугольной сетки довольно легко сформировать идеальное изображение сетки.
Если плоскости не параллельны, то проявляются перспективные искажения. Центральная часть реального изображения будет представлять собой четырёхугольники с перспективными искажениями сторон. В этом случае для формирования идеального изображения сетки можно воспользоваться точками схода противоположных сторон четырёхугольника и гармоническими свойствами полного четырёхугольника [36].
Векторная диаграмма дисторсии изображения
На рис. 2 приведена векторная диаграмма дис-торсии изображения, которая представляет собой совокупность векторов из узловых точек реального изображения в соответствующую ей узловую точку на идеальном изображении сетки. Как видно из рисунка, направления и модули векторов погрешностей различны для разных точек поля, причём характер изменения нелинейный.
Yv V V- т- г- = Ч 7 <
+ + 4- 1" - ■н -i -К ■Y
1
1 1
1 1 1
Г Г
1 1
1 1
1 1 1 Г
1 1
rl-ir 1
гМ- л Г 1 4
iH- н т +
я™ т «— — — -» Н т -к т "К + ■к
ft 3- ¡L г п — — — — -Ч Гг X ч S
Рис. 2. Векторная диаграмма дисторсии изображения В матричной форме функции поперечной аберрации X, У будут иметь вид:
х ° Ы, у° Ы, (3)
где Ху, у у - элементы матриц, представляющие величины погрешности в у узловой точке сетки. Размер прямоугольной матрицы - тхп, где т, п - соответственно количество горизонтальных и вертикальных линий сетки.
Для практического использования функции, представленной в матричной форме (3), необходимо выполнить аппроксимацию функций Дх1, Ду, зависящих от координат узловых точек реального изображения сетки. Как правило, регрессионный анализ данных производится методом наименьших квадратов [37].
Нелинейный регрессионный анализ исследуемой зависимости Классические методы аппроксимации и регрессии функции реализованы во всех современных матема-
тических системах (MathCAD, MatLAB, Mathematica, STATISTICA и пр.), которые имеют в своём составе универсальный аппарат выполнения таких операций, и дают пользователю возможность реализации достаточно сложных практических задач по обработке данных. В работе для оценивания параметров регрессии использован программный модуль «Нелинейное оценивание системы STATISTICA», который оставляет за пользователем выбор характера зависимости, что немаловажно для подбора наиболее оптимального вида функции [38 - 40]. Для исследуемого процесса в качестве оптимальной функции был подобран полином третьего порядка:
2 2
z = a + aix+ ^у+ a3 x+ a4y + щ xy+
+a5 x2 у + a5 xf + a5 x3 + a5 у3.
Дальнейшее увеличение степени полинома не дало ощутимого эффекта в повышении точности обработки данных.
Корректировка изображения
Корректировка изображения предмета заключается в переопределении координат точек контура реального изображения предмета и последующей за ним аппроксимации массива переопределённых точек для восстановления контура предмета. Для определения искажения координат точек изображения предмета решается задача интерполяции, т. е. вычисление значений функций Ax1, Ду в точках, отличных от узловых точек сетки, но принадлежащих контуру изображения предмета. Координаты точек контура скорректированного изображения предмета определяются по формуле:
x '0 = x'+ Ax',
0 . А (5)
Уо = у + ду.
Экспериментальные исследования
Для экспериментальных исследований применялась система технического зрения, в составе которой использована камера марки Sony DSC-H50. Были проведены несколько серий экспериментов, в которых варьировались размеры сетки, фокусное расстояние и разрешение камеры, расстояние от камеры до сетки. В ходе экспериментов выявлено, что дис-торсия цифрового изображения значимо зависит от величины фокусного расстояния и незначимо - от расстояния до предмета. Полученные результаты позволяют оптимизировать процесс формирования цифрового изображения и при необходимости переопределить координаты точек контура изображения. Более подробно результаты изменения дисторсии от условий съёмки будут изложены в другой статье.
В данной же работе важно принципиальное подтверждение результатами экспериментов корректности разработанного метода в целом. Результативность метода оценивается на примере обработки изображений объектов, приведённых на рис. 3: сетки с расстоянием в 10 мм между линиями и правильного шестиугольника с длиной стороны 80 мм.
Рис. 3. Объекты исследований
Полутоновое цифровое изображение сформировано в формате *.jpeg. Плоскости изображения и предмета параллельны. На изображении наблюдается отрицательная дисторсия (рис. 4).
/
/
/
/ \
/ \
с \
\
! \
\ \
\
\ / 7*-
т-
у-
у-
7*"
Рис. 4. Реальное изображение исследуемых объектов
По цифровому изображению вначале восстанавливается контур реального изображения, затем создаётся идеальное изображение сетки. По отклонениям координат соответствующих узловых точек реального и идеального изображений формируются два массива исходных данных X, У в матричной форме.
Аппроксимация экспериментальных данных полиномом третьей степени показала довольно высокую точность, величина достоверности 1=0,98 - 0,99. График аппроксимирующей функции АХ представлен на рис. 5.
Используемый в работе аппроксимирующий полином третьей степени не накладывает ограничения на линейное или угловое поле зрения тестируемого объектива.
В табл. 1 приведены результаты расчёта координат вершин изображения шестиугольника в пикселях.
На рис. 6 приведены контуры шестиугольника до и после корректировки с учётом суммарной дистор-сии изображения на фоне реального и идеального контуров сетки. Внутренний контур (пунктирная линия) - контур изображения шестиугольника до корректировки, причём стороны шестиугольника имеют отклонения от прямых линий. Наружный
контур (сплошная линия) - контур после корректировки представлен линейными сегментами.
Рис. 5. График аппроксимирующей функции Ах'
Таблица 1. Результаты корректировки положений вершин шестиугольника в пикселях
№ в е р 1 Реальные координаты Погрешности Скорректированные координаты
x1 829,2 У 1407,7 АХ 5,23 ДУ -11,87 Хо 823,97 Уо 1419,57
2 1484,6 1234,5 -3,69 -5,83 1488,25 1240,28
3 1653,7 582,9 -11,17 5,46 1664,82 577,40
4 1173,5 110,4 -5,14 18,10 1178,60 92,30
5 524,9 284,3 8,85 12,87 516,01 271,44
6 350,2 930,1 12,88 -4,28 337,32 934,36
г т т T T- r- Г" 1" Г" -] 1 1 Г Г T H
L -11 4- F I F 1 1
I- -t- 1- I +— I —i i 1 -+ -h
1 т f 'Г I | i ■ i i | "Г
1 1
\
1 1 I I f \ -t -1
I* i 1 1
j ' \
и 1 J_ I H- .4 — .1 1 -1
h н- I I— I \ ■ 1 r\ 1 1 i .4 1 4
l 1 ' (_ [ 1 1
L -1- l __| J
h -t- __ F 1 -1- -1
I- т T 4— __ _L ■J- -+ -h 1 -t- 1 7
L J. -L. — __ L- _ — __i —4 _ _L J- J
Рис. 6. Контуры шестиугольника до и после корректировки
Результаты расчёта длин сторон шестиугольника до и после корректировки с учётом суммарной дис-торсии цифрового изображения в мм приведены в табл. 2.
Из табл. 2 видно, что при восстановлении контура правильного шестиугольника по реальному изображению максимальная погрешность при определении длин сторон составляет ~2,07 мм, после корректировки изображения она уменьшается до ~0,10 мм.
Это свидетельствует о значительном повышении точности восстановления контура шестиугольника после корректировки изображения.
Таблица 2. Результаты расчёта длин сторон шестиугольника в мм
Заключение
В задаче калибровки камеры при определении её внутренних параметров необходимо учитывать суммарную дисторсию, вызванную разными типами погрешностей изготовления и сборки деталей объектива, погрешностями установки объектива при контроле дисторсии, погрешностями измерений, погрешностями установки ПЗС матрицы камеры. При этом суммарные погрешности искажений изображения надо рассматривать как систематическую погрешность объектива камеры.
Использование квадратной сетки в качестве тестового объекта по сравнению с другими объектами (геометрическими фигурами) имеет большое преимущество, так как более полно иллюстрирует нелинейный характер погрешностей оптических аберраций по всему полю изображения. Уменьшая расстояние между линиями сетки, можно увеличить точность аппроксимации функции погрешностей А Я, Ду.
Создание идеального изображения по реальному изображению тестового объекта позволяет формировать статистические данные для аппроксимации функций погрешностей Ах1, Ау.
Для разработки математической модели зависимости погрешностей оптических аберраций от координат узловых точек изображения можно использовать математические системы, имеющие в своём составе универсальный аппарат для решения задач аппроксимации и регрессионного анализа.
Использованный в работе модуль «Нелинейное оценивание системы STATISTICA» показал довольно высокие результаты аппроксимации исследуемой функции полиномом третьего порядка. Величина достоверности аппроксимации R = 0,98 - 0,99.
Результаты исследований свидетельствуют о значительном увеличении точности восстановления контуров плоских деталей после корректировки изображения с учётом дисторсии объектива камеры.
Проблема повышения точности восстановления контура исследуемых объектов имеет важное прикладное значение. Например, позволит увеличить точность калибровки систем технического зрения, что немаловажно для интеллектуальных робототех-нических комплексов, особенно при решении измерительных задач.
№ Длина стороны шестиугольника, мм
до корректировки после корректировки
1 78,91 80,09
2 78,42 79,91
3 78,44 79,98
4 78,17 79,90
5 77,93 79,98
6 78,77 80,02
Рассматриваемая задача важна также и для фундаментальных исследований при решении задач сегментации изображений и распознавания образов.
Благодарности Работа выполнена при поддержке Министерства образования и науки в рамках НИР «Исследование геометрических вопросов интеллектуальных робо-тотехнических комплексов» по Единому заказ-наряду, научный руководитель д.т.н., проф., проректор Восточно-Сибирского государственного технологического университета Найханов В.В.
Литература
1. Хорн, Б.К.П. Зрение роботов: / Б.К.П. Хорн, пер. с англ. - М.: Мир, 1989. - 487 с.
2. Мошкин, В.И. Техническое зрение роботов / В.И. Мошкин, В.С. Титов, Ю.Г. Якушенков; под общ ред. Ю.Г. Якушенкова. - М.: Машиностроение, 1990. -272 с.
3. ГОСТ 20825-75. Объективы съёмочные. Методы измерений дисторсии. - Введ. 01.07.76. - М.: Издательство стандартов. - 12 с.
4. Родионов, С.А Автоматизация проектирования оптических систем / С.А. Родионов - Л.: Машиностроение, 1982.
5. Вычислительная оптика: Справ. / под общей ред. М.М. Русинова. - Л.: Машиностроение, 1984.
6. Родионов, С.А. Обработка результатов измерения дисторсии проекционных объективов / С.А. Родионов, Н.Б. Вознесенский, Э.М. Щекольян // Изв. вузов. Приборостроение, 1991. - Т. XXXIV, № 7. - С. 61-68.
7. Tatian, B. Aberration balancing in rotationally symmetric lenses / B. Tatian // JOSA. - 1974. - Vol. 64. - P. 1083.
8. Гонин, Г.Б. Экспериментальная цифровая аэросъёмка / Г.Б. Гонин, В.И. Микеров // Геодезия и картография, 1997. - № 9. - С. 39-42.
9. Шоломицкий, А.А. Цифровые съёмки открытых горных разработок / А.А. Шоломицкий // Физико-технические проблемы горного производства: сб. н. тр., -Вып. 2. - Донецк, 1999. - С. 127-133.
10. Гельман, Н.Р. Опыт использования и калибровки цифровых камер при совместной аэрофотосъёмке с АФА / Н.Р. Гельман, Н.Ю. Никитин, А.Л. Дунц // Геодезия и картография. - 2001. - № 6. - С. 25-31.
11. Гельман, Н.Р. Ещё об оценке точности одиночной фотограмметрической модели / Н.Р. Гельман // Геодезия и картография. - 2000. - № 12. - С. 19-23.
12. Курков, В.М. Методы учёта систематических искажения аэроснимка. Самокалибровка / В.М. Курков // Изв. вузов. Геодезия и аэрофотосъёмка. - 1980. - № 6. - С. 75-79.
13. Дубиновский, В.Б. Калибровка снимков / В.Б. Дуби-новский - М.: Недра, 1982. - 224 с.
14. Шоломицкий, А.А. Создание испытательного стенда для калибровки цифровых камер / А.А. Шоломицкий, А.С. Шатохин, [Электронный ресурс]: - Режим доступа: http://info.dgtu.done1sk.ua/el izdan/geolog/sborniki/sbornik
ggf N45 2002/M4.pdf
15. Шоломицкий, А.А. Новая технология измерения стенда для калибровки цифровых камер / А.А. Шоломицкий, А.А. Лунев // Проблеми прського тиску. - 2005. -№ 1. - С. 333.
16. Лунев, АА Выбор оптимальных параметров калибровки цифровой камеры / А.А. Лунев [Электронный ресурс]: Ре-
жим доступа: http://www.info.donntueduua/el izdan/geolog/ sborniki/ggfl 11t2.pdf
17. Matsuoka, R. A new calibration system of a non_metric digital camera / R. Matsuoka [et al.] // Procs. 6th Conference on Optical 3D Measurement Techniques, Zurich, Switzerland, September 22-25, 2003. - P. 130-137.
18. Технология фотограмметрической калибровки цифровых камер. - 4 мая 2008 [Электронный ресурс]: Режим доступа: http ://photo grammetria.ru/17-tekhno lo gija-foto-grammetrichesko j j -kalibro vki. html.
19. Комиссаров, Д.В. Методика калибровки цифровых неметрических камер для наземных лазерных сканеров / Д.В. Комиссаров, А.В. Комиссаров, [Электронный ресурс]: 2006. - Режим доступа: http://www.geoprofi.ru/de-fault.aspx?mode=binary&id=671.
20. Юрченко, В.И. Способ аналитической обработки неметрических снимков / В.И. Юрченко // Геодезия и картография. - 2000. - № 11. - С. 23-30.
21. Гельман, Р.Н. Лабораторная калибровка камер с большой дисторсией / Р.Н. Гельман, А. Л. Дунц // Геодезия и картография. - 2002. - № 7. - С. 23-31.
22. Grenzdoerfer, G. Konzeption, Entwicklung und Erprobung eines digitalen integrierten flugzeuggetragenen Fernerkundungssystems fuer Precisin Farming (PFIFF). DGK Reiche C, Heft 552: Universitaet Muenchen / G. Grenzdoerfer, 2002. - P. 142.
23. Глотов, В. Особливоси визначення фокусно! ввдсташ цифрових фототеодолггних камер / В. Глотов // Геоде-зш, картографш i аерофотозшмання. - 2003 - № 63. -С. 122-127.
24. Zhang, Z. A Flexible New Technique for Camera Calibration. Technical Report MSRTR- 98-71, Microsoft Research, December 1998. Available together with the software at http://research.microsoftcom/zhang/Calib/ [Электронный ресурс]: режим доступа: http://research.microsoft.com/en-us/um/people/zhang/papers/TR98-71.pdf.
25. Труфанов, М.И. Адаптивное оптико-электронное устройство автоматической трёхмерной калибровки на основе нечёткой логики / М.И. Труфанов. - Дис. ... канд. техн. наук: 05.13.05. - Курск, 2005. - 148 с.
26. Remondino, F. Digital camera calibration methods: considerations and comparisons / F. Remondino, C. Fraser // ISPRS Commission V Symposium «Image Engineering and Vision Metrology». Volume XXXVI, Part 5, Dresden: 25-27 September 2006. - P. 266-272.
27. Ежова, К.В. Математическое моделирование фотограмметрической дисторсии / К.В. Ежова // Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики. - 2006. - Вып. 26. - С. 235-239.
28. Pollefeys, M. Stratified self-calibration with the modulus constraint / M. Pollefeys, L.V. Gool // IEEE Transactions on pattern analysis and machene intelligence. - August 1999. - Vol. 21, N 8. - P. 707-724.
29. Овчинников, А.М. Разработка алгоритмических и программных средств реализации и визуализации локальных гиперспектральных данных: Автореф. дис...канд. физ.-мат. наук: 05.13.11 / А.М. Овчинников. - М., 2009. - 18 с.
30. Работа с GPS [Электронный ресурс]: Режим доступа: http://gis.talka2000.ru/gisgps.htm
31. Сайт http ://nnm.ru/blo gs/SHAH/ptlens/
32. Родионов, С.А. Основы оптики: конспект лекций / С.А. Родионов. - СПб: СПб ГИТМО (ТУ), 2000. - 167 с. [Электронный ресурс]: Режим доступа: http://aco.if-mo.ru/el_books/basics_optics/index.html.
33. Заказнов, Н.П. Теория оптических систем: учеб. пособие для вузов / Н.П. Заказнов, С.И. Кирюшин, В.И. Кузичев - СПб.: Лань, 2008. - 446 с.
34. Гонсалес, Р. Цифровая обработка изображений / Р. Гон-салес, Р. Вудс. - Москва: Техносфера, 2005. - 1072 с.
35. Карнаухов, А.В. Многофункциональная цифровая модель системы искажения и восстановления изображений / А.В. Карнаухов, Н.С. Мерзляков, О.П. Милюкова // Компьютерная оптика. - 2000. - Т. 20. - С. 118-121.
36. Четверухин, Н.Ф. Проективная геометрия: учебник для пед. ин-тов. / Н.Ф. Четверухин, изд. 8-е. - М.: Просвещение, 1969.
37. Херн, Д. Компьютерная графика и стандарт OpenGL, 3-е изд.: / Дональд Херн, М. Паулин Бейкер. - пер. с англ. - М.: Издательский дом «Вильямс», 2005. -1168 с.
38. Боровиков, В.П. Прогнозирование в системе STATI-STICA в среде Windows. Основы теории и интенсивная практика на компьютере: учеб. пособие / В.П. Боровиков, Г.И. Ивченко. - М.: Финансы и статистика, 2000. - 384 с.
39. Боровиков, В.П. STATISTICA. Искусство анализа данных на компьютере: Для профессионалов. 2-е изд. (+CD) / В.П. Боровиков. - СПб.: Питер, 2003. - 688 с.
40. Халафян, А.А. STATISTICA 6. Статистический анализ данных: учебник, / А. А. Халафян - 3-е изд. - М.: ООО «Бином-Пресс», 2007. - 512 с.: ил.
References
1. Horh, B.K.P. Robot vision / B.K.P. Horh. - Moscow: "Mir" Publisher, 2005. - 487 p. - (in Russian).
2. Moshkin, V.I. Machine vision of robots / V.I. Moshkin, V.S. Titov, Yu.G. Yakushenkov; Under the general. ed. Yu.G. Yakushenkov. - Moscow: "Mashinostroenie" Publisher, 1990. - 272 p. - (in Russian).
3. GOST 20825-75. Lenses shooting. Methods of measurement distortion. - Intr. 07/01/1976. - Moscow: Publishing the standards. - 12 p. - (in Russian).
4. Rodionov, S.A Computer-aided design of optical systems. / S.A. Rodionov - Leningrad: "Mashinostroenie" Publisher, 1982. - (in Russian).
5. Computational optics: reference / Ed. M.M. Rusinov. - Leningrad: "Mashinostroenie" Publisher, 1984. - (in Russian).
6. Rodionov, S.A. Processing the measurement results of distortion of the projection lenses / S.A. Rodionov, N.B. Voznesensky, E.M. Schekolyan // Proceedings of the universities. Instrumentation. - 1991. - T. XXXIV, N 7. -P. 61-68. - (in Russian).
7. Tatian, B. Aberration balancing in rotationally symmetric lenses / B. Tatian // JOSA. - 1974. - Vol. 64. - P. 1083.
8. Gonin, G.B. An experimental digital aerial / G.B. Gonin, V.I. Mikerov // Geodesy and Cartography. - 1997. - N 9. - P. 39-42. - (in Russian).
9. Sholomitskii, A.A. Digital photography opencast mining / Physical and technical problems of mining industry // Proceedings of Issue 2. Donetsk. - 1999. - P. 127-133. - (in Russian).
10. Gelman, N.R. Experience in the use and calibration of digital cameras for aerial photography together with the AFA / N.R. Gelman, N.J. Nikitin, A.L. Dunts // Geodesy and Cartography. - 2001. - N 4. - P. 25-31. - (in Russian).
11. Gelman, N.R. More on assessing the accuracy of a single photogrammetric model / N.R. Gelman // Geodesy and Cartography. - 2000. - N 12. - P. 19-23. - (in Russian).
12. Kurkov, V.M. Methods to account for the systematic distortion of aerial photographs. Self-calibration / V.M. Kur-
kov // Proceedings of the universities. Surveying and aerial photographs. - 1980. - N 6. - P. 75-79. - (in Russian).
13. Dubinovsky, V.B. Calibration images / V.B. Dubinovsky.
- Moscow: "Nedra" Publisher, 1982. - 224 p. - (in Russian).
14. Sholomitskii, A.A. Creating a test bench for calibration of digital cameras / A.A. Sholomitskii, A.S. Shatokhin [Electronic resource]: - Mode of access: http://info.dgtu.do-netsk.ua/el_izdan/geolog/sborniki/sbornik_ggf_N45_2002 /M4.pdf. - (in Russian).
15. Sholomitskii, A.A. New techniques to measure the stand for the calibration of digital cameras / A.A. Sholomitskii, A.A. Lunev // Problems of Mining delagirskogo gripe. -2005. - N 1. - P. 333. - (in Russian).
16. Lunev, A.A. Selection of optimal calibration parameters digital camera / A.A. Lunev [Electronic resource]: Mode of access: http://www.info.donntu.edu.ua/el izdan/geo-log/sborniki/ggf111t2.pdf. - (in Russian).
17. Matsuoka, R. A new calibration system of a non_metric digital camera / R. Matsuoka [et al.] // Procs. 6th Conference on Optical 3D Measurement Techniques, Zurich, Switzerland, September 22-25, 2003. - P. 130-137.
18. Technology photogrammetric calibration of digital cameras. - May 4, 2008 [Electronic resource]: Mode of access: http ://photogrammetria.ru/17-tekhnologija-fotogram-metricheskojj-kalibrovki.html. - (in Russian).
19. Komissarov, D.V. Nonmetric technique for calibrating digital cameras for terrestrial laser scanners [Electronic resource] / D.V. Komissarov, A.V. Komissarov. - 2006. -Mode of access: http://www.geoprofi.ru/default.aspx?mo-de=binary&id=671. - (in Russian).
20. Yurchenko, V.I. Method analytical processing nonmetric images / V.I. Yurchenko // Geodesy and Cartography. -2000. - N 11. - P. 23-30. - (in Russian).
21. Gelman, R.N. Laboratory calibration of cameras with a large distortion / R.N. Gelman, A.L. Dunts // Geodesy and Cartography. - 2002. - N 7. - P. 23-31. - (in Russian).
22. Grenzdoerfer, G. Konzeption, Entwicklung und Erprobung eines digitalen integrierten flugzeuggetragenen Fernerkundungssystems fuer Precisin Farming (PFIFF). DGK Reiche C, Heft 552: Universitaet Muenchen / G. Grenzdoerfer, 2002. - P. 142.
23. Glotov, V. Features defining the focal distance of digital cameras fototeodolitnyh / V. Glotov // Geodesy, cartography and aerial photography. - 2003 - N 63. - P. 122-127.
- (in Russian).
24. Zhang, Z. A Flexible New Technique for Camera Calibration. Technical Report MSRTR- 98-71, Microsoft Research, December 1998. [Электронный ресурс]. - Available together with the software at htp://research.micro-soft.com/zhang/Calib/. Mode of access: htp://research.mi-crosoftcom/en-us/um/people/zhang/papers/TR98-71.pdf.
25. Trufanov, M.I Adaptive optical эlektronnoe a device for automatic three-dimensional calibration based on fuzzy logic / M.I. Trufanov [Electronic resource]: Thesis for the degree of candidate of technical sciences: 05.13.05. -Kursk, 2005. - 148 p. - (in Russian).
26. Remondino, F. Digital camera calibration methods: considerations and comparisons / F. Remondino, C. Fraser, // ISPRS Commission V Symposium «Image Engineering and Vision Metrology». - Volume XXXVI, Part 5, Dresden: 25-27 September 2006. - P. 266-272.
27. Ezhov, K.V. Mathematical modeling of photogrammetric distortion / K.V. Ezhov // Scientific and Technical Bulletin, St. Petersburg State University of Information Technologies, Mechanics and Optics. - 2006. - Vol. 26. -P. 235-239. - (in Russian).
28.
29.
30.
31.
32.
33.
34.
35.
Pollefeys, M. Stratified self-calibration with the modulus constraint / M. Pollefeys, L.V. Gool // IEEE Transactions on pattern analysis and machene intelligence. - August 1999. - Vol. 21, N 8. - P. 707-724.
Ovchinnikov, A.M. The development of algorithmic and software implementation and visualization of local hyper-spectral data: Summary of thesis for degree of Candidate of Physical and Mathematical Sciences: 05.13.11 / A.M. Ovchinnikov. - Moscow, 2009. - 18 p. - (in Russian). Working with GPS [Electronic resource]: Mode of access: http://gis.talka2000.ru/gisgps.htm. - (in Russian). Site Web http://nnm.ru/blogs/SHAH/ptlens/ - (in Russian). Rodionov, S.A. Optics bases. Lecture notes [Electronic resource] / S.A. Rodionov. - St.-Petersburg, 2000 - 167 p. - Mode of access: http://aco.ifmo.ru/el books/basics op-tics/index.html . - (in Russian).
Zakaznov, N.P. The theory of optical systems: The Manual for high schools / N.P. Zakaznov, S.I. Kiryashin, V.I. Ku-zichev. - St-Petersburg: Lan, 2008. - 446 p. - (in Russian). Gonzalez, R.C. Digital Image Processing / Rafael C. Gonzalez, Richard E. Woods. - Moscow: "Technosphere" Publisher, 2005. - 1072 p. - (in Russian). Karnaukhov, A.V. Multipurpose digital model of system of distortion and restoration of images / A.V. Karnaukhov,
N.S. Merzlyakov, O.P. Milyukova // Computer Optics. -2000. - Vol. 22. - P. 118-121. - (in Russian).
36. Chetverukhin, N.F. Projective geometry. The edition 8. The textbook for teacher training colleges / N.F. Chetverukhin. - Moscow: "Prosveshenie" Publisher, 1969. - (in Russian).
37. Hearn, D. Computer Graphics with OpenGL, 3rd ed.: Trans. from English / Donald Hearn, M. Pauline Baker. -Moscow: The publishing house «Vilyams», 2005. - 1168 p. - (in Russian).
38. Borovikov, V.P. Forecasting in system STATISTICA in Windows environment. Bases of the theory and intensive practice on the computer: the Manual / V.P. Borovikov, G.I. Ivchenko. - Moscow: "Finanse i statistics" Publisher, 2000. - 384 p. - (in Russian).
39. Borovikov, V.P. STATISTICA. Art analysis of the data on the computer: For professionals. 2 edition. (+CD) / V.P. Borovikov. - Saint-Petersburg: "Peter" Publisher, 2003. -688 p. - (in Russian).
40. Khalafyan, A.A STATISTICA 6. The statistical analysis of the data 3 edition, the Textbook / A.A. Khalafyan. - Moscow: "Binom-Press" Publisher, 2007. - 512 p. - (in Russian).
METHOD OF DETERMINING THE TOTAL DISTORTION DIGITAL IMAGES
L.D. Zhimbueva The East-Siberian State Technological University
Abstract
The article describes the method for determining the total distortion of digital images. The algorithm for the solution of the problem is developed and consists of the following steps: (i) creating the ideal image of test object according to its real image, (ii) construction of the image distortion vector diagram, (iii) nonlinear regression analysis of the approximable functions, (iv) solution of the interpolation problem, object image correction and object plane contour reconstruction. The square grid is used as a test object. The mathematical model of the image distortion is built using the Nonlinear estimation package of the STATISTICA system. The reliability value of the investigated function approximation is 0.98-0.99 with cubic polynomial applied. The developed method makes it possible to improve the object plane contour reconstruction accuracy. The method may be used in COMPUTER VISION SYSTEMS in intellectual robot-technical systems, pattern recognition etc.
Key words: image aberration, distortion, nonlinear regression analysis, image correction, contour reconstruction accuracy.
Сведения об авторе
Жимбуева Любовь Дамбиевна, кандидат технических наук, доцент кафедры инженерной и компьютерной графики Восточно-Сибирского государственного технологического университета. Область научных интересов: техническое зрение роботов, обработка изображений, калибровка камеры.
E-mail: Zhimbueva@mail.ru .
Lubov Dambievna Zhimbueva, candidate of engineering sciences, docent of Academic department «Engineering and computer graphics», East-Siberian State Technological University. Interests: robot vision, processing of images, camera calibration.
Поступила в редакцию 2 апреля 2010 г.
