医学成像术及图像处理.doc

《医学成像术及图像处理.doc》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、医学成像技术及图像处理实验报告几何均值滤波器学生：马亚红学号：3105037051考核形式：考查实验日期：2006.2几何均值滤波器一实验任务几何均值滤波器计算输入图像中大小为m*m的正方形区域的几何均值作为正方形区域中心点的滤波输出。将这一过程遍历图像的每一像素位置，即可得到输出图像。几何均值：一个长度为N的序列各元素的乘积的N次方根。本实验仅考虑m=n=2,3,4三种正方形区域。对于边缘上的点，正方形区域中超出图像边界的部分予以忽略。以实验26**中的灰度图像PCB.bmp作为实验图像，在matlab中利用函数imnoise在图像中添加不

2、同类型的噪声，试验滤波效果。PCB.bmp原始图像二设计程序实验模块框图：在MATLAB里面调用imread函数读入PCB原图像调用imnoise函数加入各类噪声将加入噪声的图像读入VC++进行滤波得到滤波处理后的图像三源程序1)matlab的加噪声程序clearall;%变量清除I=imread('PCB.bmp');%读入图片I1=imnoise(I,'gaussian');%加入高斯噪声v=zeros(300,420);%创建一个和原始图像大小相等的300*420的零矩阵v=v+0.02;%给矩阵的每个元素加上因子0.02I4=imno

3、ise(I,'salt&pepper');%加入椒盐噪声I5=imnoise(I,'speckle');%加入speckle噪声figure(1),imshow(I1);title('加入高斯噪声的图像');figure(4),imshow(I4);title('加入焦盐噪声的图像');figure(5),imshow(I5);title('加入speckle噪声的图像');2)VC++实现几何均值滤波器的源程序特别说明：该实验是在本科毕业设计已经做好的框架photostar上添加的代码。//Filter.cpp:implementation

4、oftheCFilterclass.///////////////////////////////////////////////////////////////////////////////////////////////////////////////////*****文件名：几何均值滤波器*************************************//*****文件功能：计算输入图像中大小为m*m的正方形区域的几何均值作为//*****正方形区域中心点的滤波输出。将这一过程遍历图像的每一像素位//*****置，即可得到输

5、出图像。***********************************//*****创建日期：2006.2.01**************************************//*****作者：马亚红**********************************************/#include"stdafx.h"#include"photostar.h"#include"Filter.h"#include"math.h"#ifdef_DEBUG#undefTHIS_FILEstaticcharTHIS_F

6、ILE[]=__FILE__;#definenewDEBUG_NEW#endif////////////////////////////////////////////////////////////////////////Construction/Destruction//////////////////////////////////////////////////////////////////////CFilter::~CFilter(){}CFilter::CFilter(CImageObject*pImageObject,inti

7、WinSize=3){ASSERT(pImageObject!=NULL);//声明m_pImageObject不得为空m_iWinSize=iWinSize;m_pImageObject=pImageObject;}BOOLCFilter::GeometryFilter(){if(m_pImageObject==NULL){returnFALSE;}intiWidth=m_pImageObject->GetWidth();//获取图像的宽度intiHeight=m_pImageObject->GetHeight();//获取图像的高度/**

8、*******************原始图像*******************/unsignedchar*pOldBuffer,*pNewBuffer,*pO