X=imread('I.bmp');
%X=rgb2gray(X);
figure,imshow(X);
title('原始图像');
Y=histeq(X);
figure,imshow(Y);
title('全局均衡化图像');
E=6.0;k0=0.3;k1=0.2;k2=0.3;
[m,n]=size(X);
X1=X;
Mg=sum(X1(:))/(m*n);
Vg=sum(sum((X1-Mg).^2))/(m*n);
for i=2:m-1
for j=2:n-1
X2=X1(i-1:i+1,j-1:j+1);
Ms=sum(X2(:))/9;
Vs=sum(sum((X2-Ms).^2))/9;
if(Ms<=k0*Mg&&k1*Vg<=Vs&&Vs<k2*Vg)
X1(i,j)=E*X1(i,j);
end
end
end
imwrite(X1,'x.bmp');
figure,imshow(X1);
title('直方图统计局部增强图像');
直方图统计局部增强
需积分: 50 115 浏览量
2018-05-10
15:54:41
上传
评论 3
收藏 434B RAR 举报
DIDADIDADD
- 粉丝: 16
- 资源: 7