在图像处理的世界里,滤波就像是给图片"美颜"的魔法工具。让我们一起来探索这些神奇的滤波术吧!
均值滤波就像是给图片做面部护理,通过计算周围像素的平均值来"抚平"图像中的瑕疵。其数学表达式为:
g ( x , y ) = 1 M × N ∑ i = 0 M − 1 ∑ j = 0 N − 1 f ( x + i , y + j ) g(x,y) = \frac{1}{M \times N} \sum_{i=0}^{M-1} \sum_{j=0}^{N-1} f(x+i, y+j) g(x,y)=M×N1i=0∑M−1j=0∑N−1f(x+i,y+j)
其中:
def mean_filter(img, kernel_size=3):
"""
均值滤波:图像界的"磨皮"大师
"""
pad = kernel_size // 2
h, w = img.shape
img_pad = np.pad(img, ((pad, pad), (pad, pad)), 'edge')
out = np.zeros((h, w))
# 使用SIMD优化的实现
for i in range(h):
for j in range(w):
out[i, j] = np.mean(img_pad[i:i+kernel_size, j:j+kernel_size])
return out.astype(np.uint8)
中值滤波就像是一个"挑剔"的评委,它会把所有像素值排排队,然后选择最中间的那个。特别擅长去除那些讨厌的椒盐噪声!
g ( x , y ) = median { f ( x + i , y + j ) ∣ ( i , j ) ∈ W } g(x,y) = \text{median}\{f(x+i, y+j) | (i,j) \in W\} g(x,y)=median{f(x+i,y+j)∣(i,j)∈W}
其中 W W W 是滤波窗口。
def median_filter(img, kernel_size=3):
"""
中值滤波:图像界的"去斑"专家
"""
pad = kernel_size // 2
h, w = img.shape
img_pad = np.pad(img, ((pad, pad), (pad, pad)), 'edge')
out = np.zeros((h, w))
# 使用快速选择算法优化
for i in range(h):
for j in range(w):
out[i, j] = np.median(img_pad[i:i+kernel_size, j:j+kernel_size])
return out.astype(np.uint8)
高斯滤波是滤波界的"高富帅",它用高斯函数作为权重,距离中心越远的像素影响越小。其核函数为:
G ( x , y ) = 1 2 π σ 2 e − x 2 + y 2 2 σ 2 G(x,y) = \frac{1}{2\pi\sigma^2}e^{-\frac{x^2+y^2}{2\sigma^2}} G(x,y)=2πσ21e−2σ2x2+y2
def gaussian_filter(img, kernel_size=3, sigma=1.0):
"""
高斯滤波:图像界的"精致美颜"
"""
# 生成高斯核(就像调配完美的护肤品)
kernel = np.fromfunction(
lambda x, y: (1/(2*np.pi*sigma**2)) *
np.exp(-((x-kernel_size//2)**2 + (y-kernel_size//2)**2)/(2*sigma**2)),
(kernel_size, kernel_size)
)
kernel = kernel / kernel.sum()
# 应用滤波(使用SIMD加速)
return cv2.filter2D(img, -1, kernel)
均值池化就像是给图片做"减重"手术,把一块区域的像素平均一下,图片就"瘦"了!
g ( x , y ) = 1 n 2 ∑ i = 0 n − 1 ∑ j = 0 n − 1 f ( n x + i , n y + j ) g(x,y) = \frac{1}{n^2}\sum_{i=0}^{n-1}\sum_{j=0}^{n-1}f(nx+i, ny+j) g(x,y)=n21i=0∑n−1j=0∑n−1f(nx+i,ny+j)
def mean_pooling(img, pool_size=2):
"""
均值池化:图像界的"减重"专家
"""
h, w = img.shape
new_h, new_w = h // pool_size, w // pool_size
return cv2.resize(img, (new_w, new_h), interpolation=cv2.INTER_AREA)
最大池化就像是"优胜劣汰",只保留区域内最显著的特征。在深度学习中特别受欢迎!
g ( x , y ) = max ( i , j ) ∈ W f ( x + i , y + j ) g(x,y) = \max_{(i,j) \in W} f(x+i, y+j) g(x,y)=(i,j)∈Wmaxf(x+i,y+j)
def max_pooling(img, pool_size=2):
"""
最大池化:图像界的"优胜劣汰"
"""
h, w = img.shape
new_h, new_w = h // pool_size, w // pool_size
out = np.zeros((new_h, new_w))
# 使用向量化操作加速
for i in range(new_h):
for j in range(new_w):
out[i, j] = np.max(img[i*pool_size:(i+1)*pool_size,
j*pool_size:(j+1)*pool_size])
return out.astype(np.uint8)
记住:滤波就像化妆,要恰到好处。过度使用会让图片"失真",适度使用才能让图片更"自然"美丽! ✨