基于灰度变换的骨骼图像增强设计.docVIP

燕 山 大 学 计算机视觉(图像处理技术) 课程设计报告 设计题目：基于灰度变换的骨 骼扫描图像增强设计 学院（系）： 里仁学院 年级专业： 12级检测二班 学生姓名： 曹洪乐 是否答辩： 否 指导教师： 刘斌 完成日期： 2015年1月1日  电气工程学院《课程设计》任务书 课程名称： 计算机视觉 基层教学单位：仪器科学与工程系 指导教师： 刘斌 学号 121203021028 学生姓名 曹洪乐 （专业）班级 检测二班 设计题目 基于灰度变换的骨骼扫描图像增强设计 设 计 技 术 参 数 利用灰度变换实现骨骼扫描图像的增强处理，自行设计软件流程，配置相关参数，对比不同参数的增强效果。 设 计 要 求 (1)利用Matalb软件编写相关程序； (2)撰写课程设计报告，报告内容应包含算法原理、程序流程和本人的主要工作；报告其他内容详见模板。 参 考 资 料 课程教材 《数字图像处理 Matlab版》Rafael C. Gonzalez著，阮秋琦等译 《精通Matlab GUI设计 第三版》陈垚光 等著 MATLAB软件相关参考书 周次 前半周 后半周 应 完 成 内 容 搜集资料，总体设计，编写程序； 调试程序、撰写报告，准备答辩用PPT； 指导教 师签字 基层教学单位主任签字 说明：1、此表一式四份，系、指导教师、学生各一份，报送院教务科一份。 2、学生那份任务书要求装订到课程设计报告前面。 电气工程学院 教务科 理论基础或算法原理设计流程 matlab实现软件编写： （1）读入原图“guge.tif” （2）拉普拉斯变换 （3）锐化原图 （4）soble算子运算 （5）均值平滑 （6）乘机图像 （7）再次锐化 （8）伽马变换 h1=[0,-1,0;-1,4,-1;0,-1,0]; 离散拉普拉斯所用的滤波器掩膜,用h1与原图进行卷积，得到中心像素点与邻域像素点的差值。 M1=imfilter(I,h1);%掩膜得到的拉普拉斯变换图M1 M2=M1+I; %与原图相加得到M2，锐化原图 一幅图像的边缘是通过一阶和二阶数字导数来检测得到的。边缘的宽度取决于从初始灰度级跃变到最终灰度级的斜坡的长度。这个长度又取决于斜度，而斜度又取决于模糊程度。所以，我们可知，模糊的边缘使其变粗，而清晰的边缘使其变得较细。 在实践中计算数字梯度时最常用的是Prewitt算子和Sobel算子，虽然Prewitt模版的实现比Sobel模版更简单，但是Sobel模版在抑制噪声上更胜一筹，考虑到图片中的噪声，所以在此处用Sobel模版。 模版是用于求梯度分量的，其中Sobel算子为 H1= -1 -2 -1 0 0 0 1 2 1 H2= 1 0 -1 2 0 -2 1 0 -1 实现图像边缘检测的程序代码（Matlab）为： h2=[-1,-2,-1;0,0,0;1,2,1]; h3=[-1,0,1;-2,0,2;-1,0,1]; Gx=imfilter(I,h2); Gy=imfilter(I,h3); for i=1:m for j=1:n M3(i,j)=abs(Gx(i,j))+abs(Gy(i,j)); end end % sobel算子运算结果得到M3，保留边缘去除噪声 根据最后运行的结果显示，Sobel图像的边缘要比拉普拉斯图像的边缘突出很多。拉普拉斯变换作为一种二阶微分算子，在图像细节的增强处理方面有明显的优点，但拉普拉斯变换与梯度变换相比较会产生更多的噪声。其中位于平滑区域的噪声非常显眼，梯度变换在灰度变换的区域的相应要比拉普拉斯更强烈，而梯度变换对噪声和小细节的响应要比拉普拉斯变换弱，而且可以通过均值滤波器对其进行平滑处理而进一步的降低，对梯度图像进行平滑处理并用拉普拉斯图像与他相乘。处理后的结果在灰度变化强的区域仍然保留细节，而在灰度变化相对平坦的区域则减少噪声，这种处理可以看做是将拉普拉斯变换与梯度变