声振论坛

 找回密码
 我要加入

QQ登录

只需一步,快速开始

查看: 3057|回复: 4

[综合讨论] matlab做频谱分析时如何实现图像去噪?

[复制链接]
发表于 2017-2-17 09:11 | 显示全部楼层 |阅读模式

马上注册,结交更多好友,享用更多功能,让你轻松玩转社区。

您需要 登录 才可以下载或查看,没有账号?我要加入

x
如何用matlab实现对图像的频谱分析判断噪声的类别,然后根据对应的去噪方法实现去噪?  先不考虑多噪声的情况。
回复
分享到:

使用道具 举报

发表于 2017-2-17 09:11 | 显示全部楼层
拿图像做个二维的FFT或DCT,把高频部分的幅值置0,再反变换回来,就是个简单的例子。matlab里Demo有例子。
发表于 2017-2-20 09:00 | 显示全部楼层
《基于MATLAB的频谱分析与信号去噪》参考一下
发表于 2017-2-22 09:03 | 显示全部楼层
  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. c=medfilt2(II1,[3 3]);
  29. figure;imshow(c);title('对椒盐噪声的中值滤波图像');
  30. %---------对高斯噪声做中值滤波---------%
  31. c=medfilt2(II2,[3 3]);
  32. figure;imshow(c);title('对高斯噪声的中值滤波图像');

  33. B=rgb2gray(I1);  
  34. C=B;
  35. xsize=size(B);
  36. for k=2:(xsize(1)-1)
  37.     for j=2:(xsize(2)-1)
  38.         t=B(k-1:k+1,j-1:j+1);
  39.         C(k,j)=median(t(1:9));
  40.     end;
  41. end;
  42. imshow(C);
  43. title('中值处理后的图');

  44. %采用高斯滤波
  45. a=1/1151.*[1  4  7  10 7  4   1   
  46.          4  12 26 33 26 12  4
  47.          7  26 55 71 55 26  7
  48.          10 33 71 91 71 33 10
  49.          7  26 55 71 55 26  7
  50.          4  12 26 33 26 12  4
  51.          1  4  7  10 7  4   1 ];
  52. b=conv2(a,II2);
  53. figure;imshow(b,[0 255]);title('对椒盐噪声的高斯滤波图像');
  54. h=fspecial('gaussian',[3 3],0.5);
  55. I4=imfilter(I3,h);
  56. figure;imshow(I4);
  57. title('高斯滤波后的图像');
复制代码


发表于 2017-6-11 18:56 | 显示全部楼层
学习了
您需要登录后才可以回帖 登录 | 我要加入

本版积分规则

QQ|小黑屋|Archiver|手机版|联系我们|声振论坛

GMT+8, 2024-11-28 13:34 , Processed in 0.106632 second(s), 21 queries , Gzip On.

Powered by Discuz! X3.4

Copyright © 2001-2021, Tencent Cloud.

快速回复 返回顶部 返回列表