%改進的全局直方圖均衡化
close?all;
clear;
clc;
tic;

img_name?=?‘T2.JPG‘;
Img?=?double（imread（img_name））/255;????????%?加載圖像
figure;
imshow（Img）;?
title（‘原圖像‘）;

figure;
I=rgb2gray（Img）;
imhist（I）%顯示原始圖像直方圖
ylim（‘auto‘）;
title（‘輸入圖像直方圖‘）

hsvImg?=?rgb2hsv（Img）;??
V=hsvImg（::3）;??%取圖像的第三個二維矩陣，即第三個通道V
[heightwidth]=size（V）;??
??
V?=?uint8（V*255）;??
NumPixel?=?zeros（1256）;??
for?i?=?1:height??
????for?j?=?1:?width??
????NumPixel（V（ij）?+?1）?=?NumPixel（V（ij）?+?1）?+?1;??
????end??
end??
??
??
ProbPixel?=?zeros（1256）;??
for?i?=?1:256??
????ProbPixel（i）?=?NumPixel（i）?/?（height?*?width?*?1.0）;??
end??
??
CumuPixel?=?cumsum（ProbPixel）;??
CumuPixel?=?uint8（255?.*?CumuPixel?+?0.5）;??
??
for?i?=?1:height??
????for?j?=?1:?width??
????????V（ij）?=?CumuPixel（V（ij））;??
????end??
end??
???
V?=?im2double（V）;?%把圖像數據類型轉換為double類型?
hsvImg（::3）?=?V;??
outputImg?=?hsv2rgb（hsvImg）;??
figure;
imshow（outputImg）;??
title（‘輸出圖像‘）;

figure;
II=rgb2gray（outputImg）;
imhist（II）%顯示原始圖像直方圖
ylim（‘auto‘）;
title（‘均衡后圖像直方圖‘）



toc;

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????文件?????????398??2009-08-12?22:43??MeanSquareError.m
?????文件????????1590??2017-05-31?12:34??MSEPSNR.m
?????文件?????????500??2009-08-12?22:44??PeakSignaltoNoiseRatio.m
?????文件????????1152??2017-06-01?17:47??AzHSV.m
?????文件????????1231??2017-06-01?17:47??AZZhifangtu.m
?????文件??????233946??2017-04-29?16:05??T4.bmp
?????文件??????199545??2011-08-19?14:08??T8.bmp
?????文件??????314868??2012-12-23?19:47??W8.jpg
?????文件??????595854??2012-12-10?15:50??W9.bmp
?????文件???????16065??2017-05-31?12:33??ZZ13.jpg
?????文件???????23423??2017-05-31?12:33??ZZ14.jpg
?????文件??????289174??2017-04-24?17:48??T0.bmp
?????文件??????279166??2017-04-24?17:48??T1.bmp
?????文件???????13244??2017-04-21?19:20??T2.jpg