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空間
騰訊微博
騰訊朋友
收藏
淘帖
頂
踩
