一. 拉普拉斯滤波器简介:
我们知道:
拉普拉斯算子 ↑
x方向上二阶偏导数的数值近似计算 ↑
y方向上二阶偏导数的数值近似计算 ↑
拉普拉斯算子在平面内的数值近似 ↑
拉普拉斯滤波器卷积核表示 ↑
二. 3*3的laplacian滤波器实现
# laplacian filter
def laplacian_filter(img, K_size=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] = img.copy().astype(np.float)
tmp = out.copy()
# laplacian kernle
K = [[0., 1., 0.],[1., -4., 1.], [0., 1., 0.]]
# filtering
for y in range(H):
for x in range(W):
ou