4. 图像预处理与增强
细胞形态学成像分析的第一步不是分割，也不是测量，而是让图像“可被正确测量”。显微镜成像链路上的噪声、背景漂移、光照不均、低对比度和边缘模糊，会直接改变阈值分割的边界位置、连通域面积和周长，最终影响核质比、圆形度、纹理等形态参数的数值。因此，预处理与增强的目标不是把图像“修得好看”，而是在尽量不引入伪结构的前提下，恢复灰度与真实光强之间的对应关系，并让目标边界稳定可重复。
一、先判断“该不该处理”：成像缺陷的识别
拿到一张细胞图像，先做缺陷诊断，再决定处理链。常见缺陷与典型表现如下：
- 随机噪声：暗背景上出现细密颗粒，高斯噪声表现为灰度小幅波动，椒盐噪声表现为孤立黑白点。
- 背景不均：视野中央亮、四周暗（渐晕），或一侧亮一侧暗（照明斜射），导致同一细胞在不同位置灰度不同。
- 低对比度：细胞与背景灰度差小，直方图集中在窄区间，阈值分割时边界抖动。
- 边缘模糊：离焦或衍射导致边界过渡带变宽，面积测量偏大或偏小取决于阈值取在过渡带何处。
诊断可用直方图和背景剖面：沿无细胞的空白区域取一条水平线，观察灰度随位置的变化曲线。若曲线呈明显抛物线或斜坡，说明存在光照不均，必须先校正再分割。
二、去噪：在平滑与保边之间取舍
去噪的核心矛盾是：平滑噪声的同时会模糊边界。对细胞形态分析，边界保真优先于噪声彻底消除。
高斯滤波用二维高斯核对图像卷积，标准差 σ 控制平滑尺度。σ 越大噪声越少，但边界过渡带越宽。示例：对 8 位灰度图，σ=1.0 通常可抑制传感器噪声而不过度模糊；σ=2.0 以上会明显扩大边界过渡带，导致面积测量偏差增大。
中值滤波取邻域灰度中位数，对椒盐噪声特别有效，且比同尺寸高斯核更好地保留阶跃边缘。3×3 中值适合细颗粒噪声，5×5 适合较粗的孤立坏点，但会抹掉小于核尺寸的细胞突起。
非局部均值与双边滤波在平滑时按灰度相似性加权，能更好保边，但计算量大。实际工程中常用双边滤波：空间核 σ_s 控制邻域范围，灰度核 σ_r 控制“多大灰度差算同一边界”。σ_r 过大会把边界两侧混在一起，过小则去噪不足。
三、背景校正与光照不均校正
背景校正的目标是让无细胞区域的灰度趋于一致，通常做法是估计背景场并从原图中减去或除以它。
形态学背景估计：对图像做较大结构元的灰度开运算（先腐蚀后膨胀），结构元尺寸需大于最大细胞直径，得到只含背景的低频分量 B(x,y)，再用 I_corr = I − B + C 校正，C 为常数偏置（如取背景均值），避免出现负值。示例：细胞直径约 20 px 时，结构元取 35～45 px 的圆盘。
多项式拟合校正：在空白区域采样若干背景点，拟合二次曲面 B(x,y)=a x²+b y²+c xy+d x+e y+f，再逐像素减去。适合渐晕明显的明场图像，但采样点必须避开细胞。
白参考校正：若成像系统可采集无样本的空白视野 I_white 和暗场 I_dark，则校正公式为 I_corr=(I−I_dark)/(I_white−I_dark)×mean(I_white−I_dark)。这是最物理、最可重复的光照不均校正方式，适合荧光和明场定量。
校正后应检查：空白区域灰度标准差是否下降，细胞与背景的灰度差是否在视野各处保持一致。若校正后出现环状伪影，通常是结构元尺寸接近细胞尺寸所致。
四、对比度增强与锐化
对比度增强常用直方图拉伸和 CLAHE（限制对比度自适应直方图均衡）。全局拉伸把 [p_low, p_high] 映射到 [0,255]，p_low、p_high 可取 1% 和 99% 分位数，避免极端像素拉偏。CLAHE 把图像分块（如 8×8），每块做受限均衡，clipLimit 控制对比度放大上限，防止噪声被过度放大。对细胞图像，clipLimit 取 2.0～3.0、块大小 8×8 是常见起点。
锐化用非锐化掩模：I_sharp=I+k(I−G_σ∗I)，k 为锐化强度，σ 为高斯尺度。k 过大在边界两侧产生过冲（halo），会让阈值分割把边界向外推，面积偏大。建议 k≤1.0，σ 取 1～2 px，并在锐化后复查边界过渡带宽度。
五、处理链顺序与参数起点
综合上述分析，推荐的处理顺序为：光照/背景校正 → 去噪 → 对比度增强 → 锐化（可选）→ 分割测量。理由如下：光照校正在低频层面操作，若放在去噪之后，平滑核会改变背景梯度的估计；对比度增强会放大噪声，因此应在去噪之后；锐化会改变边界梯度，若放在分割前必须谨慎，且应在对比度增强之后，避免增强过程把锐化产生的过冲进一步放大。
一组可参考的示例参数（针对 8 位、细胞直径约 20 px 的明场图像）：形态学开运算结构元取 40 px 圆盘；高斯去噪 σ=1.0；CLAHE 块大小 8×8、clipLimit=2.5；非锐化掩模 k=0.8、σ=1.5 px。这些数值是合理起点，实际需根据图像分辨率和细胞尺寸调整。
六、如何判断预处理是否“合格”
预处理没有绝对最优，只有可重复、可验证。建议用以下三个检查点：
- 背景一致性：在视野四角和中央各取一块无细胞区域，计算灰度均值与标准差。校正后各区域均值差应小于背景标准差的 1 倍。
- 边界稳定性：对同一细胞，用阈值上下浮动 ±5 灰度级分割两次，面积变化应小于 3%。若变化过大，说明边界过渡带太宽或对比度不足。
- 参数可重复性：对同一批图像重复整条处理链，关键形态参数（如平均面积、圆形度）的批间变异系数应控制在 5% 以内。
最后需要强调：预处理不是独立的美化步骤，而是形态测量链的“定标”环节。任何一步引入的灰度偏移或边界展宽，都会以面积、周长、圆形度等参数的形式传递到最终统计结果。因此，每完成一步处理，都应回到空白区灰度一致性、边界过渡带宽度和直方图形状这三个可观测指标上做复查，确认没有引入伪结构，再进入分割与测量。