close?all
clear?all
clc
%%本程序：實現《數字圖像處理》p71?圖3.15的對比度增強，通過直方圖灰度值的改變進行調節
%%通過冪次變換s=cr^gamma來實現亮度調節
%%s是變換后的新灰度值；r是原圖像的灰度值；c為正常數本程序默認為1；gamma本程序為變量，需設定
%%讀取圖像
I=imread（‘dz.jpg‘）;???????????????????????%讀取圖像
if?length（size（I））>2??????????????????????%判斷圖像是否為灰度圖
????I=rgb2gray（I）;
end
subplot（221）imshow（I）
title（‘灰度化后的圖像‘）

%%繪制直方圖
[mn]=size（I）;????????????????????????????%測量圖像尺寸參數
P=zeros（1256）;???????????????????????????%預創建存放灰度的向量
for?k=0:255
????P（k+1）=length（find（I==k））;????????????%循環獲取每個灰度值的像素個數，將其存入P中相應位置
end
subplot（222）
bar（0:255P）??????????????????????????????%繪制直方圖
title（‘原圖像直方圖‘）
xlabel（‘灰度值‘）
ylabel（‘像素個數‘）

%%對直方圖進行冪次變換
S1=zeros（