2-D обратное дискретное косинусное преобразование
Прочитайте изображение в рабочей области, а затем преобразуйте его в градации серого.
RGB = imread('autumn.tif');
I = im2gray(RGB);Выполните 2-D DCT изображения в градациях серого с помощью dct2 функция.
J = dct2(I);
Отображение преобразованного изображения с использованием логарифмического масштаба. Обратите внимание, что большая часть энергии находится в левом верхнем углу.
imshow(log(abs(J)),[])
colormap parula
colorbar
Установите значения меньше величины 10 в DCT матрице в нуль.
J(abs(J) < 10) = 0;
Реконструировать изображение с помощью функции обратного DCT idct2. Масштабировать значения до диапазона [0, 1], ожидаемого для изображений типа данных double.
K = idct2(J); K = rescale(K);
Отображение исходного изображения в градациях серого рядом с обработанным изображением. Обработанное изображение имеет меньше высокочастотных деталей, таких как текстура деревьев.
montage({I,K})
title('Original Grayscale Image (Left) and Processed Image (Right)');
Для любой матрицы A, idct2(dct2(A)) равняется A в пределах ошибки округления.
idct2 вычисляет двумерный обратный DCT, используя:
0≤m≤M−10≤n≤N−1,
где
и
[1] Джейн, А. К., Основы обработки цифровых изображений, Энглвуд Клиффс, Нью-Джерси, Прентис Холл, 1989, стр. 150-153.
[2] Пеннебакер, В. Б., и Дж. Л. Митчелл, JPEG: Still Image Data Compression Standard, Нью-Йорк, Ван Ностранд Рейнхольд, 1993.