CC BY
45
ОБРАБОТКА ИЗОБРАЖЕНИИ
ИДЕНТИФИКАЦИЯ МОДЕЛЕЙ СИСТЕМ ФОРМИРОВАНИЯ ИЗОБРАЖЕНИЙ В КЛАССЕ ФИЛЬТРОВ С БЕСКОНЕЧНОЙ ИМПУЛЬСНОЙ ХАРАКТЕРИСТИКОЙ
В. А. Фурсов, С. И. Парфенов Институт систем обработки изображений РАН
В работе рассматриваются особенности задачи идентификации линейных моделей искажающих систем в классе фильтров с бесконечной импульсной характеристикой при отсутствии априорной информации о распределении ошибок измерений. Развивается алгебраический подход, основанный на непосредственном уточнении вектора ошибок измерений по небольшому числу его компонент, задаваемых с использованием дополнительной информации, содержащейся в невязках. Приводятся результаты экспериментов, показывающие принципиальную возможность существенного повышения точности оценок по сравнению с обычным методом наименьших квадратов (МНК).
1. Постановка задачи
В работе [ 1 ] рассматривалась задача идентификации моделей систем формирования изображений в классе фильтров с конечной импульсной характеристикой (КИХ-фильтров). В частности, обсуждались вопросы оценки информативности малых фрагментов и была показана возможность непосредственной идентификации восстанавливающих КИХ-фильтров, обладающих различными регуляризую-щими свойствами, путем простых линейных преобразований данных типа взвешивания. В настоящей работе эти результаты развиваются применительно к моделям типа фильтров с бесконечной импульсной характеристикой (БИХ-фильтров).
Модель БИХ-фильтра с опорными областями -А (на входном) и Д (на выходном) изображениях в виде прямоугольников Ь,хЬ2 и К,хК2 соответственно описывается выражением вида [2]:
Л«ИЯ2)= X Ха(/1»/2)*(я1 -4**2 -/2>-/, =0/2=0
А',-1 *->-1
X £А(*1>*2 )>’(«! ,«2 -А2) + ^(Я1,Л2 ),(!)
*,=0 А,=0
(*„*,) *(0,0),
где а(/|,/2),&(£),£2) - значения импульсной характеристики, х(п].п2), у(п],п2) - отсчеты поля яркости в точках опорных областей Ох и Ц соответственно, а дискретная последовательность, определяемая как
к, ЧК2-1
£(И|,И2) = £(я1*я2)+ X ХЗД.*2)х
А-, =0 А ,=0
х£(п, - к1, п2 - к2 )>(к\, к2 ) * (О.О).
(2)
Здесь Цпх,щ) - ошибки измерения отсчетов поля яркости выходного изображения (предполагается, что входное изображение-эталон задано точно).
При реализации БИХ-фильтров (1) в системах обработки изображений возникают принципиальные трудности, связанные с тем, что для вычисления вы-
ходных отсчетов у(пх,п2) используются отсчеты
у(я, —к\,п2-к2), к, = 0,(к, -*), к7 = 0,(к2-1) из той же опорной области. Поэтому для построения физически реализуемых фильтров необходимо согласовывать направления “обхода” изображения и форму опорной области,
Указанное ограничение не имеет значения при решении задачи идентификации. Действительно, для N отсчетов на фрагменте изображения в соответствии с (1) можно записать уравнение типа линейной регрессии [3]:
У = &+$, (3)
где вектор с формируется как
с =
а матрица в имеет вид
8 = [Хж*ж1=
(здесь для компактности дальнейших записей изменена система индексации входных х(п} - , п2 - /2),
/1 =0, , /2 =0,Ь2 и выходных у(щ -кьп2 -к2),
к1 = 0,(/Ст -1), к2 = 0,(/С2 -1) отсчетов, в частности, введены новые обозначения для размерностей описывающих их массивов: 1=£;х£2, К=(К,хКг)-1, М-Ь+К. Компоненты Мх1- векторов у и ^ формируются в соответствии с соотношениями (1) и (2) по N отсчетам используемого для идентификации фрагмента изображения.
С точки зрения задачи оценивания не имеет значения в каком порядке взяты строки матрицы 8 в (3),(4). Поэтому, если все отсчеты искаженного изображения уже сформированы, задача идентификации модели БИХ-фильтра может решаться точно также как и при использовании модели КИХ-филь-тра [1]. Принципиальное отличие заключается в том, что размерность вектора искомых параметров зависит от порядка передаточной функции системы, т.е. от размеров опорных областей, как на входном, так
1 »«г- [с\>с2> ...,см г.
*1,1 Л1 У\,к
*2,1 х2,1. У2,К У2,К (4)
**,1 ХН,1 Ук,\ Уиж_
и на выходном изображении. Кроме того, компоненты вектора ошибок £ в силу (2) всегда коррелированны.
По аналогии с [1] модель инверсного БИХ-фильтра запишем в виде (предполагается, что она существует)
К,-1К2-1
х(п1,п2)= £ ЕК^ЛМ^-кьп2-к2)~ к} =0 к2=0
Л|—1Х-2-1
~ Е 2^а(игМи1 ~1\>пг-^)+¥{п\,»2\ (5)
1,=01г=0
Сл»^2 ) ^ (о,о).
Для N - отсчетов на выбранном фрагменте в соответствии с (5) можно записать матричное соотношение:
х = Исг + у, (6)
где х, V}/- Мх1-векторы, составленные из отсчетов неискаженного изображения и ошибок измерений соответственно, с, - А/х /-вектор параметров инверсной модели, К - [Уж- :Х№] - ЛхА/ -матрица (здесь используются те же, что и в (4) обозначения для матриц и Хд^ , т.к. по-прежнему М-Ь+К
(хотя теперь Ь^(Ь1хЬ2)-\ и наоборот К=К]Х.К2).
Соответствующий инверсной модели (5) восстанавливающий БИХ-фильтр:
-1*2-] _
*(«!/ »*) = ЕЕ Ь{к1' кг)у(п\ - к\> пг - кг) -
^=0*2=0
-ЕЁ 5(А* - 1Л
А= о /5 = и
'Л * (0,0)
или в матричном виде
X = йсг! (8)
где х, И , с г - являются оценками фигурирующих в уравнении регрессии (6) х, К, с г и имеют те же размерности.
В работах [4,5] обсуждалась обшая схема алгоритма оценивания по малым фрагментам изображений. Кратко она сводится к последовательности двух основных этапов; предварительный анализ информативности данных и собственно вычисление оценок с использованием преобразований типа взвешивания и/или непосредственной корректировки данных.
В настоящей работе рассматриваются особенности реализации указанных этапов при решении задач оценивания х , с и/или сг, связанные с особенностями структуры БИХ-фильтра, а также исследуется новый способ построения модифицированных МНК-оценок, основанный на непосредственном уточнении вектора ошибок измерений по небольшому числу его компонент, задаваемых с использованием дополнительной информации, содержащейся в невязках.
2. Обоснование общей схемы
построения модифицированных МНК-оценок
Традиционно для построения алгоритмов идентификации используется статистическая теория оценивания [3], в рамках которой критерий качества задается на основе априорной информации о распределении ошибок. Такой подход неприемлем при решении задачи оценивания по малым фрагментам изображений. Обсудим эти особенности подробнее. Поскольку они в равной степени справедливы для моделей (3) и (6). далее для матриц Б и Я используется одинаковое (традиционное для задач линейной регрессии) обозначение - X.
Во-первых, вместо обычного предположения Яапк(Х)=М мы вынуждены констатировать более реалистичное: Яапк(Х)^ї, поскольку всегда существует опасность попадания малого фрагмента изображения на участок фона.
Во-вторых, мы должны отказаться от двух, наиболее важных с теоретической точки зрения, предположений: А//£/=0 и соу/£/=а*Еу, где ст~- дисперсия отклонений, соу{-} - ЛЧУ - матрица ковариаций отклонений, а Еу - единичная Ух У -матрица. Связано это с тем, что свойство устойчивости статистических характеристик сигналов при малом объеме данных не проявляется в полной мере.
Отказ от указанных предположений классической регрессии имеет существенное значение для обоснования общей схемы решения задачи. Известно, что в случае гауссовых помех с ковариационной
матрицей ошибок измерений К = М ^Г} оптимальным (в срсднеквадратическом смысле) является обобщенный метод наименьших квадратов (ОМНК).
ДО = |гСг| (9)
(здесь специально из соображений удобства сопоставления с предлагаемым далее методом введено обозначение С~=К1).
Соответствующая этому критерию оценка парамегров моделей типа (3), (6):
у = Хс+$ , (10)
определяется как
с = [х7С2х] 'хгС2у , (11)
а ошибка идентификации
Дс = с-с = [хгСгх]~‘хгС2*. (12)
Если статистические характеристики шумов не известны, выбор матрицы С‘ тру дно обосновать. В то же время из (11),(12) видно, что при заданных X, у изменяя матрицу С2 можно существенным
образом влиять на ошибку идентификации. Другой, очевидный, способ изменения ошибки идентификации (12)- непосредственная корректировка вектора у в уравнении регрессии (10).
Исходя из приведенных соображений предлагается следующая модифицированная МНК-оценка
с = |хгС2х] 'хгС2(у+ Ау), (13)
где в2 - так называемая, весовая матрица размерности МсЛГ, а Ду - корректирующий Л'х 1 -вектор. При С2=Е и Ду=0 оценка (13) является обычной МНК-оценкой.
В работе (1] было показано, что варьируя параметрами преобразований взвешивания можно непосредственно управлять качеством восстановления на основе визуальной оценки изображения. В частности, для выбора матрицы в2 использовались реализации вектора невязок:
£ = у-у = Хс-у. (14)
Поскольку вектор невязок ''появляется” лишь в ходе оценивания, схема получения оценки (13) должна быть итерационной.
На первом шаге строится МНК-оценка:
с = [хгх]'‘хгу. (15)
Затем вычисляется вектор невязок (14) и анализируется содержащаяся в нем дополнительная информация об ошибках в исходных данных. С использованием этой информации определяются параметры взвешивания и/или корректировки данных (С и Ду). С их использованием осуществляется преобразование исходного набора данных, вновь вычисляется МНК-оценка (15) и т.д.
В описанной схеме важное значение имеет надежность начальной оценки. Поэтому необходимым является этап предварительной обработки данных, заключающийся в проверке их информативности на заданном фрагменте изображения [3].
5. Оценка информативности данных
В работе [3.4] предложено в качестве меры информативности данных на фрагментах использовать показатель диагонального преобладания
ф = (ггА) //гА; информационной матрицы А=ХТХ.
В данном случае особенностью задачи является то, что матрица X (8 и/или И) в соответствии с (4), (6) состоит из двух подматриц - УЛХ и \м , составленных из соответственных отсчетов входного и выходного изображений. Если размеры и форма опорных областей на входном и выходном изображениях одинаковы, то фигурирующие в этих матрицах отсчеты повторяются (за исключением столбцов, соответствующих центральным отсчетам в опорных областях Д. и Д для моделей (3) и (6) соответственно). Если порядок реальной искажающей или восстанав-ливаюшей системы соответствует априори заданно-му порядку модели (I) и/или (5), то информативность данных, из которых составлена матрица \УК , является естественным следствием информативности данных, фигурирующих в матрице Хлх . С учетом сказанного, для предварительного анализа ин-
формативности данных предлагается использовать лишь одну из указанных матриц: \ык или Х№ . При этом алгоритм контроля информативности упрощается, а надежность принимаемых решений увеличивается.
4. Анализ точности идентификации при преобразованиях взвешивания
Введение в (13) весовой матрицы С2 равносильно тому, что исходная матрица X и вектор у подверглись линейному преобразованию вида
у = Су, Х = вХ (16)
При этом вектор ошибок £ преобразуется в
с, = . а ошибка (12) - в Дс = |хг х| Х£ .
В работе [5] показана возможность повышения точности путем преобразований (16) в классе вещественных матриц 0= Расширим
этот класс.
В [5] показано, что неравенство
К>Н£. <17)
где ||ЛсЦ^, |Лс||‘. - квадраты евклидовых норм век-
торов ошибок идентификации, эквивалентно следующему
нтЛт? + тЛт?нт0т^ +
+ т0т0тнт0т0т|<о,
где Н = СТУС V = х[гх]"‘хт.
На основании (18) можно сделать следующие выводы. Если вектор ошибок £ принадлежит нуль-пространству матрицы X (Т?г £ = 0), ника ким выбором весовой матрицы нельзя улучшить точность оценок. Действительно, квадратичная форма
с,1 Т0 Тпг С1УСТ0 Т0' ^ обращается в нуль лишь при
С=аЕЛ, где а любое число. Для всякой др угой весовой матрицы С собственные значения матрицы
Т0ТП‘ в1 УСТ0Т0т неотрицательны и одновременно не обращаются в нуль, поэтому неравенство (18) (для отличающегося от нулевого вектора ошибок 4) невозможно. Этот вывод не является неожиданным. Если вектор ошибок принадлежит нуль-пространству, ошибка оценивания равна нулю, а всякое изменение его ориентации относительно пространства матрицы X может лишь ухудшить ситуацию.
Если вектор ^ принадлежит пространству матрицы X, изменением весовой матрицы С также ничего сделать нельзя, т.к. при Т^ = 0 все слагаемые в левой части (18) обращаются в нуль. Во всех остальных случаях (Т] £ ф 0, Т0Г 2, ф 0) в к лассе диагональных матриц элементы которых могут быть, как вещественными, так и комплексными (без вещественной части), всегда можно подобрпть, и при-
том не единственную, матрицу для которой в (17), (18) будет иметь место строгое неравенство. Более того, всегда можно подобрать диагональную матрицу в такую, что:
Хг| = ХгС2? = 0. (19)
Действительно, поскольку матрица С2 диагональная, в равенстве (19) ничего не изменится, если мы представим его в виде
ХгЕ-г = 0, (20)
где Е = ,$,2..^Лг) - диагональная Л-хА' -матрица, а г = [^, Ґ - ЛГх 1 - вектор, компонен-
ты которого составлены из элементов диагональной матрицы Є2
Из недоопределенной системы (20) (АГ>А 1) всегда может быть найдено, и притом не единственное, решение для искомого вектора / Поскольку
2і = ёр, / = 1, Аг, отрицательным компонентам этого вектора будут соответствовать комплексные числа (элементы матрицы О с нулевой вещественной частью.
Если для вычисления оценок с используется система нормальных уравнений, комплексность весовой матрицы Є не создает дополнительных проблем. Из соотношения (11) видно, что при формировании этой системы используется не исходная матрица в (возможно с “чисто” мнимыми элементами), а ее квадрат: Є2.
5. Построение процедур преобразования данных
Из соотношений (19),(20) видно, что для определения весовой матрицы. обеспечивающей нулевую ошибку идентификации (как и следовало ожидать) необходимо знание вектора который, конечно, никогда не известен. Однако мы можем приближенно оценить его из следующих соображений. Можно показать, что
т,Ч = т0т£, (21)
где Т0 - Мх(М-М) матрица, составленная из нормированных собственных векторов, соответствующих нулевым собственным значениям матрицы ХХТ Равенство (21) утверждает, что проекции векторов
ошибок £ и невязок Ё, на базис нуль-пространства
матрицы X совпадают.
Из (21), в частности, видно, что если вектор ошибок принадлежит нуль-пространству, он совпадает с вектором невязок. Ясно, что существует множество (не нулевой меры) векторов ошибок. ДІЯ которых имеет место близость к соответствующим векторам невязок. Признаком принадлежности вектора ошибок к этому множеству может служить сравнительно большое значение нормы вектора невязок. Более того, есть основания ожидать, что выделяющиеся компоненты вектора невязок будут
близки к соответствующим компонентам вектора ошибок.
В связи с этим возникает задача оценки тех компонент вектора ошибок, которые не могут быть предсказаны непосредственно по невязкам или с использованием других источников дополнительной информации.
Введем в рассмотрение М<1 -вектор
представляющий собой прямую
сумму векторов \м и . Предположим, что вектор задан так, что А/ его компонент равны (возможно, приблизительно) соответствующим компонентам вектора ошибок: £м = £ м . Задача заключается в определении (Л-АУ)-вектора - составленного из оставшихся компонент вектора ошибок
I
Будем искать этот вектор из условия
Т„г(г -|)= 0. (22)
Основанием для этого служит соотношение (21). Действительно, если В>=В), то вектор “поправок”
= ^ = удовлетворяющий (22), ортого-
нален нуль-пространству. Ортогональность сохраняется и в случае, когда это равенство приблизительное: \ =£ . Таким образом, реализация условия (22) позволит оценить составляющую вектора ошибок, принадтежащую пространству матрицы X.
Указанное свойство оказывается особенно полезным при использовании БИХ-фильтров. Связано это с тем, что ошибки измерений, фигурирующие в модели (3), коррелированны в силу (2), поэтому угол между вектором ошибок и пространством матрицы X, как правило, оказывается действительно малым.
В соответствии с 'Щ) сформируем стлсйу уравнений
То.у-дДл'-.и = >/> (23)
где Т0у_л, - матрица, составленная из Лг-Л/ столбцов матрицы Тс! , (Л'-А/)х1 -вектор Хом определен как *0Л/ =Т0^-Т0ГД„, а Т1М - (Л'-Л/)хЛ/ матрица, составлена из М столбцов той же матрицы Т07 .
Так как КапкТ0=Л-А/, существует единственное решение системы (23):
~ -М *<Ш • (24)
Это решение удовлетворяет равенству (23). Поэтому, если Л/ компонент заданы равными истинным
значениям компонент вектора ошибок. %м =£Л,. остальные Л-А/ компонент, вычисленные в соответствии с (24), также будут совпадать с соответствующими истинными ошибками: м - . По-
видимому, при задании значений компонент близкими к истинным: Ё,Л/ * Ъ)М есть основания ожидать
близкое к решение:
Из (22) следует также, что при задании компонент вектора равными соответствующим компонентам вектора невязок: в соответствии
с (24) будут получены оставшиеся компоненты вектора невязок: Э70 не удивительно,
т.к. задавая компоненты подобным образом мы не вносим никакого нового знания Для того чтобы улучшить результат оценивания вектора ошибок,
необходимо при задании компонент вектора £,Л, использовать дополнительную информацию (текущую измерительную или априорную). В качестве такой информации может использоваться, например, обычно имеющая место при фиксированном полезном сигнале устойчивая связь между выделяющимися компонентами реализаций £ и
Найденный описанным способом вектор \ далее может использоваться для построения весовой матрицы в или непосредственной коррекции вектора у (вычитания из него вектора поправок Ау = ).
Ясно, что успех реализации описанного метода будет зависеть от того, насколько близкими к истинным окажу тся задаваемые компоненты вектора с, и .
6. Результаты экспериментов
Для моделирования искаженного изображения использовался БИХ-фильтр порядка (2.2). который описывается передаточной функцией
я(*1,22 )= я(г,,72 )/>!(*, ,*2),
где
В(г„г,) = йоо + +&_1021~1 + Ь012п +Ь0_,2-[ +
+ Ь[\2\г2 + 27 + + Ь_^.21 2г
А(?1927) = ] + а1а2} +ЙГ_102Г1 +йг0122 +а0_12'7' +
а.лгх2г +а +й1_,г,2“| + аиг^г2
Переменная 2\ соответствует горизонтальному направлению на изображении, а г2- вертикальном}-.
Параметры передаточной функции задавались так. чтобы обеспечивалась радиальная симметрия искажений:
*10 =*01 ~*-1Э ” VI' *11 -*1-1 “*-11 =*-1 1 а10 = а01 -°-\0 - аС-1 » а11 =а\-\ =ЙП = Л-1-1‘
При вьшолнении указанных условий полиномы числителя и знаменателя представляются в следующем компактном виде
в[2.,2г) = Л0[1 +А (г, + 2~1 + 2, -*-22'1| +
\%)1
Л(2х,22)=\ + а[(г] +г, 1 +г2 + г2’)+
+ я2(г,г2 +2[122 +21221 +г,'1г2)
где
Ь0 = 0.02834. ^ = 1.2,63 =0.27, я, = -1.2,я2 =0.27.
Числовые значения параметров специально подобраны так, чтобы обеспечивалась устойчивость, как искажающего, так и соответствующего ему восстанавливающего (инверсного) фильтра.
Шум с(«!,»2) моделировался в виде смеси гауссовского с параметрами (Мх= 0, Ох= 16) и с вероятностью Р=0.03 равномерного в интервале [0-64] импульсного шума. При моделировании использовался итерационный метод реализации БИХ-фильт-ров [2]. Исходное тестовое и полученное описанным способом искаженное изображения приведены на рис.1, (а и б).
По характеру искажений видно, что для их реализации с помощью КИХ-фильтра потребовалась бы маска значительных размеров. В данном случае размеры опорных областей на входном и выходном изображениях, соответствующие указанной передаточной функции, весьма малые: 3x3.
По сформированным тестовым изображениям осуществлялась идентификация 5-ти параметров БИХ-фильтра на фрагменте, содержащем 24 наблюдения. Задача оценивания решалась итерационно. Вначале строилась МНК-оценка и вычислялся вектор невязок. Затем с использованием полученного вектора невязок по соотношению (24) оценивался вектор ошибок, осуществлялось преобразование данных и вновь вычислялась МНК-оценка.
Обработка изображ* связана с решение»/ к их задач, в которы: входные, и выходны ные являются изобр ниями. Одним из пр ров служат системы редачи изображений </)оазоаботчики сталю б)
+*:(7!г2 + 2\'2'1 + 2\27 +г
Рис. I. Тестовые изображения: а) исходное, б) искаженное
Цель эксперимента заключалась в том, чтобы показать, принципиальную возможность увеличения точности идентификации по мере более точного задания небольшого числа компонент вектора ошибок. Для этого была проведена серия экспериментов, в которых М (в данном случае 5) компонент
оценок вектора ошибок задавались с возрастающей точностью в соответствии с формулой
1м = + (1 - 1.4>*|д/.
Для сравнения точности идентификации в ходе каждого эксперимента вычислялась скалярная характеристика относительной погрешности в виде
м/и , где М и Н - евклвдовы нормы векторов ошибок оценивания и истинных параметров модели. Результаты приведены в табл. 1.
Из таблицы видно, что при задании компонент
вектора £>м равными невязкам (а=0) повышения
ожидать, не происходит. Однако по мере более точного задания этих компонент происходи! увеличение точности и при (а=0.71) погрешность
составляет всего лишь около одного процента. Заметим, что погрешность МНК-оценок при этом составляет почти 20%,
точности по сравнению с МНК. как и следовало
Таблица 1
Параметры М/И
ЬО ы Ь2 а1 а2
Истинные 0.028340 0.034008 0.007652 -1.200000 0.270000
МНК-оценки 0.022962 0.057481 0.006793 -1.029265 0.114629 0.189
© II я 0.023002 0.057663 0.006808 -1.026394 0.111991 0.191
а=0.14 0.024483 0.050960 0.008005 -1.067358 0.149081 0.147
а=0.28 0.026047 0.044293 0.009344 -1.105419 0.183520 0.104
а=0.42 0.027685 0.037637 0.010813 -1.141070 0.215756 0.065
а=0.56 0.029445 0.031336 0.012390 -1.169678 0.241525 0.034
а=0.71 0.031265 0.025105 0.014 061 -1.195811 0.265037 0.011
а=0.86 0.033172 0.019205 0.015788 -1.216209 0.283283 0.022
а=1 0.035153 0.013615 0.017550 -1.231559 0.296894 0.039
В данном случае уточнению подвергались 19 компонент вектора ошибок по заданным 5-ти компонентам. Для получения дополнительной информации об этих компонентах могут использоваться различные источники: сопоставление текущих измерений, полученных от других приборов, предварительное моделирование различных вариантов реализаций векторов ошибок на заданном тестовом изображении и др. Оказывается, что такая устойчивая дополнительная информация содержится непосредственно в невязках.
В таблице 2 приведены результаты экспериментов по сопоставлению одноименных компонент реализаций векторов ошибок и невязок. В эксперименте на фиксированном информативном фрагменте изображения задавалось множество векторов ошибок одинаковой длины (это достигалось соответствующей нормировкой), имеющих случайное направление относительно пространства матрицы X. Для одноименных компонент векторов ошибок и невязок подсчитывалось число совпадений знаков, усредненное по 100 реализациям, отношение одноименных компонент и дисперсия этих отношений. Номера компонент упорядочены по убыванию их абсолютных значений.
Из таблицы видно, что выделяющимся компонентам вектора ошибок соответствуют выделяющиеся компоненты (с тем же знаком) вектора невязок, а коэффициент связи между ними достаточно устойчив.
В ходе экспериментов исследовалась также возможность непосредственной идентификации вос-
станавливающего БИХ-фильтра. Оценивались параметры передаточной функции
2 г) = I' ^2 )'
где
А(ги гг) = а0[ 1 + л (г, + 2Х 1 + 4-2;
+ а2(г,22 + ^-1-1 , _ _-1 21 2 2 12 + г\'гг\
В(г1Г 2г) = 1 + А, (г, + 1 + + 72')
+ Ь2(г^2 + г;1 + + 2\
Точные параметры этого фильтра, соответствующие приведенной выше модели искажений, соответственно равны. ас = 35.2857, а1 = -1.2,
а2 = 0.27. Ь, =1.2, Ь3 = 0.27 . Путем непосредственной идентификации инверсной модели были получены следующие оценки параметров: а0 = 34.761, а „а, = -41.05279,
а0а2 = 8.775, Ь, = 1.216, Ь, = 0.2669
С использованием полученных оценок осуществлялась обработка изображения “часы”, подвергшегося тем же искажениям, что и тестовое изображение “текст” (рис. 1) Эти изображения приведены на рисунке 2: а) - без шума б) - с добавлением указанной выше смеси гауссовых и импульсных помех. Полученные в результате обработки изображения приведены на рисунках 2 б) и г) соответственно.
а) искаженное (vet тупа)
іj) lioccmanait.ieiiiwc (fwt шума)
в)иск(Ністш>е (с шумом)
сосетаномтнос (с шумам)
Рис. 2 Результаты восстановления изображения “часы”
Таблица 2.
N П / n Относительное число совпадений знаков Отношение т\ Дисперсия отношений M/|i|
1 0.992 0.994 0.073
2 0.975 0.987 0.177
3 0.983 1.0480 0.210
4 0.933 1.0167 0.288
5 0.933 1.0689 0.414
6 0.942 1.0962 0.483
7 0.908 1.0250 0.516
8 0.949 1.0065 0.561
9 0.883 1.1810 0.811
10 0.891 1.1005 0.776
11 0.850 1.2369 1.162
12 0.808 1.2023 1.519
13 0.816 1.2514 1.436
14 0.841 1.3291 1.549
15 0.750 1.3395 2.198
16 0.791 1.5177 3.362
17 0.699 1.4881 3.915
18 0.675 1.9899 8.199
19 0.741 2.0329 6.805
20 0.649 2.3555 24.24
21 0.733 3.0389 52.83
22 0.592 3.9789 87.83
23 0.467 5.4769 276.9
24 0.433 12.956 64387
Заключение
Несмотря на некоторые трудности, связанные с реализацией БИХ-фильтров. вычислительная сложность сквозной технологии оценивания и восстановления при большой интенсивности искажений может быть меньшей, по сравнению с использованием моделей КИХ-фильгров. Если задача минимизации вычислительных затрат не ставится, то при использовании моделей БИХ-фильтра можно реализовать более высокое качество восстановления.
Рассмотренный в работе метод построения преобразований данных открывает новые принципиальные возможности для существенного повышения точности оценивания по малому числу наблюдений. В частности, метод оценки “недостающих’' компонент вектора ошибок по небольшой части заданных позволит непосредственно использовать априорную информацию об устойчивых связях между выделяющимися по абсолютной величине компонентами реализаций векторов ошибок и невязок. Эксперименты подтверждают существование такой устойчивой связи на фиксированных информативных фрагментах изображения.
Описанная в работе итерационная схема построения модифицированных МНК-оценок с промежуточным преобразованием данных и предложенный метод их корректировки развивают подход, основанный на высказывавшейся в [5J точке зрения, что задача идентификации по малому числу наблюдений. по существу, является алгебраической.
Благодарность
Работа выполнена при поддержке Российского фонда фундаментальных исследований (грант № 96-J 5-96026)
Литература
J. В. А. Фурсов. Восстановление изображений КИХ-фильтрзми, построенными путем непосредственной идентификации инверсного тракта. Компьютерная оптика. Вып. 16, 1996, с. - 103-108.
2. Dan Е. Dudgeon, Russel М. Mersereau. Multidimensional digital signal processing. Prentice - Hall, Inc., Englewood Cliffs. 1984.488 c.
3. Дсмиденко Е.З. Линейная и нелинейная регрессии -М.: Финансы и статистика, 1981, 303с.
4. Fursov V.A. Identification of optical distorting systems with selecting image informative fragments. Workshop on Digital Image Processing and Computer Graphics. Proceedings SPIE. - 1994. 2363.
5. Fursov; Vladimir A. Constructing unified identification algorithms using a small number of observations for adaptive control and navigation systems. Journal:Proc. SPIE Vol. 3087, p. 34-44, Nav igation and Control Technologies for Unmanned Systems II, Scott A. Speigle; Ed., 1997.
