数字图像处理指导-测控10版本.doc
《数字图像处理指导-测控10版本.doc》由会员分享,可在线阅读,更多相关《数字图像处理指导-测控10版本.doc(37页珍藏版)》请在咨信网上搜索。
1、数字图像处理指导-测控10版本 作者: 日期:37 个人收集整理 勿做商业用途 实验一:图像预处理实验(验证性实验)实验目的:要求学生掌握图像读取、显示和保存的方法,掌握空间域增强方法和频率域增强方法.实验要求:上机运行,调试通过。实验内容:1) 对附录一中的练习111、练习13,进行上机运行和调试.2) 对上述练习,改变一些函数的参数,观察运行结果的变化。提交试验报告的内容:1、 实现对一幅图(比如:pout.tif)进行直方图均衡化的源代码;2、 显示直方图均衡化前后的图像及其直方图;3、 针对处理前后的图像和对应的直方图进行分析,得出直方图均衡化的作用。实验二、图像的傅立叶变换(验证性实
2、验)实验目的1了解图像变换的意义和手段;2熟悉傅里叶变换的基本性质;3熟练掌握FFT的方法与应用;4理解频率域高低频与时域图像的对应关系;5通过本实验掌握利用MATLAB编程实现数字图像的傅立叶变换。实验原理1、应用傅立叶变换进行图像处理傅里叶变换是线性系统分析的一个有力工具,它能够定量地分析诸如数字化系统、采样点、电子放大器、卷积滤波器、噪音和显示点等的作用。通过实验培养这项技能,将有助于解决大多数图像处理问题。对任何想在工作中有效应用数字图像处理技术的人来说,把时间用在学习和掌握博里叶变换上是很有必要的。2、傅立叶(Fourier)变换的定义对于二维信号,二维Fourier变换定义为:二维
3、离散傅立叶变换为: 图像的傅立叶变换与一维信号的傅立叶变换一样,有快速算法,具体参见参考书目,有关傅立叶变换的快速算法的程序不难找到.实际上,现在有实现傅立叶变换的芯片,可以实时实现傅立叶变换.3、用matlab编程观察图像频域中高低频率与时域图像的对应关系:-% 本代码主要演示图像fft2变换的频谱图中高频和低频分别对应 图像时域图中的哪些部分,加深对频域的理解-clear;close all;A=imread(lena.jpg); %读入图像A=rgb2gray(A);figure(1);imshow(A); %显示图像A=double(A);%- 对图像进行fft变换,看看频谱图和iff
4、t变换的结果 -F0=fft2(A); 计算二维傅里叶变换F1=fftshift(F0); %将直流分量移到频谱图的中心figure(2), imshow(log(abs(F1)+1),0 10);显示变换后的频谱图F0=fftshift(F1);R0=ifft2(F0); %进行反变换figure(3);imshow(uint8(R0)); %将变换的结果转换为uint8数据类型,并显示,比较反变换的结果与原图。%- 把频谱中直流和低频系数变为0,只保留高频系数,然后反变化,% 看看高频系数对应图像时域中的哪一部分?%-m,n=size(R0); 得到图像大小i0=fix(m/2); %图像
5、中心点的坐标j0=fix(n/2);F2=F1;r=50; 设置系数改变区域的半径。%把中心点附近2W2W的低频全部设置为0;for i=1:m for j=1:n d=sqrt(i-i0)2+(j-j0)2); if(d=r) F2(i,j)=0; end; end;end;F3=fftshift(F2); 将直流分量移到频谱图的四个角R1=ifft2(F3); %进行反变换figure(4), imshow(log(abs(F2)+1),0 10);显示去除低频的频谱图figure(5);imshow(uint8(R1)); 将变换的结果转换为uint8数据类型,并显示,比较反变换的结果与
6、原图。%- 把频谱中高频系数变为0,只保留低频系数,然后反变化,% 看看低频系数对应图像时域中的哪一部分?-m,n=size(R0); 得到图像大小i0=fix(m/2); %图像中心点的坐标j0=fix(n/2);F5=F1;r=50; %设置系数不改变区域的半径。把中心点附近2W*2W的低频全部设置为0;for i=1:m for j=1:n d=sqrt(i-i0)2+(j-j0)2); if(dr) F5(i,j)=0; end; end;end;F6=fftshift(F5); 将直流分量移到频谱图的四个角R2=ifft2(F6); %进行反变换figure(6), imshow(l
7、og(abs(F5)+1),0 10);显示去除低频的频谱图figure(7);imshow(uint8(R2)); 将变换的结果转换为uint8数据类型,并显示,比较反变换的结果与原图。%- 把高频和低频的反变换结果相加,与原图相比较,看看结果如何?%-R3=R1+R2;figure(8);imshow(uint8(R3));实验内容1打开计算机,安装和启动MATLAB程序;程序组中“work”文件夹中应有待处理的图像文件;2分析、运行其中的代码,观察结果。3 a)。对三幅图像内容区别较大图做fft并显示其频谱; b) 讨论时域图像内容与FFT频谱之间的对应关系。4记录和整理实验报告。提交试
8、验报告的内容:1、对图像内容区别较大的三幅图分别做fft并显示其频谱.2、讨论不同的图像内容的FFT频谱图之间的区别。实验三、数字水印(设计性实验)实验要求1了解数字水印的基本技术;2设计一个基于LSB的数字水印嵌入,提取、攻击方案;3编写代码实现基于LSB的数字水印嵌入、提取功能;参考资料:1、 matlab编程资料.2、 黄仿元,基于LSB的数字水印算法及MATLAB实现。提交的报告内容:1、 数字水印的基本知识介绍。2、 水印嵌入、提取设计方案3、 LSB数字水印的嵌入、提取代码及结果图。4、水印图像被攻击后提取的水印结果图。 附录一:MATLAB数字图像处理1 MATLAB简介11 主
9、要用途及特点主要用途:算法研究主要特点:语句功能强大;不能生成可执行文件.12 MATLAB工作环境 121 Matlab桌面 桌面包括5个子窗口:命令窗口、工作空间窗口、当前目录窗口、历史命令窗口、一个或多个图形窗口(仅在用户显示图形式出现)。命令窗口是用户在提示符()处键入MATLAB命令和表达式的地方,也是显示那些命令输出的地方。工作空间窗口显示对话中创建的变量和它们的某些信息。当前目录窗口显示当前目录的内容(即路径)。历史命令窗口包含用户已在命令窗口中输入的命令的纪录. 122 使用MATLAB编辑器创建M文件 123 获得帮助(1)help 函数名(2)look 关键词2 数字图像的
10、读取、显示、保存、数据类型和图像类型、数据类型间的转换、图像类型间的转换21 图像的读取语法:imread ( filename ) 说明:读取图像格式名称 描述 可识别扩展符TIFF加标记的图像文件格式.tif, 。tiffJPEG联合图像专家组。jpg, 。jpegGIF图形交换格式.gifBMPWindows位图.bmpPNG可移植网络图形.pngXWDX Window 转储。xwd 语法:M, N=size ( filename ) 说明:给出一幅图像的行数和列数22 图像的显示语法:imshow ( f, G) imshow (f, low high) imshow (f, )说明:
11、G是显示该图像的灰度级数; 小于或等于low的值都显示为黑色,大于或等于high的值都显示为白色. 自动将变量low设置为f的最小值,将high设置为f的最大值.23 图像的保存语法:imwrite ( f, filename) 说明:练习1 f = imread ( saturn。tif); M, N = size( f ); g = imread ( trees.tif ); imshow ( f ); figure, imshow (g); 显示另一幅图像 imwrite (f, s2.jpg)24 数据类型名称 描述double双精度浮点数,范围为10exp(308) 10exp(30
12、8),8字节uint8无符号8比特整数,1字节uint16无符号16比特整数,2字节uint32无符号32比特整数,4字节int8有符号8比特整数,1字节int16有符号16比特整数,2字节int32有符号32比特整数,4字节single单精度浮点数,范围为10exp(38) 10exp(38) ,4字节char字符logical值为0或1四种常用类型:double,uint8,char,logical。25 图像类型名称 描述亮度图像uint8类范围为0 255、double类归一化为0 1二值图像图像取值只有0和1的逻辑数组索引图像索引图像RGB图像彩色图像2 6 数据类与图像类型间的转换
13、261 数据类间的转换语法:B = data_class_name ( A )举例:若A是8位图像,则B = double ( A )转换为双精度图像.262 图像类型间的转换函数名称将输入转换为有效的输入图像数据类im2uint8uint8Logical, uint8, uint16和doubleim2uint16uint16Logical, uint8, uint16和doublemat2graydoubledoubleim2doubledoubleLogical, uint8, uint16和doubleim2bwlogicaluint8, uint16和double 练习2f=1 2;
14、 3 4;g=mat2gray(f);gb=im2bw(g, 0.6) %0。6 is a threshold3 亮度变换与空间滤波31 亮度变换函数311 基本亮度变换函数 语法:g=imadjust(f, low_in high_in, low_out high_out,gamma)说明:将图像f中的亮度值影响到g中的新值,即将low_in至high_in之间的值映射到low_out至high_out之间的值, low_in以下的值映射为low_out, high_in以上的值映射为high_out.参数gamma指定了映射曲线的形状。练习3f = imread ( pout。tif);
15、imshow(f);g1=imadjust(f, 0 1, 1 0); figure, imshow(g1);g2=imadjust(f, 0。5 0.75, 1 0, 0。5); figure, imshow(g2);312 对比度拉伸变换对数变换方法:g=c*log(1+double(f))说明:对数变换的一项主要应用是压缩动态范围,常用于傅里叶频谱显示。练习4f = imread ( pout.tif); F=fft2(f); Fourier TransformFC=fftshift(F); %将变换原点移到频率矩形的中心。imshow(abs(FC), );S2= log(1+ abs
16、(FC));figure,imshow(S2, );对比度拉伸方法:g=1。/(1+(m。/(double(f)+eps).E)说明:eps是浮点数的相对精度,可避免f出现0值时的溢出现象。32 直方图处理与函数绘图321 生成并绘制图像的直方图 语法:h=imhist(f, b) p= imhist(f, b)/num1(f) 说明:b是用于形成直方图的灰度级的个数.P是归一化直方图.练习5f = imread ( pout。tif); imhist(f);322 直方图均衡化语法:g=histeq(f, nlev)说明:nlev是为输出图像制定的灰度级数。练习6f = imread ( t
17、ire。tif); imshow(f); figure,imhist(f); ylim(auto); %自动设定y轴坐标范围和刻度 g=histeq(f,256); figure, imshow(g); figure,imhist(g); ylim(auto);33 空间滤波331 线性空间滤波 语法:g=imfilter(f, w, filtering_mode, boundary_options, size_options)说明:w为滤波掩模选项描述滤波类型corr滤波器通过使用相关来完成.该值是默认值。conv滤波器通过使用卷积来完成边界选项P输入图像的边界通过用值P来扩展。P的默认值为
18、0。replicate图像大小通过复制外边界的值来扩展。symmetric图像大小通过反射其边界来扩展。circular图像大小通过将图像看成是一个二维周期函数的一个周期来扩展大小选项full输出图像的大小与被扩展图像的大小相同same输出图像的大小与输入图像的大小相同练习7f = imread ( saturn。tif); w=ones(31); 单位矩阵掩模gd=imfilter(f,w);imshow(gd,);332 非线性空间滤波略.34 图像处理工具箱的标准空间滤波器341 线性空间滤波器语法:w=fspecial(type,parameter)说明:见下表练习8f = imrea
19、d ( saturn。tif); w4=fspecial(laplacian, 0)w8=1 1 1; 1 8 1; 1 1 1;g4=fimfilter(f, w4, replicate);g8=fimfilter(f, w8, replicate);imshow(f);figure, imshow(g4); figure, imshow(g8); 342 非线性空间滤波器语法:g=ordfilt2(f, order, domain)说明:使用邻域的一组排序元素中的第oder个元素来代替f中的每个元素,而该邻域则由domain中的非零元素指定.语法:g=medfilt2(f, m n,pad
- 配套讲稿:
如PPT文件的首页显示word图标,表示该PPT已包含配套word讲稿。双击word图标可打开word文档。
- 特殊限制:
部分文档作品中含有的国旗、国徽等图片,仅作为作品整体效果示例展示,禁止商用。设计者仅对作品中独创性部分享有著作权。
- 关 键 词:
- 数字图像 处理 指导 测控 10 版本
1、咨信平台为文档C2C交易模式,即用户上传的文档直接被用户下载,收益归上传人(含作者)所有;本站仅是提供信息存储空间和展示预览,仅对用户上传内容的表现方式做保护处理,对上载内容不做任何修改或编辑。所展示的作品文档包括内容和图片全部来源于网络用户和作者上传投稿,我们不确定上传用户享有完全著作权,根据《信息网络传播权保护条例》,如果侵犯了您的版权、权益或隐私,请联系我们,核实后会尽快下架及时删除,并可随时和客服了解处理情况,尊重保护知识产权我们共同努力。
2、文档的总页数、文档格式和文档大小以系统显示为准(内容中显示的页数不一定正确),网站客服只以系统显示的页数、文件格式、文档大小作为仲裁依据,平台无法对文档的真实性、完整性、权威性、准确性、专业性及其观点立场做任何保证或承诺,下载前须认真查看,确认无误后再购买,务必慎重购买;若有违法违纪将进行移交司法处理,若涉侵权平台将进行基本处罚并下架。
3、本站所有内容均由用户上传,付费前请自行鉴别,如您付费,意味着您已接受本站规则且自行承担风险,本站不进行额外附加服务,虚拟产品一经售出概不退款(未进行购买下载可退充值款),文档一经付费(服务费)、不意味着购买了该文档的版权,仅供个人/单位学习、研究之用,不得用于商业用途,未经授权,严禁复制、发行、汇编、翻译或者网络传播等,侵权必究。
4、如你看到网页展示的文档有www.zixin.com.cn水印,是因预览和防盗链等技术需要对页面进行转换压缩成图而已,我们并不对上传的文档进行任何编辑或修改,文档下载后都不会有水印标识(原文档上传前个别存留的除外),下载后原文更清晰;试题试卷类文档,如果标题没有明确说明有答案则都视为没有答案,请知晓;PPT和DOC文档可被视为“模板”,允许上传人保留章节、目录结构的情况下删减部份的内容;PDF文档不管是原文档转换或图片扫描而得,本站不作要求视为允许,下载前自行私信或留言给上传者【天****】。
5、本文档所展示的图片、画像、字体、音乐的版权可能需版权方额外授权,请谨慎使用;网站提供的党政主题相关内容(国旗、国徽、党徽--等)目的在于配合国家政策宣传,仅限个人学习分享使用,禁止用于任何广告和商用目的。
6、文档遇到问题,请及时私信或留言给本站上传会员【天****】,需本站解决可联系【 微信客服】、【 QQ客服】,若有其他问题请点击或扫码反馈【 服务填表】;文档侵犯商业秘密、侵犯著作权、侵犯人身权等,请点击“【 版权申诉】”(推荐),意见反馈和侵权处理邮箱:1219186828@qq.com;也可以拔打客服电话:4008-655-100;投诉/维权电话:4009-655-100。