直方图均衡化技术实现图像的灰度归一化

I=imread('D:\project1.jpg');
I=im2double(I);
[M,N]=size(I);
figure(1);
imshow(I);
title('source');
figure(2);
[H,x]=imhist(I,32);
stem(x,(H/M/N),'.');
title('source1');
%增加对比度
A=2;B=-55;
O=A.*I+B/255;
figure(3);
subplot(2,2,1);
imshow(O);
figure(4);
subplot(2,2,1);
[H,x]=imhist(O,32);
stem(x,(H/M/N),'.');
figure(5);
subplot(2,2,1);
imshow(histeq(O));
figure(6);
subplot(2,2,1);
[H,x]=imhist(histeq(O,32));
stem(x,(H/M/N),'.');
%减小对比度
A=0.5;B=-55;
O=A.*I+B/255;
figure(3);
subplot(2,2,2);
imshow(O);
figure(4);
subplot(2,2,2);
[H,x]=imhist(O,32);
stem(x,(H/M/N),'.');
figure(5);
subplot(2,2,2);
imshow(histeq(O));
figure(6);
subplot(2,2,2);
[H,x]=imhist(histeq(O,32));
stem(x,(H/M/N),'.');
%增加亮度
A=1;B=55;
O=A.*I+B/255;
figure(3);
subplot(2,2,3);
imshow(O);
figure(4);
subplot(2,2,3);
[H,x]=imhist(O,32);
stem(x,(H/M/N),'.');
figure(5);
subplot(2,2,3);
imshow(histeq(O));
figure(6);
subplot(2,2,3);
[H,x]=imhist(histeq(O,32));
stem(x,(H/M/N),'.');
%反相显示
A=-1;B=255;
O=A.*I+B/255;
figure(3);
subplot(2,2,4);
imshow(O);
figure(4);
subplot(2,2,4);
[H,x]=imhist(O,32);
stem(x,(H/M/N),'.');
figure(5);
subplot(2,2,4);
imshow(histeq(O));
figure(6);
subplot(2,2,4);
[H,x]=imhist(histeq(O,32));
stem(x,(H/M/N),'.');

进行直方图均衡化处理：