一种用于扩大视野的ct图像重建方法
【技术领域】
[0001]本发明属于生物医学成像技术领域，涉及PET/CT设备成像，尤其涉及一种用于扩大视野的CT图像重建方法，可用于PET/CT设备中的图像融合与图像校正。
【背景技术】
[0002]CT是计算机体层扫描，英文(Computed Tomography)的缩写。CT是利用X射线对人体进行体层检查。PET (派特)是正电子发射计算机断层显像，英文(Positron Emiss1nTomography)的缩写，是一种进行功能代谢显像的分子影像学设备。
[0003]PET-CT的全称叫正电子发射断层及X线计算机体层摄影成像系统，能同时提供解剖显像和功能显像，是目前影像诊断技术中最为理想的结合，特别是在肿瘤的诊断、分期、疗效评估等方面发挥重要的作用。
[0004]所谓PET/CT，是将PET和CT两个设备有机地结合在一起，使用同一个检查床和同一个图像处理工作站。PET/CT同时具有PET和CT的功能，但它不是二者功能的简单叠加，而为PET与CT优势互补，形成1+1>2的功效。PET可以显示病灶病理生理特征，更容易发现病灶；CT可以精确定位病灶，显示病灶结构变化。PET/CT除了具备PET和CT各自的功能夕卜，其独有的融合图像功能将PET图像与CT图像融合，CT图像信息可用来对PET进行衰减校正，这也是目前PET/CT中最常用的衰减校正方法。
[0005]然而，PET与CT的视野口径(最大物理成像直径)是不一样的。由于几何尺寸限制，通常CT的视野口径小于PET的视野口径，于是当体胖病人的身体一部分如肩膀超出CT物理视野以外时，这时的CT图像不能合成正确的PET衰减校正系数，可能导致PET图像上出现伪影，干扰医生的诊断，尽管PET视野口径是足够大的。
[0006]由于硬件上的限制，不能再增加CT探测器的长度，即不能增大CT的物理视野口径。为解决上述问题，人们用软件方法来扩大CT有效成像口径。目前基本上有两类方法，第一类方法是利用CT不完备数据重建技术，第二类方法是目前新近研宄的利用已有的同机PET图像对缺失的CT投影数据进行补充重建的技术。现有方法中，此类基于PET图像的CT补充重建技术大多引入较多的人为假定。比如，利用PET信息求出CT缺失部分的轮廓后，这些方法需要使用合适的组织类型来填充这些截断部分，这些组织类型的数值的确定往往是先验假定的，而且是单一介质(通常为水)，如在文献[I]和[2]所记载的，这些方法能够消除大部分伪影，因为即使是用水代替人体的CT缺失组织，也比没有好，因为水与空气的物质密度是完全不同的，水与人体软组织则相差不算太大。然而，当CT缺失部分含有骨骼时，上述方法还是可能造成较大的误差。
[0007]现有技术由于不完备数据问题引起的图像缺失和图像伪影，特别不利于缺失的CT图像中的骨骼图像的恢复。现有技术在利用CT投影空间的边缘信息时，通常采用的是数值外推(一种插值技术，比如最近邻插值)技术。然而CT投影空间结构与普通空间结构是不同的，比如在CT投影空间中很近的两个亮点，在普通空间中可能相距很远，也可能相距很近，这与扫描角度和物体所处的位置有关，因此，传统插值方法会带来新的误差。
[0008]引用文献：
[0009][I]中国发明专利，200980117410. 1，使用未经衰减校正的PET发射图像补偿不完备解剖图像，皇家飞利浦电子股份有限公司，2009。
[0010][2]Lonn A, Hsieh J, Evaluation of CT Field of View Restoration for PET-CTAttenuation Correction, 2006IEEE NSS Conference Record, Ml1-347.
【发明内容】
[0011]为了克服上述现有技术的不足，本发明提供一种用于扩大视野的CT图像重建方法，利用PET图像信息对其CT视野口径（即最大物理成像直径）以外的CT图像进行补充重建，从而克服由于不完备数据问题引起的图像缺失和图像伪影，可用于PET/CT设备中的图像融合与图像校正。
[0012]本发明的原理是：本发明提供的用于扩大视野的CT图像重建方法是一种利用已经获得的PET图像信息对缺失的CT图像进行补充重建的方法，属于【背景技术】所述第二类方法，即利用已有的同机PET图像对缺失的CT投影数据进行补充重建的方法。本发明提供方法不仅能利用PET图像轮廓信息，还充分利用了 CT断层投影空间图像（称为sino图）的边缘信息，对缺失的CT图像的形态和物质属性进行很好的补充重建，较好地克服由于不完备数据问题引起的图像缺失和图像伪影；本发明利用了一种CT投影空间图像边缘信息对缺失CT的sino图像进行逐曲线赋值的方法，最大限度利用CT投影空间边缘处的数值信息，能够实现CT投影空间边缘处的数据流连续性条件，避免了先验性的单一介质假定，较好地解决了现有方法中的投影空间数据流的连续性问题等技术难点。本发明提供方法尤其有利于缺失的CT图像中的骨骼图像的恢复，可更好地用于PET/CT设备中的图像融合与图像校
IHo
[0013]本发明提供的技术方案是：
[0014]一种用于扩大视野的CT图像重建方法，利用PET图像信息对CT视野口径以外的CT图像通过逐曲线赋值的方法进行补充重建，获得扩大视野的CT图像，包括如下步骤：
[0015]I)以PET视野口径为成像区域，比较PET图像与CT图像，找出缺失CT图像；记录缺失CT图像的所有像素点；PET视野口径比CT的视野口径大；
[0016]缺失CT图像的所有像素点设为（i，j)，其中，i = 1,2,3-1, j = 1，2，3…I，I为这些像素的个数；（i，j)是全体缺失图像像素点的下标，表示该像素点在图像空间中的位置。
[0017]2)针对步骤I)的每一个空间图像像素点，通过CT前投影方法计算得到与像素点相对应的一条正弦投影曲线；该正弦投影曲线即二维投影空间图像（sino图）中的一条曲线；
[0018]针对每一个空间图像像素点（i，j)，得到对应的一条正弦投影曲线，该正弦投影曲线即数值化二维图像中的一个连续的一维点阵，这些一维点阵用二维投影空间图像（sino图）中的一维子集表示为：m = m(n);其中，m为CT探测器阵列编号，η = 1，2，3…na，na为投影角度总个数，（m，η)表示二维投影空间图像上任一点。由于以较大的PET视野口径为成像区域，对应的扩大视野的sino图（行数）也增大了，即相当于有效地在原来CT探测器两端又增加了若干探测器单元，这些正弦投影曲线在扩大视野的sino图空间中表现为一条完整的周期正弦曲线。
[0019]3)通过计算得到上述正弦投影曲线与CT图像的断层投影空间图像（即sino图）中对应于探测器阵列编码截断边界上的相交叉点的位置，并累积交叉次数N;相交叉点的位置的sino图数值（‘图像值’)为该处的线性衰减系数线积分或CT投影值；
[0020]累积交叉次数N越大，则表示该点的线积分包含了更多的缺失图像像素的贡献，相交叉点的投影空间图像像素位置用（mO，nO)表示，其中，mO为CT探测器阵列编号，nO为投影角度，该点的sino图的‘图像值’即为投影角度nO下，编号mO的探测器单元接收到的线性衰减系数的线积分；相交叉点的位置是CT的sino图截断边界所对应的那一条直线与上述正弦投影曲线的交叉点，(CT的sino图也是一种数字图像，可表示为一个二维矩阵，行m为探测器编号，列η为投影角度编号，sino图的像素用（m，η)标注，其截断边界是一条‘直线’，所在行是对应最边缘的探测器标号，所在列为所有投影角度，sino图在像素点（m，n)上的数值（‘图像值’)为m编号探测器单元在η角度下的线性衰减系数线积分或CT投影值）。
[0021]4)将步骤3)得到的相交叉点的位置上的线积分值（交叉点处的sino图的‘图像值’）被除以相应的交叉次数N，得到值K ;在投影空间中将值K对经过该相交叉点的正弦投影曲线上的所有投影空间像素点进行累积赋值，即对该正弦投影曲线上的所有点（构成的二维投影空间即sino图中的一维正弦投影曲线点阵的那些像素集合）的投影空间图像值均赋予数值K ;累积赋值指的是：若某一投影空间像素点被多次赋值，则其取值为每次赋值的值K的加和；对正弦投影曲线上的所有点进行累积赋值后，得到赋值后的正弦曲线；
[0022]5)将赋值后得到的正弦曲线（一维点阵）的超出原视野之外的部分组合为一个二维矩阵；该二维矩阵的列和行分别为投射角度和CT探测器阵列编号；
[0023]6)将步骤5)组合成的二维矩阵与CT原sino图（二维矩阵）的探测器阵列进行拼接，得到补充后的sino图；
[0024]7)对步骤6)得到的补充后的sino图重建CT图像，得到扩大视野的CT断层图像，此时成像口径为PET 口径，大于原来CT 口径。
[0025]针对上述用于扩大视野的CT图像重建方法，在本发明实施例中，步骤I)比较人体的PET图像与CT图像，找出缺失CT图像的方法具体是：首先从PET图像获取人体的未做衰减校正的重建断层图，经过边缘增强和数值图像处理中的膨胀与腐蚀方法，利用阈值将PET人体图像进行二值化，从而确定人体的边界。上述图像处理方法同样适用于进行像素尺寸刻度后的CT图像（此时得到的CT图像与PET图像同样大小），以较大的PET视野口径为成像口径，对比PET图像与CT图像，可检测到有CT图像缺失的一些区域，CT图像缺失多发生在肩部。
[0026]步骤2)中的CT前投影方法为Radon变换方法。在本发明实施例中，采用的CT前投影方法具体是在某个视角下，在连接CT光源点与某一个CT探测器单元的连线上的线性衰减系数的路径积分，构成二维投影空间（Sino图）中的一个点，其位置行列编号分别为投射视角角度与CT探测器阵列编号；对于一个图像空间中的像素，其二维投影空间中的图像表现为为一条正弦投影曲线。
[0027]与现有技术相比，本发明的有益效果是：
[0028]本发明提供一种用于扩大视野的CT图像重建方法，该方法利用已经获得的PET图像信息对缺失的CT图像进行补充重建，通过利用CT投影空间边缘信息对缺失CT的Sino图像进行逐曲线赋值的方法，最大限度利用CT投影空间边缘处的数值信息，避免了先验性的单一介质假定，较好地解决了该类算法中的投影空间数据流的连续性问题等技术难点，即不仅能利用PET图像轮廓信息，还充分利用了