在這一篇里面介紹一下CLAHE算法的第二步對(duì)比度限制。

這個(gè)過(guò)程很簡(jiǎn)單，分為下面幾個(gè)步驟。

計(jì)算出來(lái)限制的閾值

將統(tǒng)計(jì)好的直方圖數(shù)據(jù)限制在0到閾值范圍內(nèi)。也就是將大于閾值的直方圖數(shù)據(jù)減去閾值，并將差值累計(jì)起來(lái)。

將累計(jì)的差值平均分給每個(gè)灰度。

來(lái)看一下參考的Python代碼：

def cl_hist(img, clip):
    h, w = img.shape
    n = np.zeros(256, np.uint32)
    for i in range(h):
        for j in range(w):
            n[img[i][j]] = n[img[i][j]] + 1


    t = 0
    for i in range(256):
        if n[i] > clip:
            t += n[i] - clip
            n[i] = clip
    t = t / 256
    n = n + t
    pk = n
    sk = np.zeros(256, np.float32)
    sk[0] = pk[0]
    for i in range(1, 256):
        sk[i] = sk[i - 1] + pk[i]
    sk = sk / (h * w)
    sk = sk * 255
    sk = np.around(sk)
    return sk

這個(gè)和之前的直方圖統(tǒng)計(jì)就多了中間的步驟，也就是下面這一段代碼：

    t = 0
    for i in range(256):
        if n[i] > clip:
            t += n[i] - clip
            n[i] = clip
    t = t / 256
    n = n + t

就是遍歷256個(gè)灰度的統(tǒng)計(jì)值，將大于閾值的灰度統(tǒng)計(jì)值就設(shè)置為閾值，并將兩者的差值進(jìn)行累加。最后將總的差值除以256，將最后的結(jié)果累加到每一個(gè)灰度的統(tǒng)計(jì)值上面。

來(lái)看一下最后的效果：

可以看到相對(duì)于沒(méi)有進(jìn)行限制對(duì)比度的效果要好很多。

下面這幅圖是沒(méi)有限制對(duì)比度閾值的結(jié)果：

那么這個(gè)閾值要怎么計(jì)算呢：

一般來(lái)說(shuō)我們采取這樣的方法來(lái)計(jì)算閾值，也就是通過(guò)圖像的大小，分塊的多少，和clip_limit 的值來(lái)計(jì)算出來(lái)。一般圖像大小，和分塊多少是不變的，也就是說(shuō)可以通過(guò)clip_limit的值來(lái)調(diào)整效果。

h, w = img.shape
block = 8
clip_limit = 2
clip = clip_limit * (h * w / (block * block) / 256)
clip = round(clip)

最后來(lái)看一下Verilog的實(shí)現(xiàn)，可以通過(guò)下面的代碼來(lái)實(shí)現(xiàn)對(duì)比度限制的功能。

審核編輯：劉清