matlab遥感影像均衡化,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('出现次数');