CC BY
218
УДК 621.384.3
РАСЧЕТ МИНИМАЛЬНОЙ РАЗРЕШАЕМОЙ РАЗНОСТИ ТЕМПЕРАТУР
ТЕПЛОВИЗИОННЫХ ПРИБОРОВ
АхмедМ.Л. Аль-Мзирави, Колобродов В.Г., Микитенко В.И.
Национальный технический университет Украины «Киевский политехнический институт», г. Киев, Украина e-mail: [email protected]
Статья посвящена повышению точности определения минимальной разрешаемой разности температур тепловизионных систем за счет совершенствования модели визуального восприятия термографических изображений оператором. Рассмотрен алгоритм получения уравнения для расчета минимальной разрешаемой разности температур, в основе которого лежит более достоверная аппроксимация модуляционной передаточной функции зрительной системы, предложенная Д. У. Шульцем. Получено новое выражение для коэффициента ширины полосы для аппроксимации Шульца, которое не зависит от углового размера штриха миры Фуко.
Ключевые слова: тепловидение, оценка эффективности, минимальная разрешаемая разность температур.
Введение
В последние годы благодаря совершенствованию и удешевлению элементной базы тепловизионных приборов существенно возрос спрос на их использование в различных отраслях хозяйственной деятельности человека. Многие производители тепловизоров осуществляют блочную сборку устройств. Комплектующие приобретаются у различных поставщиков или частично изготавливаются самостоятельно. Для согласования отдельных блоков тепловизоров, обоснования требований к основным функциональным узлам, сравнения устройств используются различные методы анализа и синтеза [1-6]. В зависимости от назначения устройства в методы закладываются критерии оценки качества и эффективности. Например, для потребителей систем наблюдения важны дальности обнаружения и распознавания объектов и т.п. Однако во многих случаях требуются менее обобщенные характеристики качества, которые позволяют достаточно полно оценивать свойства тепловизионных приборов. Наиболее важными характеристиками тепловизоров являются пространственная и температурная разрешающая способность, взаимосвязь между которыми устанавливает минимальная разрешаемая разность температур
(Minimum Resolvable Temperature Difference, MRTD). В некоторых учебниках и монографиях [4-6] были получены уравнения для расчета MRTD, в основе которых лежит ряд допущений, предложенных Дж. Ллойдом еще в начале 1970-х гг. [7]. В то же время экспериментальные исследования, проведенные на предприятии «Арсенал» (г. Киев), свидетельствуют о различии полученных результатов теоретических расчетов. Поэтому возникает потребность в дальнейшем усовершенствовании метода расчета MRTD с целью получения более достоверных значений этой функции.
Целью статьи является усовершенствование метода расчета MRTD тепловизоров за счет использования более совершенной модели восприятия наблюдателем термографического изображения.
Уравнение расчета минимальной разрешаемой разности температур
Минимальная разрешаемая разность температур MRTD(vX) определяется минимальной разностью эффективных температур между штрихами тест-объекта заданной пространственной частоты vx и фоном, при которой наблюдатель может визуально различить штрихи изображения объекта на экране дисплея тепловизора. Тест-
объект представляет собой четырехштриховую миру Фуко, в которой полосы и промежутки имеют одинаковую ширину Ур, а их высота равна 7Ур. Это позволяет рассматривать одномерную модель, полагая, что большие вертикальные размеры тест-объекта не вызывают заметных частотных искажений в горизонтальном направлении. Предполагается, что полосы и промежутки излучают как абсолютно черное тело. Пространственная частота может быть выражена как в линейных, так и в угловых величинах (обычно в миллиметрах-1 или миллирадианах-1).
Получим уравнение для расчета ЫКТВ, используя усовершенствованный алгоритм, приведенный в [5, 6]. Одним из дискуссионных допущений Ллойда является механизм восприятия зрительной системой изображения миры Фуко на экране дисплея. Было предложено пространственное интегрирование глаза вдоль координаты сканирования х учитывать согласованной передаточной функцией, которая определяется угловым размером ^"р изображения штриха миры на экране дисплея. При этом модуляционная передаточная функция (МПФ) зрительной системы аппроксимируется функцией [7]:
MEL (V x ) = -
Sin(%^> tp V x )
V
■ = sincfè p V x ).
(1)
Исследованию МПФ зрительной системы посвящено большое количество работ, например [8-10]. На рисунке 1 показаны графики реальной МПФ зрительной системы [10] и ее аппроксимации Дж. Ллойдом (уравнение (1)) и Д.У. Шульцем (уравнение (2)) [11]:
MEs (v x ) = 2,71[exp(-0,11v x ) - exp(-0,35v x )], (2)
где пространственная частота vx имеет размерность град-1.
Анализ приведенных графиков показывает, что, во-первых, аппроксимация Ллойда (1) заметно отличается от реальной МПФ зрительной системы, особенно на низких пространственных частотах; во-вторых, аппроксимация Шульца в значительной степени соответствует реальной МПФ зрительной системы. Поэтому для получения уравнения расчета MRTD будем использовать аппроксимацию Шульца (2).
Одномерный механизм зрительного восприятия изображения миры Фуко будем учитывать через отношение сигнал/шум на выходе
эталонного фильтра, для которого справедлива пропорциональность:
SNR
f
M v^T tp,E
NETD
(3)
где АТ - температурный контраст миры Фуко;
- функция передачи контраста тепловизора; ЫЕТВ - эквивалентная шуму разность температур.
Рисунок 1 - Модуляционная передаточная функция зрительной системы человека: 1 - экспериментально полученная; 2 - аппроксимация (1), при ^ = 0,85 мрад; 3 - аппроксимация (2)
Функция передачи контраста тепловизора с учетом зрительного восприятия связана с МПФ МДух) соотношением [6]:
8
Mtp,E =— Ms .
(4)
Пространственное интегрирование в зрительной системе увеличивает воспринимаемое
отношение сигнал/шум в у / алРл =
= С,у / ал = Сл/7/аоРо раз [5, 6],
где ,- угловые размеры изображения
штриха миры Фуко на экране дисплея вдоль соответствующих осей.
С учетом этих замечаний, воспринимаемое отношение сигнал/шум на экране имеет вид [5, 6]:
SNRf (v x %tD
E %2 NETD tp1
7 ffte k&f a dp d
(5)
где- частота кадров; ¡Е - постоянная времени зрительной системы; Е,{р - угловой размер штриха миры; а0хр0 - угловые размеры пиксела приемника излучения; к^ - коэффициент ширины полосы, показывающий отличие ре-
x
альной полосы пропускания тепловизора с учетом зрительного восприятия от эффективной шумовой полосы эталонного фильтра.
Последний коэффициент вычисляют как:
кА/
1 крб (/)м2е1 (Г)Ы28 (/)М2Е
А/ '
(6)
где ЫРБ/ - спектр мощности шума приемника излучения; Ме1(/), М<(/), МЕ/ - МПФ электронного блока, дисплея и зрительной системы соответственно.
Угловой размер изображения штриха в пространстве изображений в формуле (5) связан с угловым размером штриха ^ в пространстве предметов соотношением:
/о Ре/
Я
(7)
БЕ
где /0 - фокусное расстояние объектива; ре1 -электронное увеличение тепловизора; - расстояние между экраном и зрительной системой.
Угловая пространственная частота у,; в МПФ дисплея и зрительной системы у," связаны с временной частотой / соотношениями:
t "V
VX =—Г и V, =-Х,
аг
Г,
(8)
где 10 - время формирования одного элемента расписания изображения, с; а0 - угловой размер пиксела, мрад; Г, - угловое увеличение тепловизора.
Для практического применения соотношения (5) упростим его согласно рекомендациям Дж. Ллойда [7]: электронный блок и дисплей имеют достаточно широкие эффективные полосы пропускания по сравнению с полосой пропускания зрительной системы, объектива и приемника излучения, т.е. Ме1(/) « 1, М/ « 1; шумовой сигнал в пределах этой полосы является «белым», т.е. УР^(/) « 1.
Тогда коэффициент ширины полосы (6) с учетом (8) при угловом увеличении Г, = 1 будет иметь вид:
/ ■
7,34 <»
"АГ о
[ехр
-1,92
ав у
-ехр
- 5,67^/ ав у
]2 /.
(9)
Вычисление интеграла (9) дает следующий результат:
А/
7,34 а1
-0,138 *
а
*о А/
(10)
где а 0 - безразмерная величина, численно равная угловому размера пикселя в миллирадианах.
Коэффициент ширины полосы, полученный Дж. Ллойдом, определяется как:
к
а,
А/ £
гР*о А/
(11)
На рисунке 2 приведены графики функций (9) и (1 1 ) в зависимости от угловых размеров штриха миры. Анализ графиков показывает, что они совпадают между собой, когда угловой размер штриха равен Е)1р = 0,5 мрад = = 1,72 угл. мин.
Рисунок 2 - Зависимость коэффициента ширины полосы системы «тепловизор - зрительная система» от углового размера £ штриха миры Фуко, когда
а0 = 1,88 мрад и А/ = 0,5: 1 - для аппроксимации Шульца (9); 2 - для аппроксимации Ллойда (11)
Вероятность различения штрихов в изображении миры Фуко зависит от воспринимаемого отношения сигнал/шум 8КЯГ=8КЯЕ. Например, для вероятности распознавания 90 % необходимо, чтобы 8КЯЕ = 4,5. Если подставить это значение в формулу (5), то АТ будет искомой функцией МЯТО(^х) пространственной частоты V,:
2
АТ
■ МЯТБ (V х ) =-БКЯГ х
8
х КЕТВ-
1
М (V X )4
*р
к А/ а БР Б
(12)
7//*Е
*
Подставим (9) в (12) с учетом того, что
пространственная частота vx = 1/2^tp. MRTD будет определяться уравнением:
V ra п
MRTD„(vx) = 0,93-SNRr - NETD x D х
Тогда
Ms (V x )
a Dß I
(13)
- tofftE
В работах [5, 7] было получено подобное уравнение для расчета МКТВ, в основе которого лежит аппроксимация Дж. Ллойда МПФ зрительной системы (11):
MRTDl (vx) = 0,66 - SNRr - NETD
Ms (v x)
1
a Dß D
(14)
¥ - tofftE
Анализ уравнений (13) и (14) показывает, что при условии a D = 0,71 эти уравнения совпадают между собой.
Пример расчета минимальной разрешаемой разницы температур
В качестве примера рассчитаем MRTD тепловизора на базе передающей камеры Thermal Eye TSC, которая имеет следующие характеристики:
- приемник излучения - микроболометрическая матрица размером p х q = 160 х 120 пикселов;
- размер пиксела VD х WD = 30 х 30 мкм2;
- температурная чувствительность NETD = = 0,1 К;
- фокусное расстояние объектива f = = 16 мм при относительном отверстии 1:1;
- угловое поле зрения 2юох х aoy = 17 х 12;
- частота кадров ff = 30 Гц.
Для расчета функции MRTDS(vx) по формуле (13) используем методику, приведенную в учебнике [5]:
1. В исходной формуле (13) известны: SNRr = 4,5 (для вероятности различения штрихов миры 90 % по критериям Джонсона), NETD = 0,1 К; f = 30 Гц; tE = 0,2 с.
2. Угловые размеры пиксела микроболометрической матрицы равны:
a d = ß d = — = Jo
VD 30 -10
-3
16
= 1, 87 мрад.
3. При оптимальном согласовании полосы пропускания А/ и времени ¡0 формирования одного элемента разложения изображения справедливо соотношение: А/ ¡0 = 0,5.
4. МПФ тепловизора М^(уХ) с высококачественным объективом определяется размером пиксела а0 приемника излучения, т.е.:
Уп 30•10"3 , оп
а о = Р о = =- = 1, 87 мрад,
fo
16
ms (Vx ) ~ md (Vx ) = sin c(vxad ) =
sin(1,87rcv x) 1,87tcv v
После подстановки исходных параметров в уравнение (13) получим функцию для расчета МКТВ:
MRTDS (v x) = 0,93 - 4,5 - 0,1 - 1,87v x
1,87rcv
1,87rcv,
sin(1,87rcv )
1
1,87
2
0,5 - 30 - 0,2
= 0,845v,
sin(1,87rcv )
На рисунке 3 приведены графики MRTD тепловизора с передающей камерой Thermal Eye TSC, построенные по формулам (13) и (14).
Рисунок 3 - Минимальная разрешаемая разность температур тепловизора с камерой Thermal Eye TSC для аппроксимаций зрительной системы Шульцем (1) и Ллойдом (2)
Анализ этих функций показывает, что: - значения функций MRTDs(vX) и MRTDL(vx) отличаются между собой в 2,6 раза;
*
х
v
x
х
х
x
- функция МЯТО^х) зависит от параметра а О. При аО = 0,71 функции МЯТО^х) и МЯТО^,) совпадают между собой.
Заключение
Получено усовершенствованное уравнение для расчета минимальной разрешаемой разности температур тепловизора, в основе которого лежит более достоверная аппроксимация модуляционной передаточной функции зрительной системы, предложенная
Д.У. Шульцем.
Исследование уравнений для расчета минимальной разрешаемой разности температур, в которых использованы модуляционные передаточные функции зрительной системы, предложенные Дж. Ллойдом и Д.У. Шульцем, показало, что для рассматриваемой модели теп-ловизионной системы:
- расчетные значения минимальной разрешаемой разности температур по этим уравнениями существенно отличаются между собой;
- коэффициент ширины полосы, показывающий различие между реальной полосой пропускания тепловизора с учетом зрительного восприятия и эффективной шумовой полосой эталонного фильтра, полученный Дж. Ллойдом, зависит от углового размера штриха миры Фуко, а для аппроксимации Шульца - не зависит. В этом заключается принципиальное отличие полученного уравнения для расчета минимальной разрешаемой разности температур от известных уравнений.
Дальнейшие исследования следует направить на усовершенствование уравнения расчета минимальной разрешаемой разности температур с учетом влияния характеристик дисплея и условий наблюдения.
Список использованных источников
1. Vollmerhausen, R.H. Analysis and evaluation of sampled imaging systems / R. H. Vollmerhausen, D. Reago, R. G. Driggers. - Bellingham, SPIE Press, 2010. - 304 p.
2. Wishart, G.C. The accuracy of digital infrared imaging for breast cancer detection in women undergoing breast biopsy / G.C. Wishart, M. Campisi, M. Boswell [et al.] // European Journal of Surgical Oncology. - 2010. - Vol. 36 (6). - P. 535-540.
3. Leachtenauer, J.C. Surveillance and Reconnaissance Systems: Modeling and Performance Prediction / J. C. Leachtenauer, R. G. Driggers. -Artech House, Incorporated, 2001. - 419 p.
4. Тарасов, М.М. Инфракрасные системы «смотрящего» типа / М.М. Тарасов, Ю.Г. Якушен-ков. - М. : Логос, 2004. - 444 с.
5. Колобродов, В.Г. Проектирование тепловизионных и телевизионных систем наблюдения : учебник / В.Г. Колобродов, Н.И. Лихолит. -Киев : НТУУ «КПИ», 2007. - 364 с.
6. Мирошников, М.М. Теоретические основы оп-
тико-электронных приборов : учеб. пособие для приборостроит. вузов / М.М. Мирошников. - 2-е изд., перераб. и доп. - Л. : Машиностроение, 1983. - 696 с.
7. Ллойд, Дж. Системы тепловидения / Дж. Ллойд;
пер. с англ. - М. : Мир, 1978. - 416 с.
8. Blackwell, H.R. Contrast thresholds of the human
eyes / H.R. Blackwell // Journal of the Optical Society of America. - 1946. - Vol. 36, №11. -P. 624-643.
9. Campbell, F.W. Optical Quality of the Human Eye / F.W. Campbell, R.W. Gubisch // Journal of Physiology. - 1966. - Vol. 186(3). - P. 558-578.
10. Barten, P.G.J. Contrast sensitivity of the human eye and its effects on image quality / P.G.J. Barten; Bellingham: SPIE Press, 1999. - 211p.
11. Карасик, В.Е. Лазерные системы видения / В.Е. Карасик, В.М. Орлов. - М. : Издательство МГТУ им. Н. Э. Баумана, 2001. - 174 с.
ESTIMATION OF THERMAL IMAGER MINIMUM RESOLVABLE TEMPERATURE
DIFFERENCE
AhmedM.L. Al-Mzirawi, Kolobrodov V.G., Mykytenko V.I.
National Technical University of Ukraine «Kyiv Politechnic Institute», Kyiv, Ukraine
e-mail: [email protected]
Abstract. The paper is devoted to improving the accuracy of determining the minimum resolvable temperature difference (MRTD) of thermal imagers by means of more sophisticated models of thermal images visual perception. An MRTD calculating algorithm, which is based on a more reliable approximation of the visual system modulation transfer function by Schultz, is considered. New expression for the bandwidth of Schultz's approximation, which does not depend on the angular size of the Foucault pattern bar, was obtained.
Keywords: thermovision, performance evaluation, minimum resolvable temperature difference.
References
1. Vollmerhausen R.H., Reago D., Driggers R.G. Analysis and evaluation of sampled imaging systems. Belling-ham, SPIE Press, 2010, 304 p.
2. Wishart G.C, Campisi M., Boswell M., et al. The accuracy of digital infrared imaging for breast cancer detection in women undergoing breast biopsy. European Journal of Surgical Oncology, 2010, vol. 36 (6), pp. 535-540.
3. Leachtenauer J.C., Driggers R.G. Surveillance and Reconnaissance Systems: Modeling and Performance Prediction. Artech House, Incorporated, 2001, 419 p.
4. Tarasov M.M., Yakushenkov Yu.G. Infrakrasnye sistemy smotryashchego tipa [Forward looking Infrared systems]. Moscow, Logos Publ., 2004, 444 p. (in Russian).
5. Kolobrodov V.G., Liholit N.I. Proektirovaniye teplovisionnykh i televisionnykh system nablyudeniya [Development of thermovision and television systems]. Kiev, NTUU "KPI" Publ., 2007, 364 p. (in Russian).
6. Miroshnikov M.M. Teoreticheskiye osnovy optiko-elektronnykh priborov [Theoretical foundation for optoelectronic devices]. Leningrad, Mashinostroenie Publ., 1983. 696 p. (in Russian).
7. Lloyd J. Teplovisionnye sistemy [Thermal Imaging Systems]. Moscow, Mir Publ., 1978, 416 p. (in Russian).
8. Blackwell H.R. Contrast thresholds of the human eyes. Journal of the Optical Society of America, 1946, vol. 36, № 11, pp. 624-643.
9. Campbell F.W., Gubisch R.W. Optical Quality of the Human Eye. Journal of Physiology, 1966, vol. 186(3), pp. 558-578.
10. Barten P.G.J. Contrast sensitivity of the human eye and its effects on image quality. Bellingham: SPIE Press, 1999. - 211p.
11. Karasik V.E., Orlov V.M. Lazernye sistemy videniya [Laser imaging systems]. Moscow, MGTU im.Baumana N., 2001, 174 p. (in Russian).
Поступила в редакцию 11.12.2014.
