4. 图像预处理工程
显微图像在进入定量分析之前，几乎不可能直接使用。相机暗电流、荧光团自发闪烁、样品自发荧光、照明灯泡老化、载物台热漂移、盖玻片厚度差异，都会把真实信号叠加成一层“脏底”。预处理的目标不是把图像修得好看，而是让后续的阈值分割、面积测量、强度比较建立在可重复、可解释的数值基础上。本节按“噪声—背景—光照—漂移—批处理”五条链路展开，重点说明每一步的输入、操作、参数与结果判读。
一、去噪：先分清噪声类型，再选滤波器
显微图像常见噪声有三类：泊松散粒噪声（光子计数低时明显）、高斯读出噪声（相机电路）、以及椒盐/坏点噪声（传感器缺陷或宇宙射线）。泊松噪声与信号强度相关，不能简单用固定阈值去除；高斯噪声近似加性，适合线性滤波；坏点必须用中值或形态学方法剔除。
工程上常用三种策略：
- 高斯滤波：对高斯噪声有效，但会模糊边缘。示例参数：σ=1.0～1.5 px，核大小 5×5。σ 越大越平滑，但小结构（如 2 μm 以下的囊泡）会被抹掉。
- 中值滤波：对椒盐噪声和坏点最有效，核大小 3×3 通常足够；核过大（7×7 以上）会改变小颗粒的形态测量值。
- 非局部均值 / BM3D：保留边缘能力更强，适合低信噪比荧光图，但计算量大，批处理时需评估耗时。
判断去噪是否“过头”的实用方法：取同一视野的明场边缘（如细胞轮廓），比较去噪前后边缘的 10%～90% 上升宽度。若上升宽度增加超过 1.5 倍，说明结构已被明显模糊，应减小滤波强度。去噪不是越干净越好，而是“噪声下降、结构不丢”。
二、背景校正：把加性底从信号中减掉
背景校正处理的是加性成分：相机偏置、培养基自发荧光、玻片荧光。常用方法是“暗场参考法”或“形态学开运算估计背景”。以形态学为例：对图像做半径大于最大目标尺寸的开运算，得到低频背景估计 B(x,y)，再用 I_corrected = I_raw − B 得到校正图。若背景估计偏大，弱信号会被减成负值，因此通常再加一个常数偏移（如 +5 灰度级）避免截断。
示例：一幅 16-bit 荧光图，最大细胞核直径约 40 px，则开运算结构元半径取 50 px 以上。若半径取 20 px，背景估计会“吃掉”细胞核，导致核内强度被错误压低。判断依据：校正后背景区域均值应接近偏移常数，且目标区域强度不应出现大面积零值。
三、光照不均校正：乘性场的除法补偿
光照不均属于乘性误差：I_raw = I_true × S(x,y)，S 为照明场。校正方法是采集“空场”参考图（无样品、同物镜、同曝光），做平场校正：I_corrected = (I_raw − D) / (F − D) × mean(F − D)，其中 D 为暗场，F 为空场。若没有空场图，可用大尺度高斯低通估计 S，再做除法。
实操要点：空场图必须与样品图同曝光、同增益、同光源强度；若光源老化导致亮度漂移，空场图应在每批实验前后各采一次。校正后检查图像四角与中心的目标强度差异，若仍超过 10%，说明空场图不匹配或样品本身存在真实梯度。
四、漂移校正：时间序列的刚性/非刚性对齐
活细胞延时成像中，载物台热漂移和细胞运动会造成帧间位移。漂移校正分两步：先做刚性配准（平移+旋转），再做非刚性配准（局部形变）。刚性配准常用相位相关法：对相邻帧做 FFT，求互功率谱峰值位置得到平移量；旋转可用对数极坐标变换估计。非刚性配准则用光流或 B 样条自由形变。
判断配准质量：取一帧作为参考，计算配准后相邻帧的归一化互相关（NCC），NCC 应大于 0.95；若低于 0.9，说明形变过大或配准失败，需检查是否因细胞分裂、聚焦漂移导致。注意：漂移校正会引入插值模糊，因此应在去噪之后、定量测量之前进行，且不要反复重采样。
五、批处理策略：把上述步骤串成可复现流水线
单张图调参容易，批量处理难在“参数一致性”。推荐策略：先选 3～5 张代表性图像（含最暗、最亮、最不均）做参数标定，确定去噪 σ、背景半径、空场图、配准参考帧；然后将参数固化为配置文件，对整批数据按同一顺序执行：暗场/空场校正 → 去噪 → 背景减除 → 漂移配准 → 输出。顺序不可随意调换：若先去噪再平场，噪声会被除法放大；若先配准再背景校正，背景估计会随位移变化。
批处理必须记录每张图的处理日志：输入文件名、参数、输出文件名、NCC 值、背景均值。这样当某张图结果异常时，可以回溯是参数问题还是样本问题。建议每批保留一份“未处理原始图”和一份“处理后图”，便于复核。
六、常见误区与检查清单
预处理中最容易犯的错误有三类：一是顺序错误，例如先去噪再平场，导致噪声被除法放大；二是参数一刀切，不同曝光、不同荧光通道用同一组参数，弱通道被过度平滑；三是不留原始数据，校正后直接覆盖原图，无法回溯。建议每批处理前填写一份简表：样品名、物镜倍率、曝光时间、空场图文件名、去噪参数、背景半径、配准参考帧、NCC 阈值、输出目录。这份表不需要复杂，但能保证三个月后仍能复现同一批结果。
预处理不是“美化”，而是把物理成像过程引入的系统误差尽量剥离，让后续的阈值、面积、强度比值反映真实的生物学差异。去噪保结构、背景减加性、平场除乘性、配准对齐时间、批处理固化参数——这五条链路各司其职，顺序和参数一旦确定，就应作为实验记录的一部分固定下来。