《《医学成像与处理技术》实验.docx》由会员分享,可在线阅读,更多相关《《医学成像与处理技术》实验.docx(7页珍藏版)》请在taowenge.com淘文阁网|工程机械CAD图纸|机械工程制图|CAD装配图下载|SolidWorks_CaTia_CAD_UG_PROE_设计图分享下载上搜索。
1、医学图像处理实验(一)【实验目的】:1 .运用Matlab对矩阵进行运算;2 .运用Matlab对图像的基本应用。【实验内容】:-、矩阵运算1、1、4 33 7-1 1A减C,4、一59A乘以B。0-1-210-36230545670-330-11、2、3、二、图像处理 读取和显示彩色图像“lenajpg”;将彩色图像“lena.jpg”转换为灰度图像,保存为文件“lenal.bmp;将彩色图像“lena.jpg”转换为二值图像,保存为文件“lena2.bmp” ; 读取和显示索引图像tree.tif; 将索引图像tree.tif转化为彩色图像,保存为文件“treeLbmp” ;将图像“dus
2、t.bmp”所有像素的灰度值加50,保存为文件“dust50.bmp” ;【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的新图像,均复制或截图、粘贴到实 验报告上。)【报告上交时间】:下课之前医学图像处理实验(二)【实验目的】:1 .运用Matlab对图像进行代数运算;2 .运用Matlab对图像进行几何运算。【实验内容】:一、图像的代数运算2、将彩色图像“lenajpg”转为灰度图,然后显示其直方图。3、将图像“dust.bmp”的数据类型转化为double型。4、将图像“a
3、l.jpg”与图像“a2.jpg”进行加运算。5、将图像“al.jpg”与图像“a2.jpg”进行减运算。6、将图像“al.jpg”与图像“a2.jpg”进行乘运算。7、将图像“al.jpg”与图像“a2.jpg”进行除运算。二、图像的几何运算1、运用双立方插值,将图像“aaa.jpg”放大到原图的1.6。2、运用最近邻插值,将图像“doctor.jpg”顺时针旋转45度。3、将图像“al.jpg”中的美女裁剪下来。4、将图像“doctor.jpg”平移水平向右100像素,垂直向下50像素。5、将图像“bl.jpg”与图像“b2.jpg”拼接成一张完整的图像。【实验要求及上交内容】:1 .操作
4、所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘 贴到实验报告上。)【报告上交时间】:实验课下课之前医学图像处理实验(三)【实验目的】:1 .运用Matlab对图像进行线性与非线性变换;2 .运用Matlab对图像进行直方图的均衡化和规定化。【实验内容】:运用Matlab编写程序,实现下述图像处理操作:图像的线性与非线性变换1 .对图像“dust.bmp”进行线性变换:2 .对图像 “lenna.jpg” 进行非线性变换:g(x, j) = 21og(l + /(x, j)3 .对图像
5、“boy.tif”进行非线性变换:g(x,y) = 2/(x,y严二.图像的直方图均衡化和规定化。1.对图像“dust.bmp”进行直方图均衡化2,将图像“boy.tif”直方图规定化为图像“moon.tif”的直方图【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘 贴到实验报告上。)【报告上交时间工下周实验课之前医学图像处理实验(四)【实验目的】:1 .运用Matlab对医学图像进行平滑;2 .运用Matlab对医学图像进行锐化。【实验内容】:运用M
6、atlab编写程序,实现下述图像处理操作:二.医学图像的平滑.对图像“lennal.jpg”进行均值滤波,分别用3*3滤波器和5*5滤波器。 并比拟两个滤波器的滤波效果,说明差异在哪里?为什么?1 .对图像“Ienna2.jpg”进行中值滤波,分别用3*3滤波器和5*5滤波器。 并比拟两个滤波器的滤波效果,说明差异在哪里?为什么?三.医学图像的锐化1 .对图像“dust.bmp”采用Sobel梯度获取梯度图像。2 .对图像“lenna.jpg”采用Prewitt算子获取梯度图像,并将原图像与梯度图 像相叠加。3 .对图像“lennajpg”采用sobel算子计算梯度,并用阀值0.03进行二值化
7、 提取轮廓。4 .对图像“dust.bmp”采用8邻域的高提升(常数系数k=3)拉普拉斯滤波 器(中心系数负)进行滤波增强。【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘 贴到实验报告上。)【报告上交时间工下周实验课之前医学图像处理实验(五)【实验目的】:1 .运用Matlab对图像进行加噪声与噪声估计;2 .运用Matlab对图像进行均值滤波。【实验内容工运用Matlab编写程序,实现下述图像处理操作:1 .自己选取一张图像,对其进行高斯噪声污染。
8、2 .对含噪声图像进行噪声估计,分析估计所得参数的准确程度。3 .对含噪声图像进行滤波(自选滤波方法),对滤波结果进行分析,判断该 滤波方法是否最正确的图像恢复方法。【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘 贴到实验报告上。)【报告上交时间】:下周实验课之前医学图像处理实验(六)【实验目的】:1 .运用Matlab对图像进行噪声估计、逆滤波;.运用Matlab对图像进行有参维纳滤波、基于相关函数的维纳滤波。【实验内容】:运用Matlab编写程序
9、,实现下述图像处理操作:1 .对图像“noise.jpg”进行噪声估计,判断出噪声类型、均值、方差。f=imread(,noise.jpg,);b=imcrop(f);imhist(b)p=imhist(b);p=p./numel(b);z=0:255;z=z./255;u=sum(z.*p)0.1371 v=sum(z-u).八2).*p)0.0176根据一下直方图,推测该噪声类型是高斯噪声。.运动模糊图像carmove.jpg运动模糊的方向和距离为同一个值。请通 过实验确定该值,并进行逆滤波。提示:估计值可以通过反复试验,每次试验 的参数以步长为10进行尝试。f=imread(,carmo
10、ve.jpg,);psf二fspecial(motion, 1 ()(), 1 ()();imshow(deconvwnr(f,psf).对运动模糊带噪图像“carmovenoise.jpg”进行有参维纳滤波,所需条件请参 考如下图中的噪声:噪声图像“noise.jpg的噪声。原图像:用第二题滤波的结果。退化函数PSF:与第二题相同。f=imread(,carmove.jpg,);psf二fspecial(motionl 00,100);g=deconvwnr(f,psf);sf=abs(fft2(g).A2;fa=sum(sf(:)/numel(g);noise=imread(noise.j
11、pg);nf=abs(fft2(noise).A2;na=sum(nf(:)/numel(noise);R=na/fa;f=imread(,carmovenoise.jpg);k=deconvwnr(f,psf,R);imshow(k)2 .对运动模糊带噪图像“carmovenoise.jpg”进行基于相关函数的维纳滤波,所 需条件请参考如下。图中的噪声:噪声图像“noise.jpg的噪声。原图像:用第二题滤波的结果。退化函数PSF:与第二题相同。f=imread(,carmove.jpg,);psf=fspecial(motion 100,100);g=deconvwnr(f,psf);sf
12、=abs(fft2(g).A2;fcorr=fftshift(real(ifft2(sf);noise二imread(noise.jpg);nf=abs(fft2(noise).A2;ncorr=fftshift(real(ifft2(nf);f=imread(,carmovenoise.jpg,);k=deconvwnr(f,psf,ncorr,fcorr);imshow(k)【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;.原图像及处理后的图像或结果;2 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘贴到 实验报告上。)【报告上交
13、时间】:下周实验课之前医学图像处理实验(七)【实验目的】:1 .运用Matlab对图像进行全局阈值分割与局部阈值分割;2 .运用Matlab对图像进行区域生长法分割。【实验内容】:运用Matlab编写程序,实现下述图像处理操作:4 .对图像septagon.tif进行全局阈值分割,要求分别使用Ostu法和迭代法, 并得到两种方法计算所得的阈值。5 .对图像septagon2.tif进行局部阈值分割,并尽可能提高分割的效果,最后 对所采用方法进行评估。6 .对图像weld.tif进行区域生长法分割,把中间位置的裂缝分割出来。【实验要求及上交内容】:1 .操作所用的命令、步骤、源程序及代码;2 .原图像及处理后的图像或结果;3 .对程序及结果的说明。(操作过程用到的命令,代码,显示或生成的结果图像,均复制或截图、粘 贴到实验报告上。)【报告上交时间】:下周实验课之前