MATLAB中将灰度图像均衡化

GrayImage = imread('E:/image/matlab/MouseGray.png');%原始的灰度图像[rows , cols] = size(GrayImage);%得到图像矩阵的行、列数ZhiFang = zeros(1 , 256);%初始化一个矩阵，用来存储原来灰度图像的直方图NewZhiFang = zeros(1 , 256);%初始化一个矩阵，用来存储变化之后的灰度图像的直方图Formula = zeros(1 , 256);%初始化一个矩阵，存储经过变化之后的新的像素值信息NewGrayImage = uint8(zeros(rows , cols));%初始化一个矩阵，用来存储新的灰度图像for i = 0:255 ZhiFang(1 , i+1) = length(find(GrayImage == i)) / (rows * cols);%计算矩阵GrayImage中值与i相等的元素的个数endfor j = 0:255 for k = 0:j Formula(1 , j+1) = ZhiFang(1 , k+1) + Formula(1 , j+1);%计算新的像素值信息 endendfor i = 1:rows for j = 1:cols NewGrayImage(i , j) = round(Formula(1 , (GrayImage(i , j) + 1)) * 255);%得到新的图像信息 endendfor i = 0:255 NewZhiFang(1 , i+1) = length(find(NewGrayImage == i)) / (rows * cols);%计算矩阵NewGrayImage中值与i相等的元素的个数endimwrite(NewGrayImage , 'E:/image/matlab/ZhiFangMouseGray.png' , 'png');%绘制初始灰度图像figure(1);imshow(GrayImage);%绘制初始图像的直方图figure(2);bar(0:255 , ZhiFang , 'grouped');%根据bar绘制直方图xlabel('灰度值');ylabel('出现次数');%绘制变换后的灰度图像figure(3);imshow(NewGrayImage);%绘制变换后的图像的直方图figure(4);bar(0:255 , NewZhiFang , 'grouped');%根据bar绘制直方图xlabel('灰度值');ylabel('出现次数');