Основы цифровой обработки сигналов

4.4.3.    Влияние неидеальности восстанавливающих фильтров

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

,                         (4.86)

где  – двумерная интерполяционная функция.

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

,                       (4.87)

где  обозначает оптимальную интерполяционную функцию, например заданную формулой (4.81) или (4.83).

В этом случае ошибка восстановления для произвольной точки изображения равна:

.     (4.88)

Величина ошибки восстановления зависит от двух факторов:

· в реальной системе интерполяционная функция  может отличаться от идеальной интерполяционной функции ;

· восстановленное изображение имеет конечные размеры и поэтому интерполяционные функции отсекаются на его границах.

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

Влияние неидеальности восстанавливающего фильтра удобно исследовать, рассматривая структуру спектра восстановленного изображения, описываемого соотношением (4.77):

.               (4.89)

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

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

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

.                                                 (4.90)

Здесь интеграл

                              (4.91)

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

                                    (4.92)

описывает энергию таких составляющих идеально восстановленного изображения.

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

,                                                        (4.93)

где

                                (4.94)

есть полная энергия восстановленного изображения, а

                                                     (4.95)

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

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

                         (4.96)

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

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

Таблица 4.2 Ошибка восстановления и коэффициент потери резкости для различных разделимых двумерных интерполяционных функций

Функция

Коэффициент потери резкости , %

Ошибка восстановления ,%

sinc (x)

0,0

0,0

Прямоугольная

26,9

15,7

Треугольная

44,0

3,7

Колоколообразная

55,4

1,1

Кубический В-сплайн

63,2

0,3

Гауссова,

38
,6

10,3

Гауссова,

54,6

2,0

Гауссова,

66,7

0,3