matlab做频谱分析时如何实现图像去噪？

旧言虐心

如何用matlab实现对图像的频谱分析判断噪声的类别，然后根据对应的去噪方法实现去噪？  先不考虑多噪声的情况。

sd

拿图像做个二维的FFT或DCT，把高频部分的幅值置0，再反变换回来，就是个简单的例子。matlab里Demo有例子。

怪咖先生

《基于MATLAB的频谱分析与信号去噪》参考一下

zghitboy

1. I=imread('luo.jpg');
   
2. figure;imshow(I);title('原图像');
   
3. I1=imnoise(I,'salt & pepper',0.02);
   
4. figure;imshow(abs(I1));title('加入椒盐噪声后的图像');
   
5. I2=imnoise(I,'gaussian',0,0.005);
   
6. figure;imshow(I2);title('加入高斯噪声后的图像');
   
7. %对噪声污染的图像做均值滤波：
   
8. %---------对椒盐噪声做均值滤波 ---------%
   
9. II1=rgb2gray(I1);
   
10. a=1/9.*[1   1    1
    
11.          1   1    1
    
12.          1   1    1];
    
13. b=conv2(a,II1);
    
14. figure;imshow(b,[0 255]);title('对椒盐噪声的均值滤波图像');
    
15. %---------对高斯噪声做均值滤波---------%
    
16. II2=rgb2gray(I2);
    
17. a=1/9.*[1   1    1
    
18.          1   1    1
    
19.          1   1    1];
    
20. b=conv2(a,II2);
    
21. figure;imshow(b,[0 255]);title('对高斯噪声的均值滤波图像');  
    
22. h=fspecial('average',3)/255;
    
23. I3=filter2(h,II1);
    
24. figure;imshow(I3);
    
25. title('对高斯噪声的均值滤波后的图像');
    
26. 
    
27. %对噪声污染的图像做中值滤波：
    
28. %---------对椒盐噪声做中值滤波---------%
    
29. c=medfilt2(II1,[3 3]);
    
30. figure;imshow(c);title('对椒盐噪声的中值滤波图像');
    
31. %---------对高斯噪声做中值滤波---------%
    
32. c=medfilt2(II2,[3 3]);
    
33. figure;imshow(c);title('对高斯噪声的中值滤波图像');
    
34. 
    
35. B=rgb2gray(I1);  
    
36. C=B;
    
37. xsize=size(B);
    
38. for k=2:(xsize(1)-1)
    
39.     for j=2:(xsize(2)-1)
    
40.         t=B(k-1:k+1,j-1:j+1);
    
41.         C(k,j)=median(t(1:9));
    
42.     end;
    
43. end;
    
44. imshow(C);
    
45. title('中值处理后的图');
    
46. 
    
47. %采用高斯滤波
    
48. a=1/1151.*[1  4  7  10 7  4   1   
    
49.          4  12 26 33 26 12  4
    
50.          7  26 55 71 55 26  7
    
51.          10 33 71 91 71 33 10
    
52.          7  26 55 71 55 26  7
    
53.          4  12 26 33 26 12  4
    
54.          1  4  7  10 7  4   1 ];
    
55. b=conv2(a,II2);
    
56. figure;imshow(b,[0 255]);title('对椒盐噪声的高斯滤波图像');
    
57. h=fspecial('gaussian',[3 3],0.5);
    
58. I4=imfilter(I3,h);
    
59. figure;imshow(I4);
    
60. title('高斯滤波后的图像');

复制代码

dh492510085

学习了