欢迎来到天天文库
浏览记录
ID:22767100
大小:722.81 KB
页数:18页
时间:2018-10-31
《北航数字图象处理实验报告》由会员上传分享,免费在线阅读,更多相关内容在工程资料-天天文库。
1、数字图像处理实验报告实验二图像变换实验1.实验目的学会对图像进行傅立叶等变换,在频谱上对图像进行分析,增进对图像频域上的感性认识,并用图像变换进行压缩。2.实验内容对Lena或cameraman图像进行傅立叶、离散余弦、哈达玛变换。在频域,对比他们的变换后系数矩阵的频谱情况,进一步,通过逆变换观察不同变换下的图像重建质量情况。3.实验要求实验采用获取的图像,为灰度图像,该图像每象素由8比特表示。具体要求如下:(1)输入图像采用实验1所获取的图像(Lena、Cameraman);(2)对图像进行傅立叶变换、获得变换后的系数矩阵;(3)将傅立叶变换后系数矩阵的频谱用图像输出,观察频谱;(4)
2、通过设定门限,将系数矩阵中95%的(小值)系数置为0,对图像进行反变换,获得逆变换后图像;(5)观察逆变换后图像质量,并比较原始图像与逆变后的峰值信噪比(PSNR)。(6)对输入图像进行离散余弦、哈达玛变换,重复步骤1-5;(7)比较三种变换的频谱情况、以及逆变换后图像的质量(PSNR)。4.实验结果1.DFT的源程序及结果J=imreadC10021033.bmp’);P=fft2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P,2)+j)=P(i+1,j);endQ=sort(G);fori=1:size(Q,2)if(i3、.2)*0.95&&i+1>=size(Q,2)*0.95)t=Q⑴;endendG(abs(G)4、2)对图像进行傅立叶变换、获得变换后的系数矩阵;(3)将傅立叶变换后系数矩阵的频谱用图像输出,观察频谱;(4)通过设定门限,将系数矩阵中95%的(小值)系数置为0,对图像进行反变换,获得逆变换后图像;(5)观察逆变换后图像质量,并比较原始图像与逆变后的峰值信噪比(PSNR)。(6)对输入图像进行离散余弦、哈达玛变换,重复步骤1-5;(7)比较三种变换的频谱情况、以及逆变换后图像的质量(PSNR)。4.实验结果1.DFT的源程序及结果J=imreadC10021033.bmp’);P=fft2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P,5、2)+j)=P(i+1,j);endQ=sort(G);fori=1:size(Q,2)if(i=size(Q,2)*0.95)t=Q⑴;endendG(abs(G)6、erexp2part实辁2PSNR32.4671尿囝交挨后复SS•耳里叶变挨DCT哈达玛变疾返回1.DCT的源程序及结果J=imreadC10021033.bmp’);P二dct2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P.2)+j)=P(i+1,j);endend0=sort(G);fori=1:size(Q,2)if(i=size(Q.2)*0.95)t=Q(i);endendG(abs(G)7、G(n*size(P,2)+m):endendf2=idct2(W);f3=uint8(f2);axes(handIes.axes2);imshow(f3)axes(handIes.axesl);imshow(J)psnr1=psnr(J,f3);set(handles.text3,’string’,psnr1);Qotherexp2part$葢2PSNR33.9157叶变換DCT咍达玛变換边囝1.哈达玛变换的源程序及结果J=imr
3、.2)*0.95&&i+1>=size(Q,2)*0.95)t=Q⑴;endendG(abs(G)4、2)对图像进行傅立叶变换、获得变换后的系数矩阵;(3)将傅立叶变换后系数矩阵的频谱用图像输出,观察频谱;(4)通过设定门限,将系数矩阵中95%的(小值)系数置为0,对图像进行反变换,获得逆变换后图像;(5)观察逆变换后图像质量,并比较原始图像与逆变后的峰值信噪比(PSNR)。(6)对输入图像进行离散余弦、哈达玛变换,重复步骤1-5;(7)比较三种变换的频谱情况、以及逆变换后图像的质量(PSNR)。4.实验结果1.DFT的源程序及结果J=imreadC10021033.bmp’);P=fft2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P,5、2)+j)=P(i+1,j);endQ=sort(G);fori=1:size(Q,2)if(i=size(Q,2)*0.95)t=Q⑴;endendG(abs(G)6、erexp2part实辁2PSNR32.4671尿囝交挨后复SS•耳里叶变挨DCT哈达玛变疾返回1.DCT的源程序及结果J=imreadC10021033.bmp’);P二dct2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P.2)+j)=P(i+1,j);endend0=sort(G);fori=1:size(Q,2)if(i=size(Q.2)*0.95)t=Q(i);endendG(abs(G)7、G(n*size(P,2)+m):endendf2=idct2(W);f3=uint8(f2);axes(handIes.axes2);imshow(f3)axes(handIes.axesl);imshow(J)psnr1=psnr(J,f3);set(handles.text3,’string’,psnr1);Qotherexp2part$葢2PSNR33.9157叶变換DCT咍达玛变換边囝1.哈达玛变换的源程序及结果J=imr
4、2)对图像进行傅立叶变换、获得变换后的系数矩阵;(3)将傅立叶变换后系数矩阵的频谱用图像输出,观察频谱;(4)通过设定门限,将系数矩阵中95%的(小值)系数置为0,对图像进行反变换,获得逆变换后图像;(5)观察逆变换后图像质量,并比较原始图像与逆变后的峰值信噪比(PSNR)。(6)对输入图像进行离散余弦、哈达玛变换,重复步骤1-5;(7)比较三种变换的频谱情况、以及逆变换后图像的质量(PSNR)。4.实验结果1.DFT的源程序及结果J=imreadC10021033.bmp’);P=fft2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P,
5、2)+j)=P(i+1,j);endQ=sort(G);fori=1:size(Q,2)if(i=size(Q,2)*0.95)t=Q⑴;endendG(abs(G)6、erexp2part实辁2PSNR32.4671尿囝交挨后复SS•耳里叶变挨DCT哈达玛变疾返回1.DCT的源程序及结果J=imreadC10021033.bmp’);P二dct2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P.2)+j)=P(i+1,j);endend0=sort(G);fori=1:size(Q,2)if(i=size(Q.2)*0.95)t=Q(i);endendG(abs(G)7、G(n*size(P,2)+m):endendf2=idct2(W);f3=uint8(f2);axes(handIes.axes2);imshow(f3)axes(handIes.axesl);imshow(J)psnr1=psnr(J,f3);set(handles.text3,’string’,psnr1);Qotherexp2part$葢2PSNR33.9157叶变換DCT咍达玛变換边囝1.哈达玛变换的源程序及结果J=imr
6、erexp2part实辁2PSNR32.4671尿囝交挨后复SS•耳里叶变挨DCT哈达玛变疾返回1.DCT的源程序及结果J=imreadC10021033.bmp’);P二dct2(J);fori=0:size(P,1)-1forj=1:size(P,2)G(i*size(P.2)+j)=P(i+1,j);endend0=sort(G);fori=1:size(Q,2)if(i=size(Q.2)*0.95)t=Q(i);endendG(abs(G)7、G(n*size(P,2)+m):endendf2=idct2(W);f3=uint8(f2);axes(handIes.axes2);imshow(f3)axes(handIes.axesl);imshow(J)psnr1=psnr(J,f3);set(handles.text3,’string’,psnr1);Qotherexp2part$葢2PSNR33.9157叶变換DCT咍达玛变換边囝1.哈达玛变换的源程序及结果J=imr
7、G(n*size(P,2)+m):endendf2=idct2(W);f3=uint8(f2);axes(handIes.axes2);imshow(f3)axes(handIes.axesl);imshow(J)psnr1=psnr(J,f3);set(handles.text3,’string’,psnr1);Qotherexp2part$葢2PSNR33.9157叶变換DCT咍达玛变換边囝1.哈达玛变换的源程序及结果J=imr
此文档下载收益归作者所有