РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ ОПТИКО-ЭЛЕКТРОННОЙ СИСТЕМЫ НАБЛЮДЕНИЯ ЗА ТОЧЕЧНЫМИ ОБЪЕКТАМИ
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ …
УДК 004.942:519.876.5
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ ОПТИКО-ЭЛЕКТРОННОЙ СИСТЕМЫ НАБЛЮДЕНИЯ
ЗА ТОЧЕЧНЫМИ ОБЪЕКТАМИ
А.Л. Андреев, Т.П. Лбова
Дается описание структуры модели, используемой для синтеза, анализа и прогнозирования работы оптико-электронной системы (ОЭС) телевизионного типа, решающей задачи наблюдения за подвижными объектами в автоматическом режиме. В качестве исходных данных задаются параметры входных воздействий объектов и фоновых помех, параметры оптической системы, фотоприемника, дестабилизирующих факторов, алгоритмов и др. Ключевые слова: компьютерное имитационное моделирование систем, телевизионные системы наблюдения, обнаружение и измерение координат точечных объектов.
Введение
Целью настоящей статьи является описание принципов построения комплексной математической модели, используемой для синтеза, анализа, прогнозирования работы телевизионной системы наблюдения за малоразмерными подвижными объектами и обоснования требований к отдельным звеньям системы с учетом необходимых эксплуатационных параметров в реальных условиях функционирования.
К числу важных промежуточных задач, которые нужно решить для достижения поставленной цели, относятся, в частности: − исследование особенностей физических процессов, протекающих в многоэлемент-
ных структурах фотоприемных устройств с учетом физических свойств различных фоточувствительных материалов; − исследование физической природы источников помех, действующих в различных звеньях оптико-электронной системы и др.; − исследование различных алгоритмов обработки сигналов, используемых для автоматического обнаружения и измерения параметров объектов.
Универсальность модели предполагает возможность ее использования для различных условий функционирования подобных систем – с использованием наземного, космического или иного базирования, различных дистанций наблюдения, энергетических и динамических параметрах объектов и др.
Структура модели
Структура математической модели ОЭС представлена на рис. 1. Задача моделирования ОЭС разбивается на следующие этапы: моделирование сигнально-фоновой ситуации; моделирование различных видов помех; моделирование процедур предварительной обработки сигнала (с учетом искажений, связанных с квантованием сигналов при аналого-цифровом преобразовании); моделирование различных алгоритмов целевой обработки сигналов (связанных, в общем случае, с решением задач обнаружения, распознавания и оценки параметров исследуемых объектов, а также оценки метрологических параметров проектируемой системы с учетом заданных доверительных интервалов). Ввод исходных данных осуществляется в интерактивном режиме. При этом вводимые переменные можно подразделить на следующие группы: − параметры оптической системы; − параметры фотоприемного устройства; − параметры алгоритма первичной обработки сигнала;
10 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
− условия наблюдения; − параметры, определяющие режим статистической обработки полученных реализа-
ций. Предварительная обработка исходных данных предполагает преобразование вво-
димых величин и формирование на их основе переменных и констант, используемых при моделировании.
Рис. 1. Структурная схема модели
В общем случае фоновая составляющая сигнала представляет собой случайное поле яркости. Исходным материалом для компьютерного моделирования случайного поля яркости являются независимые числа, полученные генератором случайных чисел ЭВМ. Их совокупность можно рассматривать как случайное δ-поле – обобщение понятия дискретного белого шума на случай нескольких переменных. Моделирование δ-поля осуществляется следующим образом: пространственно-временной координате (xi , y j ,tk ) ставятся в соответствие выборочные значения из датчика нормальных слу-
чайных чисел с параметрами (0, 1). Для моделирования дискретных реализаций однородного стационарного случай-
ного поля, заданного в конечном пространстве и времени, предлагается использовать метод канонического разложения поля в пространственно-временной ряд Фурье. Неограниченные дискретные реализации однородного стационарного случайного поля можно получить с помощью алгоритма пространственно-временного скользящего суммирования δ-поля [1].
При формировании массива темновых зарядов учитывается аддитивная составляющая «геометрического» шума. Среднее значение темнового заряда, накапливаемого в ячейке, определяется выражением [2]
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
11
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ …
QT = iT AэлTн ,
где iT – средняя плотность темнового тока накопления заряда, характерная для данного полупроводникового материала при рабочей температуре кристалла; Аэл – площадь электрода, под которым осуществляется накопление.
Формирование массива сигнальных зарядов осуществляется с учетом случайного распределения темновых токов и чувствительности по кристаллу ПЗС с заданным интервалом корреляции «геометрического» шума. Если, например, распределение освещенности от точечного объекта по фоточувствительной площадке аппроксимируется гауссоидой вращения, то распределение сигнальных зарядов определяется в соответствии с выражением
E(xi ,
yj
)
=
Ф 2πR2
exp
⎢⎡− ⎣
(x
−
x0 )2 + (y 2R2
−
y0 )2
⎤ ⎥ ⎦
,
где Ф – поток оптического излучения, создающего изображение точечного объекта; R –
– радиус кружка рассеяния объектива, определяемый на уровне 0,606 Emax; x0, y0 – ко-
ординаты центра пятна рассеяния.
Динамику перемещения наблюдаемых объектов можно задать в параметрическом
виде с помощью коэффициентов А, В, С (первой, второй и третьей степени) полинома,
аппроксимирующего изменение во времени модуля угловых перемещений объектов
S(t) = A ⋅ t + B ⋅ t 2 + C ⋅ t3 .
Направление перемещения задается посредством угла α между направлением оси X и проекцией вектора скорости в плоскости анализа изображения.
При моделировании массива шумовых составляющих в ячейках ПЗС учитывается совокупность сигнальной, фоновой и темновой составляющих зарядов в каждой ячейке. При этом шумовая составляющая заряда в каждой ячейке определяется путем многократного задания флуктуаций зарядов по случайному закону с заданным распределением с учетом влияния всех вышеперечисленных факторов.
Флуктуация числа зарядов подчиняется закону Пуассона. Однако, поскольку чис-
ло накапливаемых зарядов достаточно велико ( nq > 100 ), распределение Пуассона
можно аппроксимировать нормальным законом распределения с дисперсией, равной среднему значению. Тогда среднеквадратическое значение шума на выходе ФПУ равно
σш =
σ2 ш.т
+
σ2 ш.ф
+
σ2 ш.с
+
σ2 ш.г
+
σ2 ш.п
+
σ2 ш.в.у.
,
где σш.т – среднеквадратическое число зарядов в i-ой ячейке ФПУ, характеризующее флуктуацию числа темновых зарядов; σш.ф – среднеквадратическое значение флуктуа-
ции числа зарядов, обусловленных фоновым излучением; σш.с – среднеквадратическое
значение флуктуации числа сигнальных зарядов, накапливаемых в ячейках ФПЗС, ха-
рактеризующее фотонный шум; σш.г – среднеквадратическое отклонение среднего
числа темновых зарядов в различных ячейках ФПУ, характеризующее геометрический
шум; σш.п – среднеквадратическое значение флуктуаций заряда при переносе (для
ФПЗС); σш.в.у – среднеквадратическое значение шума выходного устройства, пересчи-
танное к числу зарядов. В результате выполнения описанных выше процедур формируется массив сум-
марных зарядов, представляющий собой реализацию двумерного нестационарного случайного процесса, результат моделирования которого иллюстрирует рис. 2.
12 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
Рис. 2. Распределение освещенности на фоточувствительной площадке ФПУ (по осям X, Y отложены номера пикселей по горизонтали и вертикали)
Квантование сигнала в блоке АЦП осуществляется с учетом заданной разрядно-
сти и необходимости согласования динамического диапазона входных уровней АЦП с
реальным динамическим диапазоном изменения сигнала на выходе ПЗС. При этом
принимается во внимание возможность «привязки» нижнего уровня динамического
диапазона АЦП к уровню сигнала, соответствующего минимальному темновому заря-
( )ду. Максимальный сигнал Qmax, поступающий на АЦП, должен соответствовать 2N −1
уровням квантования, где N – разрядность АЦП. Тогда произвольному сигналу Q соот-
ветствует k уровней квантования:
( )k
=
INT ⎛⎝⎜⎜
Q
⋅ 2N −1 Qmax
⎟⎟⎞⎠
.
На этапе предварительной обработки сигнала из исходной последовательности
чисел, поступающих с выхода АЦП, выделяются и запоминаются подмассивы чисел,
соответствующие областям полезного сигнала. В выделенных подмассивах осуществля-
ется компенсация сигнала от фона, который в общем случае может быть неоднород-
ным. Возможные алгоритмы, используемые на стадии предварительной обработки сиг-
нала, подробно описаны в [3].
Оценка координат объектов заключается в вычислении координат характерных то-
чек, связанных с наблюдаемыми объектами, например, энергетических центров их изобра-
жений или координат точек, соответствующих максимумам освещенности в изображениях
объектов [3]. Если стоит задача нахождения угловых координат объекта ψ и χ (выраженных в радианах), то пересчет координат изображения точечного объекта x и y (выраженных в числе пространственных периодов матричной структуры) осуществляется с учетом следующих соотношений:
tgψ = x ⋅ d /f ; tgχ = y ⋅ d′/f ,
где d и d′ – пространственный шаг элементов матричной структуры по горизонтали и вертикали соответственно, f – расстояние от объектива до плоскости фокусировки изо-
бражения. Очевидно, что подобное соотношение устанавливает связь между средне-
квадратическими погрешностями измерения координат σx и σy. Статистическая обработка результатов моделирования включает в себя подсчет числа
ложных выбросов, числа пропусков объектов, оценку вероятности пропусков и среднеквадратических погрешностей измерения их координат.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
13
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ …
Определение вероятности выделения ложного выброса Pл.в. осуществляется методом статистических испытаний.
В режиме «поиск» расчетно-прогнозируемая вероятность выделения ложного выбро-
са в кадре может быть определена по формуле
N
∏ [ ]Pлв.к. = 1 − 1 − Pлв.(i) ,
(1)
i =1
где Pлв.(i) – вероятность ложного выброса в i-м элементе кадра. Поскольку (как указывалось
ранее) число зарядов, накапливаемых в ячейках значительно больше 10, то флуктуации за-
рядов в ячейках можно приблизительно описывать нормальным законом распределения с
дисперсией, равной среднему значению накопленных зарядов. При этом вычитание посто-
янной составляющей, соответствующей среднему числу зарядов, обеспечивается в резуль-
тате выполнения процедуры «селекция» на этапе предварительной обработки сигнала. Та-
ким образом,
∫Pлв.(i) =
1 2π
∞ ξ′
exp
⎣⎡
x2
/2⎤⎦
dx ,
где ξ′ – относительный порог в разностном сигнале. Учитывая также, что величина от-
носительного порога ξ′ для каждого элемента устанавливается одинаковой (благодаря
учету «карты» распределения темновых токов и чувствительности по кристаллу на эта-
пе калибровки), можно считать, что Pлв.(i) не зависит от i. Тогда
[ ]Pлв.к. = 1 − 1 − Pлв. (i) N .
(2)
В режиме «слежение» при вычислении Pлв.к. в выражения (1) или (2) следует вместо числа элементов в кадре N подставлять число элементов N′ в стробируемых окнах.
Обработка результатов серии косвенных измерений угловой координаты наблюдаемого объекта ψ сводится первоначально к нахождению результатов измерения ~x , оцен-
ки среднего квадратического отклонения (СКО) результата измерения S(~x ) , довери-
тельной границы случайной погрешности результата измерения εx и доверительной границы неисключенной систематической погрешности (НСП) результата измерения θx , полученных при обработке прямых измерений [4]. Далее определяется оценка результата косвенного измерения ψ(x) и рассчитывается доверительная граница погреш-
ности результата прямого измерения Δ x при доверительной вероятности 0,95:
⎧ ⎪εx ⎪
при
θx
S( x%)
< 0,8;
ψ% = arctg ( x ⋅ d
/
f
),
Δx
= ⎨⎪⎪θx ⎪
при
θx
S( x%)
< 8;
⎪
⎪ ⎩⎪
KSΣ
при
0,8 <
θx
S( x%)
< 8;
вычисляются значения частной производной
∂ψ ∂x
при найденной величине
~x
и опреде-
ляется доверительная граница погрешности результата измерения угловых координат
при доверительной вероятности α = 0,95 :
Δψ
=
∂ψ ∂x
⋅Δx
.
14 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
Итог измерения угловой координаты записывается в виде ψ = ψ% ± Δψ ; α = 0,95 . Аналогичные расчеты производятся для угловой координаты χ .
Заключение
На базе предложенной методики комплексного моделирования возможно на ранней стадии проектирования проводить анализ степени влияния каждого из физических факторов на обнаружительные, метрологические или иные характеристики ОЭС в целом. Модульная структура позволит проводить дальнейшее совершенствование модели путем дополнения ее необходимыми модулями (с учетом реальных потребностей и условий функционирования ОЭС) или путем замены отдельных модулей на альтернативные варианты без изменения остальных частей и структуры в целом.
Литература
1. Андреев А.Л. Автоматизированные телевизионные системы наблюдения. Часть I. Аппаратные средства и элементная база: Учебное пособие для курсового и дипломного проектирования. – СПб: СПбГУИТМО, 2005.
2. Андреев А.Л. Автоматизированные телевизионные системы наблюдения. Часть II. Арифметико–логические основы и алгоритмы: Учебное пособие для курсового и дипломного проектирования. – СПб: СПбГУИТМО, 2005.
3. Андреев А.Л., Ярышев С.Н. Методы моделирования ОЭС с многоэлементными анализаторами изображения: Методические указания к лабораторным работам. – СПб: СПбГУ ИТМО, 2006.
4. Бегункова А.Ф., Парфенов В.Г. Обработка результатов наблюдений при выполнении лабораторных работ физического практикума: Учебное пособие. – Л.: ЛИТМО, 1980.
5. Ишанин Г.Г., Панков Э.Д., Челибанов В.П. Приемники излучения. – СПб: Папирус, 2003.
6. Ферапонтов М.М., Крицына Н.А., Деев Д.Л. Моделирование случайных воздействий на ЭВМ: Учебное пособие. – М.: МИФИ, 1995.
Андреев Андрей Леонидович Лбова Татьяна Павловна
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат технических наук, доцент, and-andr@yandex.ru
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, студентка, tlbova@gmail.com
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
15
УДК 004.942:519.876.5
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ ОПТИКО-ЭЛЕКТРОННОЙ СИСТЕМЫ НАБЛЮДЕНИЯ
ЗА ТОЧЕЧНЫМИ ОБЪЕКТАМИ
А.Л. Андреев, Т.П. Лбова
Дается описание структуры модели, используемой для синтеза, анализа и прогнозирования работы оптико-электронной системы (ОЭС) телевизионного типа, решающей задачи наблюдения за подвижными объектами в автоматическом режиме. В качестве исходных данных задаются параметры входных воздействий объектов и фоновых помех, параметры оптической системы, фотоприемника, дестабилизирующих факторов, алгоритмов и др. Ключевые слова: компьютерное имитационное моделирование систем, телевизионные системы наблюдения, обнаружение и измерение координат точечных объектов.
Введение
Целью настоящей статьи является описание принципов построения комплексной математической модели, используемой для синтеза, анализа, прогнозирования работы телевизионной системы наблюдения за малоразмерными подвижными объектами и обоснования требований к отдельным звеньям системы с учетом необходимых эксплуатационных параметров в реальных условиях функционирования.
К числу важных промежуточных задач, которые нужно решить для достижения поставленной цели, относятся, в частности: − исследование особенностей физических процессов, протекающих в многоэлемент-
ных структурах фотоприемных устройств с учетом физических свойств различных фоточувствительных материалов; − исследование физической природы источников помех, действующих в различных звеньях оптико-электронной системы и др.; − исследование различных алгоритмов обработки сигналов, используемых для автоматического обнаружения и измерения параметров объектов.
Универсальность модели предполагает возможность ее использования для различных условий функционирования подобных систем – с использованием наземного, космического или иного базирования, различных дистанций наблюдения, энергетических и динамических параметрах объектов и др.
Структура модели
Структура математической модели ОЭС представлена на рис. 1. Задача моделирования ОЭС разбивается на следующие этапы: моделирование сигнально-фоновой ситуации; моделирование различных видов помех; моделирование процедур предварительной обработки сигнала (с учетом искажений, связанных с квантованием сигналов при аналого-цифровом преобразовании); моделирование различных алгоритмов целевой обработки сигналов (связанных, в общем случае, с решением задач обнаружения, распознавания и оценки параметров исследуемых объектов, а также оценки метрологических параметров проектируемой системы с учетом заданных доверительных интервалов). Ввод исходных данных осуществляется в интерактивном режиме. При этом вводимые переменные можно подразделить на следующие группы: − параметры оптической системы; − параметры фотоприемного устройства; − параметры алгоритма первичной обработки сигнала;
10 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
− условия наблюдения; − параметры, определяющие режим статистической обработки полученных реализа-
ций. Предварительная обработка исходных данных предполагает преобразование вво-
димых величин и формирование на их основе переменных и констант, используемых при моделировании.
Рис. 1. Структурная схема модели
В общем случае фоновая составляющая сигнала представляет собой случайное поле яркости. Исходным материалом для компьютерного моделирования случайного поля яркости являются независимые числа, полученные генератором случайных чисел ЭВМ. Их совокупность можно рассматривать как случайное δ-поле – обобщение понятия дискретного белого шума на случай нескольких переменных. Моделирование δ-поля осуществляется следующим образом: пространственно-временной координате (xi , y j ,tk ) ставятся в соответствие выборочные значения из датчика нормальных слу-
чайных чисел с параметрами (0, 1). Для моделирования дискретных реализаций однородного стационарного случай-
ного поля, заданного в конечном пространстве и времени, предлагается использовать метод канонического разложения поля в пространственно-временной ряд Фурье. Неограниченные дискретные реализации однородного стационарного случайного поля можно получить с помощью алгоритма пространственно-временного скользящего суммирования δ-поля [1].
При формировании массива темновых зарядов учитывается аддитивная составляющая «геометрического» шума. Среднее значение темнового заряда, накапливаемого в ячейке, определяется выражением [2]
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
11
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ …
QT = iT AэлTн ,
где iT – средняя плотность темнового тока накопления заряда, характерная для данного полупроводникового материала при рабочей температуре кристалла; Аэл – площадь электрода, под которым осуществляется накопление.
Формирование массива сигнальных зарядов осуществляется с учетом случайного распределения темновых токов и чувствительности по кристаллу ПЗС с заданным интервалом корреляции «геометрического» шума. Если, например, распределение освещенности от точечного объекта по фоточувствительной площадке аппроксимируется гауссоидой вращения, то распределение сигнальных зарядов определяется в соответствии с выражением
E(xi ,
yj
)
=
Ф 2πR2
exp
⎢⎡− ⎣
(x
−
x0 )2 + (y 2R2
−
y0 )2
⎤ ⎥ ⎦
,
где Ф – поток оптического излучения, создающего изображение точечного объекта; R –
– радиус кружка рассеяния объектива, определяемый на уровне 0,606 Emax; x0, y0 – ко-
ординаты центра пятна рассеяния.
Динамику перемещения наблюдаемых объектов можно задать в параметрическом
виде с помощью коэффициентов А, В, С (первой, второй и третьей степени) полинома,
аппроксимирующего изменение во времени модуля угловых перемещений объектов
S(t) = A ⋅ t + B ⋅ t 2 + C ⋅ t3 .
Направление перемещения задается посредством угла α между направлением оси X и проекцией вектора скорости в плоскости анализа изображения.
При моделировании массива шумовых составляющих в ячейках ПЗС учитывается совокупность сигнальной, фоновой и темновой составляющих зарядов в каждой ячейке. При этом шумовая составляющая заряда в каждой ячейке определяется путем многократного задания флуктуаций зарядов по случайному закону с заданным распределением с учетом влияния всех вышеперечисленных факторов.
Флуктуация числа зарядов подчиняется закону Пуассона. Однако, поскольку чис-
ло накапливаемых зарядов достаточно велико ( nq > 100 ), распределение Пуассона
можно аппроксимировать нормальным законом распределения с дисперсией, равной среднему значению. Тогда среднеквадратическое значение шума на выходе ФПУ равно
σш =
σ2 ш.т
+
σ2 ш.ф
+
σ2 ш.с
+
σ2 ш.г
+
σ2 ш.п
+
σ2 ш.в.у.
,
где σш.т – среднеквадратическое число зарядов в i-ой ячейке ФПУ, характеризующее флуктуацию числа темновых зарядов; σш.ф – среднеквадратическое значение флуктуа-
ции числа зарядов, обусловленных фоновым излучением; σш.с – среднеквадратическое
значение флуктуации числа сигнальных зарядов, накапливаемых в ячейках ФПЗС, ха-
рактеризующее фотонный шум; σш.г – среднеквадратическое отклонение среднего
числа темновых зарядов в различных ячейках ФПУ, характеризующее геометрический
шум; σш.п – среднеквадратическое значение флуктуаций заряда при переносе (для
ФПЗС); σш.в.у – среднеквадратическое значение шума выходного устройства, пересчи-
танное к числу зарядов. В результате выполнения описанных выше процедур формируется массив сум-
марных зарядов, представляющий собой реализацию двумерного нестационарного случайного процесса, результат моделирования которого иллюстрирует рис. 2.
12 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
Рис. 2. Распределение освещенности на фоточувствительной площадке ФПУ (по осям X, Y отложены номера пикселей по горизонтали и вертикали)
Квантование сигнала в блоке АЦП осуществляется с учетом заданной разрядно-
сти и необходимости согласования динамического диапазона входных уровней АЦП с
реальным динамическим диапазоном изменения сигнала на выходе ПЗС. При этом
принимается во внимание возможность «привязки» нижнего уровня динамического
диапазона АЦП к уровню сигнала, соответствующего минимальному темновому заря-
( )ду. Максимальный сигнал Qmax, поступающий на АЦП, должен соответствовать 2N −1
уровням квантования, где N – разрядность АЦП. Тогда произвольному сигналу Q соот-
ветствует k уровней квантования:
( )k
=
INT ⎛⎝⎜⎜
Q
⋅ 2N −1 Qmax
⎟⎟⎞⎠
.
На этапе предварительной обработки сигнала из исходной последовательности
чисел, поступающих с выхода АЦП, выделяются и запоминаются подмассивы чисел,
соответствующие областям полезного сигнала. В выделенных подмассивах осуществля-
ется компенсация сигнала от фона, который в общем случае может быть неоднород-
ным. Возможные алгоритмы, используемые на стадии предварительной обработки сиг-
нала, подробно описаны в [3].
Оценка координат объектов заключается в вычислении координат характерных то-
чек, связанных с наблюдаемыми объектами, например, энергетических центров их изобра-
жений или координат точек, соответствующих максимумам освещенности в изображениях
объектов [3]. Если стоит задача нахождения угловых координат объекта ψ и χ (выраженных в радианах), то пересчет координат изображения точечного объекта x и y (выраженных в числе пространственных периодов матричной структуры) осуществляется с учетом следующих соотношений:
tgψ = x ⋅ d /f ; tgχ = y ⋅ d′/f ,
где d и d′ – пространственный шаг элементов матричной структуры по горизонтали и вертикали соответственно, f – расстояние от объектива до плоскости фокусировки изо-
бражения. Очевидно, что подобное соотношение устанавливает связь между средне-
квадратическими погрешностями измерения координат σx и σy. Статистическая обработка результатов моделирования включает в себя подсчет числа
ложных выбросов, числа пропусков объектов, оценку вероятности пропусков и среднеквадратических погрешностей измерения их координат.
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
13
РАЗРАБОТКА СТРУКТУРЫ КОМПЛЕКСНОЙ МОДЕЛИ …
Определение вероятности выделения ложного выброса Pл.в. осуществляется методом статистических испытаний.
В режиме «поиск» расчетно-прогнозируемая вероятность выделения ложного выбро-
са в кадре может быть определена по формуле
N
∏ [ ]Pлв.к. = 1 − 1 − Pлв.(i) ,
(1)
i =1
где Pлв.(i) – вероятность ложного выброса в i-м элементе кадра. Поскольку (как указывалось
ранее) число зарядов, накапливаемых в ячейках значительно больше 10, то флуктуации за-
рядов в ячейках можно приблизительно описывать нормальным законом распределения с
дисперсией, равной среднему значению накопленных зарядов. При этом вычитание посто-
янной составляющей, соответствующей среднему числу зарядов, обеспечивается в резуль-
тате выполнения процедуры «селекция» на этапе предварительной обработки сигнала. Та-
ким образом,
∫Pлв.(i) =
1 2π
∞ ξ′
exp
⎣⎡
x2
/2⎤⎦
dx ,
где ξ′ – относительный порог в разностном сигнале. Учитывая также, что величина от-
носительного порога ξ′ для каждого элемента устанавливается одинаковой (благодаря
учету «карты» распределения темновых токов и чувствительности по кристаллу на эта-
пе калибровки), можно считать, что Pлв.(i) не зависит от i. Тогда
[ ]Pлв.к. = 1 − 1 − Pлв. (i) N .
(2)
В режиме «слежение» при вычислении Pлв.к. в выражения (1) или (2) следует вместо числа элементов в кадре N подставлять число элементов N′ в стробируемых окнах.
Обработка результатов серии косвенных измерений угловой координаты наблюдаемого объекта ψ сводится первоначально к нахождению результатов измерения ~x , оцен-
ки среднего квадратического отклонения (СКО) результата измерения S(~x ) , довери-
тельной границы случайной погрешности результата измерения εx и доверительной границы неисключенной систематической погрешности (НСП) результата измерения θx , полученных при обработке прямых измерений [4]. Далее определяется оценка результата косвенного измерения ψ(x) и рассчитывается доверительная граница погреш-
ности результата прямого измерения Δ x при доверительной вероятности 0,95:
⎧ ⎪εx ⎪
при
θx
S( x%)
< 0,8;
ψ% = arctg ( x ⋅ d
/
f
),
Δx
= ⎨⎪⎪θx ⎪
при
θx
S( x%)
< 8;
⎪
⎪ ⎩⎪
KSΣ
при
0,8 <
θx
S( x%)
< 8;
вычисляются значения частной производной
∂ψ ∂x
при найденной величине
~x
и опреде-
ляется доверительная граница погрешности результата измерения угловых координат
при доверительной вероятности α = 0,95 :
Δψ
=
∂ψ ∂x
⋅Δx
.
14 Научно-технический вестник Санкт-Петербургского государственного университета
информационных технологий, механики и оптики, 2009, № 5(63)
А.Л. Андреев, Т.П. Лбова
Итог измерения угловой координаты записывается в виде ψ = ψ% ± Δψ ; α = 0,95 . Аналогичные расчеты производятся для угловой координаты χ .
Заключение
На базе предложенной методики комплексного моделирования возможно на ранней стадии проектирования проводить анализ степени влияния каждого из физических факторов на обнаружительные, метрологические или иные характеристики ОЭС в целом. Модульная структура позволит проводить дальнейшее совершенствование модели путем дополнения ее необходимыми модулями (с учетом реальных потребностей и условий функционирования ОЭС) или путем замены отдельных модулей на альтернативные варианты без изменения остальных частей и структуры в целом.
Литература
1. Андреев А.Л. Автоматизированные телевизионные системы наблюдения. Часть I. Аппаратные средства и элементная база: Учебное пособие для курсового и дипломного проектирования. – СПб: СПбГУИТМО, 2005.
2. Андреев А.Л. Автоматизированные телевизионные системы наблюдения. Часть II. Арифметико–логические основы и алгоритмы: Учебное пособие для курсового и дипломного проектирования. – СПб: СПбГУИТМО, 2005.
3. Андреев А.Л., Ярышев С.Н. Методы моделирования ОЭС с многоэлементными анализаторами изображения: Методические указания к лабораторным работам. – СПб: СПбГУ ИТМО, 2006.
4. Бегункова А.Ф., Парфенов В.Г. Обработка результатов наблюдений при выполнении лабораторных работ физического практикума: Учебное пособие. – Л.: ЛИТМО, 1980.
5. Ишанин Г.Г., Панков Э.Д., Челибанов В.П. Приемники излучения. – СПб: Папирус, 2003.
6. Ферапонтов М.М., Крицына Н.А., Деев Д.Л. Моделирование случайных воздействий на ЭВМ: Учебное пособие. – М.: МИФИ, 1995.
Андреев Андрей Леонидович Лбова Татьяна Павловна
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, кандидат технических наук, доцент, and-andr@yandex.ru
– Санкт-Петербургский государственный университет информационных технологий, механики и оптики, студентка, tlbova@gmail.com
Научно-технический вестник Санкт-Петербургского государственного университета информационных технологий, механики и оптики, 2009, № 5(63)
15