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₀实时测一遍。

避坑指南: 我曾经因为X射线管电压不稳定,导致重建图像出现条状伪影。后来在DAS前端加了实时剂量监测,才彻底解决。记住:CT图像质量,一半靠硬件,一半靠校准。

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),先滤波再反投影。

核心要点: Radon变换是CT的「编码」过程,反投影是「解码」过程。没有数学,就没有CT图像。

我记得刚学这部分时,总觉得反投影很抽象。后来我用一个简单例子:假设一个2x2的像素矩阵,从0°和90°分别投影,再反投影回来——手动算一遍,就全明白了。建议你也试试。

1.3 CT系统整体架构:三大件

一台CT机,说白了就三部分:扫描架、检查床、控制台。

部件 功能 嵌入式要点
扫描架(Gantry) 承载X射线球管、探测器、旋转电机 高速旋转控制(0.3s/圈)、滑环数据传输
检查床(Couch) 承载患者,精确进床 步进电机控制、位置反馈(编码器)
控制台(Console) 操作界面、图像重建、存储 实时操作系统、GPU加速、DICOM通信

我参与过一个项目,扫描架的滑环磨损导致数据传输丢包。嗯,这个问题很隐蔽——图像偶尔出现错位,但重启后又好了。后来我们在滑环通信协议里加了CRC校验和重传机制,才算彻底搞定。

注意: 扫描架的旋转精度直接影响图像质量。我曾经见过一个案例,因为编码器安装松动,导致重建图像出现几何失真。所以,机械装配和校准流程,一定要严格。

1.4 数据采集系统(DAS)概述

DAS是CT的「耳朵」和「嘴巴」。它负责:

  1. 接收信号:探测器把X射线转换成电信号
  2. 放大与整形:信号太弱,需要低噪声放大器
  3. 模数转换:把模拟信号变成数字信号
  4. 数据传输:通过滑环或光纤传到重建计算机

DAS的设计难点在于:

  • 速度:每秒钟要采集几千个通道的数据
  • 精度:16位甚至20位的ADC,噪声要极低
  • 同步:所有通道必须同时采样,不能有相位差

我个人习惯在DAS前端加一个「暗电流校准」步骤。为什么?因为探测器在没有X射线时也会有微小输出,这个底噪必须扣除。我曾经因为忽略这一步,导致低对比度图像出现假阳性病灶——差点误诊。

经验之谈: DAS的采样时钟抖动,是图像噪声的隐形杀手。建议用低相噪的晶振,并在FPGA里做时钟域同步。别问我怎么知道的——都是踩坑踩出来的。

好了,这一章的内容就到这里。CT成像原理,说白了就是「用X射线穿透,用数学重建,用硬件实现」。下一章我们会深入讲探测器阵列和信号链设计,到时候再聊。