实验2 图像增强(一)

《实验2 图像增强(一)》由会员上传分享，免费在线阅读，更多相关内容在行业资料-天天文库。

1、实验二图像增强（一）实验目的：1.熟悉亮度变换函数imadjust的原理和使用；2.会分析图像的直方图，理解直方图均衡化的基本原理，熟练掌握直方图显示函数imhist、均衡化函数histeq；3.练习自己写函数。实验内容：1、做出四幅心脏图像（heart1-heart4）的直方图，指出哪幅图的效果最好，结合直方图分析质量不好的图像出现的问题和可能原因。试用我们学过的方法将三幅质量不好的图片变好。图像偏暗图像偏亮图像动态范围小图像动态范围适中>>A=imread('heart1.bmp');>>B=imread(

2、'heart2.bmp');>>C=imread('heart3.bmp');>>D=imread('heart4.bmp');>>imshow(A),figure,imhist(A)>>imshow(B),figure,imhist(B)>>imshow(C),figure,imhist(C)>>imshow(D),figure,imhist(D)修改方案：>>f1=imadjust(A,[0,0.2],[0,1]);>>imshow(f1)>>figure,imhist(f1)>>f2=imadjust(B

3、,[0.6,1],[0,1]);>>imshow(f2)>>figure,imhist(f2)>>f3=imadjust(C,[0.3,0.6],[0,1]);>>imshow(f3)>>figure,imhist(f3)2、分别做出图片heart1和moon的直方图，并对其进行均衡化。将得到的结果图片和原图片比较，效果如何？为什么对不同的图片进行均衡化会有如此差别？根据直方图总结什么样的图片不适合用直方图均衡化方法进行增强。通过下图效果对比发现，如果直方图偏向一边即直方条多数在一个灰度区段而且是最小灰度级的，

4、那么直方图均衡化后的效果就没那么好。>>A=imread('heart1.bmp');>>imshow(A),figure,imhist(A)>>g1=histeq(A);>>imshow(g1),figure,imhist(g1)>>E=imread('moon.tif');>>imshow(E),figure,imhist(E)>>g5=histeq(E);>>imshow(g5),figure,imhist(g5)3、编写一个m文件，完成以下功能：输入参数：图片文件处理：（1）计算图像灰度的最大、最小值和

5、平均值。（2）将图像做反色变换，并画出变换函数的映射图；输出：将每次处理的结果显示出来用breast.tif图片进行测试，看一下你编写的函数是否可以达到课本48页图（b）的效果。M文件function[max,min,ave,g]=test(f)max=f(1,1);min=f(1,1);g=imadjust(f,[0,1],[1,0]);[m,n]=size(f);ave=sum(sum(f))/(m*n);fori=1:mforj=1:niff(i,j)>=maxmax=f(i,j);endiff(i,j)

6、<=minmin=f(i,j);endendend>>f=imread('breast.tif');>>[max,min,ave,g]=test(f);>>imshow(g)>>maxmax=255>>minmin=21>>aveave=66.7099>>plot(f,g);4、对于三段的分段线性变换，是否能用imadjust函数实现？考虑三段灰度变换如何编写m文件im3adjust实现，输入为图片、a，b和c，d的值（见PPT），输出为经过变换后的图片。用pout.tif图片进行测试，并画出变换函数（仿照课本

7、59页）。我认为能用imadjust实现三段的分段线性变换。function[g]=im3adjust(f,a,b,c,d)[m,n]=size(f);g=zeros(m,n);fori=1:mforj=1:niff(i,j)>0&f(i,j)<=ag(i,j)=(c/a)*f(i,j);elseiff(i,j)>a&f(i,j)<=bg(i,j)=((d-c)/(b-a))*(f(i,j)-a)+c;elseiff(i,j)<=255g(i,j)=((255-d)/(255-b))*(f(i,j)-b)+d;

8、endendendendendg=mat2gray(g);>>imshow(f);>>imhist(f);>>[g]=im3adjust(f,70,160,0,255);>>imshow(g);>>imhist(g);>>plot(f,g);5、将图片lung.bmp进行旋转90度，缩小为原图一半的变换。>>I=imread('lung.bmp');>>Lung1=imrotate(I