彩色图片直方图均衡化

 

1,统计每个像素出现的次数

2,计算每个像素出现的概率

3,累计求和

4,统一量纲

5,数值替换

6,呈现效果

%读取一个照片
I = imread('t01a44596706ed343dd.jpg');
subplot(2,3,1);
imshow(I),title('正常图像');
J = I/2;
subplot(2,3,2);
imshow(J),title('原图像');
[ROW,COL,H] = size(J);
%统计每个像素出现的次数
count = zeros(H,256);
for k = 1:H
    for i = 1:ROW
        for m = 1:COL
            count(k,J(i,m,k)+1) = count(k,J(i,m,k)+1)+1;
        end
    end
end
count = double(count);
f = zeros(H,256);
%每个像素出现的概率
for k = 1:H
    for i = 1:256
        f(k,i) = count(k,i)/(ROW*COL);
    end
end
 
%累计求和
for k = 1:H
    for i = 2:256
        f(k,i) = f(k,i-1) + f(k,i);
    end
end
%统一量纲
f = 255*f;
 plot(f),title('f');

%数值替换
for k = 1:H
    for i = 1:ROW
        for m = 1:COL
             J(i,m,k) = f(k,J(i,m,k)+1);
        end
    end
end
Q= uint8(J);
subplot(2,3,5);
imshow(Q),title('增强后的图片');
            
 
 

评论
添加红包

请填写红包祝福语或标题

红包个数最小为10个

红包金额最低5元

当前余额3.43前往充值 >
需支付:10.00
成就一亿技术人!
领取后你会自动成为博主和红包主的粉丝 规则
hope_wisdom
发出的红包
实付
使用余额支付
点击重新获取
扫码支付
钱包余额 0

抵扣说明:

1.余额是钱包充值的虚拟货币,按照1:1的比例进行支付金额的抵扣。
2.余额无法直接购买下载,可以购买VIP、付费专栏及课程。

余额充值