DIP实验三.docVIP

DIP实验三 Image Enhancement Using Intensity Transformations 班级：08计算 学号：08 姓名： 采用技术及主要计算公式： 软件环境：MATLAB 7.0 通过对数变换对图像灰度进行变换，使一窄带低灰度输入图像值映射为一个宽带输出值，计算公式为：s = c log(1+r) 3) 通过幂次变换对图像灰度进行变换，计算公式为：s = c r ^γ 其中，当γ1时，输入的窄带暗值将映射到宽带高输出值。 当γ1时，输入的窄带高值将映射到宽带暗输出值。 结果: (b) 图 3.1 (a)原始人的脊椎骨折的核磁共振图 (b)对式 s = c log(1+r) (取c=1)，应用 对数变换的结果。 (b) (c) (d) (e) (f) 图 3.2 (a) 原始人的脊椎骨折的核磁共振图 (b)~(f) 应用式s = c r ^γ并且c=1，γ分别等于0.6，0.4，0.3，0.2, 1, 2时的变换结果。 实验小结： 对数变换和幂次是进行动态范围处理的基本工具，在MATLAB 7.0 中，对数变换和幂次变换分别通过如下表达式实现：(f为输入图像矩阵) 对数变换：g = c*log(1 + double(f))，其中c是一个常数，在实验中取c = 1。 幂次变换：g= c*(double(f) + eps).^ γ,其中c和γ为常数，在实验中取c=1，γ取一系列整数值。 对数变换、幂次变换共同的特点是使用非线性变换关系进行图像变换。 对数变换实现了图像灰度扩展和压缩的功能。它扩展低灰度值而压缩高灰度值，让图像的灰度分布更加符合人的视觉特征。 幂次图像可以应用在对比度的操作中，对图3.8(a)进行幂次变换函数处理，其中，图3.8(b)到图3.8(f)的相应伽马值分别为0.6，0.4，0.3，1，2，从变换后的图像可以看出，当伽马值从0.6到0.4时，可看见更多细节；当伽马值进一步减少到0.3时，背景中的细节进一步增强，但图像看上去开始有轻微“冲淡”外观的那一点的对比度开始减小。当伽马值增加到大于1是，图像变暗。比较所有结果，可辨别细节的最好增强效果在γ=0.4时。 附录： % 图像增强 function imEnhancement f=imread(Fig3.08(a).jpg); imshow(f); c=255/log(256); x=0:255; y=c*log(1+x); figure,plot(x,y); [m,n]=size(f); f1=double(f); for i=1:m for j=1:n g(i,j)=c*log(f1(i,j)+1); end end q=mat2gray(g); figure, imshow(q); imwrite(q,Fig3.08(b).jpg) %此为指数的变换 f2=im2uint8(mat2gray((double(f1)+eps).^0.1)); figure,imshow(f2); imwrite(f2,Fig3.08(0.1).jpg) f3=im2uint8(mat2gray((double(f1)+eps).^0.4)); figure,imshow(f3); imwrite(f3,Fig3.08(0.4).jpg) f4=im2uint8(mat2gray((double(f1)+eps).^0.6)); figure,imshow(f4); imwrite(f4,Fig3.08(0.6).jpg) f5=im2uint8(mat2gray((double(f1)+eps).^1)); figure,imshow(f5); imwrite(f5,Fig3.08(1).jpg) f6=im2uint8(mat2gray((double(f1)+eps).^8)); figure,imshow(f6); imwrite(f6,Fig3.08(8).jpg)