cvMan
工程师
工程师
  • 最后登录2016-12-14
  • 发帖数91
  • 金币1093枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
阅读:3251回复:10

傅里叶变换学习笔记(冈萨雷斯数字图像处理)

楼主#
更多 发布于:2014-03-10 15:06
1,二维离散傅里叶变换公式的理解,DFT公式中,我们可以将指数项扩展为正弦项和余弦项的形式,
所以从公式的角度可以理解F(0,0)即傅里叶变换的直流分量是最大值,这也就是为什么居中二维傅里叶变换频谱图都是中心较亮。(图a为原始图像,图b为对数变换后增强的DFT图像)

图片:2.bmp


                                                         (图 a)

图片:2-fft.bmp

                                                             (图 b)
2,图像的平移对傅里叶频谱是没有影响的,那是因为离散傅里叶变换的计算过程是累积计算的过程,在周期不变的情况下,对于图像的平移是不会影响最后累积计算的结果的。(图c是对图a的平移)

图片:1-offset-equal2.bmp

                                                       (图 c)
       

图片:2-fft.bmp

                                                          (图 d)
3.对频域滤波器的理解,低通滤波(包含理想低通滤波、n阶巴特沃兹低通滤波、高斯低通滤波等),高通滤波(包含包含理想高通滤波、n阶巴特沃兹高通滤波、高斯高通滤波等),低通滤波是保持傅里叶变换的低频,削弱高频会使图像变得模糊,相反的高通滤波是保持傅里叶变换的高频,削弱低频,会使图像变得清晰。但从滤波器的定义公式可以简单理解设置的截止频率D0的作用,就是将大于截止频率或低于截止频率的部分*0或*1来实现的,但要了解为什么会消弱了低频,保持高频就能将图像变得清晰,还要理解从原始图像到傅里叶变换的过程,简单地讲一下傅里叶变换的原理,对于一维傅里叶变换,就是将一个一维空间的函数曲线分解出一个周期内的不同的正弦、余弦曲线,而这些曲线有不同的频率和幅度,以频率为横轴,幅度为纵轴,就可以得到这个一维空间函数对应的频谱图,扩展为二维空间,就是两个频率方向(U,V)幅度的叠加而成的频谱图。理解了这些,我们可以知道对于图像中的平坦部分,就是梯度值很小的部分所产生的曲线为接近于DC的直流低频部分,而对于图像的轮廓附近即梯度较大的部分就会产生频率较高的高频部分,也就可以理解为什么保持高频,削弱低频图像变清晰,削弱高频,保持低频图像变模糊了。为了更明白这个道理,我还将低通、高通滤波后的图像再做一次DFT变换,观察频谱图中低频及高频部分的变换是与前面讲到的逻辑是一致的。
下面3幅图依次为原始图像、高通滤波后效果图、低通滤波后效果图
               

图片:origin.JPG

图片:hp.jpg

图片:lp.jpg

下面3幅图依次为原始图、高通滤波、低通滤波图对应的频谱图

图片:originFFT0.jpg

图片:hp-fft.jpg

图片:lp-fft.jpg

4.Matlab 中频域滤波dftfilt步骤:
step1.使用函数paddedsize获得填充参数;
step2.得到使用填充的傅里叶变换
step3.生成一个大小为PQ(1)*PQ(2)的滤波函数H(注意滤波与变换的居中fftshit,要同步)
step4.将变换乘以滤波函数
step5.获得傅里逆变换的实部
step6.将左上部矩形修剪为原始大小
[cvMan于2014-03-20 17:04编辑了帖子]
cvMan
工程师
工程师
  • 最后登录2016-12-14
  • 发帖数91
  • 金币1093枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
沙发#
发布于:2014-03-10 16:01
后面怎么会多三张图?
danruochenxi
工程师
工程师
  • 最后登录2016-08-20
  • 发帖数138
  • 金币799枚
  • 贡献勋章0枚
  • 社区居民
板凳#
发布于:2014-03-11 11:43
cvMan:后面怎么会多三张图?回到原帖
后面的三张图片应该是验证带通滤波器的图片
cvMan
工程师
工程师
  • 最后登录2016-12-14
  • 发帖数91
  • 金币1093枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
地板#
发布于:2014-03-11 12:17
danruochenxi:后面的三张图片应该是验证带通滤波器的图片回到原帖
这个是我自己写贴的时候,粘图误操作,后面的三张图片没有用,帖子内容应该到4.Matlab 中频域滤波dftfilt6个步骤后面结束
cvMan
工程师
工程师
  • 最后登录2016-12-14
  • 发帖数91
  • 金币1093枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
4楼#
发布于:2014-03-11 12:20
cvMan:这个是我自己写贴的时候,粘图误操作,后面的三张图片没有用,帖子内容应该到4.Matlab 中频域滤波dftfilt6个步骤后面结束回到原帖
哦,后来编辑的时候删不了了,哈哈
st_leger
初学者
初学者
  • 最后登录2014-06-19
  • 发帖数25
  • 金币68枚
  • 贡献勋章0枚
  • 社区居民
5楼#
发布于:2014-03-11 15:14
日后慢慢学习  
cvMan
工程师
工程师
  • 最后登录2016-12-14
  • 发帖数91
  • 金币1093枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
6楼#
发布于:2014-03-11 16:10
st_leger:日后慢慢学习回到原帖
自己的一点心得,表达得可能不是很好,一起学习吧!哈哈哈
nian16
高级工程师
高级工程师
  • 最后登录2016-11-05
  • 发帖数158
  • 金币1935枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
7楼#
发布于:2014-03-14 09:50
好好学习吧~~~~~~~~~~
zkdhu1989
技术员
技术员
  • 最后登录2015-01-13
  • 发帖数24
  • 金币179枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
8楼#
发布于:2014-06-12 12:41
最近做的东西也涉及到fft  学习一下
manchao91
工程师
工程师
  • 最后登录2016-05-09
  • 发帖数22
  • 金币822枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
9楼#
发布于:2014-06-12 13:05
可得好好学习学习了~~~
MVision007
助理工程师
助理工程师
  • 最后登录2016-12-20
  • 发帖数181
  • 金币523枚
  • 贡献勋章0枚
  • 社区居民
  • 忠实会员
10楼#
发布于:2016-12-04 20:34
傅里叶变换学习笔记(冈萨雷斯数字图像处理)
游客

返回顶部