Ct扫描机模式、图像重建方法及高分辨ct扫描机的制作方法
【专利摘要】本发明涉及一种CT机扫描模式、图像重建方法及高分辨CT扫描机，所述扫描模式的扫描轨迹是在普通轴扫的同时，使辐射源和探测器单元与被扫描物体间沿旋转中心轴线发生定值的相对位移，旋转一周对应的轴向相对位移的大小等于1?2个所述探测器单元的宽度，所述图像重建方法的核心是采用公式进行反向投影插值，以考虑轴向相对位移。所述CT扫描机是用于执行所述扫描模式并设置了所述图像重建方法的CT扫描机。本发明能在基本上不改变现有硬件设备的前提下显著提高重建图像的纵向空间分辨率。
【专利说明】
CT扫描机模式、图像重建方法及高分辨CT扫描机
技术领域
[0001] 本发明涉及医疗器械技术领域，具体涉及一种CT机的新的扫描模式、与该扫描模 式相配套的扫描图像重建方法以及采用上述扫描模式和图像重建方法的CT扫描机。
【背景技术】
[0002] 目前CT机中要想提高轴扫模式得到的图像的纵向分辨率通常可以使用（1)高分辨 率探测器单元，简单说就是将每一个探测器单元在Z方向的尺寸做小一点；（2)飞焦点 (Flying Focal Spot)设计方法，这种方法来提高纵向空间分辨率对X-光源球管的设计有 较高的要求，从而导致制造成本的提高。事实上，这两种方法都会要求对CT机硬件设备有较 大的变动，成本会变高，而且也有其他方面各自的不足之处。
【发明内容】
[0003] 本发明的目的是提供一种CT机扫描模式、图像重建方法及高分辨CT扫描机，能在 不改变现有硬件设备的前提下，例如不提高对X-光源球管的设计要求的情况下显著提高重 建图像的纵向空间分辨率。
[0004] 本发明的主要技术方案是： 一种CT机扫描模式：以极小螺距(也可以称为微螺距)的螺旋运动路线为扫描轨迹对被 扫描物体进行断层信息记录(但只扫描一周），所述螺旋运动是指辐射源和探测器单元围绕 所述被扫描物体同步做定轴相对旋转的同时所述辐射源和探测器单元二者与所述被扫描 物体间轴向相对移动（即沿旋转中心轴线相对位移，旋转360度对应的轴向相对位移值等于 一个螺距），所述极小螺距为不大于2个所述探测器单元的宽度的螺距。本方案可以看成用 极小螺距的螺旋扫描(但只扫描一周）来近似普通轴扫。为描述方便起见，本文图中使用了 平板探测器来代替实际医用CT机中的柱面探测器，并将探测器投影到了旋转中心z轴上。
[0005] 所述辐射源的中心到所述探测器单元的垂线与所述旋转中心轴线垂直相交。
[0006] 所述相对旋转和所述轴向相对移动可以是连续的或者是步进的，且两者同步进 行。
[0007] 所述360度螺旋运动的实现方式可以是:所述辐射源和探测器单元旋转同时沿所 述旋转中心轴线平移，所述被扫描物体静止。
[0008] 所述360度螺旋运动也可以是:所述辐射源和探测器单元只旋转不平移，所述被扫 描物体只沿所述旋转中心轴线平移不旋转。
[0009] 所述极小螺距的大小优选等于1-2个所述探测器单元的宽度，这里讲的探测器单 元宽度指的是探测器单元投影到旋转轴z上的长度。
[0010] 所述极小螺距的大小优选等于1个所述探测器单元的宽度、1.5个所述探测器单元 的宽度或2个所述探测器单元的宽度。
[0011] -种CT机图像重建方法，用于前述任意一种所述的CT机扫描模式的图像重建，图 像中待重建的点在空间中的位置用笛卡尔坐标来标识，辐射源和探测器单元以z 轴为旋转中心轴线旋转，所述辐射源的中心到所述探测器单元的垂线与所述旋转中心轴线 的交点为所述笛卡尔坐标系的原点〇,（r, Φ,ζ)是与笛卡尔坐标二）相对应的柱坐标， 即尤《『〇^#，:3^^11(_，其特征在于包括反向投影步骤，采用以下公式计算反向投 影中的所需的坐标？：
其中，是辐射源和探测器单元绕Z轴旋转的旋转角； Η为螺距（即所述的极小螺距）； D是辐射源中心到ζ轴的距离； 所述CT机图像重建方法优选采用公式(1)获得重建图像：
其中，函数的)是CT机采集到的投影数据，甩t 表示，函数，(:)是将投影数据 芦0加权后再进行滤波得到的加权滤波后的投影数据；K，€和K的计算方法分别由公式 (2-4)给出。，
[0012] -种高分辨CT扫描机，包括辐射源、探测器单元和平台，还包括数据采集模块和图 像重建模块，所述平台为用于承载被扫描物体的可直线移动载体，所述辐射源和探测器单 元与所述平台之间通过采用所述的CT机扫描模式所对应的运动轨迹完成扫描，所述数据采 集模块采集所述探测器单元的投影数据，所述图像重建模块根据所述投影数据进行图像重 建，获得被扫描物体的断层平面图像，图像重建的方法采用所述的CT机图像重建方法。
[0013] 本发明的有益效果为： 和普通轴扫模式相比，相同之处在于都只扫描一周，不同之处在于本发明增加了一个 自由度，即在普通轴扫的基础上辐射源和探测器单元相对被扫描物体沿旋转中心轴线发生 一个精确的微小的直线位移（即所述微螺距H)，从而在数据中引入了更多的信息量，由此为 提高重建出来的图像的纵向分辨率提供了可能。也因此，本发明能够在基本上不改变现有 硬件设备的前提下，以及在不使用飞焦点设计的前提下，显著提高CT机重建出来的图像的 纵向空间分辨率。
[0014]本发明将会在需要较高纵向分辨率的案例中（如支架造影stent imaging,内耳 造影inner ear等等)显著提高图像的质量。
[0015]图像重建过程中，本发明采用含有微螺距Η的公式（3)来做反向投影 (Backprojection)中的插值，不仅不会导致重建图像产生伪影，还能显著提高纵向空间分 辨率。
【附图说明】
[0016] 图1是图像空间中一点的笛卡尔坐标与柱坐标也怎:)的对应关系图 示； 图2是本发明的扫描模式的扫描轨迹的几何坐标关系示意图； 图3是普通轴扫得到的模拟数据的重建图像； 图4是采用本发明的扫描模式且微螺距为2个探测器单元宽度时得到的模拟数据的重 建图像； 图5采用本发明的扫描模式且微螺距为1.5个探测器单元宽度时得到的模拟数据的重 建图像； 图6采用本发明的扫描模式且微螺距为1个探测器单元宽度时得到的模拟数据的重建 图像。
【具体实施方式】
[0017] 本发明公开了一种C T机扫描模式:C T机以螺旋运动路线为扫描轨迹对被扫描物体 进行断层信息记录，所述螺旋运动是指福射源（Source)和探测器单元(detector cell)围 绕所述被扫描物体同步做定轴相对旋转，在旋转的同时所述辐射源和探测器单元二者与所 述被扫描物体间沿旋转中心轴线发生精确的定值的相对位移。为了方便描述，以下简称本 发明的CT机扫描模式为微螺旋轴扫。
[0018] 所述辐射源可以包括X射线、单光子、正电子、超声、微波等。
[0019] 上述对于所述辐射源的描述进一步优选为对所述辐射源的中心的描述。
[0020] 所述辐射源的中心到所述探测器单元的垂线优选与所述旋转中心轴线垂直相交。
[0021] 所述相对旋转和所述轴向相对移动可以是连续的，也可以是步进的。要实现本发 明所述的CT机扫描模式，可以使所述辐射源和探测器单元旋转同时沿所述旋转中心轴线平 移，而保持所述被扫描物体静止;也可以使所述辐射源和探测器单元只旋转不平移，所述被 扫描物体只沿所述旋转中心轴线平移不旋转。由于所述被扫描物体通常是静置于CT机的平 台（c 〇uch)(或称床)上的，因此对所述床的运动控制就等于实现了对所述被扫描物体的运 动控制。
[0022] 作为进一步优化的技术方案，旋转360度对应的轴向相对位移的大小(微螺距的大 小)为1-2个所述探测器单元的宽度，例如，优选为1个、1.5个或者2个所述探测器单元的宽 度，以方便图像重建的运算。
[0023] 所述微螺旋轴扫与普通的CT机轴扫不同的是本发明的扫描模式可以理解为在普 通轴扫的同时还伴有沿旋转中心轴线的相对位移，比普通轴扫多一个自由度，这是本发明 能够在普通的CT机(无飞焦点功能设计的CT机)上提高纵向分辨率的基础。由于扫描模式的 上述不同，普通轴扫的重建方法中不体现螺距，如果简单地使用普通轴扫的重建方法（如 FDK方法)来处理本发明所产生的数据，重建的图像中将会出现伪影(artifacts)，因此普通 轴扫模式的图像重建方法不能用于本发明的所述微螺旋轴扫的图像重建。
[0024] 所述微螺旋轴扫与通常的螺旋扫描不同，包括:（1)所述微螺旋轴扫只需扫描一圈 (旋转360度），通常的螺旋扫描需要扫描多圈，这在相应重建算法公式中体现为前者的积分 限是2π，后者的积分限通常是z的函数(z轴是指与所述旋转中心轴线重合的笛卡尔坐标系 的坐标轴，本文所称的纵向也指z轴方向）。（2)所述微螺旋轴扫采用微螺距扫描，由于探测 器单元的宽度通常比较微小，至少与CT机上用于承载被扫描物体的床(或称平台）的人容易 感知到的宏观平移距离相比是极其微小的（由此，相应螺旋运动的螺距被称为极小螺距或 微螺距），因此说所述微螺旋轴扫的螺距比通常的螺旋扫描的螺距小很多。由于只扫描一周 且螺距极小，通常的螺旋CT的重建方法无法重建出普通轴扫模式中所能重建的图像层数。 综上，普通螺旋扫描的图像重建算法也不能应用于本发明所述微螺旋轴扫进行图像重建。
[0025] 为此，本发明还公开了一种与本发明任意一种扫描模式相配套的CT机图像重建方 法。/'是要重建的图像函数，图像中待重建的点在空间中的位置（即被扫描物体上某体素点 的空间位置)用笛卡尔坐标(為来标识，辐射源和探测器单元以z轴为旋转中心轴线旋 转。所述辐射源的中心到所述探测器单元的垂线与所述旋转中心轴线的交点为所述笛卡尔 坐标系的原点0。/(x,y.Z)在CT图像中指的是在坐标++诊处的Hounsfield Unit(HU) 值。在物理上，针对经过（Λ乂F, Z)点的每条来自辐射源的射线，被扫描物体上（Λ%乃Z)处 体素点都有一个对应的线性衰减系数(linear attenuation coefficient)，求解断层平面 (被扫描物体上垂直于z轴的截面）内某一体素点的就相当于求解对应于该断层 平面内所有经过该体素点的射线的线性衰减系数值。CT机图像重建的最终目的是得到
[0026] (r,也0是与笛卡尔坐标：_,:3以.)相对应的柱坐标，即r._cos(_，= ， 二者间关系参见图1。
[0027]本发明的重建方法优选基于FDK类型的解析算法，具体为：
其中，择是辐射源和探测器单元绕z轴旋转的旋转角（参见图2)，函数ρ.β是CT机采集 到的投影数据(即探测器单元上测得的数据），是(F,?)的函数。从整个空间上来讲，ΡΗ也是 菸的函数，因此也用祕約表示。投影数据的获得可以用图2简单表示。函数 PT)是加权滤波后的投影数据，即将投影数据Pi )加权后再进行滤波所得。
[0028] 参见图2，χ光透过人体后的投影数据被探测器Β记录下来。图2表示由光源发出的 经过扇面Α上的X光透过人体后的投影数据被探测器Β上的一行所记录，图中带箭头的线表 示扇面A中的一束X光线对应于探测器B中一个探测器单元的数据。探测器B上的所有行记录 了 X光在转角时的锥束投影数据。探测器B和光源同时绕z轴旋转就可以得到对应于一族 转角#的锥束投影数据。图像重建的任务就是要从这一族锥束投影数据中获得目标函数 /(? ，也就是人体的断层图像。
[0029] 在重建的过程中，对于每一个给定的图像空间坐标.?為爽游，需要在中找到 相应的每一个转角择时的?κ,Ο位置处的值，然后对所有的转角梦进行加和，即积分，这个过 程就是反投影(Backprojection)过程。
[0030] 给定(私次,破>，按照公式(2)、（3)确定(fJ):
其中，D是辐射源中心到z轴的距离(或者说辐射源中心到原点0的距离），H是本发明中 的所述螺距； 本发明的图像重建方法是针对本发明的扫描模式开发的新算法，其中采用考虑了螺距 Η的公式(3)做反向投影中的坐标?是本发明的图像重建方法的关键，这也正是采用本重建 算法不仅不会产生伪影反而可以提高纵向分辨率的主要原因。
[0031] Μζ ?2 ^τ2 2 r? slap -- φ) ； (4) 公式（1)-(4)所表示的算法是基于反投影重建算法的新算法，也属于FDK类型的解析算 法。事实上，可以采用任何的普通轴扫中其他FDK类型的方法（公式（1)只是其中的一个例 子)来处理本发明所产生的数据。需要注意的是无论用哪种类型的FDK算法，上面提到的反 投影过程中需使用(3)式来计算在旋转轴上的坐标?，这体现了本发明的算法的核心思想。 另外，在本发明的图像算法的核心思想框架下，目前所有的迭代算法都可以拿来套用。 [0032]为了验证本发明的创新思路，特进行了若干组模拟实验，图3-6所示为模拟实验数 据的结果，这些图像均为冠状图像(coronal image)，线条宽度和缝隙宽度均为1个探测器 单元宽度。通过上述模拟数据的重建图像结果可以看出本发明（图4-6分别对应了螺距为2、 1.5和1个探测器单元宽度的重建结果)相比于普通轴扫（图3)可以明显提高纵向分辨率，由 此验证了本发明中的思想和方法是行之有效的。
[0033] 因此，本发明的扫描模式和图像重建方法可以在不改变现有硬件设备的前提下， 提高重建出来的图像的纵向（即旋转中心轴线方向）分辨率。
[0034] 本发明还公开了一种高分辨CT扫描机，包括辐射源、探测器单元和平台，还包括数 据采集模块和图像重建模块，所述平台为用于承载被扫描物体的可直线移动载体，所述辐 射源和探测器单元与所述平台之间通过采用前述的任意一种CT机扫描模式，以相应的运动 轨迹完成扫描，所述数据采集模块采集所述探测器单元的投影数据，所述图像重建模块根 据所述投影数据进行图像重建，获得被扫描物体的断层平面图像，图像重建的方法采用前 述的任意一种CT机图像重建方法。
[0035] 该CT扫描机的辐射源和探测器单元围绕所述平台做同步的定轴旋转运动，转动范 围是一圈，二者可固定在同一装置上，也可以分别固定在不同的装置上，随各自所在装置做 旋转运动，在此过程中，所述平台沿旋转轴做同步的直线运动，所移动距离微小但精确，控 制在1-2个探测器单元的宽度范围内，形成所述螺旋运动。这些是实现高的纵向分辨率的扫 描图像的硬件基础。而内置了所述CT机图像重建方法所对应软件程序的图像重建模块是实 现高的纵向分辨率的扫描图像的必备软件条件。
[0036] 本发明就是要在不使用飞焦点的前提下，辐射源和探测器单元仍然遵循传统轴扫 模式下的轨迹运行，只要配合CT机的床沿z轴方向做一个微小的螺距位移以及采用配套的 图像重建方法进行图像重建，即可得到纵向空间分辨率明显高于CT机采用普通轴扫方式的 断层图像。
【主权项】
1. 一种CT机扫描模式，其特征在于w极小螺距的螺旋运动路线为扫描轨迹对被扫描物 体进行断层信息记录，所述螺旋运动是指福射源和探测器单元围绕所述被扫描物体同步做 定轴相对旋转的同时所述福射源和探测器单元二者与所述被扫描物体之间轴向相对移动， 所述极小螺距为不大于2个所述探测器单元的宽度的螺距。2. 如权利要求1所述的CT机扫描模式，其特征在于所述福射源的中屯、到所述探测器单 元的垂线与所述旋转中屯、轴线垂直相交。3. 如权利要求2所述的CT机扫描模式，其特征在于所述相对旋转和所述轴向相对移动 是连续的或者是步进的，且两者同步进行。4. 如权利要求3所述的CT机扫描模式，其特征在于所述福射源和探测器单元旋转同时 沿所述旋转中屯、轴线平移，所述被扫描物体静止。5. 如权利要求3所述的CT机扫描模式，其特征在于所述福射源和探测器单元只旋转不 平移，所述被扫描物体只沿所述旋转中屯、轴线平移不旋转。6. 如权利要求1、2、3、4或5所述的CT机扫描模式，其特征在于所述极小螺距的大小等于 1-2个所述探测器单元的宽度。7. 如权利要求6所述的CT机扫描模式，其特征在于所述极小螺距的大小等于1个所述探 测器单元的宽度、1.5个所述探测器单元的宽度或2个所述探测器单元的宽度。8. -种CT机图像重建方法，用于权利要求1-7中任意一项权利要求所述的CT机扫描模 式的图像重建，图像中待重建的点在空间中的位置用笛卡尔坐标（.尤躬句来标识，福射源 和探测器单元Wz轴为旋转中屯、轴线旋转，所述福射源的中屯、到所述探测器单元的垂线与 所述旋转中必轴线的交点为所述笛卡尔坐标系的原点〇，:紅,球，装I是与笛卡尔坐标 載句相对应的柱坐标，即龙CCM妹），?' sin轉)，其特征在于包括反向投影步 骤，采用W下公式做反向投影中的插值计算经过点(礙少V思;I的射线与平面B的交点的Z坐标 f .其中，巧是福射源和探测器单元绕Z轴旋转的旋转角； Η为螺距； D是福射源中屯、到Ζ轴的距离； 平面Β是过Ζ轴的与探测器单元平行的平面。9. 如权利要求8所述的CT机图像重建方法，其特征在于采用公式(1)获得重建图像：其中，函数狄，提CT机采集到的投影数据，用約祇欺表示，函数孩"巧是将投影数据 Ρ。)加权后再进行滤波得到的加权滤波后的投影数据；F，《和Κ的计算方法分别由公式2- 4给出；10. -种高分辨CT扫描机，包括福射源、探测器单元和平台，其特征在于还包括数据采 集模块和图像重建模块，所述平台为用于承载被扫描物体的可直线移动载体，所述福射源 和探测器单元与所述平台之间通过采用如权利要求1-7中任意一项权利要求所述的CT机扫 描模式所对应的运动轨迹完成扫描，所述数据采集模块采集所述探测器单元的投影数据， 所述图像重建模块根据所述投影数据进行图像重建，获得被扫描物体的断层平面图像，图 像重建的方法采用如权利要求8或9所述的CT机图像重建方法。
【文档编号】G06T7/00GK106097411SQ201610491671
【公开日】2016年11月9日
【申请日】2016年6月29日
【发明人】施大新, 邹宇
【申请人】沈阳开普医疗影像技术有限公司