Научная статья на тему 'Сравнение по точности алгоритмов определения координат центров изображений в оптико-электронных приборах'

Сравнение по точности алгоритмов определения координат центров изображений в оптико-электронных приборах Текст научной статьи по специальности «Медицинские технологии»

CC BY
656
137
i Надоели баннеры? Вы всегда можете отключить рекламу.
Ключевые слова
аппроксимация окружностью / коллиматор / определение контуров изображений / оптико-электронный прибор / фотоприемник / центр тяжести / circle approximation / collimator / determination of image contours / optoelectronic device / photodetector / centre of gravity

Аннотация научной статьи по медицинским технологиям, автор научной работы — Н О. Старосотников, Р В. Фёдорцев

В метрологических оптико-электронных системах, таких как автоколлиматоры, звездные датчики, датчики Шака–Гартмана, в схемах геометрической калибровки цифровых фотокамер для аэрои космической съемки, в различных системах следящего типа, основанных на определении координат изображений простых форм, одним из важных и существенных параметров является точность определения этих координат. В статье приводится описание математической модели измерительного стенда на базе коллиматора, проецирующего тест-объект на фотоприемник оптико-электронного прибора. Математическая модель учитывает характерные шумы для фотоприемников: дробовый полезного сигнала (фотонный) и дробовый темнового сигнала, считывания и пространственную неоднородность элементов ПЗС-матрицы. С целью снижения влияния шумов предлагается применять фильтр Винера для сглаживания изображения и однозначной его идентификации, а также вводить порог по уровню яркости. Приводится сравнение двух алгоритмов определения координат по энергетическому центру тяжести и по контуру. При определении контура тест-объекта использовались детекторы Собела, Прюитт, Робертса, лапласиана гауссиана, Канни. Сущность алгоритма определения координат заключается в поиске контура изображения в виде окружности с последующей его аппроксимацией и определением центра изображения. Произведен расчет погрешности определения координат центра тяжести для тест-объектов различных диаметров (5, 10, 20, 30, 40, 50 пикс) фотоприемника, а также значений отношения сигнал/шум 200, 100, 70, 20, 10. Отношение сигнал/шум рассчитывалось как разница максимальной интенсивности изображения тест-объекта и фона, деленная на среднеквадратическое отклонение фона. С увеличением отношения сигнал/шум улучшается точность определения координат на 0,5–1 порядок. Улучшение точности с увеличением диаметра тест-объекта характерно для больших отношений сигнал/шум: 70 и более. На основании проведенных исследований установлено, что алгоритм определения координат по энергетическому центру тяжести является более точным по сравнению контурными методами и требует меньших вычислительных мощностей (для программного пакета MatLab), что связано с дискретностью при определении контура.

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

Похожие темы научных работ по медицинским технологиям , автор научной работы — Н О. Старосотников, Р В. Фёдорцев

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

Accuracy Comparison of Algorithms for Determination of Image Center Coordinates in Optoelectronic Devices

Accuracy in determination of coordinates for image having simple shapes is considered as one of important and significant parameters in metrological optoelectronic systems such as autocollimators, stellar sensors, Shack-Hartmann sensors, schemes for geometric calibration of digital cameras for aerial and space imagery, various tracking systems. The paper describes a mathematical model for a measuring stand based on a collimator which projects a test-object onto a photodetector of an optoelectronic device. The mathematical model takes into account characteristic noises for photodetectors: a shot noise of the desired signal (photon) and a shot noise of a dark signal, readout and spatial heterogeneity of CCD (charge-coupled device) matrix elements. In order to reduce noise effect it is proposed to apply the Wiener filter for smoothing an image and its unambiguous identification and also enter a threshold according to brightness level. The paper contains a comparison of two algorithms for determination of coordinates in accordance with energy gravity center and contour. Sobel, Pruitt, Roberts, Laplacian Gaussian, Canni detectors have been used for determination of the test-object contour. The essence of the algorithm for determination of coordinates lies in search for an image contour in the form of a circle with its subsequent approximation and determination of the image center. An error calculation has been made while determining coordinates of a gravity center for test-objects of various diameters: 5, 10, 20, 30, 40, 50 pixels of a photodetector and also signaltonoise ratio values: 200, 100, 70, 20, 10. Signal-to-noise ratio has been calculated as a difference between maximum image intensity of the test-object and the background which is divided by mean-square deviation of the background. The accuracy for determination of coordinates has been improved by 0.5-1 order in case when there was an increase in a signal-to-noise ratio. Accuracy improvement due to increase of a diameter in a test-object is typical for large signal-to-noise ratios: 70 or more. The conducted investigations have made it possible to establish that the algorithm for determination of coordinates of the energy gravity center is more accurate in comparison with contour methods and requires less computing power (for the MatLab software package), which is related to discreteness while determining a contour.

Текст научной работы на тему «Сравнение по точности алгоритмов определения координат центров изображений в оптико-электронных приборах»

DOI: 10.21122/2227-1031 -2018-17-1-79-86 УДК 681.723.078, 681.775.078, 681.777.078

Сравнение по точности алгоритмов определения координат центров изображений в оптико-электронных приборах

Асп. Н. О. Старосотников1*, канд. техн. наук, доц. Р. В. Фёдорцев1*

^Белорусский национальный технический университет (Минск, Республика Беларусь)

© Белорусский национальный технический университет, 2018 Belarusian National Technical University, 2018

Реферат. В метрологических оптико-электронных системах, таких как автоколлиматоры, звездные датчики, датчики Шака-Гартмана, в схемах геометрической калибровки цифровых фотокамер для аэро- и космической съемки, в различных системах следящего типа, основанных на определении координат изображений простых форм, одним из важных и существенных параметров является точность определения этих координат. В статье приводится описание математической модели измерительного стенда на базе коллиматора, проецирующего тест-объект на фотоприемник оптико-электронного прибора. Математическая модель учитывает характерные шумы для фотоприемников: дро-бовый полезного сигнала (фотонный) и дробовый темнового сигнала, считывания и пространственную неоднородность элементов ПЗС-матрицы. С целью снижения влияния шумов предлагается применять фильтр Винера для сглаживания изображения и однозначной его идентификации, а также вводить порог по уровню яркости. Приводится сравнение двух алгоритмов определения координат по энергетическому центру тяжести и по контуру. При определении контура тест-объекта использовались детекторы Собела, Прюитт, Робертса, лапласиана гауссиана, Канни. Сущность алгоритма определения координат заключается в поиске контура изображения в виде окружности с последующей его аппроксимацией и определением центра изображения. Произведен расчет погрешности определения координат центра тяжести для тест-объектов различных диаметров (5, 10, 20, 30, 40, 50 пикс) фотоприемника, а также значений отношения сигнал/шум 200, 100, 70, 20, 10. Отношение сигнал/шум рассчитывалось как разница максимальной интенсивности изображения тест-объекта и фона, деленная на среднеквадратическое отклонение фона. С увеличением отношения сигнал/шум улучшается точность определения координат на 0,5-1 порядок. Улучшение точности с увеличением диаметра тест-объекта характерно для больших отношений сигнал/шум: 70 и более. На основании проведенных исследований установлено, что алгоритм определения координат по энергетическому центру тяжести является более точным по сравнению контурными методами и требует меньших вычислительных мощностей (для программного пакета MatLab), что связано с дискретностью при определении контура.

Ключевые слова: аппроксимация окружностью, коллиматор, определение контуров изображений, оптико-электронный прибор, фотоприемник, центр тяжести

Для цитирования: Старосотников, Н. О. Сравнение по точности алгоритмов определения координат центров изображений в оптико-электронных приборах / Н. О. Старосотников, Р. В. Фёдорцев // Наука и техника. 2018. Т. 17, № 1. С. 79-86. DOI: 10.21122/2227-1031-2018-17-1-79-86

Accuracy Comparison of Algorithms for Determination of Image Center Coordinates in Optoelectronic Devices

M. A. Starasotnikau1*, R. V. Feodortsau1*

^Belarusian National Technical University (Minsk, Republic of Belarus)

Abstract. Accuracy in determination of coordinates for image having simple shapes is considered as one of important and significant parameters in metrological optoelectronic systems such as autocollimators, stellar sensors, Shack-Hartmann

Адрес для переписки

Фёдорцев Ростислав Валерьевич Белорусский национальный технический ул. Я. Коласа, 22,

220013, г. Минск, Республика Беларусь

Тел.: +375 17 292-62-86

[email protected]

Address for correspondence

Feodortsau Rostislav V.

Belarusian National Technical University

22 Ya. Kolasa str.,

220013, Minsk, Republic of Belarus

Tel: +375 17 292-62-86

[email protected]

H Наука

итехника. Т. 17, № 1 (2018)

sensors, schemes for geometric calibration of digital cameras for aerial and space imagery, various tracking systems. The paper describes a mathematical model for a measuring stand based on a collimator which projects a test-object onto a photo-detector of an optoelectronic device. The mathematical model takes into account characteristic noises for photodetectors: a shot noise of the desired signal (photon) and a shot noise of a dark signal, readout and spatial heterogeneity of CCD (charge-coupled device) matrix elements. In order to reduce noise effect it is proposed to apply the Wiener filter for smoothing an image and its unambiguous identification and also enter a threshold according to brightness level. The paper contains a comparison of two algorithms for determination of coordinates in accordance with energy gravity center and contour. Sobel, Pruitt, Roberts, Laplacian Gaussian, Canni detectors have been used for determination of the test-object contour. The essence of the algorithm for determination of coordinates lies in search for an image contour in the form of a circle with its subsequent approximation and determination of the image center. An error calculation has been made while determining coordinates of a gravity center for test-objects of various diameters: 5, 10, 20, 30, 40, 50 pixels of a photodetector and also signal-to-noise ratio values: 200, 100, 70, 20, 10. Signal-to-noise ratio has been calculated as a difference between maximum image intensity of the test-object and the background which is divided by mean-square deviation of the background. The accuracy for determination of coordinates has been improved by 0.5-1 order in case when there was an increase in a signal-to-noise ratio. Accuracy improvement due to increase of a diameter in a test-object is typical for large signal-to-noise ratios: 70 or more. The conducted investigations have made it possible to establish that the algorithm for determination of coordinates of the energy gravity center is more accurate in comparison with contour methods and requires less computing power (for the MatLab software package), which is related to discreteness while determining a contour.

Keywords: circle approximation, collimator, determination of image contours, optoelectronic device, photodetector, centre of gravity

For citation: Starasotnikau M. A., Feodortsau R. V. (2018) Accuracy Comparison of Algorithms for Determination of Image Center Coordinates in Optoelectronic Devices. Science and Technique. 17 (1), 79-86. DOI: 10.21122/2227-1031-201817-1-79-86 (in Russian)

Введение

В основе математической обработки данных автоколлиматоров [1], звездных датчиков [2], датчиков Шака-Гартмана [3], схем геометрической калибровки оптико-электронных приборов [4], различных следящиих систем лежит определение координат изображения, спроецированного на фотоприемник оптико-электронного прибора. Для повышения точности измерений необходимо исследовать влияние различных параметров на точность определения координат изображения. В [5] исследованы алгоритм определения координат изображения по центру тяжести и способы повышения его точности. В [6] для калибровки цифровой камеры предложено предварительно определять контур проецируемого массива точек, а в последующем - их центры.

В данной работе исследуется точность алгоритма определения координат, заключающегося в определении контура изображения в виде окружности, аппроксимации данного контура и нахождении его центра, который в свою очередь совпадает с центром изображения. Определение контура выполнялось при помощи детекторов Собела, Прюитт, Робертса, лапласиана гауссиана, Канни. Данные детекторы наиболее часто используются в задачах цифровой обработки изображений для обнаружения границ по разрыву яркости [7].

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

Описание математической модели

В основе исследований заложена математическая модель часто встречающейся схемы контроля, которая описана в [5]. Объектив коллиматора проецирует тест-объект в виде круга на фотоприемник оптико-электронного прибора. Для упрощения объективы коллиматора и оптико-электронного прибора идентичны, с фокусным расстоянием 1000 мм. Длина волны подсветки тест-объекта коллиматора 0,6 мкм. Расчет производился для тест-объектов диаметрами 5, 10, 20, 30, 40, 50 пикс фотоприемника. Размер пикселя фотоприемника 4,5 мкм. Шаг дискретизации в модели - 5 (0,9 мкм). Таким образом, на один пиксель фотоприемника

Наука

итехника. Т. 17, № 1 (2018)

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

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

В модели учитывались следующие шумы, характерные для современных промышленных цифровых камер: дробовый полезного сигнала (фотонный) и дробовый темнового сигнала, считывания и пространственная неоднородность элементов фотоприемника. Дробовый шум описывается распределением Пуассона, шум считывания - распределением Гаусса со среднеквадра-тическим отклонением (СКО) 15 электрон, тем-новой сигнал - 15 электрон. Пространственная неоднородность - 1 %. Также учитывалось квантование фотоприемником изображения тест-объекта. Шаг кантования - 10 бит.

Для снижения влияния шумов применялся фильтр Винера размером 5x5 пикс, а также вводился порог [8], значения яркости ниже которого принимались равными нулю. Величина порога была в 10 раз больше среднего значения фоновой яркости. Размер фильтра и величина порога определялись для достижения достоверного вычисления контуров, т. е. при данных значениях определяется только один замкнутый контур изображения тест-объекта. Помимо этого, такие меры повышают точность расчета центра тест-объекта по центру тяжести.

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

■ Наука

итехника. Т. 17, № 1 (2018)

новным критерием сравнения алгоритмов определения центра изображения. Погрешность является периодической и систематической [5]. Максимум погрешности приходится на 0,25 и 0,75 пикс смещения изображения тест-объекта относительно матричной структуры фотоприемника, минимум - на 0,5 и 1,0 пикс. Однако погрешность также носит случайный характер, поскольку неизвестно первоначальное расположение изображения тест-объекта относительно матричной структуры фотоприемника, и кроме как при математическом моделировании, невозможно предварительно ее исключить. Поэтому стоит задача в ее исследовании, поиске оптимального размера диаметра тест-объекта, других способов математической обработки для уменьшения погрешности.

Описание алгоритмов определения

центра изображений

Алгоритм определения координат по центру тяжести осуществляется по следующей формуле

N M

Z Z*

X Л! — Я =1 =1

л0>У0 - NM

Z Z.

n —1 m —1

n m

n ,m '

(1)

где 1пт - величина яркости в пикселе п, т; п, т - порядковые номера пикселей по столбцам и строкам матрицы яркостей; N М - общее число столбцов и строк матрицы яркостей.

С помощью (1) производится непосредственно оценка координат изображения тест-объекта, а также осуществляется интерполяция результатов измерений дискретной матричной структурой на внутренние точки интервалов дискретизации [5].

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

центром изображения на фотоприемнике. Определение контура заключается в обработке изображения с фотоприемника скользящим детектором, что осуществляется путем свертки. От детектора будут зависеть форма и достоверность контура. Детекторы Собела, Прюитт, Робертса используют следующие фильтры для численного приближения производных [9] по соответствующим направлениям: • Собела:

-1 -2 -1 -1 0 1

N 0 0 0 ? -2 0 2

1 2 1 -1 0 1

■ Прюитт:

-1 -1 -1 -1 0 1

Gx,Gy - 0 0 0 1 -1 0 1

1 1 1 -1 0 1

"-1 0" " 0 -1"

0 1 ? 1 0

■Робертса:

Gx, Gy -

Лапласиан гауссиана обнаруживает края, выполняя поиск пересечений нулевого уровня после фильтрации гауссианом. При использовании детектора Канни [10] изображение первоначально сглаживается гауссовым фильтром для уменьшения шума. В каждой точке вычисляются градиент и направление края. Точки перепада определяются как точки локального максимума градиента. Алгоритм отслеживает этот максимум и присваивает точкам вне максимума нулевые значения. Результатом является контур. Данный контур подвергается двух-пороговой обработке, в результате добавляются дополнительные точки для формирования плавного контура. Результатом обработки изображения тест-объекта (рис. 1а) будут логические массивы данных, которые представлены на рис. 1Ь-1£

b

a

c

d е f

Рис. 1. Представление изображения тест-объекта с указанием центра: a - изображение тест-объекта на фотоприемнике; b, c, d, e, f - результат обработки детекторами Канни, лапласиана гауссиана, Прюитта, Робертса и Собела соответственно

Fig. 1. Representation of test-object image with center indication: a - test-object image on photodetector; b, c, d, e, f - processing results while using Canny, Laplasian Gaussian, Pruitt, Roberts and Sobel detectors, respectively

Наука

итехника. Т. 17, № 1 (2018)

Описание алгоритма аппроксимации окружностью

Массивы контуров, полученные при помощи применения разных детекторов, аппроксимировались окружностью, чтобы затем определить ее центр. Аппроксимация является минимизацией следующей суммы квадратов путем дифференцирования 1 по х0, у0 и Я [9], где суммирование ведется по всем контурным элементам:

N , 2 2 \2

1=Ех-хо) +(у-Уо0-Я),

7=1

где х7, у7 - контурные элементы; - вес контурного элемента, 7 = 1 ... п.

Тогда х0, у0 и Я задаются формулами:

xo =

ByCx - BxCy AyCx - AxCy AxBy - AyBx ' Уо AyBx - AxBy '

r2 = W Zw (( x - xo)2+( y - Уо)2)2 ;

где

N 1 N 1 N

W=Zw; х=wwZwx у=wwZwy;

i=1

Ax=Zw ( x-х ) x ;Bx=Zw ( x- х ) y ;

г=1

i=1

1 N

Cx=2 Z w ( x- х )( xi+y2 );

2 i=1

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

N

N

Ay=Zw ( У - у ) x ; By=Zw ( У - у ) У ;

i=1

i=1

1 N

Су=1 Z у-у)(х2+у2).

2 7=1

Влияние на точность аппроксимации различных факторов, вызывающих отклонение от окружности, может быть уменьшено путем присвоения элементам контура весов В первоначальной итерации веса всех элементов контура будут равны 1. Затем все контурные точки повторно оцениваются в зависимости от их отклонения от аппроксимированной окруж-

■ Наука

итехника. Т. 17, № 1 (2018)

ности, и вычисление х0, у0 и Я повторяется. Наиболее естественным выбором для весов является обратное квадрату расстояние от отдельных элементов контура х7, у7 к ранее определенной окружности, т. е. = 1/й 7,

где й 7 =^(х7 - х0)2 +(у 7 - у0)2 - Я - невязки.

Помимо этого, исключаются элементы массива, которые явно не относятся к истинной окружности. Веса определяются

w =

= Jl- d i К 2a) для |d |< 2a;

0 иначе,

где a2 - СКО

N

a2 =Z d2 ( n-1).

7=1

Процедура определения весов и пересчета параметров аппроксимирующей окружности выполняется итерационно. Критерием остановки итераций является прекращение уменьшения а2. Обычно необходимо выполнить три итерации.

Центр окружности и будет являться центром изображения тест-объекта на фотоприемнике.

Сравнение двух алгоритмов

Численные исследования проводились на разработанной математической модели. Произведен расчет погрешности определения координат центра тяжести различными алгоритмами - энергетическим и контурными - для тест-объектов диаметров 5, 10, 20, 30, 40, 50 пикс фотоприемника, а также значений отношения сигнал/шум 200, 100, 70, 20, 10. Отношение сигнал/шум рассчитывалось как разница максимальной интенсивности изображения тест-объекта и фона, деленная на СКО фона. Расчет производился в МаЛаЬ для выборки, состоящей из 100 изображений, для каждого диаметра тест-объекта. Изображения отличались друг от друга из-за случайного характера шумов фотоприемника. По каждому изображению определялся центр тяжести. Задавалось смещение тест-объекта на 0,25 пикс и вычислялась погрешность как разница между заданным и фактическим смещениями тест-объекта отно-

сительно матричной структуры фотоприемника. Затем значения погрешности усреднялись по всей выборке с учетом весов. Полученные средневзвешенные результаты зависимости по-

грешности определения центра изображения тест-объекта на фотоприемнике от диаметра тест-объекта для различных алгоритмов представлены на рис. 2.

0,15 j 0,14 0,13 0,12

I

с 0,10

S3

X 0,07

В- 0,0!

С

0,03 0,02 0,01 О

0,15 0,14 0,13

0 0,12

1 С.11 с 0,10 ¡5 0,09

0 сев

X 01,07

1

£ ='.Jr О 0,04 С 0.0 J 0,02 0,01 О

0,15 0,14 0,13 0,12 0,11 0,10 0,09 0,08 0,07 0,06 0,05 0,04 0,03 0,02 0,01 О

0,15 0,14 0,13 0,12 0,11 0,10 0,09 , 0,08 0,07 1 0,06 0,05 0,04 0,03 0,02 ' 0,01 1 О

Отношение сигнал/шум 200 (signal-to-noise ratio 200)

\

\

x\ \ \

\ \ N

*.\ 1 N

V. N -\- m

4 Ч-

Vi .............♦-t-^ Sv.................. ----ж ----

10 15 20 25 30 35

Диаметр тест-объекта, пикс

40 45

V Отношение сигнал/шум 100 (signal-to-noise ratio 100)

\ - -*

iL

I V ч

\

£ \ \

• Ч v N

v * •Л \ ---------1

\ V .t \ _ _ - ✓

4 ___- N. •^««HU«—j ___^--4-«'............ У ✓

\ Ж....... -------

20 25 30 35

Диаметр тест-объекта, пикс

Отношение сигнал/шум 70 (signal-to-noise ratio 70)

Ш _

/ N

К ' 4 \

\ ' V A4 \

N

/ \

1 V'__ \_ __ _ —ЯЯ. .__

Т» N ■---- -'i- - ---- N - М» —

^^ \m............. i^- Ш...... N .............я........_ - __ — -rr...... ---- - -я

х

X

0,006 л I-

0,0 05 О х

О'Сч § 0,011.1 О

с

0,002 0,001 О

0,010 0,009

0,008 о

0,00.' 1

о ооо ¡£ о

0,00= о

' 0.004 |

I ö 1". в-0,0.10? | 0,001 о

0,010

0,009

0,000

ООО/ х X

■0,001=. £ I-

0,005 g

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

0 00:4 3

(D

0,7.0=

О

0 002 П. 0,001 о

Диаметр тест-объекта, пикс Отношение сигнал/шум 20 (signal-to-noise ratio 20)

----«----_

---1

_ — X ßr'' Г -Ч

\ / \ v N N __ _ -

M - 4 V >4 ^— — ~ ~~

'.......Щ......... —1 ■■•■■•■■•■■•■•■■•■ ft ■ • W — — — — .........f............. —..........J _ ■

20 25 30 35

Диаметр тест-объекта, пикс

Рис. 2. Погрешность (ось ординат) определения центра изображения тест-объекта на фотоприемнике различными алгоритмами для разных диаметров (ось абсцисс)

Fig. 2. Error (ordinate axis) for determination of test-object image center on photodetector while using various algorithms for different diameters (abscissa axis)

Наука

итехника. Т. 17, № 1 (2018)

0,15 0,14 0,13 0,12 0,11 0,10 0,09 0,08 1 0,07 0,06 0,05 J 0,04 0,03 0,02 0,01 0

Отношение сигнал/шум 10 (signal-to-noise ratio 10)

20 25 30 35

Диаметр тест-объекта, пикс

— • — Канни / Canny

— ■ — Лапсласиан гауссианэ/LOG

Прюитт/ Prewitt Роберте/ Roberts Собел /Sobel

• Энергетический центр тяжести (правая ось) / Energy center of gravity (right axis) Рис. 2. Окончание Fig. 2. Termination

Как видно из рис. 2, точность определения координат по энергетическому центру тяжести превосходит контурные алгоритмы для всех рассмотренных диаметров тест-объекта и отношений сигнал/шум примерно на порядок.

С увеличением отношения сигнал/шум улучшается точность определения координат на 0,5-1 порядок. Улучшение точности с увеличением диаметра тест-объекта наиболее характерно для больших отношений сигнал/шум: 70 и более. В этом случае для детектора Канни характерное улучшение точности заметно для тест-объектов диаметром 20 пикс и более, для остальных - 10 пикс и более.

Для тест-объектов диаметром 10 пикс и менее в детекторах Канни и лапласиан гауссиана погрешность оказывается постоянной и, следовательно, СКО равняется нулю при выборке 100, в отличие от других детекторов.

Уступают контурные алгоритмы по величине вычислительных мощностей (табл. 1), на каж-

Таблица 1

Сравнение алгоритмов по величине вычислительных мощностей

Comparison of algorithms according to computing power values

Алгоритм

Энергетический Контурный

Канни Лапласиан гауссиана Прюитт Робертс Собел

x 14,4x 14,1x 12,5x 12,5x 12,2x

■ Наука

итехника. Т. 17, № 1 (2018)

дый из которых необходимо, по сравнению с энергетическим алгоритмом, затратить больше времени, в частности за счет дополнительной операции аппроксимации окружностью и проведения нескольких итераций во время расчета. Значения необходимых вычислительных мощностей определялись средствами Ма1ЬаЬ по затрачиваемому времени вычислений на каждый алгоритм.

ВЫВОДЫ

1. Математически смоделирована оптическая схема коллиматора, проецирующего тест-объект на фотоприемник оптико-электронного прибора с параметрами, которые часто встречаются в реальных схемах контроля. Для сравнения алгоритмов определения координат изображений проведен численный расчет для отношений сигнал/шум 200, 100, 70, 20, 10 и диаметров тест-объекта 5, 10, 20, 30, 40, 50 пикс с учетом шумов фотоприемника. Установка порога, значения яркости ниже которого принимались нулю, и применение фильтра Винера с подобранными параметрами позволяют сгладить изображение и устранить влияние шумов на автоматическое определение контура. Это дает возможность однозначно его идентифицировать, т. е. исключить определение ложных контуров, а также повысить точность алгоритма определения координат по центру тяжести от 1 до 2 порядков.

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

3. Для определения координат центров изображений в оптико-электронных приборах целесообразнее использовать алгоритм определения по энергетическому центру тяжести по сравнению с контурными алгоритмами.

ЛИТЕРАТУРА

1. Старосотников, Н. О. Высокоточный цифровой автоколлиматор для измерения малых углов / Н. О. Старо-сотников // Новые направления развития приборостроения: материалы 7-й Междунар. студ. науч.-техн. конф., 23-25 апр. 2014 г. Минск: БНТУ, 2014. С. 244.

2. Jinyun, Yan. Dynamic Imaging Model and Parameter Optimization for a Star Tracker / Yan Jinyun, Jiang Jie, Zhang Guangjun // Optics Express. 2016. Vol. 24, Iss. 6. Р. 5961-5983. DOI: 10.1364/ oe.24.005961.

3. Adaptive Thresholding and Dynamic Windowing Method for Automatic Centroid Detection of Digital Shack-Hartmann Wavefront Sensor / Xiaoming Yin [et al.] // Applied Optics. 2009. Vol. 48, Iss. 32, P. 6088-6098. DOI: 10.1364/ao.48.006088.

4. Фотограмметрические параметры оптико-электронной аппаратуры / С. А. Архипов [и др.] // Вестник МГТУ имени Н. Э. Баумана. Сер. Приборостроение. 2008. № 4. С. 105-115.

5. Старосотников, Н. О. Оценка точности определения координат энергетического центра тяжести тест-объекта коллиматора в схемах контроля оптико-электронных приборов с матричными фотоприемниками / Н. О. Ста-росотников, Р. В. Фёдорцев // Наука и техника. 2015. № 5. С. 71-76.

6. Accurate Projector Calibration Based on a New Point-to-Point Mapping Relationship Between the Camera and Projector Images / Huang Zhengrong [et al.] // Applied Optics. 2015. Vol. 54, Iss. 3. Р. 347-356. DOI: 10.1364/ ao.54.000347.

7. Старосотников, Н. О. Метод снижения влияния шумов фонового сигнала при определении координат энергетического центра тяжести изображения в ОЭП / Н. О. Старосотников, Р. В. Фёдорцев // Приборострое-ние-2016: материалы 9-й Междунар. науч.-техн. конф., 23-25 нояб. 2016 г. Минск: БНТУ, 2016. С. 133-135.

8. Гонсалес, Р. Цифровая обработка изображений в среде MatLab; пер. с англ. / Р. Гонсалес, Р. Вудс, С. Эддинс. М.: Изд-во «Техносфера», 2006. 616 с.

9. Determination of Ocular Torsion by Means of Automatic Pattern Recognition / Е. Groe [et al.] // IEEE Transactions on Biomedical Engineering. 1996. Vol. 43, Iss. 5. Р. 471-479. DOI: 10.1109/10.488795.

10. Canny, J. A Computational Approach to Edge Detection / J. Canny // IEEE Transactions on Pattern Analysis and Machine Intelligence. 1986. Vol. PAMI-8, Iss. 6, Nov. P. 679-698. DOI: 10.1109/ tpami.1986.4767851. DOI: 10. 1109/tpami.1986.4767851.

Поступила 14.08.2017 Подписана в печать 10.11.2017 Опубликована онлайн 29.01.2018

REFERENCES

1. Starosotnikau N. O. (2014) Highly-Accurate Digital Auto-collimator for Measuring Small Angles. Novye Naprav-leniya Razvitiya Priborostroeniya: Materialy 7-i Mezh-dunarodnoi Studencheskoi Nauch.-Tekhn. Konferentsii, 23-25 Aprelya 2014 g. [New Directions in Instrumentation Development: Materials of the 7th International Students' Scientific and Technical Conference]. Minsk, Bela-rusian National Technical University, 244 (in Russian).

2. Jinyun Yan, Jie Jiang, Guangjun Zhang (2016) Dynamic Imaging Model and Parameter Optimization for a Star Tracker. Optics Express, 24 (6), 5961-5983. DOI: 10.1364/ oe.24.005961.

3. Xiaoming Yin, Xiang Li, Liping Zhao, Zhongping Fang (2009) Adaptive Thresholding and Dynamic Windowing Method for Automatic Centroid Detection of Digital Shack-Hartmann Wavefront Sensor. Applied Optics, 48 (32), 6088-6098. DOI: 10.1364/ao.48.006088.

4. Arkhipov S. A., Gasitch G. V., Zavarzin V. I., Morozov S. A. (2008) Photogrammetric Parameters of Optical and Electronic Apparatus. VestnikMoskovskogo Gosudarstvennogo Tekhnicheskogo Universiteta imeni N. E. Baumana = Herald of the Bauman Moscow State Technical University. Series Instrument Engineering, (4), 105-115 (in Russian).

5. Starosotnikau N. O., Feodortsau R. V. (2015) Estimation of Accurate Determination for Coordinates of Gravity Energy Center in Collimator Test-Object in Respect of Control Schemes for Optoelectronic Devices with Matrix Photodetectors. Nauka i Tekhnika = Science and Technique, (5), 71-76 (in Russian).

6. Zhengrong Huang, Jiangtao Xi, Yanguang Yu, Qing-hua Guo (2015) Accurate Projector Calibration Based on a New Point-to-Point Mapping Relationship Between the Camera and Projector Images. Applied Optics, 54 (3), 347-356. DOI: 10.1364/ao.54.000347.

7. Starosotnikau N. O., Feodortsau R. V. (2016) Method for Decreasing Influence of Background Signal Noise while Determining Energy Gravity Centre Coordinates for Images in Electrooptical Devices. Priborostroenie-2016: Materialy 9-i Mezhdunarodnoi Nauch.-Tekhn. Konferentsii, 23-25 Noyabrya 2016 g. [Instrumentation-2016: Materials of the 9th International Scientific and Technical Conference]. Minsk, Belarusian National Technical University, 133-135 (in Russian).

8. Gonzalez R., Woods R., Eddins S. (2004) Digital Image Processing Using MatLab. New Jersey, Prentice Hall.

9. Groen E., Bos J. E., Nacken P.F.M., de Graaf B. (1996) Determination of Ocular Torsion by Means of Automatic Pattern Recognition. IEEE Transactions on Biomedical Engineering, 43 (5), 471-479. DOI: 10.1109/10.488795.

10. Canny J. (1986) A Computational Approach to Edge Detection. IEEE Transactions on Pattern Analysis and Machine Intelligence, PAMI-8 (6), 679-698. DOI: 10.1109/ tpami.1986.4767851.

Received: 14.08.2017 Accepted: 10.11.2017 Published online: 29.01.2018

Наука

итехника. Т. 17, № 1 (2018)

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