荧光补偿算法——液相芯片多色荧光检测
来自液相芯片项目的荧光补偿算法文档,阐述软件全矩阵补偿(Full Matrix Compensation)的数学原理和计算方法。 相关:流式细胞仪详细设计说明书-CLFCS-V1.0、半自动双激光流式细胞仪串口通信协议-液路上样单元命令全集
一、核心原理
多色荧光检测中,由于荧光发射光谱重叠,一个荧光通道的信号会”泄漏”到其他通道(串色/串扰)。荧光补偿的目的就是从实测的混合信号中,还原出每种荧光的真实强度。
二、软件全矩阵补偿步骤
- 在每个 PMT 上对于每一个荧光计算**“叠加荧光素”的比率**(即串扰系数/泄漏比)
- 通过矩阵运算去除或加上交叉干扰
三、数学模型
3.1 实测信号公式
PMT = K × FL + A
| 符号 | 含义 | 说明 |
|---|---|---|
| PMT | 实测荧光信号向量 | 每个通道探测器实际读到的信号值(含串扰) |
| FL | 真实荧光信号向量 | 每种荧光染料的真实发射强度(待求) |
| K | 荧光泄露矩阵 | K_ij 表示荧光 FL_i 泄漏到探测器 PMT_j 的比例 |
| A | 噪声/背景信号向量 | 暗电流、环境背景光等固有偏移 |
3.2 荧光泄漏的物理过程
- 每个荧光源 FL_i 的信号,大部分进入自己对应的探测器 PMT_i
- 小部分信号会泄漏到其他通道(K_ij 即泄漏比例)
- 所有通道的分配比例之和为 1
3.3 补偿(解串扰)公式
从实测信号反推真实荧光,即求逆运算:
FL = K⁻¹ × (PMT - A)
其中 K⁻¹ 是荧光泄漏矩阵 K 的逆矩阵,也称为补偿矩阵。
运算步骤:
- 实测信号 PMT 减去背景噪声 A,得到纯荧光贡献的观测值
- 乘以串扰矩阵的逆矩阵 K⁻¹,解混得到各通道真实荧光强度 FL
四、两色荧光补偿示例(FITC / PE)
以 FITC(FL1)和 PE(FL2)两色荧光为例:
4.1 实测信号
PMT1 = FITC + x% × PE (PE 有 x% 泄漏到 PMT1)
PMT2 = PE + y% × FITC (FITC 有 y% 泄漏到 PMT2)
4.2 补偿矩阵(K 的逆矩阵)
1 [ 1 -x% ]
补偿矩阵 = ———— ×
1 - x%y% [ -y% 1 ]
4.3 补偿后的真实荧光
x% × PMT2 - PMT1
FITC = ————————————————
x% × y% - 1
y% × PMT1 - PMT2
PE = ————————————————
x% × y% - 1
五、关键概念定义
| 名词 | 定义 |
|---|---|
| PMT | 实测荧光信号(探测器读数) |
| FL | 待求的真实荧光矩阵 |
| A | 噪声/背景信号 |
| K | 荧光泄露矩阵,k_ij 是荧光信号 FL_i 泄露到荧光检测 PMT_j 的荧光泄露比 |
| K⁻¹ | 补偿矩阵(泄漏矩阵的逆矩阵) |
六、工程实现要点
- 串扰系数标定:需要用单荧光素对照样品(compensation beads/单染管)预先测定每个荧光对各通道的泄漏比例,构建 K 矩阵
- 背景扣除:补偿前必须先扣除探测器暗电流和环境背景 A
- 矩阵求逆:N 色荧光需要 N×N 矩阵求逆,两色用 2×2 解析解,多色需要数值求逆
- 扩展到多色:两色补偿可直接推广到 N 色(N 个荧光通道对应 N×N 的串扰矩阵),只需保证矩阵可逆(即各荧光光谱不完全线性相关)
关联
- 流式细胞仪详细设计说明书-CLFCS-V1.0 — 流式细胞仪软件整体设计,荧光检测是核心模块
- OLYMPUS FV1000共聚焦显微镜硬件用户手册 — 共聚焦中也有多通道荧光串扰问题
- 信号与系统习题解答手册-Oppenheim-第二版 — 线性系统理论基础