数图实验报告总汇 (8000字) - 实验报告 - 书业网.doc

数图实验报告总汇 (8000字) - 实验报告 - 书业网 实验二：数字图像增强 实验（一）为必做，写入实验报告；在时间宽裕的情况下继续做实验（二），但不写入报告。 实验目的：验证数字图像增强算法 实验设备与软件：硬件：P4微机； 软件：Matlab 实验原理与实验内容： （一）采用对比度线性展宽算法增强人物图像。 对比度线性展宽，原图像像素灰度 f(i,j) 与处理后图像像素灰度 g(i,j) 的映射关系可用下图表示。 其中，K1、K2、K3为三段折线的斜率。用计算公式表达像素的灰度映射关系为： I(i,j)?fa?K1?I(i,j) ? J(i,j)??K2?(I(i,j)?fa)?ga fa?I(i,j)?fb ?K3?(I(i,j)?f)?gI(i,j)?fb bb? (i?1,2,...,m;j?1,2,...,n) K1?ga/fa K2?(gb?ga)/(fb?fa);K3?(255?gb)/(255?fb); 给定灰度图像girl.bmp，采用matlab编程，请对其进行线性对比度展宽处理。原图中重 要景物灰度分布在[fa, fb]范围，希望处理后图像的重要景物灰度分布范围为[ga, gb]。相关参数如下：fa=30, fb=140, ga=60, gb=200。实现上述算法，最后显示处理后的图像，并将其取名为new-girl，在当前路径下保存为bmp格式。 编程思路提示： 1) 对fa, fb, ga, gb进行赋值； 2) 读取图像girl.bmp的数据，存放到矩阵I； 3) 读取矩阵I的行与列的大小，赋给变量m与n 4) 把I的数据类型从uint8转换为double 5) 计算三段折线的斜率k1、k2、k3; 6) 用循环语句，根据对比度线性展宽的计算公式，计算对每个像素的新灰度值，赋值给 新矩阵J； 7) 把I与J的数据类型从double转换为uint8； 8) 将矩阵J保存为新图像“new_girl.bmp”。 9) 在同时显示原始图像girl.bmp与新图像new_girl.bmp。并分别用“原始图像”与“新 图像”为标题。 对比度线性展宽程序： 1. clc,clear 2. fa=70; 3. fb=180; 4. ga=40; 5. gb=220; 6. I=imread(#39;girl.bmp#39;); 7. [m,n]=size(I); 8. k1=4/7; 9. k2=18/11; 10. k3=7/15; 11. I=double(I); 12. for i=1:m 13. for j=1:n 14. if I(i,j)lt;fa 15. J(i,j)=k1*I(i,j); 16. elseif I(i,j)lt;fb 17. J(i,j)=k2*(I(i,j)-fa)+ga; 18. else I(i,j)gt;=fb 19. J(i,j)=k3*(I(i,j)-fb)+gb; 20. end 21. end 22. end 23. I=uint8(I); 24. J=uint8(J); 25. imwrite(I, #39;girl.bmp#39;) 26. imwrite(J, #39;new girl.bmp#39;) 27. subplot(1,2,1),imshow(I),title(#39;原始图像#39;); 28. subplot(1,2,2),imshow(J),title(#39;新图像#39;); 10) （二）采用灰级窗算法增强CT图像中的肺部区域。 灰级窗算法只显示指定灰度级范围内的信息，并充分其增强对比度 本实验根据给定的CT图像，通过灰级窗算法，充分增强肺部区域，尽量不显示其他人体组织。用计算公式灰度映射关系请同学们自行推导。已知肺部区域的图像灰度范围大致为：fa=45；fb=160。采用matlab编程实现上述算法，最后显示处理后的图像，并将其取名为lung_window，在当前路径下保存为bmp格式。 实验三：数字图像的几何变换 实验（一）为必做；实验时间宽裕的同学继续做实验（二）。 实验目的：验证数字图像的几何变换算法 实验设备与软件：硬件：P4微机； 软件：Matlab 实验原理： 图像的几何变换是通过将图像中所有像素按要求进行移动来实现的。通过像素坐标变换公式，将原