УДК 519.248, 621.384.3 doi:10.15217/issnl684-8853.2017.2.13
НЕРАВНОМЕРНОСТЬ ЧУВСТВИТЕЛЬНОСТИ ФОТОПРИЕМНИКА КАК МУЛЬТИПЛИКАТИВНАЯ ПОМЕХА ПРИ МЕЖКАДРОВОЙ ОБРАБОТКЕ В ПРЕЦИЗИОННЫХ ТЕЛЕВИЗИОННЫХ СИСТЕМАХ
В. А. Рогачева, канд. техн. наук, старший научный сотрудник
A. А. Закутаева, старший научный сотрудник, начальник лаборатории
B. Д. Лиференкоа, доктор техн. наук, профессор М. О. Колбанев6, доктор техн. наук, профессор
аВоенно-космическая академия им. А. Ф. Можайского, Санкт-Петербург, РФ 6Санкт-Петербургский государственный экономический университет, Санкт-Петербург, РФ
Постановка проблемы: в телевизионных системах визуализации достаточно широко используется компенсация неравномерности чувствительности фотоприемника. Однако, несмотря на видимые улучшения, компенсация неравномерности чувствительности фотоприемника не изменяет пороговой чувствительности. Цель: оценка влияния неравномерности чувствительности фотоприемника при межкадровой обработке в прецизионных телевизионных системах обнаружения. Результаты: получены выражения параметров распределений решающих статистик при влиянии неравномерности чувствительности фотоприемника. Показано, что компенсация неравномерности чувствительности фотоприемника приводит к перераспределению влияния и не обеспечивает полного устранения неравномерности чувствительности фотоприемника. Практическая значимость: полученные результаты позволяют оценить степень влияния неравномерности чувствительности фотоприемника для основных режимов работы прецизионной телевизионной системы обнаружения, а также показывают, что в режиме ограничения фоновыми шумами при малых фоновых шумах параметр распределения решающей статистики не зависит от неравномерности чувствительности фотоприемника.
Ключевые слова — межкадровая обработка, обнаружение сигналов, режимы работы фотоприемника, неравномерность чувствительности фотоприемника, компенсация неравномерности чувствительности фотоприемника.
Введение
Прецизионные телевизионные системы (ПТС) характеризуются, помимо прочих особенностей, существенной неравномерностью выходного сигнала фотоприемника — «геометрическим шумом» [1].
Эта неравномерность определяется двумя параметрами: неравномерностью темнового тока, проявляющейся при отсутствии внешней освещенности, и неравномерностью чувствительности фотоприемника (НЧФ). Причиной этих не-равномерностей, как правило, является разброс технологических параметров при изготовлении фотоприемника.
Неравномерность чувствительности фотоприемника представляет собой помеху, которая при равномерной засветке поля фотоприемника дает различные значения сигнала для различных элементов [1].
Относительная неравномерность выходного сигнала фотоприемника согласно ГОСТ 28953-91 [2] оценивается следующим образом: Д = (Umax -- U . )/(U + U ), где U и U . — макси-
mi^ max mm" H max min
мальное и минимальное значения сигналов отдельных элементов фотоприемника.
Компенсация НЧФ как мультипликативной помехи производится умножением сигнала
каждого элемента фотоприемника на величину, обратную его чувствительности (в одной или нескольких точках), и является гораздо более сложной операцией, чем компенсация темнового тока, представляющего собой аддитивную помеху [3-7].
В телевизионных системах визуализации неравномерность выходного сигнала достаточно хорошо корректируется (компенсируется) [3-7]. Однако при обнаружении сигналов компенсация НЧФ не улучшает пороговой чувствительности [8].
Ранее [9, 10] были рассмотрены задачи обнаружения в ПТС без учета НЧФ. В данной работе рассмотрим особенности влияния НЧФ на характеристики обнаружения в различных режимах работы ПТС с компенсацией и без компенсации НЧФ.
Ограничимся случаем межкадровой обработки сигнала фотоприемника, когда производится поэлементное сравнение соседних кадров и можно сравнить элементы фотоприемника с одинаковой чувствительностью. Оценки влияния НЧФ на характеристики обнаружения ПТС уже рассматривались, но только для случая внутрикад-ровой обработки [11].
Анализ проведем на основе модели с двумя источниками шумов, применяемой в видеоинформационных (телевизионных) системах [12].
Такая модель приводит к появлению многих режимов работы ПТС, обусловленных взаимным соотношением токов и шумов на выходе фотоприемника [13].
Световой отклик каждого элемента фотоприемника, как правило, различен, тем самым и образуется НЧФ. Введем относительный коэффициент чувствительности для ¿-го элемента фотоприемника, имеющий различные значения для различных элементов: = /иш1п. Значение относительного коэффициента чувствительности будет лежать в диапазоне 1 < < ишах/иш1п.
Рассмотрим влияние НЧФ раздельно для каждого режима, в котором может работать система. При этом полагаем, что размеры проекции обнаруживаемого объекта могут быть любыми в пределах матрицы фотоприемника. Для обработки имеется несколько кадров М, в которых нет полезного сигнала («опорных» кадров), и несколько кадров Ы, в которых сигнал присутствует. Полезный сигнал занимает некоторое количество элементов, его амплитуда может изменяться от элемента к элементу, ее значение неизвестно и стационарно от кадра к кадру. Характеристики НЧФ (среднее значение чувствительности и ее разброс) считаем стационарными по полю фотоприемника.
Для определения «контраста» между кадрами вычисляется некоторая функция (статистика), зависящая от режима, в котором работает ПТС [9, 11, 13].
Режим ограничения внутренним шумом
В этом случае наличие НЧФ приводит к изменению параметров модели с нормальным распределением по сравнению со случаем отсутствия НЧФ.
Проверяемые гипотезы о наличии Н1 или отсутствии Н0 сигнала будут иметь вид:
Н,
Нл
: х е N(d, ст2 j, у е N(d, ст2 j; : х е N(d, ст2 у е N(d + ys, ст2 j,
где d — темновой ток; ст2 — дисперсия внутреннего шума; у — относительный коэффициент чувствительности данного элемента фотоприемника; в — сигнал.
После вычитания темнового тока (полагаем, что дисперсия оценки темнового тока пренебрежимо мала) модель изменится и будет иметь следующий вид:
Щ
: х е N ^0, ст2 е N (о, ст2 ); : х е N(о, ст2 у е N(ys, ст2
Для принятия решения о наличии или отсутствии сигнала для этого режима необходимо вычислить статистику Стьюдента [9, 11, 13]:
M
К* У >1
/ (м )>
где в числителе вычисляется среднее по сигнальной области с числом элементов M, а в знаменателе — оценка дисперсии по фоновой области окна обработки с числом элементов N.
Данная статистика имеет нецентральное распределение Стьюдента с параметром нецентральности (pt = s/a (при отсутствии неравномерности чувствительности) [9-13]. Этот параметр нецентральности пропорционален отношению сигнал/ шум.
Определим, каким образом влияет НЧФ на параметр нецентральности распределения Стьюдента.
Вследствие наличия НЧФ числитель параметра нецентральности изменится и примет значение
Ф? (0 = (*У г Р )= (® / у I)/
где I — номер элемента фотоприемника.
Таким образом, наличие НЧФ приводит к изменению параметра нецентральности и, следовательно, к изменению мощности (вероятности правильного обнаружения) алгоритма Стьюдента.
После компенсации НЧФ модель будет иметь следующий вид:
Н,
0 : х е N (О, а2 ), у е N (о, а2 );
: х е N ^0, ст2 ), у е N (s, ст2 / у2 ).
Параметр нецентральности распределения Стьюдента не изменится:
ф'/У (0 = (*)/>/а2 / у? = («у;)/ а.
Таким образом, компенсация НЧФ не устраняет влияние неравномерности чувствительности в режиме ограничения внутренним шумом.
Режим обнаружения случайного сигнала
Для этого случая НЧФ будет проявляться следующим образом:
н.
0 : х е N ^(f + уЪ, ст2 +уab j, г/ е N (d + yb, ст2 +yab j; Щ : х е N (d + yb, ст2 +у abj,
у еЫ+ уЬ, ст2 + уаЪ + уав|,
где аЬ — дисперсия фонового шума; Ь — фоновый ток.
После вычитания темнового тока модель изменится и будет иметь следующий вид:
Н.
Нл
: х е N{уЪ, ст2 + уаЪг/ е N{уЪ, ст2 + уаЪ : х е N{уЪ, ст2 + уаЪг/ е N{уЪ, ст2 + уаЬ + уав
Для принятия решения о наличии или отсутствии полезного сигнала необходимо вычислить статистику Фишера [9-13]:
г =
N
- у)2 / 1)
7=1
М _ ^
-х) / (М-1)
.1=1
-\2
В этой статистике используются оценки средних и дисперсии в «опорных» и «сигнальных» кадрах.
Данная статистика имеет центральное распределение Фишера с параметром, пропорциональным отношению сигнал/шум ф4 = (ст2 + аЬ + ав)/ (ст2 + аЬ) (при отсутствии неравномерности чувствительности) [9-13].
С учетом влияния НЧФ данный параметр будет равен
ф! (0 =
о о
ст + у+ у;ав / ст + у;аЬ
Таким образом, наличие НЧФ приводит к изменению мощности алгоритма Фишера.
После компенсации НЧФ модель будет иметь следующий вид:
Н0:х е N(ъ, ст2 / у2 +аЬ / у),
у £ N{ъ, ст2 / у2 +аЬ/ у);
Я^х£N(ъ, ст2 /у2 +аЬ/у),
у £ N(ь, ст2 / у2 +аЬ / у + / у),
поэтому параметр нецентральности распределения Фишера также изменится:
ф:
^ (О -
о о
ст / у;+аЬ + ав / ст / у;+аЬ
Таким образом, компенсация НЧФ не устраняет влияние неравномерности чувствительности в режиме обнаружения случайного сигнала.
Необходимо отметить, что если внутренние шумы малы, т. е. ст2 << аЬ + ав, то параметр нецентральности не зависит от НЧФ.
Режим ограничения фоновым шумом
Для этого случая НЧФ будет проявляться следующим образом:
Н0: х £ N(уЬ, уаЪ), -у £ N(уЬ, уаЪ); Н^. хе N(yЬ, уаЪ), уе Ы(уЬ + у8, уаЬ + уав),
где аЬ — дисперсия фонового шума; ав — дисперсия сигнального шума.
Для принятия решения о наличии или отсутствии полезного сигнала необходимо вычислить статистику [8, 10]:
Ч =
N
) /м
.}=1
м
)2 / м
г=1
Данная статистика имеет дважды нецентральное распределение Фишера с параметром распределения фд = [(Ь + в)2 + а(Ь + в)]/(Ь2 + в) (при отсутствии неравномерности чувствительности) [9, 10].
В случае наличия НЧФ параметр распределения, пропорциональный отношению сигнал/ шум, изменится:
Фа (0 =
у? (Ъ + в)2 +у;аЬ + у;£
У^2 +УьаЬ
Таким образом, НЧФ приводит к изменению мощности алгоритма.
После компенсации НЧФ модель будет иметь следующий вид:
Н0: х £ N{Ъ, аЪ / у), -у £ N(Ь, аЬ / у); : х £ N(Ь, аЬ / у), у £ N(Ь + в, аЬ / у + ав / у),
поэтому параметр нецентральности дважды нецентрального распределения Фишера примет следующую форму:
Фд/У (0 = (ь+в Г + аЪ / У г + «в / Уг
Ь2+ аЬ / у;
Таким образом, компенсация НЧФ не устраняет влияние неравномерности чувствительности в режиме ограничения фоновым шумом.
Следует отметить, что если шумы, вызываемые фоном, малы, т. е. Ь2 >> аЬ, то параметр нецентральности не зависит от НЧФ в этом режиме.
Режим общий
Для этого случая НЧФ будет проявляться следующим образом:
Н,
х £ N^ + уЬ, ст2 + уаЪ^, у £ N^ + уЪ, ст2 + уаЬ
Щ: х е N(d + yb, ст2 + yabj, у е N|d + yb + ys, ст2 + yab + yas j.
После вычитания темнового тока модель изменится и будет иметь следующий вид:
Н,
х е N(yb, ст2 + yabj, г/ е N{уЪ, ст2 + yabj; Н^. х е N[уЪ, ст2 + yab j, у е N^yb + ys, ст2 + yab + yas j.
Для принятия решения о наличии или отсутствии полезного сигнала необходимо вычислить модифицированную статистику Фишера [9, 10]:
г =
N
-х)2 /N /=1
М
£(х; -X) /(М-1)
í=1
-\2
Данная статистика имеет нецентральное распределение Фишера с параметром фг = в2 / (ст2 + а& + ав) (при отсутствии НЧФ) [9, 10].
Для этого режима при наличии НЧФ параметр распределения, пропорциональный отношению сигнал/шум, будет равен
Фг (¿) =
2 2 Y;s
ст + y^ab + у ¿as
Таким образом, НЧФ приводит к изменению мощности алгоритма.
После компенсации НЧФ модель будет иметь следующий вид:
Нп
Нл
у. хе N(ь, ст2 /у2 +аЬ/у), у е ^(Ь, ст2 /у2 +аЬ/у); ^ хеN(ь, ст2 /у2 +аЬ/у), у е N+ ст2 /у2 +аЬ/у + /у).
Параметр нецентральности дважды нецентрального распределения Фишера будет иметь вид
Ф1/7 (I) = [в2 ] / [ст2 / у2 + аЪ / у; + ав / у;
Таким образом, компенсация НЧФ не устраняет влияние неравномерности чувствительности в общем режиме.
Для каждого из режимов работы ПТС получены выражения параметров распределений решающих статистик при отсутствии НЧФ, учитывающие влияние НЧФ, а также при компенсации НЧФ (таблица). Параметры для каждого из рассмотренных режимов пропорциональны отношению сигнал/шум и зависят от НЧФ.
Наличие НЧФ изменяет постоянную и случайную составляющие как сигнала, так и фона и не затрагивает внутренние шумы. Компенсация НЧФ перераспределяет (переносит) влияние НЧФ и не приводит к полному устранению влияния НЧФ.
Таким образом, наличие НЧФ вызывает изменение параметров распределений решающих статистик, что в свою очередь приводит к изменению вероятностей правильного обнаружения.
Влияние неравномерности чувствительности на параметры функций вероятностей правильного обнаружения
Режим Параметр статистики при условии
отсутствия НЧФ наличия НЧФ компенсации НЧФ
Ограничения внутренним шумом Ф^ К ) ф? « - Ы )/у/) Ф^ (í) = (s)/Jст2 / y¿
Обнаружения случайного сигнала Ф^ = [ст2 +ab + as] / [ст2 + ab] Ф1 (¿) = [°2 + !iab + у ¿as] / [ст2 + у tab ] Ф1/^ (í) = [ст2 / Y¿ + ab + as] / [ст2 / y¿ + ab]
Ограничения фоновым шумом Ф^= (b + s)2 + a(b + s) / [b2 + ab ] Фд(¿) = T¿ (& + s)2 + ytab + ytas / [yfb2+ytab] Фд/У (¿) = [(^ + s)2 + ab / y ¿+as / Y¿]/ [b2 +ab / y¿ ]
Общий = s2 / (ст2 + ab + as) Фl (í) = [ifs2 ] / [ст2 + y¿ab + y¿as ] Ф1/7 (í) = [s2 ] / [ст2 / y2 +ab/у г + as/Yi ]
Следует отметить, что в режиме ограничения фоновыми шумами при малых фоновых шумах параметр распределения решающей статистики не зависит от НЧФ. Аналогичная ситуация складывается и в режиме обнаружения случайного сигнала при малых внутренних шумах.
Заключение
В работе получены выражения параметров распределений решающих статистик при влиянии неравномерности чувствительности фотоприемника.
Литература
1. Solid State Imaging/Ed. by P. Jespers, F. Van de Wiele, M. White. — Noordhoff-Leyden, 1976. — 752 p.
2. ГОСТ 28953-91. Приборы фоточувствительные с переносом заряда. Методы измерения параметров. — М.: Изд-во стандартов, 2004. — 31 с.
3. Брондз Д. С., Харитонова Е. Н. Коррекция геометрического шума МФПУ с помощью аппроксимации методом наименьших квадратов передаточных характеристик матрицы полиномом Т-порядка // Журнал радиоэлектроники. 2008. № 11. С. 1-29.
4. Сергунов А. А. Методы коррекции неравномерности чувствительности инфракрасных матричных фотоприемников // Изв. вузов. Приборостроение. 2009. Вып. 8. С. 38-42.
5. Scribner D. A., Kruer M. R., and Griedly C. J. Phisical Limitations to Nonuniformity Correction in Focal Plane Arrays//Proc. SPIE. 1987. Vol. 865. P. 185201.
6. Wang R., Chen P., and Tsien P. An Improved Nonuniformity Correction Algorithm for Infrared Focal Plane Arrays which is Easy to Implement// Infrared Physics and Technology. 1998. Vol. 39. P. 15-21.
7. Emerson G. P., Little S. J. Flat-Fielding for CCDs in AAVSO Observations//AAVSO. 1999. Vol. 27. P. 49-54.
Показано, что операция компенсации не обеспечивает полного устранения неравномерности чувствительности фотоприемника.
Полученные результаты позволяют оценить степень влияния неравномерности чувствительности фотоприемника для основных режимов работы прецизионной телевизионной системы обнаружения, а также показывают, что в режиме ограничения фоновыми шумами при малых фоновых шумах параметр распределения решающей статистики не зависит от неравномерности чувствительности фотоприемника.
8. Рогачев В. А., Фантиков О. И. Влияние неравномерности чувствительности фотоприемника на пороговую чувствительность ТВ-систем // Техника средств связи. Сер. Техника телевидения. 1987. № 2. С. 48-55.
9. Колбанев М. О., Рогачев В. А. Оптимизация выделения полезного сигнала в многорежимных информационных системах // Вопросы радиоэлектроники. Сер. РЛТ. 2010. Вып. 1. С. 92-101.
10. Колбанев М. О., Рогачев В. А. Анализ проблемы обнаружения в инфракрасных системах // Информационно-управляющие системы. 2010. № 5. С. 51-54.
11. Колбанев М. О., Рогачев В. А., Закутаев А. А. Неравномерность чувствительности фотоприемника и ее компенсация при обнаружении в инфракрасных системах // Вопросы радиоэлектроники. Сер. Техника телевидения. 2015. Вып. 5. С. 79-88.
12. Цыцулин А. К., Адамов Д. Ю., Манцветов А. А., Зубакин И. А. Твердотельные телекамеры: Накопление качества информации. — СПб.: Изд-во СПбГЭТУ «ЛЭТИ», 2014. — 272 с.
13.Рогачев В. А. Особенности обработки сигнала в прецизионных телевизионных системах // Вопросы радиоэлектроники. Сер. Техника телевидения. 2011. Вып. 2. С. 18-26.
14. Lehmann E. L., Romano J. P. Testing Statistical Hypotheses. — Springer, 2005. — 800 p.
UDC 519.248, 621.384.3 doi:10.15217/issn1684-8853.2017.2.13
Photodetector Sensitivity Non-Uniformity as a Multiple Interference in High-Precision Television Systems with Interframe Processing
Rogachev V. A.a, PhD, Tech., Senior Researcher, [email protected]
Zakutaev A. A.a, Senior Researcher, Head of a Research Laboratory, [email protected]
Liferenko V. D.a, Dr. Sc., Tech., Professor, [email protected]
Kolbanev M. O.b, Dr. Sc., Tech., Professor, [email protected]
aA. F. Mozhaiskii Military Space Academy, 13, Zhdanovskaia St., 197198, Saint-Petersburg, Russian Federation bState University of Economics, 21, Sadovaya St., 191023, Saint-Petersburg, Russian Federation
Introduction: Television visualization systems often compensate the non-uniformity of photodetector sensitivity. However, in spite of the apparent improvement, this compensation does not alter the threshold sensitivity. Purpose: The goal of this research is assessing the influence of non-uniform photodetector sensitivity over inter-frame processing in high-precision TV detection systems. Results: We
have obtained the distribution parameters for decisive statistics when the photodetector sensitivity is non-uniform. The compensation of this non-uniformity redistributes its influence and does not provide its complete elimination. Practical relevance: The obtained results allow you to assess the degree of influence of photodetector sensitivity non-uniformity for the main operation modes of a high-precision TV detection system. In the background limitation mode, when the background noise is low, the distribution parameter of the decisive statistics does not depend on the non-uniformity.
Keywords — Interframe Processing, Signal Detection, Photodetector Operation Modes, Photodetector Sensitivity Non-Uniformity, Photodetector Sensitivity Non-Uniformity Compensation.
References
1. Solid State Imaging. Ed. by P. Jespers, F. Van de Wiele, M. White. Noordhoff-Leyden, 1976. 752 p.
2. State Standard 28953-91. Methods of Parameters Measurement. Photosensitive Devices with Charge Transfer. Moscow, Standartov Publ., 2004. 31 p. (In Russian).
3. Brondz D. S., Haritonova E. N. Nonuniformity Correction of Matrix Photodetector with Least Squares Approximation of the Transfer Characteristics of the Matrix using T-oder Polynom. Zhurnal radioelektroniki [Journal of Radio Electronics], 2008, no. 11, pp. 1-29 (In Russian).
4. Sergunov A. A. Methods of Nonuniformity Sensitivity Correction of Matrix Infrared Photodetectors. Izvestiia vuzov. Priborostroenie [Journal of Instrument Engineering], 2009, no. 8, pp. 38-42 (In Russian).
5. Scribner D. A., Kruer M. R., and Griedly C. J. Phisical Limitations to Nonuniformity Correction in Focal Plane Arrays. Proc. SPIE, 1987, vol. 865, pp. 185-201.
6. Wang R., Chen P., and Tsien P. An Improved Nonuniform-ity Correction Algorithm for Infrared Focal Plane Arrays which is Easy to Implement. Infrared Physics and Technology, 1998, vol. 39, pp. 15-21.
7. Emerson G. P., Little S. J. Flat-Fielding for CCDs in AAVSO Observations. AAVSO, 1999, vol. 27, pp. 49-54.
8. Rogachev V. A., Fantikov O. I. The Influence of Nonuni-formity Sensitivity of the Photodetector on the Threshold Sensitivity of TV Systems. Tekhnika sredstv sviazi. Ser. Tekhnika televideniia [Communications Equipment Tech-
nology. Ser. Television Technology], 1987, no. 2, pp. 48-55 (In Russian).
9. Kolbanev M. O., Rogachev V. A. Useful Signal Detection Optimization at the Multimode Information Systems. Voprosy radioelektroniki. Ser. Radiolokatsionnaia tekhnika [Questions of Radio-electronics. Ser. Radio Locating Technique], 2010, iss. 1, pp. 92-101 (In Russian).
10. Kolbanev M. O., Rogachev V. A. An Analysis of the Detection Problem in Infrared Systems. Informatsionno-upravli-aiushchie sistemy [Information and Control Systems], 2010, no. 5, pp. 51-54 (In Russian).
11. Kolbanev M. O., Rogachev V. A. Zakutaev A. A. Nonuni-formity Sensitivity of the Photodetector and his Compensation at the Infrared Detection Systems. Voprosy radioelek-troniki. Ser. Tekhnika televideniia [Questions of Radio-electronics. Ser. Television Technology], 2015, iss. 5, pp. 79-88 (In Russian).
12. Tsitsulin A. K., Adamov D. U., Mantsvetov A. A., Zubakin I. A. Tverdotel'nye telekamery: Nakoplenie kachestva informatsii [Solid State TV Cameras: the Accumulation of Information]. Saint-Petersburg, LETI Publ., 2014, 272 p. (In Russian).
13. Rogachev V. A. Signal Processing Features in Precision Television Systems. Voprosy radioelektroniki. Ser. Tekhnika televideniia [Questions of Radio-electronics. Ser. Television Technology], 2011, iss. 2, pp. 18-26 (In Russian).
14. Lehmann E. L., Romano J. P. Testing Statistical Hypotheses. Springer, 2005. 800 p.