直方圖均衡化算法分為三個(gè)步驟,第一步是統(tǒng)計(jì)直方圖每個(gè)灰度級(jí)出現(xiàn)的次數(shù),第二步是累計(jì)歸一化的直方圖,第三步是計(jì)算新的像素值。
第一步:
for(i=0;i<height;i++){
for(j=0;j<width;j++){ n[s[i][j]]++; } } for(i=0;i<L;i++){
p[i]=n[i]/(width*height);
}
這里,n[i]表示的是灰度級(jí)為i的像素的個(gè)數(shù),L表示的是最大灰度級(jí),width和height分別表示的是原始圖像的寬度和高度,所以,p[i]表示的就是灰度級(jí)為i的像素在整幅圖像中出現(xiàn)的概率(其實(shí)就是p[]這個(gè)數(shù)組存儲(chǔ)的就是這幅圖像的歸一化之后的直方圖)。
第二步:
for(i=0;i<=L;i++){
for(j=0;j<=i;j++){ c[i]+=p[j]; } } c[]這個(gè)數(shù)組存儲(chǔ)的就是累計(jì)的歸一化直方圖。
第三步:
max=min=s[0][0];
for(i=0;i<height;i++){ for(j=0;j<width;j++){ if(max<s[i][j]){ max=s[i][j];
}else if(min>s[i][j]){
min=s[i][j];
}
} } 找出像素的最大值和最小值。
for(i=0;i<height;i++){
for(j=0;j<width;j++){ t[i][j]=c[s[i][j]]*(max-min)+min; } } t[][]就是最終直方圖均衡化之后的結(jié)果。
對(duì)于彩色的圖片來說,直方圖均衡化一般不能直接對(duì)R、G、B三個(gè)分量分別進(jìn)行上述的操作,而要將RGB轉(zhuǎn)換成HSV來對(duì)V分量進(jìn)行直方圖均衡化的操作。
|
|