第3讲 空间域图像增强.ppt

C、实例及结论 图4-24 噪声平滑实验图像 （a） Lena原图； （b） 高斯噪声； （c） 椒盐噪声； (d) 对(c)平均平滑；  (e) 对(b)平均平滑； (f) 对(b)5×5中值滤波； (g) 对(c)5×5中值滤波 a b c d e f g 4 直方图匹配（规定化） 直方图规定化 有时人们希望增强后的图像,其灰度级的分布不是均匀的， 而是具有规定形状的直方图,这样可突出感兴趣的灰度范 围。 令 和 分别为原始图像和期望图像的灰度概 率函数，对两者均作直方图均衡化处理，应有: 经过上述处理后的原图像概率密度函数 及理想图像 的概率密度函数 是相等的.于是可用变换后的原始 图像灰度级 代替式 中的 . 综上,组合变换函数为: 扩展到离散情况: 直方图均衡与匹配比较 直方图均衡与匹配比较 5 局部增强 6 在图像增强中使用直方图统计学 图像平均灰度值： 图像的方差： 图像局部平均值： 图像局部方差： 3.4 用算术/逻辑操作增强 1 算术操作 （1）加法操作 C(x,y) = A(x,y) + B(x,y) A、图像叠加（特效处理） B、图像平均处理（去除噪声） （2）减法操作 C(x,y) = A(x,y) - B(x,y) 减法的最主要作用是突出两幅图像的差异，常用于医学影像中的变化监测，或固定场景中的运动监测。 （3）乘法操作（乘以常数、模板操作等） 除法操作（一幅图像乘以另一幅图像的取反） （4）MATLAB实现 K1=imadd(I, J); %两幅图像加 K2=imsubtract(rice,50); %图像减一个常数 K3=immultiply(I, 2); %图像乘一个常数 K4=imdivide(I,2); %图像除一个常数 减法例子： rice=imread(‘rice.tif’); background=imopen(rice,strel(‘disk’,15)); rice2=imsubtract(rice,background); subplot(1,2,1),imshow(rice); subplot(1,2,2),imshow(rice2); 加法例子： I=imread(‘rice.tif’); J=imrea(‘cameraman.tif’); k=imadd(I,J); imshow(K); 乘法例子： I=imread(‘moon.tif’); J=immultiply(I,1.2); subplot(1,2,1), imshow(I); subplot(1,2,1), imshow(J); 除法例子： rice=imread(‘rice.tif’); I=double(rice); J=I*0.43+90; rice2=uint8(J); Ip=imdivide(rice,rice2); imshow(Ip,[ ]); 四则运算例子： I=imread(‘rice.tif’); I2=imread(‘cameraman.tif’); K=imdivide(imadd(I,I2),2); imshow(K); 2 逻辑操作（把灰度值作为二进制串） 随堂练习：39的“非”，和25的“与”及“或”。 （1）与、或 可用于从一幅图像中提取子图像。 （2）非 可以实现图像取反。 （3）异或