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 с помощью:
где
и
[1] Jain, A. K., Fundamentals of Digital Image Processing, Englewood Cliffs, NJ, Prentice Hall, 1989, pp. 150-153.
[2] Pennebaker, W. B., and J. L. Mitchell, JPEG: Still Image Data Compression Standard, New York, Van Nostrand Reinhold, 1993.