Я реализовал свой собственный 2D DCT, используя 1D DCT и IDCT. Мои результаты DCT совпадают с реализацией Matlab, но IDCT дает другой результат. Разница не совсем исчезла, как видно на реконструированном изображении.
Проблема восстановления изображения с использованием обратного DCT в Matlab
Ответы (1)
Я попробовал это - ваши уравнения DCT / IDCT мне показались не совсем правильными. Я использовал формулы DCT-2 и DCT-3 из документации SciPy здесь.
original_img = imread('nggyu.jpeg');
transformed_img = permute(dct1d(permute( ...
dct1d(double(original_img)), ...
[2,1,3])), [2,1,3]);
recovered_img = uint8(permute(idct1d(permute( ...
idct1d(transformed_img), ...
[2,1,3])), [2,1,3]));
figure('position', [0, 0, 600, 200])
subplot(1,3,1), imshow(original_img), title 'Original'
subplot(1,3,2), imshow(log(abs(transformed_img)),[]), title 'DCT'
subplot(1,3,3), imshow(recovered_img), title 'IDCT'
function y = dct1d(x)
% Compute normalized DCT-2 over the first dimension of the input.
N = size(x, 1);
y = zeros(size(x));
n = (1:N)';
for k = 1:N
if k == 1
scale = sqrt(1/(4*N));
else
scale = sqrt(1/(2*N));
end
y(k,:,:) = scale * 2 * sum(x(n,:,:) .* cos((pi/(2*N)) * (2*n-1) * (k-1)), 1);
end
end
function x = idct1d(y)
% Compute normalized DCT-3 over the first dimension of the input.
N = size(y, 1);
x = zeros(size(y));
k = (2:N)';
for n = 1:N
x(n,:,:) = y(1,:,:)/sqrt(N) + sqrt(2/N) * sum(y(k,:,:) .* cos((pi/(2*N)) * (2*n-1) * (k-1)), 1);
end
end
person
aosborne
schedule
27.02.2021
Чтобы быть немного более конкретным - проблема в вашем коде, которая привлекла мое внимание, заключалась в индексации уравнения DCT/IDCT. MATLAB имеет индекс 1, но DCT/IDCT обычно записываются с индексом 0. Ваш код должен учитывать это.
- person aosborne; 27.02.2021