НАУЧНО-ТЕХНИЧЕСКИИ ВЕСТНИК ИНФОРМАЦИОННЫХ ТЕХНОЛОГИИ, МЕХАНИКИ И ОПТИКИ май-июнь 2015 Том 15 № 3 ISSN 2226-1494 http://ntv.ifmo.ru/
SCIENTIFIC AND TECHNICAL JOURNAL OF INFORMATION TECHNOLOGIES, MECHANICS AND OPTICS
May-June 2015
Vol. 15 No 3 ISSN 2226-1494
http://ntv.ifmo.ru/en
УДК 535:631.373.826
ВЫСОКОТОЧНОЕ ОПРЕДЕЛЕНИЕ УГЛОВОГО ПОЛОЖЕНИЯ ТОЧЕЧНОГО ИСТОЧНИКА ИЗЛУЧЕНИЯ С ПЗС-ЛИНЕЙ КАМИ Е.Г. Лебедькоа, Е.Н. Зверева3, Нгуен Ву Тунгь
a Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация
b Государственный технический университет имени Ле Куй Дона, Ханой, 236 Хоанг Куок Вьет, Кау Заи, Вьетнам Адрес для переписки: [email protected] Информация о статье
Поступила в редакцию 22.01.15, принята к печати 01.04.15
doi:10.17586/2226-1494-2015-15-3-398-404
Язык статьи - русский
Ссылка для цитирования: Лебедько Е.Г., Зверева Е.Н., Нгуен Ву Тунг. Высокоточное определение углового положения точечного источника излучения с ПЗС-линейками // Научно-технический вестник информационных технологий, механики и оптики. 2015. Т. 15. № 3. С. 398-404.
Аннотация
Предмет исследования. Предложен метод измерения углового положения точечного источника в системе с ПЗС-линейками посредством преобразования «угол - время - код» в условиях малой освещенности и при повышенных соотношениях сигнал/шум. Проведена оценка потенциальной точности метода, которая определяется неустранимыми аппаратурными случайными погрешностями измерения в условиях оптимальной обработки поступающей информации при однократном считывании ее с ПЗС-линеек.
Метод. Предлагается схема оптико-электронной системы измерения углового положения с ПЗС-линейками с вытянутыми чувствительными элементами и опорной точкой отсчета. При этом на линейки проецируются изображения как опорного точечного источника, так и измеряемого точечного источника, угловое положение которого следует определять с высокой точностью. С выхода ПЗС-линейки сигналы поступают на линейный оптимальный (или близкий к нему) фильтр, а с него - на схему фиксации положения максимума сигнала. Схема обеспечивает минимальную погрешность, обусловленную влиянием шумов. Фронты импульсов, соответствующие максимумам сигналов, формируют временной интервал, который заполняется счетными импульсами высокой частоты. Число импульсов в этом интервале будет соответствовать измеряемому угловому положению точечного источника.
Основные результаты. Проведен анализ случайных погрешностей с позиции теории статистических решений. Анализ учитывает спектральную функцию сигналов, определяемую передаточными функциями оптической системы и ПЗС-линейки как анализатора изображения. Проведено исследование точности измерения в зависимости от тактовой частоты считывания информации с ПЗС-линеек при различных значениях отношения сигнала к шуму. Показано, что даже при однократном считывании ПЗС-линейки по предлагаемой оптико-электронной схеме погрешность измерения не превышает 0,01 размера чувствительного элемента.
Практическая значимость. Результаты работы могут найти применение в высокоточных измерительных оптико-электронных системах звездных датчиков для определения углового положения маломощных точечных источников излучения. Ключевые слова
угловое положение точечного источника, потенциальная точность определения углового положения точечного источника, спектральная функция сигнала на выходе приемника.
HIGH-PRECISION DETERMINATION OF THE ANGULAR POSITION FOR POINT
LIGHT SOURCE WITH CCD-ARRAYS E.G. Lebedko3, E. N. Zvereva3, Nguyen Vu Tungb
a ITMO University, Saint Petersburg, 197101, Russian Federation b Le Quy Don Technical University, Hanoi, 236 Hoang Quoc Viet, Cau Giay, Viet Nam Corresponding author: [email protected] Article info
Received 22.01.15, accepted 01.04.15 doi:10.17586/2226-1494-2015-15-3-398-404 Article in Russian
For citation: Lebedko E.G., Zvereva E.N., Nguyen Vu Tung. ffigh-precision determination of the angular position for point light source with CCD-arrays. Scientific and Technical Journal of Information Technologies, Mechanics and Optics, 2015, vol. 15, no. 3, pp. 398-404.
Abstract
Subject of study. A method of angular position measuring for point light source in the system with CCD arrays is proposed by converting of angle - time - code in low light conditions and at high signal / noise ratio. Assessment of the method
potential accuracy is given. It is determined by instrumental irremovable random errors of measurement in terms of optimal processing of incoming information with a single reading it with CCD - arrays.
Method. This work introduces an optoelectronic system circuit with CCD arrays with stretched sensing elements and a point of reference for angular position measuring. In this case the arrays have images projections of both the reference point source and the target point source, which angular position is measured with high precision. From the CCD array output the signals arrive at an optimal (or apt) linear filter, and then to the signal peak position detection circuit. The scheme provides minimum error due to the influence of noise. Pulse edges, corresponding to the signals maximum, make a time interval filled with high-frequency counting pulses. The number of pulses in this interval will correspond to the measured angular position of the target point source.
Main results. Analysis of random errors has been carried out in terms of statistical decision theory. Analysis takes into account the spectral function of the signals defined by transfer functions of the optical system and a CCD array as an image analyzer. Research of measurement accuracy is carried out depending on the clock frequency of reading information from CCD arrays for different values of signal-to-noise ratio. It has been shown that even with a single readout with CCD array by the proposed opto-electronic circuit measuring error does not exceed 0.01 of the sensing element size.
Practical significance. The results are usable in high-precision measuring opto-electronic systems of star sensors for
determining the angular position of the low-power point light sources.
Keywords
point source angular position, potential accuracy of point source angular position measurement, signal spectral function at the receiver output
Введение
Высокоточное определение углового положения точечных источников излучения широко используется в измерительной технике. Распространение получили измерители на ПЗС-матрицах с последующей цифровой обработкой для нахождения центра тяжести пятна рассеяния на чувствительных элементах [1-11]. При этом для получения высокой точности измерений используется осреднение за N считываний матрицы. Наряду с таким методом измерения авторами предложен высокоточный метод измерения, не требующий компьютерной обработки при однократном считывании.
Целью настоящей работы являлось рассмотрение метода измерения углового положения точечного источника относительно реперной марки в заданном поле посредством преобразования «линейный размер - временной интервал - код» и определение случайной погрешности измерения, приходящейся на один светочувствительный элемент ПЗС-линейки, в условиях оптимальной линейной фильтрации сигналов на выходе фотоприемника при однократном считывании информации.
Метод и устройство измерения углового положения точечного источника излучения
Предлагаемый метод высокоточного определения углового положения точечного источника излучения в пространстве основан на использовании двух ПЗС-линеек с оптимальной линейной фильтрацией сигналов, снимаемых с них, и применением для дальнейшей обработки преобразования «линейный размер - временной интервал - код» (ЛРВК).
Устройство, реализующее этот метод, состоит из оптико-электронной приемной системы, которая включает в себя приемный объектив (ПО), реперный точечный источник излучения (РИИ), полупрозрачное зеркало (ППЗ), две ПЗС-линейки с вытянутыми чувствительными элементами и схемой управления линейками (СУ), а также измерительные схемы (ИС) (рис. 1).
РИИ
Рис. 1. Оптико-электронная приемная система (пояснения в тексте)
Структурная схема измерительного устройства приведена на рис. 2. Измерительная схема по каждой координате (X и У) содержит линейный согласованный фильтр (СФ), позволяющий минимизировать погрешность измерения, обусловленную шумами, схему выделения максимума сигнала (СВМ), которая состоит из дифференцирующей цепи (ДЦ), двух пороговых устройств ПУ1 и ПУ2, схемы совпадений (СС1), и схемы преобразования «линейный размер - временной интервал - код». Это простейшая схема, включающая в себя триггер (Тр), генератор счетных импульсов (ГСИ), формирующий последовательность коротких импульсов высокой частоты, и схемы совпадения (СС2) и счетчика импульсов (СИ).
Рис. 2. Структурная схема измерительного устройства
и 1
2 *
9
Рис. 3. Временные диаграммы, поясняющие работу измерительного устройства.
Номера диаграмм совпадают с номерами сигналов на структурной схеме (рис. 2).
УСПУ1 и УСПУ2 - уровни срабатывания соответственно пороговых устройств ПУ1 и ПУ2
В схеме выделения максимума сигнала ПУ1 и СС1 представляют собой временной селектор, который пропускает только импульс, фронт которого соответствует максимуму сигнала, так как уровень срабатывания ПУ2 близок к нулю. Работа схемы выделения максимума сигнала подробно описана в [12].
Работа устройства осуществляется следующим образом. Слегка расфокусированное изображение реперного и измеряемого точечных источников проецируется на ПЗС-линейки (соответственно по координатам X и У). При считывании сначала на выходе линеек появляется сигнал от реперного источника, после оптимальной фильтрации согласованным фильтром сигнал поступает на дифференцирующую цепочку ДЦ и пороговое устройство ПУ1 схемы выделения максимума. Продифференцированный ДЦ сигнал подается на пороговое устройство ПУ2, уровень срабатывания которого близок к нулю, вследствие
3
X
4
X
5
X
6
7
8
чего ПУ2 срабатывает от любого воздействия, в том числе и при переходе продифференцированного сигнала через ноль. Однако благодаря временной селекции, осуществляемой ПУ1 и схемой совпадения СС1, импульс с порогового устройства ПУ1 открывает схему совпадений СС1 только в момент перехода продифференцированного сигнала через ноль, и на выходе СС1 появится импульс, фронт которого соответствует положению максимума сигнала, снимаемого с согласованного фильтра. Этим фронтом импульса опрокидывается триггер Тр, открывается схема СС2, и счетные импульсы от ГСИ поступают в счетчик СИ. С появлением измеряемого сигнала аналогичным образом фронт импульса, соответствующего максимуму сигнала, снимаемого с согласованного фильтра, приводит Тр в исходное состояние, схема СС2 закрывается, и счетные импульсы от ГСИ прекращают поступать в счетчик.
Таким образом, в счетчике будет записано число импульсов, которое соответствует линейному положению изображения точечного источника относительно реперной марки, а следовательно, и угловое положение точечного источника в заданном поле. При этом измерение проводится при однократном считывании информации с линеек. Работу схемы поясняют временные диаграммы (рис. 3), номера диаграмм совпадают с номерами сигналов на структурной схеме (рис. 2).
Так как реперный и измеряемый сигналы проходят через один и тот же приемно-измерительный тракт, то почти полностью исключаются любые погрешности, вызванные частотно-фазовыми искажениями под действием дестабилизирующих факторов, и обеспечивается стабильность измерений [13].
Потенциальная точность измерительного устройства
Потенциальные возможности рассмотренного метода измерения по точности, определяемые случайными составляющими погрешностей измерения с различными законами распределения, можно оценить по формуле
_ _±А. ±А2 ± 3а. ± 3а2
|р=0,997 т ' ^
с
где Дх- максимальная случайная погрешность измерения временного интервала; Д2 - максимальная случайная погрешность, обусловленная дестабилизирующими факторами при измерении временного интервала; и с2 - средние квадратичные погрешности определения положения максимума реперного и измеряемого сигналов соответственно; Тс=И¥ - период считывания каждого светочувствительного элемента линейки (р - частота считывания элемента).
Формула (1) определяет величину случайной погрешности с вероятностью р=0,997, приходящейся на один чувствительный элемент (пиксель) фотоприемника.
Для удобства расчетов будем исходить из предположения, что максимальная случайная погрешность Д2 составляет в худшем случае 20% от погрешности Д1=+Г0 (Т0 - период следования счетных импульсов ГСИ), а средняя квадратичная погрешность фиксации реперного источника в пять раз меньше соответствующей погрешности фиксации измеряемого сигнала. В этом случае формула (1) принимает вид
5 ±1,2То ±3,6а
|р_0,997 гт ' (2)
с
Естественно, что основной составляющей погрешности измерения может оказаться погрешность определения положения максимума реперного и измеряемого сигналов, которая при оптимальной линейной фильтрации определяется соотношением [14]
а _
2 2 |ю2 (;ю)| ёю
пО 0
(3)
Здесь энергетический спектр О шумов на входе линейного фильтра будем считать равномерным, а спектральную функцию сигнала 5ф (/ю) можно будет представить в виде
^ф (Л®) _ {{ {_„ ^ф ( , 2 ) еХР ^^) + У 2^) _ Ю^ ] ё У 2 Л , (4)
где V! и у2 - пространственные частоты,
^ф (( , ^2 ) = ^ (( , ^2 ) КА (( , ГV2 ) = ^ (, Л'РУ 2 ) К0 (Л V, >2 ) КА (ЛЧ, ^2 )
- спектральная пространственная функция лучистого потока на выходе анализатора изображения, которым является многоэлементный приемник, 5б(]'у1,]\2) - спектральная пространственная функция освещенности на чувствительном слое приемника, КАЦ\ 1,}\-2) - передаточная функция анализатора изображения, К0Ц\1,}\-2) - передаточная функция оптической системы, 5ь(/Ру1,УРу2) - спектральная функция яркости источника излучения, которую для точечного источника можно представить в виде
^ ( , 2 ) = т1 еХР [_Л ( *0 + V2У0 )] , (5)
где т1] = /0Ж(/0)^/0 - среднее значение энергетической силы света /0; - плотность вероятно-
стей энергетической силы света.
Весовая функция характеризует оптическую систему как по способности концентрировать поток излучения, так и по способности распределения энергии в пространстве. При этом эти свойства независимы друг от друга. Исходя из этого, весовую функцию можно представить в виде двух сомножителей:
£о(X У) = £о • £о(^ У) ■
Здесь д0 (х, у) - импульсная характеристика объектива, а т(Х) Д02
g° = D + 4 f0
(6)
где £>0 - диаметр входного зрачка объектива; т(Х) - спектральное пропускание оптической системы. Будем считать, что т(Х) = т0, т.е. не зависит от длины волны в рабочем диапазоне частот.
Для простоты расчетов объективов со сферическими аберрациями удобно весовую функцию объектива аппроксимировать гауссовой функцией [15]
g°( X У) = g°-7exP Po
( xo + y2 Л Po
(7)
где р0 - радиус пятна рассеяния примерно на уровне 0,5 от максимального значения. В этом случае передаточную функцию объектива с учетом (6) и (7) можно представить в виде
K°( j'Vl. jV 0) = g° eXP
- Рл (v0+v 0)
(8)
В матричном фотоприемнике можно считать, что передаточная функция анализатора представляет собой передаточную функцию анализатора с линейным сканированием вида
j\ —i
KA (jV1 , jV0) = KA° (jV1 , jV0 )Z eXP [-j ( + V0 Л* )) =
(9)
где КА0(/У1,уу2) - передаточная функция элементарной ячейки матрицы.
Исходя из предположения, что при М чувствительных элементах линейки, которые имеют форму прямоугольника со сторонами а по оси £ и Ь по оси п, расстояния между чувствительными элементами Д ничтожно малы, передаточную функцию анализатора изображения КА(_можно представить в виде
KAc (jV1' jV0) = abSa
av1
л (bv лsin
. M(a + A)v1
Sa
0
. (a + A)v1 sin-L
-exp
- j
(M - 1)(a + A)v1 0
= abSa
(av Л
Sa
(bv Л sin
. aMv1
-exp
sin
J
a(M - 1)v1
(aM v1 Л ( bMv0 Л
= abMSa Sa exp
l 0 V l 0 V
J
a(M - 1)v1
,(1°)
где Sa(x) = (sin x) / x - функция отсчетов. Так как считывание с элементов линейки осуществляется с
a 1
частотой fc, то функцию ^(t) в (4) можно представить в виде |(t) = V.t = — t, где Tc = —. В этом случае
f
формула (4) принимает вид
^(Jra) = S* Jv0) ^ Jv0)exp J
0л-
a tdt
v1 — ю
[ 1T _
d v1d v 0.
(11)
Интеграл в скобках - интегральное представление дельта-функции, следовательно, используя
фильтрующее свойство дельта-функции и формулу 5 (ах) = —-5(х), соотношение (11) можно записать как
1а1
T fi ( юT |
^°'Ю) = 2л í-™1 J~V,Jv0 Г0.
(10)
Тогда спектральная функция сигнала на выходе приемника (12) с учетом (5), (8), (9), (10) равна
ад
1ТС Г ( аТс ч ( аТс ч / юГс ч
1ТсМЬд0т1]
2п
( ElZkl^l^ ■ (
fТсрх0 мтМ-1
-J(ü |-+ Т,
)
<í Sa
lbv0 ^
exp
l 0 .V
- ^ v 0 + J>pv0 4л
d v0
(13)
k =°
0
где l - чувствительность фотоприемника. Вычисление интеграла приводит к следующему соотношению:
¥
exp
Р2
4П v2 +i>oPv2
2л
т
exp
р2
ch
р2
erf
йл/л ~2р"
0 У
Учитывая, что линейное увеличение в «1, соотношение (14) приводится к виду HL erf f tH },
Р I 2Ро
(14)
(15)
2
где erf(x) =—=•[ exp(-t2)dt - функция Крампа. ч/л 0
Таким образом, выражение для спектральной функции сигнала на выходе приемника учетом (13) и (15) можно представить в виде
S. (» = PSa \MTT^ I exp
.P0T>2 ^
4ла2
erf
2р.
exp
о У
TPXíL + Т
M -1
(16)
где Р = 1ТсМд0т1].
Таким образом, с учетом (16) формула (3) для средней квадратичной погрешности фиксации временного положения максимума сигнала принимает вид
о,бРо7;
ст = -
( ла2M2 I erf ¿%/л
1 2р2 У 1 2Р0 У
(17)
где д - отношение максимума сигнала к среднему квадратичному значению шума на выходе оптимального фильтра.
На графиках рис. 4 приведены, рассчитанные по формуле (2) с учетом (17), кривые изменения погрешности 5|Р=о,997 в зависимости от частоты считывания чувствительного элемента при различных значениях отношения сигнал/шум д и при Тс = 10"8 с (т.е. при частоте импульсов ГСИ - 100 МГц) и отношении пятна рассеяния р0 к размеру чувствительного элемента линейки а, равному десяти. Из рис. 4 видно, что, несмотря на то, что при снижении частоты считывания средняя квадратичная погрешность фиксации временного положения сигнала возрастает, погрешность измерения, приходящаяся на один пиксель падает. Такая высокая точность измерения имеет место при однократном считывании.
1 ,4 1 ,2
1
0,8 0,6
0,4
0 ,2 0
ц=10 ц=50
ц=100
ц=1000
105 106 107 108 f, Гц
Рис. 4. Погрешность измерения в зависимости от частоты считывания светочувствительного элемента
Следует заметить, что вместо линеек с вытянутыми чувствительными элементами можно использовать обычные ПЗС-линейки в сочетании с цилиндрическими линзами. Однако в этом случае будем
иметь энергетический проигрыш в л/2 раз.
Заключение
Проведенный анализ погрешностей рассмотренного метода измерения углового положения точечного источника показал, что потенциальная точность определения углового положения точечного источника излучения зависит, главным образом, от частоты считывания чувствительных элементов линеек и величины отношения сигнал/шум. Показано, например, что при тактовой частоте считывания 105 Гц потенциальная точность может составлять о,о1 и менее размера чувствительного элемента.
Литература
1.
2.
3.
4.
5.
6.
7.
Konyakhin I.A., Turgalieva T.V. Three-coordinate digital autocollimator // Journal of Optical Technology (A Translation of Opticheskii Zhurnal). 2013. V. 80. N 12. P. 772-777. doi: 10.1364/JOT.80.000772 Konyakhin I.A., Kopylova T.V., Konyakhin A.I. Optic-electronic autocollimation sensor for measurement of the three-axis angular deformation of industry objects // Proceedings of SPIE - The International Society for Optical Engineering. 2012. V. 8439. Art. 84391N. doi: 10.1117/12.922096.
Kleshchenok M.A., Anisimov A.G., Lashmanov O.U., Timofeev A.N., Korotaev V.V. Alignment control optical-electronic system with duplex retroreflectors // Proceedings of SPIE - The International Society for Optical Engineering. 2014. V. 9131. Art. 91311X. doi: 10.1117/12.2052290
Konyakhin I., Timofeev A., Usik A., Zhukov D. The experimental research of the systems for measuring the angle rotations and line shifts of the large aperture radio-telescope components // Proceedings of SPIE - The International Society for Optical Engineering. 2010. V. 7544. Art. 75443P. doi: 10.1117/12.885604 Konyakhin I.A., Timofeev A.N., Konyakhin A.I. Three-axis optic-electronic autocollimation system for the inspection of large-scale objects // Proceedings of SPIE - The International Society for Optical Engineering. 2013. V. 8788. Art. 87882C. doi: 10.1117/12.2020343
Konyakhin I.A., Timofeev A.N., Yaryshev S.N. High precision angular and linear measurements using universal opto-electronic measuring modules in distributed measuring systems // Key Engineering Materials. 2010. V. 437. P. 160-164. doi: 10.4028/www.scientific.net/KEM.437.160
Konyakhin I.A., Kopylova T.V., Konyakhin A.I., Smekhov A.A. Optic-electronic systems for measurement the three-dimension angular deformation of axles at the millimeter wave range radiotelescope // Proceedings of SPIE - The International Society for Optical Engineering. 2013. V. 8759. Art. 87593E. doi: 10.1117/12.2014605
Turgalieva T.V., Konyakhin I. A. Research of autocollimating angular deformation measurement system for large-size objects control // Proceedings of SPIE - The International Society for Optical Engineering. 2013. V. 8788. Art. 878832. doi: 10.1117/12.2020861
Konyakhin I.A., Timofeev A.N., Usik A.A., Zhukov D.V. Optic-electronic systems for measuring the angle deformations and line shifts of the reflecting elements at the rotateable radio-telescope // Proceedings of SPIE - The International Society for Optical Engineering. 2011. V. 8082. Art. 80823R. doi: 10.1117/12.890059
10. Pantyushin A., Korotaev V. Control measurement system for railway track position // Proceedings of SPIE -The International Society for Optical Engineering. 2012. V. 8486. Art. 84861B. doi: 10.1117/12.930503
11. Anisimov A.G., Pantyushin A.V., Lashmanov O.U., Vasilev A.S., Timofeev A.N., Korotaev V.V., Gordeev S.V. Absolute scale-based imaging position encoder with submicron accuracy // Proceedings of SPIE - The International Society for Optical Engineering. 2013. V. 8788. Art. 87882T. doi: 10.1117/12.2021022
12. Лебедько Е.Г. Системы импульсной оптической локации. СПб.: Лань, 2014. 368 с.
13. Maraev A.A., Timofeev A.N. Energetic sensitivity of optical-electronic systems based on polychromatic optical equisignal zone // Proceedings of SPIE - The International Society for Optical Engineering. 2013. V. 8788. Art. 878836. doi: 10.1117/12.2020757
14. Лебедько Е.Г. Теоретические основы передачи информации. СПб.: Лань, 2011. 352 с.
15. Левшин В.Л. Обработка информации в оптических системах пеленгации. М.: Машиностроение, 1978. 168 с.
9.
Лебедько Евгений Георгиевич Зверева Елена Николаевна Нгуен Ву Тунг
Evgeny G. Lebedko Elena N. Zvereva Nguyen Vu Tung
доктор технических наук, профессор, профессор, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected] старший преподаватель, Университет ИТМО, Санкт-Петербург, 197101, Российская Федерация, [email protected]
кандидат технических наук, зам. заведующего кафедрой, Государственный технический университет имени Ле Куй Дона, Ханой, 236 Хоанг Куок Вьет, Кау Заи, Вьетнам, [email protected]
D.Sc., Professor, Professor, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
senior lecturer, ITMO University, Saint Petersburg, 197101, Russian Federation, [email protected]
PhD, Deputy Department head, Le Quy Don Technical University, Hanoi, 236 Hoang Quoc Viet, Cau Giay, Viet Nam, [email protected]