Повышение точности измерения температуры по тепловизионным изображениям

ПОВЫШЕНИЕ ТОЧНОСТИ ИЗМЕРЕНИЯ ТЕМПЕРАТУРЫ ПО ТЕПЛОВИЗИОННЫМ ИЗОБРАЖЕНИЯМ

Попов С.А., Доцин И.И., Емельянов Г.М.

Новгородский государственный университет им. Ярослава Мудрого

Аннотация

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

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

Введение

Для измерения температуры дистанционным методом используется микроболометрический модуль 851Я, который регистрирует тепловизионные изображения. В этом модуле используется матрица фотодетекторов типа иЬ03-162. Расчёт температуры базируется на теории теплового излучения абсолютно чёрного тела (АЧТ) [1]. На этапе градуировки матрицы фотодетекторов составляется градуиро-вочная таблица, связывающая температуру и амплитуды пикселей. На этапе измерения по заданным амплитудам пикселей по градуировочной таблице отыскивается соответствующая эффективная температура АЧТ. Недостатком такого метода является использование интерполяции градуировочной таблицы, что вносит определённые погрешности измерения температуры. Повысить точность оценивания температуры с помощью тепловизора и получить статистические характеристики этих оценок в этом случае можно путём использования многоотклико-вой модели зависимости амплитуды пикселей от температуры. Возможность использования методов регрессионного анализа в задачах обработки изображений показана в работе [2].

1. Построение многооткликовых моделей для изображений

Тепловизионное изображение является прямоугольной матрицей амплитуд пикселей И(/'ь ,2), где , = 1, т1 - номер строки изображения, ,2 = 1, т2 -

номер столбца изображения. Для целей моделирования каждого пикселя это изображение представляется вектором пикселей, в котором номер пикселя выражается в виде , = ,2 + т2 (,1 -1), а т = т1т2 -

количество пикселей в изображении [3].

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

и = Г

где и = {и1,и2,...,ит}т - вектор амплитуд для каждого пикселя при заданной температуре г, В = {Ь1,Ь2,...,Ьк}т - вектор коэффициентов модели,

Г(В, г) = {^ (В, г), ^ (В, г),... ,4 (В, г )}т - вектор функций зависимости амплитуды от эффективной температуры для каждого пикселя,

Е = {е1, е2,..., ет }т - вектор ошибок наблюдений амплитуды.

Вектор оценок коэффициентов в общем виде рассчитывается с помощью следующей итерационной процедуры [4]:

В"1 = В' +

V-1 XР

¡=1

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

Таким образом, входными данными для метода построения модели измерения эффективной температуры служат значения задаваемых температур г,,j = 1, п и соответствующий вектор измеренных для каждого

полученный

пикселя амплитуд и, = {и1,, и2,,..., ит1 }т

' ] 1 1Р 2 ¡'"''"щ* для каждой температуры.

Матрица Р(В, г) определяется как

Р (B, г ) =

ЭГ

ЭВ

Э ^

Э ^

ЭВ

ЭВ

ЭВ

В качестве оценки ковариационной матрицы УЕ использовалась матрица 8Е, которая рассчитывается по формуле [3]:

ГУ Я( и, - и)(и, - и )т,

1 ,=1

^Е = "

где а - количество точек наблюдения амплитуды, для которых температура практически не изменяется.

Ковариационная матрица оценок коэффициентов рассчитывается по формуле

Ув =

X Р (В', г,) Уе1Р ( В ', tj )т

¡=1

где В - вектор оценок коэффициентов, полученных по формуле (2).

2. Анализ ковариационной матрицы ошибок наблюдений и проверка адекватности модели 2.1. Методика анализа ковариационной матрицы Исследование ковариационной матрицы ошибок наблюдений 8е состоит в оценивании элементов этой матрицы и в проверке гипотезы о некоррелированности некоторых ошибок наблюдений. Проверяется гипотеза о значимости элементов этой матрицы (// = 1,т; / ф / ) [4]. Для этого вычисляется коэффициент корреляции в виде

V.

Е//

4

Если наблюдения некоррелированы, то величина 1 = (п - 2) / (1 - г/2) имеет г -распределение, нечувствительное при п>10 к числу степеней свободы. Незначимые элементы матрицы 8е принимались равными нулю.

2.2. Проверка на нулевое среднее Для проверки гипотезы о равенстве нулю математического ожидания остатков использовалась ^-статистика Хотеллинга [5] в виде

1

Т2 = -

«Т ^е'к/

пт - к/1

где Я = и - Р(Б, г) - вектор остатков модели, а матрица 8е рассчитывалась по формуле (3).

Если остатки имеют нулевое математическое ожидание, то статистика (6) приблизительно имеет распределение Фишера со степенями свободы / = п - к / т и /2 = а -1.

Тогда при выполнении неравенства Т2 < F

где Ff, / 1-а - квантиль распределения Фишера со степенями свободы /1 и /2 на уровне значимости а, гипотеза о равенстве нулю математического ожидания остатков принимается с доверительной вероятностью 1 - а .

3. Расчёт оценки температуры по вектору амплитуд В качестве метода расчёта температуры использовался метод Ньютона [4], который основан на линеаризации целевой функции ф(Б,и, г) = и -f (Б,г) в *

окрестности точки г в виде:

. Эф( Б, и, г*) ф(Б, и, г —-(г - г ) = 0,

Эг

Эф (Б, и, г) ^

где —1- - производная функции ф(Б,и, г) по

Эг

температуре.

Тогда итерационная процедура расчёта оценки температуры принимает следующий вид:

г*+1 = гА -ф

Эф

Эг

Выражение (8) было обобщено на случай вектора оценок температуры следующим образом:

Т^+1_

Эф

ЭТ

ф(б, т , и),

где Т = {г,,г2,...,гт}т - вектор оценок температуры по каждому пикселю.

Векторная функция Ц (Б, Т, и) включает формулы зависимости амплитуды от температуры для каждого пикселя и имеет следующий вид:

ф( Б, Т, и) = {ф1 (Б,и,г), ф2 (Б, и,г).....фт (Б,и, г )}т.

Ковариационная матрица УТ оценок температур по всем пикселям равна

УТ = ЛУБЛт + О8ЕОт,

где Л = Л1 - Л2, О = О1 - О2,

Л1 =

Эф

ЭБ

Л2 = ф

О =

Эи

О2 = ф( Б, Т, и )т

Эф

ЭТ

Эф

ЭТ

Эф

ЭТ

Эф

ЭТ

Э2 ф( Б, Т, и) ЭТЭБ

Э2 ф( Б, Т,и) ЭТЭи

Общая оценка температуры г минимизирует квадратичную форму (1г - Т)т У-1 (1г - Т), где J -вектор размерности т, элементы которого равны

единице, Т - вектор частных оценок температуры по каждому пикселю, которые рассчитываются по формуле (9). Приравнивая производную от этого выражения по г нулю, получим оценку температуры / в виде:

г = ( Jт уТ^ )-1 Jт уТ-1Т .

Дисперсия а2 оценки г (11) равна

г2 = ( J т У^ )-1,

что меньше наименьшей из дисперсий частных оценок.

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

4. Анализ результатов эксперимента Экспериментальные данные получают путём регистрации излучения от АЧТ в диапазоне температур,

максимально приближенных к рабочему (0 - 150°С). Последовательно проводились п = 26 наблюдений с различными значениями установленной температуры АЧТ г,, в процессе которых регистрировались и сохранялись тепловизионные изображения АЧТ. Каждое изображение представляется матрицей амплитуд пикселей, общее количество пикселей равно 384x288 = 110592 пикселей.

Точки для расчёта ковариационной матрицы 8Е (3) значения а < п выбирались при низких температурах, где зависимость f¡ (В, г) (, = 1, т ) на начальном участке хорошо описывается линейной функцией (рис. 1). Таким образом, матрица 8Е рассчитывалась как ковариационная матрица остатков линейной по г модели, построенной по участку низкотемпературных данных (0 - 20)°С. Структура этой матрицы имеет следующий вид: дисперсии ошибок наблюдений откликов пикселей равны 78 - 86 (в относительных единицах), а коэффициенты корреляции (5) равны 0,95 - 0,98.

Амплитуда

3000

0 20 40 60 80 100 120

Температура, °С

Рис. 1. Зависимость амплитуды пикселя от эффективной температуры для полиномиальной модели

В [1] показано, что зависимость амплитуды и пикселя, формируемого ИК камерой, от температуры АЧТ описывается формулой:

и =-^-,

ехр( ^ + 273) -1

где г - температура чёрного тела (°С), В = {Ь1, Ь2}т -вектор коэффициентов, зависящих от используемого объектива, его относительного отверстия и фильтров камеры. Кроме этой модели, исследовались полиномиальные модели: линейная, квадратичная, полиномы третьей и четвёртой степени. Расчёты оценок коэффициентов модели (1) выполнялись по формуле (2).

Для каждой модели рассчитывалась величина путём решения уравнения

Т 2 = РМгР,

где ^р - квантиль распределения Фишера для вероятности Р.

Для модели АЧТ (13) получены расчётные значения T 2 = 11,16 (6), что для числа степеней свободы f = 22 и f2 = 5 соответствует P = 0,993 (14); такая модель принимается неадекватной на уровне значимости более 0,007.

Для модели в виде полинома четвёртой степени значение статистики (6) T 2 = 2,37, что соответствует P = 0,828; это недостаточно для признания этой статистики значимой на уровне значимости менее

0.172. Следовательно, дисперсии SE (3) и VR (7) считаются одинаковыми и гипотеза о равенстве нулю математического ожидания остатков принимается, т. е. модель следует признать адекватной. Типичное положение экспериментальных точек для одного пикселя и построенные модели показаны на рис. 1.

В результате использования формул (11,12) относительная погрешность расчёта температуры st /1 составила около 0,015%.

Высокие значения ковариаций в матрице SE приводят к существенному повышению дисперсии оценки t по сравнению с диагональной матрицей. Для диагональной матрицы SE относительная погрешность составляет 0,005%.

Для сравнения были проведены расчёты оценки температуры с использованием кусочно-линейной интерполяции для каждого отдельного пикселя с последующим усреднением частных оценок. Относительная погрешность расчёта температуры в этом случае составила приблизительно 0,045%.

Выводы

Использование многооткликовых регрессионных моделей для описания зависимости амплитуды пикселей тепловизионного изображения от температуры позволяет учесть корреляционные связи ошибок наблюдений откликов пикселей, что повышает объективность и точность метода расчёта температуры.

Литература

1. Госсорг, Ж. Инфракрасная термография. Основы, техника, применение / Ж. Госсорг. - М.: Мир, 1988.— 416 с.

2. Сергеев, В.В. Сравнительный анализ методов аппроксимации функций в задачах обработки изображений /

B.В. Сергеев, В.Н. Копенков, А.В. Чернов // Компьютерная оптика — 2004. — № 26. — С. 119-122.

3. Попов, С.А. Статистический метод калибровки цифровых фотокамер для улучшения цветовоспроизведения / С.А. Попов // Автометрия. — 2003. — Т. 39, № 2. —

C. 92-99.— ISSN 0320-7102.

4. Мейндоналд, Дж. Вычислительные алгоритмы в прикладной статистике / Дж. Мейндоналд. — М.: Финансы и статистика, 1988. — 350 с.

5. Rencher, A.C. Methods of Multivariate Analysis / Alvin C. Rencher. — A John Wiley & Sons, Inc. Publication. Brigham Young University, 2002. — 708 p.

References

1. Gaussorgues, G. La Thermographie Infrarouge.Principes-Technologie-Applications. Technique et Documentation / Gilbert Gaussorgues. — Paris: Lavoisier, 1984. — 481 p.

2. Sergeev, V.V. Comparative analysis of function approximation methods in problems of image processing / V.V. Sergeev, V.N. Kopenkov, A.V. Chernov // Computer optics - 2004. - N 26. - P. 119-122. - (In Russian).

3. Popov, S.A. Statistical method of calibration of digital photocameras for improvement of color reproduction / S.A. Popov // Autometria. - 2003. - V. 39, N 2. - P. 92-99.- (In Russian).

4. Maindonald, J. Computational Algorithms in Applied Statistics / J. Maindonald. - Moscow: "Finansy i Sta-tistika" Publisher, 1988. - 352 p. - (In Russian).

5. Rencher, A.C. Methods of Multivariate Analysis / Alvin C. Rencher. - A John Wiley & Sons, Inc. Publication. Brigham Young University, 2002. - 708 p.

ACCURACY INCREASE OF TEMPERATURE MEASUREMENT ON THERMAL IMAGES

S.A. Popov, I.I. Docin, G.M. Emelyanov Yaroslav-the-Wise Novgorod State University

Abstract

Multiresponse regression model is used for dependence of thermal image pixel amplitudes from temperature to increase accuracy of remote indirect method of temperature measurement. Obtained model allows to calculate integral estimation of temperature and its dispersion on all pixels taking into account correlation relationships of amplitude observation errors.

Key words: remote temperature measurement, thermal image, multiresponse model, covariance matrix of observation errors, temperature estimate.

Сведения об авторах

Попов Станислав Алексеевич, зав. кафедрой начертательной геометрии и компьютерной графики Новгородского государственного университета им. Ярослава Мудрого, д.т.н., профессор. Основные научные интересы - построение многооткликовых регрессионных моделей процессов и систем, векторная и обратная оптимизация и разработка методов классификации изделий и технологических процессов.

Stanislav Alekseevich Popov, holder of chair of descriptive geometry and computer graphics at Yaroslav-the-Wise Novgorod State University, doctor of engineering sciences, professor. The main scientific interests are construction of multiresponse regression models of objects and processes, development of vector and reverse optimization procedures and methods of classification of products and technological processes.

Доцин Илья Игоревич, аспирант Новгородского государственного университета им. Ярослава Мудрого.

Ilya Igorevich Docin, postgraduate student of Yaroslav-the-Wise Novgorod State University.

Емельянов Геннадий Мартинович, доктор технических наук, профессор. В настоящее время - профессор кафедры информационных технологий и систем Новгородского государственного университета им. Ярослава Мудрого. Научные интересы включают построение проблемно-ориентированных вычислительных систем обработки и анализа изображений.

Gennady Martinovich Yemelyanov, doctor of engineering sciences, professor. At present is a professor of the department of information technology and systems of Yaroslav-the-Wise Novgorod State University. Scientific interests include construction of problem-oriented computer systems of image processing and analysis.

Поступила в редакцию 6 ноября 2012 г.

Дизайн: Я.Е. Тахтаров. Оформление и верстка: М.А. Вахе, С.В. Смагин и Я.Е. Тахтаров.

Подписано в печать 15.03.2013 г. Усл. печ. л. 15,48. Отпечатано в типографии ООО «Предприятие «Новая техника». Заказ № 11/1. Тираж 322 экз. Печать офсетная. Формат 62х84 1/8.