clc; clear all; close all
% 图像的预处理，读入彩色图像将其灰度化
PS = imread('C:\Users\Administrator\Desktop\11.JPG'); % 读入彩色图像文件
figure;
subplot(2, 2, 1); imshow(PS); title('输入的图像')
PS=rgb2gray(PS); % 灰度化后的数据存入数组
% 绘制直方图
[m,n] = size(PS); % 测量图像尺寸参数
GP = zeros(1,256); % 预创建存放灰度出现概率的向量
for k = 0:255
GP(k+1) = length(find(PS==k))/(m*n);
% 计算每级灰度出现的概率，将其存入GP中相应位置
end
subplot(2, 2, 2); bar(0:255, GP, 'g'); title('原图像直方图')
xlabel('灰度值'); ylabel('出现概率')
% 直方图均衡化
S1 = zeros(1,256);
for i = 1:256
for j = 1:i
S1(i) = GP(j)+S1(i);
% 计算Sk
end
end
S2 = round((S1*256)+0.5); % 将Sk归到相近级的灰度
for i=1:256
GPeq(i) = sum(GP(find(S2==i))); % 计算现有每个灰度级出现的概率
end
subplot(2, 2, 3); bar(0:255, GPeq, 'b'); title('均衡化后的直方图')
xlabel('灰度值'); ylabel('出现概率')
% 图像均衡化
PA = PS;
for i = 0:255
PA(find(PS==i)) = S2(i+1);
% 将各个像素归一化后的灰度值赋给这个像素
end
subplot(2, 2, 4); imshow(PA); title('均衡化后图像')
|
QQ好友和群
QQ空间
腾讯微博
腾讯朋友
收藏
淘帖
顶
踩
