2D обратное дискретное косинусное преобразование
Считайте изображение в рабочую область, затем преобразуйте изображение в шкалу полутонов.
RGB = imread('autumn.tif');
I = im2gray(RGB);Выполните 2D 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] Джайн, А. К., Основные принципы Цифровой обработки изображений, Englewood Cliffs, NJ, Prentice Hall, 1989, стр 150-153.
[2] Pennebaker, W. B., и Дж. Л. Митчелл, JPEG: Стилл-Имэдж-Дэта-Компрешен-Стэндард, Нью-Йорк, Ван Нострэнд Райнхольд, 1993.