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相等的元素的个数 end for j = 0:255 for k = 0:j Formula(1 , j+1) = ZhiFang(1 , k+1) + Formula(1 , j+1);%计算新的像素值信息 end end for i = 1:rows for j = 1:cols NewGrayImage(i , j) = round(Formula(1 , (GrayImage(i , j) + 1)) * 255);%得到新的图像信息 end end for i = 0:255 NewZhiFang(1 , i+1) = length(find(NewGrayImage == i)) / (rows * cols);%计算矩阵NewGrayImage中值与i相等的元素的个数 end imwrite(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('出现次数');