二维福利叶变换图像处理技术.doc
《二维福利叶变换图像处理技术.doc》由会员分享,可在线阅读,更多相关《二维福利叶变换图像处理技术.doc(20页珍藏版)》请在咨信网上搜索。
1、第一章 基于二维傅里叶变换的图像处理技术1.1二维傅里叶变化基础 假如以正方形网格采样得到的数字图像用f(x,y)表示,由一维傅里叶变换的定义可以得到f(x,y)的二维傅里叶变换为: 其反变换为:1.2傅里叶变换的性质1 分离性:2 平移性3 周期性4 共轭对称性5 旋转性6 分配性和比例性1.3 matlab中二维傅里叶变换的函数以及应用 在MATLAB中,用来实现数字图像傅里叶变换用到的函数主要有fft2函数和fftshift函数。前者用来实现数字信号的二维离散傅里叶变换,后者用来将傅里叶变换的零频率部分移到频谱的中间。fft2函数的语法格式如下:Y = fft2(X)Y = fft2(X
2、,m,n)fftshift函数的语法格式如下:Y = fftshift(X)Y = fftshift(X,dim)例:制作64*64 大小的简单黑白二值图像用Matlab 的imread 函数读取并做二维DFT 变换代码:clear%创建256*256的二值图像f = zeros(256,256);f(124:130,117:137) = 1;%原始图像显示figure(1);imshow(f);title(原始图像);%傅里叶变换F = fft2(f);%频谱中心化F2 = fftshift(abs(F);%结果显示figure(2);x=1:256;y=1:256;mesh(x,y,F2(
3、x,y); colormap(jet); colorbartitle(傅立叶变换结果);图 1-1 a) 图1-1 b)结果分析: 从变换结果图1-1 b)可以看出,图像的能量主要集中在某一频率段,如果将系数较小的部分舍去,即可实现图像的压缩。例:更改DFT 系数为整数,做IDFT 观察图像的变化; %将DFT系数的绝对值四舍五入后作为新的dft系数F1=round(abs(F);f2=ifft2(F1);figure(3);imshow(f2);实验结果:图1-2 a)实验结论: 可见图像无法还原到最初。例:更改小幅值的DFT 系数为0,做IDFT 变化观察图像的变化:%将绝对值小于a的都改
4、为零,在做IDFT变换:%取a=1;F(abs(F)1)=0;f2=ifft2(F);figure(4);imshow(f2);实验结果:图1-3 a)结果分析: 在执行二维离散傅里叶变换过程中,只取离散傅里叶变换系数大于1的系数,所示的反变换结果如图1-3 c)和原图像1-1 a)比较,通过肉眼观测,很难看出两者的区别,如果将系数小于100的值设为0,反变换得到的结果如图1-3 d)所示,图像的清晰度越来越低,图像质量越来越差。F(abs(F)100)=0;f2=ifft2(F);figure(4);imshow(f2);图1-3 d)例:选取不同类型(自然景观、人物照片、卡通图片)的实际图
5、片重复上述处理并分析结果用下面这幅图进行分析: 代码:f=imread(e:/26794054.jpg);figure(1);imshow(f);F=fft2(f);figure(2);f1=ifft2(F);imshow(f1);图1-4结果分析: 可见图像已经发生严重失真。例:将彩色图像转换为黑白二值图像之后才能进行二维傅里叶变换这里我们取阈值为0.5代码:f=imread(e:/26794054.jpg);f1=im2bw(f,0.5);figure(1);imshow(f1);title(原始图像);F=fft2(f1);%取F的实部再取整F1=ceil(real(F);f2=ifft
6、2(F1);figure(2);imshow(f2);%将DFT系数小于十的全取零F(abs(F)10)=0;f3=ifft2(F);figure(3);imshow(f3);实验结果:图1-5 a) 图1-5 b) 图1-5 c)结果分析:可见如果要压缩彩图,不可能保留其颜色信息,只有将其先转换为黑白二值图像之后才能进行压缩,这样至少图像基本形状信息不会缺失。第2章 基于离散余弦变换的图像压缩技术2.1离散余弦变换基础一维离散余弦变换和其反变换定义如下:式中二维离散余弦变换和反变换的公式为:2.2 Matlab实例操作1、dct2函数 利用该函数可以实现数字图像的二维离散余弦变换,该函数使用
7、了一个基于FFT的算法,提高了较大矩阵时的处理速度。其语法格式为:B = dct2(A)B = dct2(A,m,n)B = dct2(A,m n)例:对D:matlabgirl.jpg进行离散余弦变换:代码:clear all%读入图像并显示RGB = imread(D:matlabgirl.jpg); figure(1);imshow(RGB);title(原始图像) ;%将真彩图像转换为灰度图像A = rgb2gray(RGB); %实现离散余弦变换B = dct2(A); %结果显示figure(2);imshow(log(abs(B), ), colormap(jet(64), co
8、lorbartitle(离散余弦变换结果);实验结果:图2-1 a) 图2-1 b)图2-1 c)结果分析: 从图 2-1 b)所示的变化结果及图2-1 c)所示的变化系数可以看出,经过离散余弦变换,图像的能量主要集中在图像2-1 b)的左上角,其余大部分系数接近于0,这说明离散余弦变换也可用于图像的压缩。2. 、idct2函数MATLAB 中,可以用函数Idct2实现离散余弦反变换,其语法格式为:B = idct2(A)B = idct2(A,m,n)B = idct2(A,m n)例:对D:matlabgirl.jpg进行离散余弦变换,在进行反变换之前,将系数小于10的值设为0,再进行反变
9、换代码: clear all%读入图像并转换为灰度图像RGB = imread(D:matlabgirl.jpg);A = rgb2gray(RGB); %显示图像figure(1);imshow(A);title(原始图像灰度显示) ;%离散余弦变换并显示结果B = dct2(A); figure(2);imshow(log(abs(B),), colormap(jet(64), colorbar title(离散余弦变换结果图像);%将变换系数中小于10的值设为0B(abs(B) P1; 8 8-x ;D-P2Bdct = blkproc(I,8 8,P1*x*P2,D,D); % 系数选
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 二维 福利 变换 图像 处理 技术
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【a199****6536】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【a199****6536】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。