图像通过幂次变换增强对比度 matlab

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（