《医学成像与处理技术》实验.docx
医学图像处理实验（一）【实验目的】：1 .运用Matlab对矩阵进行运算；2 .运用Matlab对图像的基本应用。【实验内容】:-、矩阵运算1、1、4 33 7-1 1A减C,4、一59>A乘以B。0-1-210-36230545670-330-11、2、3、二、图像处理 读取和显示彩色图像“lenajpg”；将彩色图像“lena.jpg”转换为灰度图像，保存为文件“lenal.bmp"将彩色图像“lena.jpg”转换为二值图像，保存为文件“lena2.bmp” ； 读取和显示索引图像"tree.tif"； 将索引图像"tree.tif"转化为彩色图像，保存为文件“treeLbmp” ；将图像“dust.bmp”所有像素的灰度值加50,保存为文件“dust50.bmp” ;【实验要求及上交内容】：1 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。（操作过程用到的命令，代码，显示或生成的新图像，均复制或截图、粘贴到实 验报告上。）【报告上交时间】：下课之前医学图像处理实验（二）【实验目的】：1 .运用Matlab对图像进行代数运算；2 .运用Matlab对图像进行几何运算。【实验内容】：一、图像的代数运算2、将彩色图像“lenajpg”转为灰度图，然后显示其直方图。3、将图像“dust.bmp”的数据类型转化为double型。4、将图像“al.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 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。（操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘 贴到实验报告上。）【报告上交时间】：实验课下课之前医学图像处理实验（三）【实验目的】:1 .运用Matlab对图像进行线性与非线性变换；2 .运用Matlab对图像进行直方图的均衡化和规定化。【实验内容】：运用Matlab编写程序，实现下述图像处理操作:图像的线性与非线性变换1 .对图像“dust.bmp”进行线性变换：2 .对图像 “lenna.jpg” 进行非线性变换：g(x, j) = 21og(l + /(x, j)3 .对图像“boy.tif”进行非线性变换：g(x,y) = 2/(x,y严二.图像的直方图均衡化和规定化。1.对图像“dust.bmp”进行直方图均衡化2,将图像“boy.tif”直方图规定化为图像“moon.tif”的直方图【实验要求及上交内容】：1 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。(操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘 贴到实验报告上。)【报告上交时间工下周实验课之前医学图像处理实验(四)【实验目的】:1 .运用Matlab对医学图像进行平滑；2 .运用Matlab对医学图像进行锐化。【实验内容】：运用Matlab编写程序，实现下述图像处理操作：二.医学图像的平滑.对图像“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进行二值化 提取轮廓。4 .对图像“dust.bmp”采用8邻域的高提升（常数系数k=3）拉普拉斯滤波 器（中心系数负）进行滤波增强。【实验要求及上交内容】：1 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。（操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘 贴到实验报告上。）【报告上交时间工下周实验课之前医学图像处理实验（五）【实验目的】:1 .运用Matlab对图像进行加噪声与噪声估计；2 .运用Matlab对图像进行均值滤波。【实验内容工运用Matlab编写程序，实现下述图像处理操作：1 .自己选取一张图像，对其进行高斯噪声污染。2 .对含噪声图像进行噪声估计，分析估计所得参数的准确程度。3 .对含噪声图像进行滤波(自选滤波方法)，对滤波结果进行分析，判断该 滤波方法是否最正确的图像恢复方法。【实验要求及上交内容】：1 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。(操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘 贴到实验报告上。)【报告上交时间】：下周实验课之前医学图像处理实验(六)【实验目的】：1 .运用Matlab对图像进行噪声估计、逆滤波；.运用Matlab对图像进行有参维纳滤波、基于相关函数的维纳滤波。【实验内容】：运用Matlab编写程序，实现下述图像处理操作：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(,carmove.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.jpg');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=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 .对程序及结果的说明。(操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘贴到 实验报告上。）【报告上交时间】：下周实验课之前医学图像处理实验（七）【实验目的】：1 .运用Matlab对图像进行全局阈值分割与局部阈值分割；2 .运用Matlab对图像进行区域生长法分割。【实验内容】：运用Matlab编写程序，实现下述图像处理操作：4 .对图像septagon.tif进行全局阈值分割，要求分别使用Ostu法和迭代法， 并得到两种方法计算所得的阈值。5 .对图像septagon2.tif进行局部阈值分割，并尽可能提高分割的效果，最后 对所采用方法进行评估。6 .对图像weld.tif进行区域生长法分割，把中间位置的裂缝分割出来。【实验要求及上交内容】：1 .操作所用的命令、步骤、源程序及代码；2 .原图像及处理后的图像或结果；3 .对程序及结果的说明。（操作过程用到的命令，代码，显示或生成的结果图像，均复制或截图、粘 贴到实验报告上。）【报告上交时间】：下周实验课之前
