傅里叶变换在数字图像讲述.docxVIP

数学文化课程报告傅里叶变换在数字图像处理中的应用成绩：傅里叶变换在数字图像处理中的应用朱德辉中国石油大学华东地球科学与技术学院山东青岛 266580摘要：文章将“傅里叶变换”这种计算工具与图像分析有机结合起来,将“傅里叶变换”应用到空间频谱处理当中,并将结果以最直观的方式—图形表示出来;文章意在说明用一种实用的方法对数字图像进行处理。关键词:傅里叶变换数字图像处理图像滤波图像压缩图像恢复Application of Fourier Transform in Digital Image ProcessZhuDehuiChina University of Petroleum Global Scientific and Technological ShanDong Qingdao 266580Abstract: This article states the organic combination of the formula named Fourier transform and the analysis of image .The Fourier transform can be used in dealing with the 3-dimen-tion-space frequency curve chart ;meanwhile , it can also show us directly by image . The article is meant to seek a practical way being used in digital image process .Keywords: Fourier transform , digital image processing , image filter , image compression , image restoration 引言近年来,随着电子技术、图像处理方法和信号理论的迅猛发展, 数字图像处理技术得到飞速发展, 它广泛应用于几乎所有与成像有关的领域。传统的光学系统在信号处理时 ,存有它自身很难克服的不足 : 第一 ,它对空间频谱平面的处理很难,尤其在低频和甚低频时 ,即使可通过大量仪器来实现 ,但代价往往很高 ;第二 ,光学处理由于采样孔径(即传感单元)太窄而不能起到抗混叠作用 ,不能除去高频信息。而傅里叶变换和线性移不变系统有紧密联系, 它有一个很好的理论背景来指导它在图像处理中的作用 , 可以方便有效地克服上述不足 , 使其在数字图像处理中占有一席之地。方法与原理本文以计算机为工具,以C语言编程为手段,通过快速离散傅里叶变换与C语言利用傅里叶变换对空间频谱平面进行处理。在本文中,我们所用到的有关公式和原理有：二维离散傅里叶正变换为：二维离散傅里叶逆变换：其中,h(k1,k2)是一个定义在二维网格0≤k1≤N1-1 ,0 ≤k2≤N2-1 上的复函数(就像用等间距的矩形网格对一个二维连续函数采样所得的数组),H(n1 ,n2)是二维离散傅里叶变换后的复函数, N1、N2为数组的有效界长。本程序用基-2频域抽取法蝶形运算,有关此运算详细资料请参阅《数字信号处理》一书中“离散傅里叶变换”一章。程序中所用主要子程序简介:①init():将图形驱动软件装入内存,使屏幕适配器设置为图形方式,即图形初始化;②make_original_data():置数程序,以取得有关二维数组,进行傅里叶变换;③do_2D_DFT():完成二维离散傅里叶正变换,通过绘图程序将运算结果显示在屏幕上;④do_2D_IDFT():完成二维离散傅里叶逆变换,通过绘图程序将运算结果显示在屏幕上,使其与原图形相比较,以检验变换的恢复情况;⑤filter():在这个子程中,我们可以通过选择有效对空间频谱进行滤波,滤去不需要谱线,保留有效谱线,突出重要部分,削弱次要因素,达到对空间频谱处理目的。程序设计图一程序流程图主要子程序:do_2D_DFT(){ getch();fft2(xx ,yy , 1); disp_fft(200 , 360 , 6, 0 .03);gotoxy(30 ,25);printf(” 2D-FFT”);}{int imin , imax, jmin , jmax , fg=1 ; char c ;outtextxy(100, 430 ,”filter? (Y/ N):”); c=getch();delline();if(c == yc== Y ){while(fg) {gotoxy(60 , 21);printf(” 0 ==%d :” ,(size 1)); gotoxy(60 , 22);printf(”input imin :”);scanf(” %d”