
- •Методичні вказівки до вивчення дисципліни
- •«Цифрова обробка сигналів та зображень»
- •Для студентів заочної форми навчання
- •За напрямом 6.050102 «Комп’ютерна інженерія»
- •1 Зміст навчальної дисципліни
- •2 Методичні вказівки до виконання розрахунковографічної роботи та варіанти завдань
- •2.1. Можливості цифрової обробки зображень в Matlab
- •Типи зображень
- •2.1.1. Геометричні перетворення зображень.
- •2.1.2. Аналіз зображень.
- •2.1.3. Поліпшення зображень.
- •2.1.4. Фільтрація зображень.
- •2.1.5. Сегментація зображень.
- •2.2. Алгоритми усунення розмитості цифрових зображень
- •2.2.1 Усунення розмитості зображення за допомогою алгоритму сліпої деконволюції.
- •Зчитування зображення.
- •2.2.2. Моделювання Blur.
- •2.2.3.Відновлення розмитого зображення.
- •2.2.4. Усунення розмитості зображення з допомогою фільтру Вінера
- •2.2.5. Імітація (Motion Blur) розмитості рухом.
- •2.2.6. Відновлення розмитого зображення.
- •3 Методика виконання розрахунково-графічної роботи
- •3.1 Цілі та завдання ргр
- •3.2. Порядок виконання ргр
- •4 Запитання для підготовки до екзамену
- •Рекомендована література
- •Приклад оформлення розрахунковографічної роботи
2.1.2. Аналіз зображень.
Для роботи з окремими елементами зображень використовуються такі функції як imhist, impixel, mean2, corr2 та інші.
Однією з найбільш важливих характеристик є гістограма розподілу значень інтенсивностей пікселів зображення, яку можна побудувати за допомогою функції imhist (рис. 4.).
Синтаксис:
L=imread('cameraman.tif');
figure, imshow(L);
figure, imhist(L);
Рис. 4. Оригінальне зображення та його гістограма
Досить часто при проведенні аналізу зображень виникає необхідність визначити значення інтенсивностей деяких пікселів. Для цього необхідно в інтерактивному режимі використовувати функцію impixel.
Синтаксис:
>> impixel
ans =
170 170 170
164 164 164
158 158 158
70 70 70
12 12 12
Рис. 5. Зображення із позначеними на ньому крапками інтенсивності
2.1.3. Поліпшення зображень.
Серед вбудованих функцій, які реалізуютья найбільш відомими методами поліпшення зображень, розглянемо наступні – histeq, imadjust та imfilter (fspecial).
Як вже зазначалося раніше, гістограма зображення є однією з найбільш інформативних характеристик. На основі аналізу гістограми можна судити про яскравості зображення спотворене, тобто сказати про те, чи є зображення затемненим або засвітленим. Відомо, що в ідеальному випадку на цифровому зображенні в рівній кількості повинні бути присутніми пікселі з усіма значеннями яскравостей, тобто гістограма повинна бути рівномірною. Перерозподіл яскравості пікселів на зображенні з метою отримання рівномірної гістограми виконує метод еквалізації, який в системі Matlab реалізовано за допомогою функції histeq (вирівнювання гістограми).
Синтаксис:
L=imread('cameraman.tif');
figure, imshow(L);
L1=histeq(L);
figure, imshow(L1);
Рис. 6. Перерозподіл яскравостей
Досить часто при формуванні зображень не використовується весь діапазон значень інтенсивностей, що негативно відбивається на якості візуальних даних. Для корекції динамічного діапазону сформованих зображень використовується функція imadjust (контрастування з гамма-корегуванням).
Синтаксис:
L=imread('cameraman.tif');
figure, imshow(L);
L1=imadjust(L);
figure, imshow(L1);
figure, imhist(L);
figure, imhist(L1)
Рис. 7. Гама-корегування за гістограмою
Також під час вирішення завдань поліпшення зображень використовується функція imfilter в парі з функцією fspecial. Функція fspecial дозволяє задавати різні типи масок фільтру. Розглянемо приклад використання маски фільтру, що підвищує різкість зображення.
Рис. 8. Приклад використання маски наперед зазначеного фільтра
2.1.4. Фільтрація зображень.
Пакет Image Processing Toolbox володіє потужним інструментарієм з фільтрації зображень. Серед вбудованих функцій, які вирішують задачі фільтрації зображень, особливу увагу необхідно приділити функціям fspecial, ordfilt2 (рангова фільтрація), medfilt2 (медіанна фільтрація).
Fspecial є функцією завдання маски зумовленого фільтра. Ця функція дозволяє формувати маски:
високочастотного фільтру Лапласа
h = fspecial('laplacian', alpha);
фільтру, аналогічного послідовному застосуванню фільтрів Гауса і Лапласа, так званого лапласіана-гауссіана
h = fspecial('log', hsize, sigma);
усереднюючий низькочастотний фільтр
h = fspecial('average', hsize);
Розглянемо приклади застосування названих вище фільтрів:
Синтаксис:
L=imread('cameraman.tif');
figure, imshow(L);
h=fspecial('laplacian',0.5);
L1 = imfilter(L,h,'replicate');
figure, imshow(L1);
Рис. 9. Наперед зазначений високочастотний фільтр Лапласа
Синтаксис:
h=fspecial('log', 3, 0.5);
L1 = imfilter(L,h,'replicate');
figure, imshow(L1);
Рис. 10. Наперед зазначений фільтр лапласіана-гауссіана
Синтаксис:
h=fspecial('average', 3);
L1 = imfilter(L,h,'replicate');
figure, imshow(L1);
Рис. 11. Наперед зазначений низькочастотний фільтр