Устройство адаптивной фильтрации видеосигналов

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

 

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

При решении задач сглаживания изображений с целью ослабления действующих аддитивных помех рассматривается модель статического изображения, которая представляет собой двумерную дискретную последовательность Yi,j, i = 1, N ¯ , j = 1, M ¯ вида:

Y i , j = S i , j + η i , j , i = 1, N ¯ , j = 1, M ¯ , ( 1 )

где Si,j - полезная двумерная составляющая (исходное неискаженное изображение), ηi,j - аддитивная шумовая составляющая, N - количество строк, М - количество столбцов двумерного массива изображения.

Изображение можно рассматривать как прямоугольную матрицу Y={Yi,j}, со строками i и столбцами j, где N и М - определяют размер матрицы изображения в пикселях.

В случае динамического изображения рассматривается модель, которая представляет собой двумерную дискретную последовательность Y i , j ( k ) , i = 1, N ¯ , j = 1, M ¯ , k = 1, K ¯ , вида:

Y i , j ( k ) = S i , j ( k ) , i = 1, N ¯ , j = 1, M ¯ , k = 1, K ¯ , ( 2 )

где S i , j ( k ) - полезная двумерная составляющая (исходное неискаженное изображение), S i , j ( k ) - аддитивная шумовая составляющая, N - количество строк, М - количество столбцов двумерного массива изображения, К - количество кадров двумерного массива динамического изображения.

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

Подобная задача может возникнуть: 1) в работе приемо-передающих устройств дальней или космической связи; 2) в радиотехнике при обработке сигналов; 3) в системах цифровой обработки динамических изображений; 4) в метеорологии и экономике при обработке результатов измерений.

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

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

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

Известен способ скользящего среднего [Гонсалес Р. Цифровая обработка изображений / Р.Гонсалес, Р.Вудс. - М.: Техносфера. - 2005. - 1072 с.]. Для его использования достаточно одной реализации Yi,j, i = 1, N ¯ , j = 1, M ¯ исходного процесса.

Для исходного изображения определяется размер маски сглаживающего фильтра m, т.е. натуральное число m<N. Способ скользящего среднего предполагает запоминание исходного изображения Yi,j, i = 1, N ¯ , j = 1, M ¯ , определение размера маски фильтра m (ширины «скользящего окна»), для которого производится вычисление среднего арифметического, S ¯ i , j = 1 m m i = 1 m j = 1 m Y i , j , замену центрального из значений Yi,j, i = 1, N ¯ , j = 1, M ¯ найденным средним S ¯ i , j , сдвиг «скользящего окна» на одно значение вправо, вычисление среднего арифметического выбранных значений реализации, и так до тех пор, пока маска фильтра не переместится по всему изображению.

Ширину "окна" выбирают нечетной, т.к. сглаженное значение рассчитывается для центрального значения.

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

Недостатками известного устройства являются:

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

- способ скользящего среднего вызывает автокорреляцию остатков, даже если она отсутствовала в исходной полезной составляющей (эффект Слуцкого-Юла).

Причины, препятствующие достижению требуемого технического результата, заключаются в следующем:

- если ширина "окна" сглаживания равна 2р+1, то первые p и последние p значений исходной реализации результатов измерений не подвергаются обработке;

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

В качестве нелинейных фильтров используются фильтры, основанные на порядковых статистиках [Гонсалес Р. Цифровая обработка изображений / Р.Гонсалес, Р.Вудс. - М.: Техносфера. - 2005. - 1072 с.]. Отклик такого фильтра определяется предварительным упорядочиванием (ранжированием) значений пикселей, покрываемых маской фильтра, и последующим выбором значения, находящегося на определенной позиции упорядоченной последовательности (т.е. имеющего определенный ранг). Фильтрация сводится к замещению исходного значения (в центре маски) на полученное значение отклика фильтра. Наиболее известен медианный фильтр, который заменяет центральное значение маски фильтра на значение медианы распределения всех значений результатов измерений, принадлежащих области маски фильтра. Чтобы выполнить медианную фильтрацию для элемента изображения, необходимо сначала упорядочить по возрастанию значения пикселей внутри маски, затем найти значение медианы и присвоить полученное значение обрабатываемому элементу.

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

Недостатками известного устройства являются:

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

- вследствие нелинейности метода обработки нельзя строго разграничить влияние медианной фильтрации на сигнал и шум;

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

Причины, препятствующие достижению требуемого технического результата, заключаются в следующем:

- медианная фильтрация является нелинейным методом обработки;

- зависимость эффективности сглаживания результатов измерений от формы полезной и шумовой составляющей.

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

Известен способ наименьших квадратов и устройство для кусочно-линейной аппроксимации [Бендат Дж., Пирсол А. Прикладной анализ случайных данных: Пер. с англ. - М.: Мир, 1989. - 540 с., авторское свидетельство №1624479]. Для использования данного способа достаточно одной одномерной реализации Y1, Y2, …, YN исходного процесса.

Способ наименьших квадратов позволяет для результатов измерений Yl, Y2, …, YN исходного процесса получить оценку, k = 1, N ¯ , минимизируя целевую функцию вида:

k = 1 N ( Y k S ¯ k ) 2 min

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

Дифференцируя выражение (3) по а и b и приравнивания к нулю, получаем систему линейных уравнений:

Решением системы является:

b = k = k = 1 N k Y k 1 N k = 1 N k k = 1 N Y k k = 1 N k 2 1 N ( k = 1 N k 2 )

При оценке сумма квадратов отклонений значений оценки от значений реализации измерений является минимальной.

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

Недостатками известного способа являются:

- при использовании данного способа необходима априорная информация о функции полезного сигнала;

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

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

- ограниченность способа наименьших квадратов к распараллеливанию и построению системы многоканальной обработки.

Причины, препятствующие достижению требуемого технического результата, заключается в следующем:

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

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

Известен способ и устройство фильтрования шума видеосигналов (патент №2364937, МПК G06K 9/40). Изобретение относится к области обработки видеоинформации, и в частности, к способам пространственно-временного устранения шумов в видеоизображениях. Техническим результатом является повышение скорости фильтрования шумов видеоизображения с обеспечением высокого качества шумоподавления без внесения искажений. Способ фильтрования шума видеосигналов устойчив к движению и основан на определении локальной структуры и на нелокальном усреднении в соответствии с определенной структурой кадра и с использованием данных о текущем кадре и соседних кадрах, при котором определяют локальную структуру изображения путем последовательного свертывания предопределенных шаблонов с соседними пикселями и путем отбора образцового шаблона RPc, который выдает наименьшую ошибку после свертывания, выполняют пространственно-временную фильтрацию с помощью взвешенного усреднения значений пикселя в поисковом окне от опорного кадра, причем результат пространственно-временной фильтрации вычисляют следующим образом:

Re s P x l = 1 N o r m 2 T / 2 T / 2 K / 2 K / 2 L / 2 L / 2 w f p i x r ( x r + p , y r + s , t r )

где pixr(xr, yr, tr) обозначает опорный пиксель от опорного кадра tr=tc+t с координатами хr и yr, индекс p изменяется от L 2 до L 2 в направлении X, индекс s изменяется от K 2 до K 2 в направлении Y, индекс t изменяется от T 2 до T 2 в области времени, Т - число опорных кадров, вовлеченных в процесс фильтрования, wf - вес опорного пикселя pixr(xr, yr, tr), нормирующий множитель вычисляют как:

N o r m 2 = T / 2 T / 2 K / 2 K / 2 L / 2 L / 2 w f

вес wf каждого опорного пикселя определяют следующим образом

w f = { w ( x r , y r , σ r , σ f ) е с л и R P c = R P r 0 в п р о т и в н о м с л у ч а е ,

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

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

Недостатками известного способа являются:

- при использовании данного способа необходима априорная информация о функции полезного сигнала;

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

- ограниченность способа наименьших квадратов к распараллеливанию и построению системы многоканальной обработки;

- ресурсоемкие вычислительные затраты.

Причины, препятствующие достижению требуемого технического результата, заключается в следующем:

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

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

Наиболее близким к изобретению является адаптивный двумерный способ размножения оценок и устройство, его реализующее (патент №2461874, МПК G06F 17/17). Рассматриваемое устройство-прототип предполагает формирование адаптивных областей двумерного сигнала для каждого значения пикселя, в связи чем задаются восемь направлений h = 1 ,8 ¯ , в которых определяются интервалы квазистационарности (фиг.1). Условие квазистационарности проверяется с помощью вычисления случайной величины τ, равной сумме числа инверсий значений пикселей в каждом из восьми направлений двумерного сигнала Yi,j, i = 1, N ¯ , j = 1, M ¯ .

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

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

Процедура получения области Ω и вычисление оценки S ˜ i , j повторяется для каждого значения пикселя Yi,j, i = 1, N ¯ , j = 1, M ¯ , при этом формируется весовая функция Wi,j, i = 1, N ¯ , j = 1, M ¯ , значения которой равны количеству размноженных оценок для каждого пикселя.

Результирующая оценка изображения S ¯ i , j определяется как среднее арифметическое размноженных адаптивных оценок:

S ¯ i , j = t = 1 K S ˜ i , j t W i , j , i = 1, N ¯ , j = 1, M ¯

Недостатками известного устройства-прототипа являются:

- невозможность реализации известного способа в реальном масштабе времени;

- большие вычислительные затраты.

Причины, препятствующие достижению требуемого технического результата, заключаются в следующем:

- для использования способа размножения необходимо запоминать всю входную реализацию, то есть кадр видеопоследоватльности.

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

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

На основе адаптивного двумерного способа размножения оценок и устройства, его реализующего, предлагается устройство адаптивной фильтрации. На первом этапе происходит выделение кадра из видеопоследовательности. Далее для формирования адаптивных областей двумерного сигнала для каждого значения пикселя задаются восемь направлений h = 1 ,8 ¯ (фиг.1), в которых определяются интервалы квазистационарности. Условие квазистационарности проверяется с помощью вычисления случайной величины τ, равной сумме числа инверсий значений пикселей в каждом из восьми направлений (фиг.1) двумерного сигнала Yi,j, i = 1, N ¯ , j = 1, M ¯ .

Например, сумма числа инверсий для направления 5 равна:

τ d = l = 0 d 1 k = l + 1 d u ( Y i + l , j , Y i + k , j ) , d = 1, R ¯

u ( Y i + l , j , Y i + k , j ) = { 1, Y i + l , j > Y i + k , j , 0, Y i + l , j Y i + k , j .

где i=const, j=const; Yi,j - текущее значение пикселя изображения с координатами (i, j); Yi+l,j, l=i+1…d - последующие значения пикселей изображения по j-му столбцу (движение в направлении 5), d≤R, R - максимальная длина интервала квазистационарности.

Количество сочетаний, для которых вычисляется сумма инверсий, составляет:

X = d ( d 1 ) 2 .

Первая альтернатива (убывающий сигнал) принимается, если:

τ d > c 1 = τ α 2 X + X / 2 0,5

Правило для принятия второй альтернативы (возрастающий сигнал) имеет вид:

τ d > c 2 = τ 1 α 2 X + X / 2 0,5

Гипотеза о стационарности сигнала принимается, если с2≤τd≤с1, где α - априорно задаваемое значение ошибки первого рода.

По полученным границам интервалов для каждого из восьми секторов, образованных направлениями 1-2, 2-3, 3-4, 5-6, 7-8, 8-1, формируются области квазистационарности. Для этого используется линейная интерполяция границ смежных интервалов уравнением прямой, проходящей через две точки:

j j 1 j 2 j 1 = i i 1 i 2 i 1 , i = i 2 i 1 j 2 j 1 j + i 1 j 2 i 2 j 1 j 2 j 1 ,

где (i1, j1) - координаты границы направления h, (i2, j2) - координаты границы направления h+1.

Для направлений h = 1 ,8 ¯ на изображении для одномерных реализаций, полученных из значений пикселей по вертикальным и диагональным направлениям от центрального пикселя исходного изображения, определяются интервалы квазистационарности с помощью способа инверсий. Данные границы позволяют получить интервалы с монотонным изменением яркости сигнала. Далее все восемь полученных секторов объединяются в одну область Ω. Таким образом, для каждого пикселя формируется окрестность пикселей, близких по значению яркости.

Значения пикселей, попавшие между всеми направлениями и интерполирующими прямыми, проходящими через границы интервалов квазистационарности, объединяются в одну область Ω и аппроксимируются поверхностью первого порядка вида , где ωi,j, i = 1,2 R 1 ¯ , j = 1,2 R 1 ¯ значения бинарной маски, которая принимает значения, равные:

ω i , j = { 1, i , j Ω 0, i , j Ω , i = 1,2 R 1 ¯ , j = 1 ,2 R 1 ¯

Значения коэффициентов А, В и С определяются с помощью двумерного способа наименьших квадратов [Патент №2362207 Российская Федерация, С2, МПК G06F 17/17. 2007127727/09; Заяв. 19.07.2007; Опубл. 20.07.2009, Бюл. №20], для нахождения которых минимизируется целевая функция вида:

i = 1 2 R 1 j = 1 2 R 1 ( ( Y i , j S ˜ i , j ) ω i , j ) 2 min

i = 1 2 R 1 j = 1 2 R 1 ( ( Y i , j A i B j C ) ω i , j ) 2 min ( 2 )

Дифференцируя выражение (2) по А, В и С и приравнивания к нулю, получаем систему линейных уравнений:

{ ( i = 1 2 R 1 j = 1 2 R 1 ( Y i , j A i B j C ) ω i , j ) 2 / A = 0, ( i = 1 2 R 1 j = 1 2 R 1 ( Y i , j A i B j C ) ω i , j ) 2 / B = 0, ( i = 1 2 R 1 j = 1 2 R 1 ( Y i , j A i B j C ) ω i , j ) 2 / C = 0 ; ( 3 )

Решением системы (3) является:

A = i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j i ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j j 2 i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j j ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j i j + i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j j j + i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i j ) D

B = i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j i ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j j + i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j i ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j j j i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i ) D

C = i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j i ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j j j + i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i j ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j i j i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i ) D + + i = 1 2 R 1 j = 1 2 R 1 Y i , j ω i , j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i i = 1 2 R 1 j = 1 2 R 1 ω i , j i j i = 1 2 R 1 j = 1 2 R 1 ω i , j i j ) D

D = i = 1 2 R 1 j = 1 2 R 1 ω i , j j j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i + i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j i ) i = 1 2 R 1 j = 1 2 R 1 ω i , j i j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j i = 1 2 R 1 j = 1 2 R 1 ω i , j i j 2 i = 1 2 R 1 j = 1 2 R 1 ω i , j i i = 1 2 R 1 j = 1 2 R 1 ω i , j j ) i = 1 2 R 1 j = 1 2 R 1 ω i , j j ( i = 1 2 R 1 j = 1 2 R 1 ω i , j j i = 1 2 R 1 j = 1 2 R 1 ω i , j i i

Процедура получения области Ω и вычисления оценки S ˜ i , j повторяется для каждого значения пикселя Yi,j, i = 1, N ¯ , j = 1, M ¯ , при этом формируется весовая функция Wi.j, i = 1, N ¯ , j = 1, M ¯ , значения которой равны количеству размноженных оценок для каждого пикселя.

Результирующая оценка S ¯ i , j t (где k=1:К) определяется как среднее арифметическое размноженных оценок S ¯ i , j t = k = 1 K S ˜ i , j k K , i = 1, N ¯ , j = 1, M ¯ , t - текущий номер кадра видеопоследовательности. Для увеличения эффективности обработки динамических изображений предлагается использовать корреляционную связь между кадрами. Предполагается, что кадры в пределах небольшого интервала времени изменяются незначительно, при этом в качестве пространства обработки используется зависимость амплитуда - номер кадра. В связи с этим предлагается использовать накопление для L оценок кадров с последующим их усреднением и получением результирующей оценки t-го кадра видеопоследовательности.

Алгоритм состоит из следующих этапов:

- выделение кадра из последовательности;

- записываются значения входного двумерного сигнала Yi,j, i = 1, N ¯ , j = 1, M ¯ ;

- задаются параметры способа: значение вероятности ошибки первого рода α и максимальная длина интервала квазистационарности R;

- вычисляется случайная величина τ d = l = 0 d 1 k = l + 1 d u ( Y i + l , j , Y i + k , j ) , d = 1, R ¯ и пороговые значения c 1 = τ α 2 X + X / 2 0,5 , и c 2 = τ 1 α 2 X + X / 2 0,5 ;

- проверяется гипотеза о стационарности сигнала с2≤τd≤с1, при выполнении которой происходит формирование бинарной маски ωi,j, i = 1,2 R 1 ¯ , j = 1,2 R 1 ¯ ;

- вычисляется оценка S ˜ i , j , с помощью аппроксимации значений Yi,j поверхностью первого порядка вида S ˜ i , j = ( A i + B j + C ) ω i , j , где i = 1,2 R 1 ¯ , j = 1,2 R 1 ¯ и вычисляется весовая функция Wi,j, i = 1, N ¯ , j = 1, M ¯ ;

- данная процедура повторяется для всех значений пикселей исходного изображения Yi,j i = 1, N ¯ , j = 1, M ¯ , в результате чего получается набор оценок S ˜ i , j ;

- определяется результирующая оценка S ¯ i , j t (где k=1:К) как взвешенная сумма размноженных адаптивных оценок S ¯ i , j t = k = 1 K S ˜ i , j k W i , j , где i = 1, N ¯ , j = 1, M ¯ , t - текущий номер кадра видеопоследовательности.

Устройство адаптивной фильтрации видеосигналов (фиг.2) содержит блок хранения входной реализации 1, вход которого является информационным входом устройства, к выходу которого подключен вход блока выделения кадра 2, выход которого подключен к первому входу поэлементного умножителя 5 и первому входу блока определения участков квазистационарности 3, к выходу которого подключен вход блока формирования маски 4, к выходу которого подключен второй вход поэлементного умножителя 5, к выходу которого подключен вход блока аппроксимации 6, к выходу которого подключен первый вход блока хранения оценок 7, к выходу которого подключен первый вход блока усреднения оценок 8, к выходу которого подключен вход блока накопления кадров 9, к выходу которого подключен вход блока усреднения кадров 10, выход которого подключен к входу блока хранения оценки полезной составляющей 11, чей выход является информационным выходом устройства; выход счетчика текущей строки 12.1 подключен ко второму входу блока определения участков квазистационарности 3 и входу блока задержки 12.2, к выходу которого подключены четвертый вход поэлементного умножителя 5 и вход блока задержки 12.3, к выходу которого подключен второй вход блока хранения оценок 7 и первый вход генератора сигнала усреднения 12.7, к выходу которого подключен второй вход блока усреднения оценок 8; выход счетчика текущего столбца 12.4 подключен к третьему входу блока определения участков квазистационарности 3 и входу блока задержки 12.5, к выходу которого подключены третий вход поэлементного умножителя 5 и вход блока задержки 12.6, к выходу которого подключены третий вход блока хранения оценок 7 и второй вход генератора сигнала усреднения 12.7; синхронность работы устройства обеспечивается генератором тактовых импульсов 13.

Устройство адаптивной фильтрации видеосигналов реализуется следующим образом. Видеосигнал поступает в блок хранения входной реализации, далее выделяется кадр видеопоследовательности. В каждый момент времени выбирается одна из точек исходного кадра. Для нее находится область с пикселями близкой яркости с помощью способа инверсий. Для каждой полученной области производится аппроксимация значений двумерного сигнала плоскостью, описывающейся уравнением первой степени с помощью двумерного способа наименьших квадратов. Таким образом, получается множество перекрывающихся оценок полученных областей. Результирующая двумерная оценка полезной составляющей определяется как среднее арифметическое по объему полученных оценок. Для увеличения эффективности обработки динамических изображений предлагается использовать корреляционную связь между кадрами. Предлагается использовать накопление для L оценок кадров с последующим их усреднением и получением результирующей оценки t-го кадра видеопоследовательности. Полученные значения поступают на выход устройства.

Устройство адаптивной фильтрации видеосигналов работает следующим образом. В блок хранения входной реализации 1 записывается видеосигнал, далее выделяется кадр видеопоследовательности в блоке выделения кадра 2. Блоки счетчика текущей строки 12.1 и счетчика текущего столбца 12.4 формируют координаты пикселя, обрабатываемого в данный момент времени. Координаты поступают на вход блока определения участков квазистационарности 3. Происходит выборка значений пикселей, близких к выбранному из блока выделения кадра 2, и определение координат пикселей, принадлежащих области близкой яркости. Полученные координаты передаются на блок формирования маски 4, где формируется маска фиксированного размера, содержащая «0» в местах с координатами пикселей, не принадлежащих полученной области, и «1» для принадлежащих. Центр маски соответствует текущему выбранному пикселю, координаты которого являются координатами маски. Полученная маска подается на поэлементный умножитель 5, где путем поэлементного умножения пикселей исходного изображения на маску с учетом ее координат выделяется область, поступающая дальше на блок аппроксимации 6. Оценка, полученная путем аппроксимации поверхностью первого порядка способом наименьших квадратов, запоминается в блоке хранения оценок 7, в котором кроме самой оценки запоминаются и ее координаты, поступившие с счетчиков номера строки 12.1 и столбца 12.4 через блок задержки. Последовательно выполняется перебор всех пикселей изображения, полученные оценки запоминаются в блоке хранения оценок 7. По окончании перебора блок генерации сигнала усреднения 12.7 подает на блок усреднения оценок 8 сигнал, после чего сохраненные в блоке 7 оценки усредняются. Далее двумерная оценка полезной составляющей поступает на вход блока накопления кадров 9, в котором накапливаются L оценок кадров, которые поступают в блок усреднения кадров 10, в котором получается результирующая оценка кадра видеопоследовательности, в блоке хранения оценки полезной составляющей 11 запоминается полученная оценка полезной составляющей. Синхронность работы устройства обеспечивается генератором тактовых импульсов 13.

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

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



 

Похожие патенты:

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

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

Изобретения относятся к вычислительной технике и могут быть использованы для обнаружения неисправностей спутников и корректировки таких неисправностей. Техническим результатом является возможность определения типа неисправности.

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

Изобретение относится к средствам обработки информации для прогнозирования стационарных и нестационарных случайных процессов. .

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

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

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

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

Изобретение относится к автоматике и вычислительной технике и может быть использовано для прогнозирования стационарных и нестационарных случайных процессов. Технический результат заключается в возможности получения оценки второй производной по формуле численного дифференцирования для равноотстоящих узлов функции в n-й (текущей) расчетной точке предыстории входной сглаженной дискретной последовательности. Для этого в адаптивное цифровое прогнозирующее и дифференцирующее устройство введен третий субблок расчета второй производной в первой n-й (текущей) расчетной точке предыстории входного процесса, содержащий блок инверторов и сумматор, на вход первого слагаемого которого через блок инверторов заведен выход сумматора первого вычитателя, вход второго слагаемого сумматора соединен с выходом мультиплексора блока адаптации, а выход сумматора субблока является четвертым информационным выходом устройства. 5 ил.

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

Заявленное изобретение относится к области термометрии и может быть использовано для коррекции на основе квантовой теории температуры радиационного термометра. Заявленный способ позволяет обнаружить параметры, отражающие структуру энергетических уровней, посредством приспособления эффективной физической модели, чтобы калибровать систему радиационного термометра. Для этого задают систему радиационного термометра в состояние измерения температуры, чтобы измерять температуру объекта, и обнаруживают значение энергии излучения объекта через оптическую систему. Определяют параметры, отображающие структуру энергетических уровней, посредством вычисления и обработки через PC или MCU в системе радиационного термометра согласно физической модели. Используют указанные параметры, отображающие структуру энергетических уровней, для калибровки радиационного термометра. Для реализации указанного способа коррекции температурных данных предложена также система радиационного термометра, содержащая модуль измерения температуры, модуль калибровки, PC или MCU, используемый для калибровки радиационного пирометра согласно параметрам, отображающим структуру энергетических уровней. Технический результат - повышение точности результатов измерений. 2 н. и 10 з.п. ф-лы, 6 ил., 1 табл.

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

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

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

Изобретение относится к устройствам прогнозирования результатов измерений. Технический результат заключается в повышении достоверности результатов измерения за счет адаптации времени прогноза. Устройство содержит блок хранения входной реализации, выполненный в виде последовательно соединенных N линий задержки, выходы которых, а также вход первой линии задержки, соединены с входами блока аппроксимации полиномом первой степени и входами блока оценки среднего уровня сигнала. Выход блока аппроксимации полиномом первой степени соединен с входом блока определения периода упреждения, а также с первым входом блока прогноза и первым входом блока анализа. Выход блока прогноза соединен с входом блока обратного преобразования, выход которого соединен с первым входом коммутатора и третьим входом блока анализа, а второй вход коммутатора и четвертый вход блока анализа соединены с входом блока нелинейного преобразования. Выход блока оценки среднего уровня сигнала соединен с пятым входом блока анализа, выход которого соединен с управляющим входом коммутатора. Выход коммутатора является выходом всего устройства. 4 з.п. ф-лы, 10 ил.

Изобретение относится к области вычислительной техники, а именно к анализу и обработке изображений. Технический результат - обеспечение реконструкции значений пикселей динамических двумерных сигналов в условиях неполной априорной информации. Устройство редактирования видеопоследовательностей содержит: блок хранения данных, блок хранения пикселей, блок создания словаря, блок хранения словаря, блок поиска подобия, блок обработки, блок вычисления приоритета, блок заполнения изображения, блок управления, блок хранения маски, блок выбора меток случайным образом, блок поиска меток, блок задержки, блок хранения кадров, блок установки меток пользователя, блок выделения объектов с помощью альфа-канала, генератор тактовых импульсов. 4 ил.

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

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

Наверх