ch45.m.txt

灰度图象的直方图均衡化；梯度；直方图规定化；灰度图象的直方图数据；灰度直接变换:增强对比度

%灰度图象的直方图数据
clc
clear
fn='rice.tif';
I=imread (fn);
J=I;
%计算灰度图象的直方图数据
L=256;                        %灰度级  
Ps = zeros(1,L);          %统计直方图结果数据
nk=zeros(L,1);
[row,col]=size(I);
n=row*col;                        %总像素个数
for i = 1:row
    for j = 1:col
        num = double(I(i,j))+1;  %获取像素点灰度级
        nk(num,1) = nk(num,1)+1; %统计nk  
    end
end
%计算直方图概率估计
Ps=nk/n;                         %矩阵除以标量
figure(1),subplot(3,1,1);
imshow(J),title('原图');
subplot(3,1,2),plot(nk),title('统计直方图（nk）');
subplot(3,1,3),plot(Ps),title('直方图（Ps）');

% 
% figure(2),subplot(3,1,1),imhist(I),title('统计直方图');
% p=imhist(I);         %生成直方图函数,返回灰度图象直方图数据
% subplot(3,1,2),plot(p);
% title('统计直方图（imhist）');
% subplot(3,1,3),plot(p/n);
% title('直方图（imhist）');