2D обратное дискретное косинусное преобразование
B = idct2(A)B = idct2(A,m,n)B = idct2(A,[m n])Этот пример показывает, как удалить высокие частоты из изображения с помощью двумерной дискретной передачи косинуса (DCT).
Считайте изображение в рабочую область, затем преобразуйте изображение в шкалу полутонов.
RGB = imread('autumn.tif');
I = rgb2gray(RGB);Выполните 2D DCT полутонового изображения с помощью функции dct2.
J = dct2(I);
Отобразите преобразованное изображение с помощью логарифмического масштаба. Заметьте, что большая часть энергии находится в левом верхнем углу.
figure imshow(log(abs(J)),[]) colormap(gca,jet(64)) colorbar
Установите значения меньше, чем значение 10 в матрице DCT, чтобы обнулить.
J(abs(J) < 10) = 0;
Восстановите изображение с помощью обратной функции DCT idct2.
K = idct2(J);
Отобразите исходное полутоновое изображение вместе с обработанным изображением.
figure imshowpair(I,K,'montage') 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.