Матрица дискретного косинусного преобразования
Чтение изображения в рабочую область и приведение его в класс double.
A = im2double(imread('rice.png'));
imshow(A)
Вычислите дискретную матрицу косинусного преобразования.
D = dctmtx(size(A,1));
Умножение входного изображения A около D чтобы получить DCT столбцов A, и D' для получения обратного DCT столбцов A.
dct = D*A*D'; imshow(dct)
Если у вас есть nоколо-n изображение, A, то D*AA и D'*A - обратный DCT столбцов A.
Двухмерный DCT A может быть вычислено как D*A*D'. Это вычисление иногда происходит быстрее, чем использование dct2, особенно если вы вычисляете большое количество небольших DCT, потому что D необходимо определить только один раз.
Например, при сжатии JPEG вычисляется DCT каждого блока 8 на 8. Для выполнения этого вычисления используйте dctmtx определить D, а затем вычислить каждый DCT с помощью D*A*D' (где A каждый блок 8 на 8). Это быстрее, чем вызов dct2 для каждого отдельного блока.