1、CT成像原理与系统架构:X射线物理基础、CT成像数学原理(Radon变换与反投影)、CT系统整体架构(扫描架、检查床、控制台)、数据采集系统(DAS)概述
各位同学,欢迎来到《医疗CT影像链嵌入式系统设计》的第一课。
说实话,每次讲CT成像原理,我都觉得像是在讲一个「用X光给人体做切片」的故事。你想想看,我们怎么才能看到身体内部的结构?不是直接看,而是用射线穿透过去,再通过数学把信息「算」出来。嗯,这背后就是Radon变换和反投影的功劳。
1.1 X射线物理基础:穿透与衰减
CT成像的起点,是X射线。我个人习惯把X射线想象成「一束能穿透人体的光」。但光穿过人体后,能量会衰减——骨头衰减得多,软组织衰减得少,空气几乎不衰减。
这个衰减规律,遵循朗伯-比尔定律:
I = I₀ * e^(-μx)
其中:
- I₀:入射X射线强度
- I:出射X射线强度
- μ:线性衰减系数(不同组织不同)
- x:穿透路径长度
我在项目中遇到过一个问题:球管老化后,I₀会漂移。当时我们没注意,结果重建出来的图像有环形伪影。后来我学乖了,每次开机先做空气校准,把I₀实时测一遍。
1.2 CT成像数学原理:Radon变换与反投影
说白了,CT成像就是「先投影,再反投影」。
Radon变换:把物体在不同角度下的投影数据记录下来。每个角度得到一条一维投影曲线,所有角度合起来就是正弦图(Sinogram)。
数学上,Radon变换可以写成:
p(θ, s) = ∫∫ f(x,y) * δ(x cosθ + y sinθ - s) dx dy
其中:
- f(x,y):物体内部衰减系数分布
- θ:投影角度
- s:探测器位置
- δ:狄拉克函数(只取路径上的点)
反投影:把每个角度的投影数据「抹」回图像空间。但直接反投影会有模糊效应——你想想看,一个点光源反投影回去会变成星形。所以实际用的是滤波反投影(FBP),先滤波再反投影。
我记得刚学这部分时,总觉得反投影很抽象。后来我用一个简单例子:假设一个2x2的像素矩阵,从0°和90°分别投影,再反投影回来——手动算一遍,就全明白了。建议你也试试。
1.3 CT系统整体架构:三大件
一台CT机,说白了就三部分:扫描架、检查床、控制台。
| 部件 | 功能 | 嵌入式要点 |
|---|---|---|
| 扫描架(Gantry) | 承载X射线球管、探测器、旋转电机 | 高速旋转控制(0.3s/圈)、滑环数据传输 |
| 检查床(Couch) | 承载患者,精确进床 | 步进电机控制、位置反馈(编码器) |
| 控制台(Console) | 操作界面、图像重建、存储 | 实时操作系统、GPU加速、DICOM通信 |
我参与过一个项目,扫描架的滑环磨损导致数据传输丢包。嗯,这个问题很隐蔽——图像偶尔出现错位,但重启后又好了。后来我们在滑环通信协议里加了CRC校验和重传机制,才算彻底搞定。
1.4 数据采集系统(DAS)概述
DAS是CT的「耳朵」和「嘴巴」。它负责:
- 接收信号:探测器把X射线转换成电信号
- 放大与整形:信号太弱,需要低噪声放大器
- 模数转换:把模拟信号变成数字信号
- 数据传输:通过滑环或光纤传到重建计算机
DAS的设计难点在于:
- 速度:每秒钟要采集几千个通道的数据
- 精度:16位甚至20位的ADC,噪声要极低
- 同步:所有通道必须同时采样,不能有相位差
我个人习惯在DAS前端加一个「暗电流校准」步骤。为什么?因为探测器在没有X射线时也会有微小输出,这个底噪必须扣除。我曾经因为忽略这一步,导致低对比度图像出现假阳性病灶——差点误诊。
好了,这一章的内容就到这里。CT成像原理,说白了就是「用X射线穿透,用数学重建,用硬件实现」。下一章我们会深入讲探测器阵列和信号链设计,到时候再聊。