clear;
Irgb_a = imread('1.jpg'); %读入图片
[m, n, l] = size(Irgb_a);
Igray_a = rgb2gray(Irgb_a); %转化成灰度图
Igfft_a = fft2(Igray_a); %对灰度图进行dct变换
I1 = im2double(Irgb_a); %将数据转换成双精度格式
for i = 1:m
for j = 1:n
for k = 1:l
I(i,j+(k-1)*n) = I1(i,j,k);
end
end
end
Irfft = fft2(I);
for i = 1:m
for j = 1:n
for k = 1:l
Irifft_a(i,j,k) = I(i,j+(k-1)*n);
end
end
end
%imshow(Igray_a);
%title('原始灰度图');
%figure, imshow(log(abs(Igfft_a)),[]), colormap(jet(64)),colorbar;
%title('空间频谱图');
%figure,imshow(Irgb_a);
%title('原始彩图');
%figure,imshow(Irifft_a);
%title('复原彩图');
Irgb_b = imread('2.jpg'); %读入图片