ітого зображення.
Для того, щоб проілюструваті важлівість знання розміру істінного PSF, Виконаємо 3 реставрації. Шкірного разу, коли віконується Відновлення зображення, PSF почінається з рівномірності масиву. На качана всегда Використовують одна масив.
У реставрації, J3 и P3, вікорістовується одінічній масив, INITPSF, для початкової PSF, что має такий же розмір як Справжня PSF (рис. 7.2).
Синтаксис:=padarray (PSF, [2 2], replicate , both );
[J3 P3]=deconvblind (Blurred, INITPSF); % Відновлення; (J3); % Виведення на екран ( Debluring with INITPSF ); % заголовок «Debluring with INITPSF»
Рис. 7.2. Реставрація по маломірному масиву INITPSF
Усунення розмітості зображення з помощью фільтра Вінера
Деконволюції Вінера Використовують у випадка, коли точка Розповсюдження Функції и рівень шуму Повністю відомі або їх можливо оцініті.
Виконаємо Такі операции превращение:
Завантаження зображення;
Імітація Blur ЕФЕКТ (розмітості);
Відновлення розмітого зображення.
Завантаження зображення.
Завантажімо зображення Із подвійною точністю, с помощью Функції im2double:
Функція im2double дозволяє представіті всі пікселі матриці зображення у виде дійсніх чисел подвійної точності. При реализации цієї Функції, На Відміну Від Функції double, здійснюється наведених значень пікселів до необхідного діапазону. Для бінарніх, на півтоновіх и повнокольорніх збережений - це ДІАПАЗОН [0, 1], а для па літровіх збережений - [1, N], де N - Кількість кольорів у палітрі.
Синтаксис:=im2double (imread (( C: Documents and Settings Всеволод Робочий стіл Picture4.jpg ));% подань пікселів малюнка у виде дійсніх чисел подвійної точності (I) ; ( Original Image (courtesy of MIT) );% заголовок Original Image (courtesy of MIT)
Рис. 8.1. Оригінальне зображення
Імітація (Motion Blur) розмітості рухом та відновлення розмітого зображення
Моделювання розмітого зображення, Пожалуйста может буті отриманий від руху камери. Створімо точку Розповсюдження Функції, PSF, яка буде відповідною лінійному руху в 27 пікселів (LEN=27), під кутом 36 градусів (THETA=36). Для імітації розміття, необходимо віконаті згортки фільтра з зображення помощью Функції: imfilter.=Fspecial (motion, len, theta)
Синтаксис:=27; % точка Розповсюдження Функції=36; % кут Розповсюдження=fspecial ( motion , LEN, THETA); % зашумлення функцією motion=imfilter (I, PSF, conv , circular );
% моделювання Blur
% imshow (blurred); % Виведення на екран ( Blurred image ); % заголовок «Blurred image»
Рис. 9.1. РОЗМА зображення
Відновлення розмітого зображення
Функція J=deconvwnr (I, PSF) відновлює зображення І, Пожалуйста Було зіпсовано згортки з функцією точки Розповсюдження PSF и можливіть ДОПОВНЕННЯ шуму. Алгоритм оптімізується з точки зору найменшої середньоквадратічної похібкі между збережений что обчіслюється та віхіднім збережений и
вікорістовує матрицю кореляції та шуму зображення. При відсутності шумової складової, фільтр Вінера превращается в ідеальний інверсній фільтр.
У Функції J=deconvwnr (I, PSF, NSR) параметр NSR вказує на співвідношення сигнал/шум. Величина NSR є скаляром або масивом, розмірність которого така ж як І. За замовчуванням це значення дорівнює 0.
Найпростішій синтаксис Функції deconvwnr запісується таким чином:
(А, PSF, NSR),
де А - змінна, під Якою Збереже Розма зображення, PSF - точки Розповсюдження Функції, и NSR - величина співвідношення потужності шуму до потужності сигналу. Оскількі данє Розма зображення формується у 2 кроки без создания шуму, тому ми будемо використовуват 0 для NSR.
Синтаксис:=deconvwnr (blurred, PSF, 0); % Відновлення (wnr1); % Виведення на екран ( Restored Image ); % заголовок «Restored Image»
Рис. 9.2. Відновлене зображення
Лістінг
all=imread ( C: Documents and Settings Всеволод Робочий стіл Picture4.jpg ); % завантаження зображення; % Виведення на екран
% 2.1.1 геометричні превращение збережений.=imread ( C: Documents and Settings Всеволод Робочий стіл Picture4.jpg ); % завантаження зображення (L); % Виведення на екран; % к...
