数字图形处理实验频域滤波.docVIP

XXXXXXXX大学（数字图形处理）实验报告 实验名称 图像的频域处理 实验时间 年 月 日 专 业 姓 名 学 号 预 习 操 作 座 位 号 教师签名 总 评 一、实验目的： 二、实验原理： 1、快速离散傅里叶变换 对于一个有限长序列 ，它的傅里叶变换由下式表示: 令 因此，傅里叶变换对可写成下式 2、低通滤波 3、高通滤波 1）理想高通 2）巴特沃斯高通 3）指数高通 三、实验内容： void CDipView::OnFft() { // TODO: Add your command handler code here long w,h; long i,j; w=m_DibHead-biWidth; h=m_DibHead-biHeight; complexdouble *TD=new complexdouble[w*h]; complexdouble *FD=new complexdouble[w*h]; for(i=0;ih;i++) for(j=0;jw;j++) { int flag=-1; if((i+j)%2==0) flag=1; TD[i*w+j]=complexdouble(m_Image[i*w+j]*flag,0); } zjbh bh; bh.DFFTC(FD,TD,w,h); double temp,r1,i1; for(i=0;ih;i++) for(j=0;jw;j++) { r1=FD[i*w+j].real(); i1=FD[i*w+j].imag(); temp=sqrt(r1*r1+i1*i1)/w; if(temp=255) m_Image[i*w+j]=255; else m_Image[i*w+j]=(unsigned char)(temp); } Invalidate(); delete []TD; delete []FD; } 低通滤波的代码如下： void CDipView::OnDtlb() { // TODO: Add your command handler code here long w,h; long i,j; w=m_DibHead-biWidth; h=m_DibHead-biHeight; complexdouble *TD=new complexdouble[w*h]; complexdouble *FD=new complexdouble[w*h]; for(i=0;ih;i++) for(j=0;jw;j++) { int flag=-1; if((i+j)%2==0) flag=1; TD[i*w+j]=complexdouble(m_Image[i*w+j]*flag,0); } zjbh bh; bh.DFFTC(FD,TD,w,h); for(i=0;ih;i++) for(j=0;jw;j++) { double d=sqrt((i-h/2)*(i-h/2)+(j-w/2)*(j-w/2)); if(d20) FD[i*w+j]=complexdouble(0,0); } bh.DIFFT(TD,FD,w,h); for(i=0;ih;i++) for(j=0;jw;j++) { m_Image[i*w+j]=abs(TD[i*w+j].real()); } Invalidate(); delete []TD; delete []FD; } 高通滤波的代码如下： void CDipView::OnGtlb() { // TODO: Add your command handler code here long w,h; long i,