Matlab进行快速傅里叶变换和反变换

gaohaikuo1

1. % 读取图片
   
2. image = imread('image.jpg');
   
3. 
   
4. % 将图像转换为灰度图像（如果需要）
   
5. grayImage = rgb2gray(image);
   
6. 
   
7. % 进行二维快速傅里叶变换
   
8. fftImage = fft2(grayImage);
   
9. 
   
10. % 频谱中心化
    
11. fftImageShifted = fftshift(fftImage);
    
12. 
    
13. % 进行反变换
    
14. ifftImageShifted = ifftshift(fftImageShifted);
    
15. ifftImage = ifft2(ifftImageShifted);
    
16. 
    
17. % 显示原始图像、频谱图和反变换后的图像
    
18. figure;
    
19. subplot(1, 3, 1), imshow(grayImage), title('原始图像');
    
20. subplot(1, 3, 2), imshow(log(1 + abs(fftImageShifted)), []), title('中心化频谱图');
    
21. subplot(1, 3, 3), imshow(abs(ifftImage), []), title('反变换图像');

复制代码

gaohaikuo1

1. % 读取图片
   
2. image = imread('image.png');
   
3. image1 = image(50:280,50:300,:)
   
4. image2 = image(50:280,400:650,:)
   
5. % 将图像转换为灰度图像（如果需要）
   
6. grayImage1 = rgb2gray(image1);
   
7. grayImage2 = rgb2gray(image2);
   
8. % 进行二维快速傅里叶变换
   
9. fftImage1 = fft2(grayImage1);
   
10. fftImage2 = fft2(grayImage2);
    
11. % 频谱中心化
    
12. fftImageShifted1 = fftshift(fftImage1);
    
13. fftImageShifted2 = fftshift(fftImage2);
    
14. 
    
15. % 进行反变换
    
16. ifftImageShifted1 = ifftshift(fftImageShifted1);
    
17. ifftImage1 = ifft2(ifftImageShifted1);
    
18. ifftImageShifted2 = ifftshift(fftImageShifted2);
    
19. ifftImage2 = ifft2(ifftImageShifted2);
    
20. 
    
21. % 显示原始图像、频谱图和反变换后的图像
    
22. figure;
    
23. subplot(2, 3, 1), imshow(grayImage1), title('原始图像');
    
24. subplot(2, 3, 2), imshow(log(1 + abs(fftImageShifted1)), []), title('中心化频谱图');
    
25. subplot(2, 3, 3), imshow(abs(ifftImage1), []), title('反变换图像');
    
26. subplot(2, 3, 4), imshow(grayImage2), title('原始图像');
    
27. subplot(2, 3, 5), imshow(log(1 + abs(fftImageShifted2)), []), title('中心化频谱图');
    
28. subplot(2, 3, 6), imshow(abs(ifftImage2), []), title('反变换图像');

复制代码

gaohaikuo1

gaohaikuo1

1. % 读取图片
   
2. image = imread('aaa.jpg');
   
3. 
   
4. % 将图像转换为灰度图像（如果需要）
   
5. grayImage = rgb2gray(image);
   
6. 
   
7. % 进行二维快速傅里叶变换
   
8. fftImage = fft2(grayImage);
   
9. 
   
10. % 频谱中心化
    
11. fftImageShifted = fftshift(fftImage);
    
12. %加入噪声点
    
13. 
    
14. fftImageShifted(220:222,330:331)=2e6;
    
15. %fftImageShifted(180:222,428:end)=0;
    
16. % 进行反变换
    
17. ifftImageShifted = ifftshift(fftImageShifted);
    
18. ifftImage = ifft2(ifftImageShifted);
    
19. 
    
20. % 显示原始图像、频谱图和反变换后的图像
    
21. figure;
    
22. subplot(1, 3, 1), imshow(grayImage), title('原始图像');
    
23. %imwrite(grayImage,'image_noise.jpg');
    
24. subplot(1, 3, 2), imshow(log(1 + abs(fftImageShifted)), []), title('中心化频谱图');
    
25. img =double(abs(ifftImage));
    
26. img = uint8(img/max(max(img))*255);
    
27. subplot(1, 3, 3), imshow(img), title('反变换图像');
    
28. imwrite(img,'image_noise.jpg');

复制代码