способ коррекции неоднородности многоэлементных фотоприемных устройств со сканированием

Формула изобретения

1. Способ коррекции неоднородности многоэлементных фотоприемных устройств со сканированием, заключающийся в том, что сканируют сцену и регистрируют отсчеты сигналов элементов фотоприемного устройства, корректируют отсчеты сигналов сцены, отличающийся тем, что для проведения коррекции сканируют в последовательности, обеспечивающей попадание потоков одних и тех же элементов сцены на соседние элементы фотоприемного устройства и попадание потоков соседних элементов сцены на один и тот же элемент устройства, для каждой пары соседних элементов фотоприемного устройства, i-го и (i+1)-го, по множеству пар отсчетов сигналов S _i(t) и S_i+1(t), i N, N - множество элементов фотоприемного устройства и t - параметр времени, и каждая пара отсчетов соответствует одному и тому же состоянию элемента сцены, с помощью аппроксимации определяют корректирующие воздействия - функции связи сигналов соседних элементов так, чтобы эти функции имели наименьшее среднеквадратичное отклонение на множествах пар отсчетов S_i (t) и S_i+1(t) и отображали любое значение из диапазонов переменных s_i, соответствующих диапазонам сигналов S_i(t), в значения переменных s_i+1, соответствующие значениям из диапазонов сигналов S_i+1(t), корректируют отсчеты сигналов S_i(t) их заменой на где S_i(t) - отсчеты сигналов i-го элемента, скорректированные сигналы i-го элемента относительно j-го элемента, i,j N,

где F_j,i(S _i(t)) - функция коррекции сигнала S_i (t) относительно j-го элемента, которую получают последовательным применением аргументов функций связей F_i+1,i (), F_i+2,i+1(), ..., F_j,j-1 () сигналов пар соседних элементов i+1,i; i+2,i+1; ...j,j-1, составляющих путь от элемента i к элементу j, для i=j =S_i(t), где S_i (t) - аргумент первой функции, а аргументом каждой последующей функции является результат предыдущей функции, от F _i+1,i() до F_i,j-1().

2. Способ по п.1, отличающийся тем, что функции F_i+1,i () являются линейными зависимостями

имеющими наименьшее среднеквадратичное отклонение

от измеренных зависимостей S_i+1 (t) от S_i(t), корректирующие воздействия: коэффициенты линейных зависимостей A_i+1,i - отношение чувствительностей соседних элементов и В _i+1,i - относительное смещение соседних элементов определяются по формулам

где Т_i+1,i - множество пар отсчетов S_i(t) и S_i+1 (t) по одним и тем же элементам сцены, упорядоченное по времени: значению t=1 соответствует первый по времени отсчет, t=2 - следующий и так далее до Т_i+1,i-го отсчета; S _i и S_i+1 - средние значения сигналов i-го и (i+1)-го элементов на интервале Т_i+1,i , скорректированные относительно элемента j сигналы элементов i=1, 2,..., N определяются по исходным некорректированным сигналам S, (t) и коэффициентам A_i+l,i+l-1 и В_i+l,i+l-1, 1=1, 2,...,j-i, следующим образом:

Описание изобретения к патенту

Изобретение относится к оптико-электронным системам формирования и обработки инфракрасных изображений, для которых актуальна задача устранения неоднородности сигналов, обусловленной различиями в чувствительности к входному потоку и в темновом токе элементов фотоприемных устройств (ФПУ). Изобретение может использоваться в тепловизионных системах со сканированием при использовании многорядных ФПУ.

Известен способ [1, 2] коррекции неоднородности многоэлементных МФПУ в сканирующих тепловизорах (его рассматриваем в качестве аналога), заключающийся в том, что сканируют сцену, регистрируют отсчеты сигналов элементов ФПУ от наблюдаемой сцены, по отсчетам сигналов усреднением определяют строчные корректирующие коэффициенты, корректирование осуществляют умножением значений отсчетов сигналов на соответствующие элементам строчные корректирующие коэффициенты.

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

Наиболее близким по технической сущности к заявляемому изобретению является способ [3] двухточечной коррекции, который обеспечивает компенсацию по чувствительности и по смещению независимо от сцены, так как использует для определения корректирующих коэффициентов сигналы двух эталонных источников - высокой и низкой интенсивности излучения. Способ заключается в том, что сканируют сцену и регистрируют отсчеты сигналов элементов фотоприемного устройства, сканируют эталонные источники и регистрируют отсчеты их сигналов, по отсчетам сигналов эталонных источников определяют корректирующие воздействия - коэффициенты для коррекции по чувствительности и смещению элементов фотоприемного устройства, корректируют отсчеты сигналов сцены - вычитают из значений отсчетов корректирующие коэффициенты по смещению и умножают разности на корректирующие коэффициенты по чувствительности.

Сущность способа состоит в следующем. При воздействии эталонного оптического сигнала низкого уровня, одинакового для всех элементов ФПУ, средние значения электрических сигналов низкого уровня для всех элементов МФПУ равны

где - среднее значение сигнала низкого уровня от i-го элемента,

i=1,2,...,N, N - число элементов МФПУ;

Т - интервал усреднения, t - номера отсчетов сигнала;

- сигнал i-го элемента на t-м отсчете;

А _i и B_i - чувствительность и смещение i-го элемента;

Р^н - поток эталонного излучения низкого уровня;

Е_i - усредненный на интервале Т шум элемента, Т выбирают достаточно большим, чтобы можно было пренебречь влиянием этой составляющей.

При воздействии эталонного оптического сигнала высокого уровня, одинакового для всех элементов МФПУ, средние значения электрических сигналов высокого уровня равны

где Р^в - поток эталонного излучения высокого уровня.

При воздействии произвольных сигналов сцены текущие сигналы отдельных элементов S _i(t) равны

Корректирование сигналов сцены осуществляется заменой исходных сигналов сцены S_i(t) на скорректированные сигналы сцены , определяемые следующим образом:

где Р^вн=(Р^в-Р ^н) - разность эталонных потоков высокого и низкого уровней, i=1,2,...,N.

Корректирующими неоднородность коэффициентами являются - коэффициенты по смещению, и - по чувствительности. Применение этих коэффициентов к сигналам сцены - вычитание из сигналов коэффициентов по смещению и умножение разностей на коэффициенты по чувствительности - обеспечивает однородность скорректированных сигналов: одну и ту же по всем элементам чувствительность - величина 1/ P^вн, и одно и то же смещение, равное (-P^н/ P^вн).

Основным недостатком является то, что использование источников эталонных сигналов сопряжено с необходимостью их подстройки под сцену. Способ обеспечивает коррекцию при условии линейной зависимости сигналов элементов от оптического потока на их входе. Для того чтобы в пределах динамического диапазона наблюдаемой сцены отклонения от линейной зависимости были минимальными, диапазон значений эталонных потоков согласовывается с диапазоном значений сигналов сцены. Кроме того, требуется поддерживать однородность засветки элементов МФПУ эталонными источниками, так как от равномерности засветки зависит точность определения корректирующих коэффициентов. В некоторых случаях наличие эталонных источников неприемлемо по условиям размещения их в аппаратуре.

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

Цель достигается благодаря тому, что сканируют сцену и регистрируют отсчеты сигналов элементов фотоприемного устройства, определяют по отсчетам наблюдаемой сцены корректирующие воздействия, корректируют отсчеты сигналов сцены. Сцену сканируют в последовательности, обеспечивающей попадание потоков одних и тех же элементов сцены на соседние элементы фотоприемного устройства и попадание потоков соседних элементов сцены на один и тот же элемент устройства, для каждой пары соседних элементов фотоприемного устройства, i-го и i+1-го, по множеству пар отсчетов сигналов S_i(t) и S_i+1(t), i N, N - множество элементов фотоприемного устройства и t - параметр времени, где каждая пара отсчетов соответствует одному и тому же состоянию элемента сцены, с помощью аппроксимации определяют корректирующие воздействия - функции связи сигналов соседних элементов

так, чтобы эти функции имели наименьшее среднеквадратичное отклонение на множествах пар отсчетов S_i (t) и S_i+1(t) и отображали любое значение из диапазонов переменных s_i, соответствующих диапазону сигнала S_i(t), в значения переменных s_i+1, соответствующих значениям из диапазонов сигналов S_i+1(t), корректируют отсчеты сигналов сцены S_i(t) их заменой на , где S_i(t) - отсчеты сигналов i-го элемента, - скорректированные сигналы i-го элемента относительно j-го элемента, i, j N,

где F_j,i(S((t)) - корректирующее воздействие - функция коррекции сигнала S_i (t) относительно j-го элемента, которую получают последовательностью применения - начальный аргумент S_i(t), функций связей F_i+1,i(), F _i+2,i+1(),...,F_j,j-1() сигналов пар соседних элементов i+1,i; i+2,i+1;...j,j-1, составляющих путь от элемента i к элементу j, для i=j .

Сущность способа заключается в следующем.

Условие сканирования сцены - попадание потоков одних и тех же элементов сцены на соседние элементы фотоприемного устройства - позволяет из последовательностей отсчетов S _i(t) и S_i+1(t) выбрать пары отсчетов от одних и тех же элементов сцены при одинаковом их состоянии (одинаковых потоках) и по множеству таких пар определить с помощью аппроксимации корректирующие воздействия - функции связи соседних элементов F_i+1,_i(),

которые отображают значения отсчетов i-го элемента на значения отсчетов элемента i+1. Отображение означает, что значению s_i отсчета сигнала S_i (t), которому соответствует некоторое значение потока Р на его входе, ставится в соответствие значение s_i+1 отсчета сигнала S_i+1(t), которому соответствует то же самое значение потока Р на его входе.

Условие сканирования - попадание потоков соседних элементов сцены на один и тот же элемент фотоприемного устройства - обеспечивает полноту множества функций связи элементов F_i+1, _i() такую, что для всех путей, связывающих с любым элементом фотоприемного устройства все остальные элементы устройства, будут определены функции связи всех пар соседних элементов этих путей.

Выполнение процедуры корректирования элементов фотоприемного устройства относительно j-го элемента состоит в определении значений функций коррекции F_j,i(S _i(t)) сигналов S_i(t), i,j N, определение каждой осуществляется последовательностью применений функций связи (последовательностью корректирующих воздействий) на пути от i-го к j-му элементу.

Первое в последовательности применений функций связи - применение функции F_i+1,i(S_i(t)) - дает значение - скорректированное значение сигнала S _i(t) относительно элемента i+1,

Так как функция связи определяется на последовательности элементов сцены, общей для i-го и i+1-го элементов МФПУ, то подстановка в функцию аргумента S_i(t), которому соответствует некоторое значение потока элемента сцены, в качестве результата функции дает значение S_i+1(t), которому соответствует то же значение потока.

Аналогично, для последующих применений, 1=1,2,...,j-1,

которые определяют - скорректированные относительно элементов i+1 значения сигнала S_i(t), и которые, следуя определениям функций связи, равны S_i(t). Начальному аргументу применения функций - S_i(t), которому соответствует некоторое значение потока элемента сцены, функции связи ставят в соответствие значения S_i+l (t), которым соответствует то же самое значение потока, что и сигналу S_i(t).

Это значит, что скорректированные по отношению к одному из элементов фотоприемного устройства сигналы остальных элементов будут одинаковы при одинаковых сигналах сцены по всем элементам МФПУ, что соответствует цели коррекции.

В случае небольших изменений сигналов в пределах сцены зависимости S_i+1(t) от S_i(t) близки к линейным и аппроксимируются прямыми

имеющими наименьшее среднеквадратичное отклонение

от измеренных зависимостей S_i+1 (t) от S_i(t).

Корректирующие воздействия: коэффициенты линейных зависимостей A_i+1,i - отношение чувствительностей соседних элементов, и B _i+1,i - относительное смещение соседних элементов, определяются по формулам

где T_i+1,i, - множество пар отсчетов S_i(t) и S_i+1 (t) по одним и тем же элементам сцены, упорядоченное по времени: значению t=1 соответствует первый по времени отсчет, t=2 - следующий и так далее до T_i+1,i-го отсчета;

S_i и S_i+1 - средние значения сигналов i-го и i+1-го элементов на интервале Т _i+1,i. При увеличении интервала отношение A _i+1,i стремится к A_i+1/A _i.

Формулы для определения отношений чувствительностей и относительных смещений пар соседних элементов получены из линейной регрессии отсчетов сигналов S_i(t) по отсчетам S_i+1(t) и линейной регрессии отсчетов S _i+1(t) по S_i(t). Для устранения смещения величин A_i+1,i по отношению к А _i+1/А_i из-за накопления шума в формулах используется перемножение соседних, некоррелированных по шуму, отсчетов.

Скорректированные относительно элемента j сигналы элементов i=1,2,...,N определяются по исходным, некорректированным сигналам S_i(t) и коэффициентам A _i+l,i+l-1 и В_i+l,i+l-1, 1=1,2,...,j-i, следующим образом:

На фиг.1 приведена, в сопоставлении с реализацией способа-прототипа [3], схема варианта устройства, реализующего предлагаемый способ,

где

1 - эталон сигнала низкого уровня;

2 - эталон сигнала высокого уровня;

3 - сканер;

4 - эталонный сигнал низкого уровня;

5 - сигнал сцены;

6 - эталонный сигнал высокого уровня;

7 - фотоприемное устройство;

8 - аналого-цифровой преобразователь сигналов;

9 - определитель корректирующих коэффициентов;

10 - корректор сигналов сцены.

Пунктиром выделены устройства и сигналы, необходимые при использовании способа-прототипа и не нужные для заявляемого способа. В предлагаемом способе корректирующие коэффициенты A_i+1,i и В _i+1,i определяются по сигналам сцены определителем корректирующих коэффициентов. В прототипе этот блок определяет корректирующие коэффициенты по чувствительности и смещению с использованием эталонных источников. В предлагаемом способе корректор сигналов сцены определяет скорректированные относительно элемента j сигналы элементов i=1,2,...,N по исходным, некорректированным сигналам S_i(t) и коэффициентам A _i+l,i+l-1 и B_i+l,i+l-1, l=1,2,...,j-i и заменяет S_i(t) на . В прототипе в корректоре вычитают из значений отсчетов сцены корректирующие коэффициенты по смещению и умножают разности на корректирующий коэффициент по чувствительности и этим скорректированным значением заменяют исходные значения отсчетов сцены.

Условие сканирования сцены в последовательности, обеспечивающей попадание потоков одних и тех же элементов сцены на соседние элементы фотоприемного устройства и попадание потоков соседних элементов сцены на один и тот же элемент устройства, поясним с помощью фиг.2, где представлен вариант расположения фоточувствительных элементов для многорядной линейки формата 4×N. Изображение сканируется в направлении, перпендикулярном линейкам фоточувствительных элементов A, B,...F, каждый пиксель изображения сцены регистрируется последовательно элементами A_i, В_i , С_i, D_i, а соседний пиксель - элементами Е_i, F _i, G_i, Н_i. При этом сигналы элементов А_i, В _i, С_i, D_i могут быть скорректированы по отношению к сигналам любого элемента этой группы, а сигналы элементов Е_i, F _i, G_i, H_i могут быть скорректированы по отношению также любого элемента из этой группы. При втором сканировании изображение может быть сдвинуто на половину шага элементов в линейках. При этом последовательность пикселей, регистрируемых в первом сканировании элементами A _i, B_i, С_i, D_i, будет регистрироваться элементами Е _i, F_i, G_i, H_i. Таким образом, можно скорректировать сигналы элементов Е_i, F _i, G_i, H_i по отношению к элементам А_i, В _i, С_i, D_i, а также сигналы элементов А_i+1, B _i+1, C_i+1, D_i+1 по отношению к элементам E_i, F _i, G_i, Н_i и таким образом все элементы МФПУ связываются отношениями коррекции.

Источники информации

1. Патент РФ 2113065. Способ выравнивания неравномерной чувствительности фотоприемников сканирующих линеек тепловизоров, Белоконев В.М., Дегтярев Е.В., Рудый И.В., Малышев И.А., Павлова В.А., Тетерин В.В., Демеш О.В., Кабанов В.Ф. - Опубл. 10.06.98, Бюллетень №16.

2. Павлова В.А. и др. "Иконический подход к решению проблемы коррекции неоднородностей чувствительности многоэлементных МФПУ в сканирующих тепловизорах", "Оптический журнал", том 64, №2, 1997.

3. David L.Perry, Eustase L.Dereniak. Linear theory of nonuniformity correction in infrared staring sensors. - "OPTICAL ENGINEERING", August 1993, Vol.32, №8, 1854-1859.