LoG濾波器

一. LoG濾波器

        LoG 即 高斯-拉普拉斯（Laplacian of Gaussian）的縮寫。為了防止拉普拉斯濾波器計算二次微分時使圖像的噪聲更加明顯，所以先使用高斯濾波器使圖像更加平滑，再使用拉普拉斯濾波器使圖像的輪廓更加清晰。       

LoG 濾波器，s 代表高斯濾波的標准差 ↑

 

---

二. 實驗：LoG濾波器的實現及用於圖像邊緣檢測

        實驗中，我們先將BGR圖像轉換為灰度圖像，再給灰度圖像添加高斯噪聲，再使用LoG濾波器對圖像進行處理，最后輸出處理結果

import cv2
import numpy as np

# Gray scale
def BGR2GRAY(img):
    b = img[:, :, 0].copy()
    g = img[:, :, 1].copy()
    r = img[:, :, 2].copy()

    # Gray scale
    out = 0.2126 * r + 0.7152 * g + 0.0722 * b
    out = out.astype(np.uint8)

    return out

# add gaussian noise
def gasuss_noise(image, mean=0, var=0.001):
    '''
        添加高斯噪聲
        mean : 均值
        var : 方差
    '''
    image = np.array(image/255, dtype=float)
    noise = np.random.normal(mean, var ** 0.5, image.shape)
    out = image + noise
    if out.min() < 0:
        low_clip = -1.
    else:
        low_clip = 0.
    out = np.clip(out, low_clip, 1.0)
    out = np.uint8(out*255)
    return out

# LoG filter
def LoG_filter(img, K_size=5, sigma=3):
    H, W = img.shape

    # zero padding
    pad = K_size // 2
    out = np.zeros((H + pad * 2, W + pad * 2), dtype=np.float)
    out[pad: pad + H, pad: pad + W] = gray.copy().astype(np.float)
    tmp = out.copy()

    # LoG Kernel
    K = np.zeros((K_size, K_size), dtype=np.float)
    for x in range(-pad, -pad + K_size):
        for y in range(-pad, -pad + K_size):
            K[y + pad, x + pad] = (x ** 2 + y ** 2 - sigma ** 2) * np.exp( -(x ** 2 + y ** 2) / (2 * (sigma ** 2)))
    K /= (2 * np.pi * (sigma ** 6))
    K /= K.sum()

    # filtering
    for y in range(H):
        for x in range(W):
            out[pad + y, pad + x] = np.sum(K * tmp[y: y + K_size, x: x + K_size])

    out = np.clip(out, 0, 255)
    out = out[pad: pad + H, pad: pad + W].astype(np.uint8)

    return out


# Read image
img = cv2.imread("../paojie.jpg")

# grayscale
gray = BGR2GRAY(img)
# add gaussian noise
gray = gasuss_noise(gray)
# LoG filtering
out = LoG_filter(gray, K_size=5, sigma=3)

# Save result
cv2.imwrite("out_g_g.jpg",gray)
cv2.imshow("result_g_g",gray)
cv2.imwrite("out.jpg", out)
cv2.imshow("result", out)
cv2.waitKey(0)
cv2.destroyAllWindows()

 

---

三. 實驗結果：

圖1：灰度圖像添加高斯噪聲后 ↑

 

圖2：LoG濾波器對圖1圖像處理后 ↑

 

        從圖像中可以觀察到，LoG濾波器不僅有效地抑制了高斯噪聲，還一定程度上增強了圖像的邊緣。

---

四. 參考內容：

　　https://www.jianshu.com/p/c70be0cccd46