直方图均衡化

function histEqualization(srcImage)
if length(size(srcImage)) > 2 ;
srcImage = rgb2gray(srcImage);
end
tempImage = double(srcImage);
subplot(1,4,1),imshow(uint8(tempImage));
[M,N] = size(tempImage);
tempSeq = reshape(tempImage,1,M*N);
X = unique(tempSeq);
[nf,f] = hist(tempSeq,X); %列出原始直方图
subplot(1,4,2),hist(tempSeq,X);
L = size(f,2);
nf = nf/size(tempSeq,2);
gf = zeros(1,L);
gf(1) = nf(1);
for i=2:L;
gf(i) = gf(i-1) + nf(i); %计算原始累积直方图
end
g = floor((max(tempSeq)-1)*gf+0.5);
for i=1:M;
for j=1:N;
tempImage(i,j) = g(f==tempImage(i,j)); %确定映射关系
end
end
subplot(1,4,3),imshow(uint8(tempImage));
tempSeq = reshape(tempImage,1,M*N);
X = unique(tempSeq);
subplot(1,4,4),hist(tempSeq,X); %计算新直方图