图像降噪原理与实践 - 数码照片去噪完全指南

图像噪声的类型与机制 - 噪声为何产生

图像噪声是像素值偏离真实场景亮度的随机波动。理解噪声的成因和类型是选择正确降噪方法的前提。

噪声类型:

• 散粒噪声 (Shot Noise)：光子到达传感器的统计波动。服从泊松分布，与信号强度的平方根成正比。暗部更明显
• 读出噪声 (Read Noise)：传感器电路在读取信号时引入的电子噪声。与 ISO 无关，是固定的底噪
• 热噪声 (Thermal Noise)：传感器温度引起的暗电流。长曝光和高温环境下显著
• 固定模式噪声 (FPN)：像素间灵敏度差异导致的系统性噪声。可通过暗帧减除校正

ISO 与噪声的关系：提高 ISO 实际上是放大传感器信号，同时也放大了噪声。ISO 100 到 ISO 6400 噪声增加约 8 倍 (64 倍信号放大的平方根)。现代传感器的基础 ISO (如 ISO 100) 读出噪声约 2-5 电子。

噪声度量: SNR (信噪比) = 信号/噪声。通常以 dB 表示：SNR_dB = 20×log10(信号/噪声)。30dB 以上视觉上几乎无噪声，20dB 以下噪声明显可见。

补充：读出噪声的大小依赖于 ADC (模拟-数字转换器) 的精度与放大器的特性。天文摄影中会使用冷却 CCD 来抑制这类噪声。

经典降噪算法 - 滤波基础

经典降噪方法基于空间域或频率域滤波，通过平滑像素值来抑制噪声。核心挑战是在降噪和保持细节之间取得平衡。

均值滤波：用邻域像素的平均值替代中心像素。简单快速但严重模糊边缘。核大小越大降噪越强但模糊越严重。仅适用于均匀区域的轻微噪声。

高斯滤波：加权平均，权重按高斯分布衰减。比均值滤波稍好但仍模糊边缘。σ 控制平滑程度。OpenCV: cv2.GaussianBlur(img, (5,5), sigma)

中值滤波：用邻域像素的中值替代中心像素。对椒盐噪声 (脉冲噪声) 特别有效，能保持边缘锐度。OpenCV: cv2.medianBlur(img, 5)

双边滤波 (Bilateral Filter)：同时考虑空间距离和像素值差异的加权滤波。相似像素权重高，差异大的像素 (边缘) 权重低，实现边缘保持降噪。OpenCV: cv2.bilateralFilter(img, d=9, sigmaColor=75, sigmaSpace=75)

频率域方法：噪声通常集中在高频。低通滤波可去除噪声但也会模糊细节。维纳滤波 (Wiener Filter) 在频率域根据信噪比自适应调整滤波强度，理论上是 MSE 最优的线性滤波器。

补充：双边滤波作为 Photoshop “减少杂色” 滤镜的基础技术，被广泛使用。

非局部均值与 BM3D - 最先进的经典方法

非局部均值 (NLM) 和 BM3D 利用图像的自相似性，在保持细节的同时实现强力降噪，是深度学习之前的最先进方法。

非局部均值 (NLM, 2005)：对于每个像素，在整张图像中搜索相似的图像块 (patch)，用所有相似块的加权平均进行降噪。权重由块间相似度 (欧氏距离) 决定。关键参数：搜索窗口大小、块大小、滤波强度 h。OpenCV: cv2.fastNlMeansDenoisingColored(img, None, h=10, hColor=10, templateWindowSize=7, searchWindowSize=21)

BM3D (2007)：当前经典方法的巅峰。三步流程：

• 块匹配 (Block Matching)：为每个参考块找到相似块，组成 3D 数组
• 协同滤波 (Collaborative Filtering)：对 3D 数组进行 3D 变换 (通常为 DCT + Haar)，硬阈值/维纳滤波去噪，逆变换
• 聚合 (Aggregation)：将所有处理后的块加权平均回原位置

BM3D 性能：在 σ=25 高斯噪声下，PSNR 通常比双边滤波高 3-5dB。处理速度：1MP 图像约 2-5 秒 (CPU)。Python: pip install bm3d,bm3d.bm3d(noisy, sigma_psd=25/255)

局限性: NLM 和 BM3D 在高噪声 (σ>50) 下性能下降明显，且处理速度较慢。对于纹理丰富的区域，可能出现过度平滑。这些局限推动了 AI 降噪的发展。

补充：2000 年代出现的 Non-Local Means (NLM) 与 BM3D 大幅拓展了经典方法的极限，达到了 AI 出现之前的最高性能。Non-Local Means (NLM) 由 Buades 等人于 2005 年提出。BM3D (Block-Matching and 3D Filtering) 由 Dabov 等人于 2007 年提出：第 1 阶段把图像分块并将相似块分组 (Block-Matching)，第 2 阶段在三维变换域中进行协同滤波。面向彩色图像的扩展版本 CBM3D 也已存在。

AI 降噪 - 深度学习革命

自 2016 年以来，基于深度学习的降噪方法在质量上全面超越经典方法，尤其在高噪声和复杂纹理保持方面优势显著。

DnCNN (2017)：开创性的 CNN 降噪网络。学习噪声残差 (输入-输出=噪声) 而非直接学习干净图像。17 层 CNN，批归一化 + ReLU。在已知噪声水平下比 BM3D 高 0.5-1dB PSNR。

FFDNet (2018)：接受噪声水平图作为额外输入，单一模型处理不同噪声强度。支持空间变化的噪声 (如高 ISO 照片中暗部噪声更强)。

NAFNet (2022)：简化的网络架构，去除非线性激活函数，使用简单的通道注意力。在多个基准上达到最先进水平，同时计算效率高。

Restormer (2022)：基于 Transformer 的图像恢复网络。多尺度自注意力捕获长距离依赖，在降噪、去模糊、去雨等多任务上表现优异。

商业 AI 降噪工具:

• Adobe Lightroom AI 降噪: 2023 年集成，一键处理 RAW 文件噪声
• DxO PureRAW：专业 RAW 降噪，保持细节能力出色
• Topaz DeNoise AI：独立降噪软件，支持批量处理

补充：DnCNN (Denoising Convolutional Neural Network) 是 Zhang 等人于 2017 年提出的、专用于降噪的 CNN 架构，它采用残差学习 (Residual Learning)，让网络只估计噪声成分。NAFNet (Nonlinear Activation Free Network) 是 2022 年提出的架构，不使用非线性激活函数，结构十分简洁。Topaz DeNoise AI 则是可独立运行、也可作为 Photoshop 插件使用的 AI 降噪工具。

实用降噪工作流程 - 从拍摄到成片

最佳降噪效果来自拍摄阶段的噪声控制与后期降噪的结合。以下是完整的实用工作流程。

拍摄阶段的噪声控制:

• 尽可能使用低 ISO (基础 ISO)
• 充分曝光 (向右曝光/ETTR)：在不过曝的前提下尽量增加曝光，后期降低亮度。高信号=高 SNR
• 使用三脚架延长曝光时间替代提高 ISO
• 多帧堆叠：拍摄多张相同场景，平均后噪声降低 √N 倍

后期降噪顺序:

• RAW 阶段降噪 (最佳)：在 RAW 处理器中降噪，保留最多信息
• 色度降噪优先：人眼对色度噪声更敏感，可强力降噪而不明显损失细节
• 亮度降噪适度：过度降噪会丢失纹理，产生"蜡像"效果
• 降噪后再锐化：先降噪再锐化，避免锐化放大噪声

参数调整策略:

• 按区域调整：天空等平坦区域强降噪，纹理区域轻降噪
• 使用亮度蒙版：暗部强降噪，亮部轻降噪 (暗部噪声更严重)
• 100% 放大检查：在实际像素级别评估降噪效果

补充：拍摄阶段的噪声最小化：ISO 感度尽量压低，ISO 100-400 的范围最为理想；采用适度曝光 (ETTR: Expose To The Right)，让直方图偏右 (略亮) 曝光，可相对降低暗部噪声。RAW 显影阶段的降噪自由度远高于 JPEG，Lightroom、Capture One 等软件可以分别调节亮度噪声与色彩噪声。

编程实现降噪 - OpenCV 和 Python 实现

使用 OpenCV 和 Python 实现各种降噪算法的具体代码和参数调优指南。

OpenCV 降噪函数:

import cv2
import numpy as np
# 高斯滤波
gaussian = cv2.GaussianBlur(img, (5,5), 1.0)
# 双边滤波
bilateral = cv2.bilateralFilter(img, 9, 75, 75)
# 非局部均值 (彩色图像)
nlm = cv2.fastNlMeansDenoisingColored(img, None, 10, 10, 7, 21)

自适应降噪：根据局部噪声水平自动调整降噪强度。估计噪声：对 Laplacian 滤波结果取中值绝对偏差 (MAD),σ = MAD / 0.6745。然后根据 σ 设置 NLM 的 h 参数：h ≈ σ × 0.8-1.2。

批量处理流水线:

from pathlib import Path
for img_path in Path('input').glob('*.jpg'):
img = cv2.imread(str(img_path))
denoised = cv2.fastNlMeansDenoisingColored(img, None, 10, 10, 7, 21)
cv2.imwrite(f'output/{img_path.name}', denoised)

深度学习降噪部署：使用 ONNX Runtime 部署预训练模型。加载 NAFNet 或 Restormer 的 ONNX 模型，输入归一化到 [0,1] 的图像张量，输出降噪结果。GPU 推理 1080p 图像约 50-200ms。

补充：OpenCV 中的 Non-Local Means：用 cv2.fastNlMeansDenoisingColored(img, None, h=10, hForColorComponents=10, templateWindowSize=7, searchWindowSize=21) 即可应用 NLM 滤波。对 ISO 3200 左右的图像，h=10-15 是合适的范围。深度学习模型方面，可以在 PyTorch 或 TensorFlow 上使用预训练的降噪模型，只需下载 restormer 或 nafnet 的公开模型并执行推理，就能获得最先进的降噪效果；若有 GPU，即使 4K 图像也能在数秒内处理完毕。批量处理时，用 glob 获取文件列表，再以多进程 (multiprocessing.Pool) 并行处理，进度显示使用 tqdm，并保留处理日志以便质量管理。