3d图像集在不同空间之间的最佳变换的制作方法
【专利摘要】本发明提供了一种3D图像集在不同空间之间的最佳变换。多模态系统（10）包括核成像系统（12）和计算机断层摄影（CT）扫描器（14）。核系统（12）包括PET扫描器（28），该PET扫描器（28）采集由PET重建处理器（50）重建成PET团块图像的电子数据。CT扫描器（14）采集由CT重建处理器（56）重建成3DCT体素图像的扫描数据。插值处理器（62）将PET团块图像直接插值到CT体素空间内。一旦PET和CT图像在同一空间内，它们就会由合并装置（110）合并。视频处理器（66）处理接收到的合成PET－CT数据以在显示器（68）上显示。
【专利说明】3D图像集在不同空间之间的最佳变换
[0001] 本发明涉及诊断成像系统和方法。其尤其与多模态系统，如PET - CT系统联合应 用。应当意识到，本发明还可应用于SPECT、CT、超声、MRI和荧光透视检查等的各种结合。
[0002] 在多模态断层摄影系统中，采用两个或多个不同感应模态对物体空间内的不同组 分进行定位和测量。在PET - CT系统中，PET创建身体高代谢活动的图像，而不是创建周 围解剖结构的图像。CT扫描允许医生观察人体的内部结构。在进行PET - CT扫描之前，患 者接受放射性药物剂量。药物通过血液运送并集中到特定器官或区域中，从而导致从血液 和该器官或区域发射出辐射。在扫描期间，射出辐射的轨迹由系统检测，由此创建患者体内 放射性药物的分布图像。该图像显示出循环体系和/或放射性药物在各种区域或器官内的 相对吸收。PET - CT图像中来自CT扫描的解剖结构数据与来自PET扫描的代谢数据的合 并向医生给出可视化信息以确定是否存在病变、病变的位置和程度，并且跟踪病变扩散的 速度。PET - CT系统对于难于治疗区域(例如，头颈区、纵隔、术后腹部）以及患者接受放射 性治疗或化疗的治疗区定位特别有帮助。
[0003] 典型地，在PET - CT图像的重建过程中，PET图像表示从PET图像空间变换到CT 图像空间中以创建用于临床解释的合成图像。在扫描前，PET和CT系统机械地对准；然而， 在大多数情况下，在融合之前需要施加微小校正以使PET图像和CT图像准确对准。目前， 变换参数在校准期间得出，如采用在两种模态中均可成像的模型。
[0004] 过去，PET和CT图像数据都被重建成3D体素矩阵，但是比例不同（不同尺寸的体 素)。较低分辨率PET图像和较高分辨率CT图像的笛卡尔网格之间的变换是相对直接的， 但是有时增加了变换或插值误差，从而导致精度和图像质量的下降。
[0005] 现在，PET图像通常重建到具有非笛卡尔网格的blob (团块)空间内。与全部是均 匀值的体素不同，每个团块内的值是不均匀的。此外，团块重叠。目前，重建的PET团块图 像首先变换成常规的基于体素的笛卡尔网格图像。由于PET体素与CT体素仍然是不同的 尺寸，因而将PET体素图像变换成CT体素尺寸。进行两次插值增加了固有变换相关误差。
[0006] 本发明提供一种新型和改进的成像设备和方法，其克服了上述问题和其它问题。
[0007] 根据本发明的一个方面，公开一种诊断成像系统。第一扫描器获得对象感兴趣区 域的第一组电子数据。第二扫描器获得对象感兴趣区域的第二组电子数据。第一重建装置 将第一组重建成非体素图像空间内的非体素空间第一扫描器图像表示。第二重建装置将第 二组重建成第二图像空间内的第二图像表示。一装置将第一组的非体素空间第一扫描器图 像表示直接转换到第二图像空间内。
[0008] 根据本发明的另一方面，公开一种诊断成像方法。获得对象感兴趣区域的第一组 电子数据。获得对象感兴趣区域的第二组电子数据。将第一组重建成非体素图像空间内的 非体素空间第一扫描器图像表示。将第二组重建成第二图像空间内的第二图像表示。将第 一组的非体素空间第一扫描器图像表示直接转换到第二图像空间内。
[0009] 本发明的一个优点在于提高总体图像质量。
[0010] 另一个优点在于减少计算和存储需求。
[0011] 另一个优点在于更精确地融合图像。
[0012] 在阅读和理解下列优选实施例的详细说明后，对本领域技术人员而言，本发明的 其它优点和有益效果将变得更为明显。
[0013] 本发明可采用各种部件和部件设置、以及各种步骤和步骤设置的形式。附图仅用 于解释说明优选实施例的目的，而不应解释为限制本发明。
[0014] 图1是多模态诊断成像系统的示意表示； 图2示意性示出体心立方网格； 图3示意性示出投影到笛卡尔网格图像空间内的2D团块； 图4是由平滑曲线函数表示的团块内灰度强度值的示意图；和 图5是多模态诊断成像系统一部分的示意表示。
[0015] 参见图1，多模态系统10包括核成像系统12和计算机断层摄影（CT)扫描器14。 CT扫描器14包括非旋转台架16。X射线管18安装在旋转台架20上。孔22形成CT扫描 器14的检查区域24。辐射检测器阵列26布置在旋转台架20上以在X射线穿过检查区域 24后接收来自X射线管18的辐射。可替换地，检测器阵列26可放置在非旋转台架16上。
[0016] 核成像系统12优选包括安装在轨道30上的正电子发射断层摄影（PET)扫描器 28。当然，也可考虑采用SPECT和其它核成像系统。轨道30平行于对象支撑台或床32的纵 轴延伸，从而使CT扫描器14和PET扫描器28形成封闭系统。提供移动装置34,如电机和 驱动器，以移动扫描器28进出该封闭位置。检测器36设置在形成检查区域40的孔38的 周围。在PET系统中，检测器36优选设置成静止的环状，虽然也可考虑采用可旋转头。在 SPECT系统中，检测器36优选合并在各个头中，其安装成可相对于患者转动或径向移动。床 移动装置42,如电机或驱动器，提供床32在检查区域24,40内的纵向运动和垂直调整。
[0017] 如下面所更详细描述的，PET扫描器28采集电子数据，该电子数据由PET重建处 理器或处理过程50重建成PET团块图像并存储到PET团块图像存储器52内。合适的团块 重建算法在Lewitt(1990)和Lewitt(1992)中有描述。
[0018] 参见图2和3,团块在概念上可以是每个以体心立方设置的角顶点或者体心点为 中心的球体。该球体比角顶点52的间距大，例如，直径等于角顶点与角顶点之间间距的2 - 1/2倍。在每个团块中，灰度(或强度值)在中心处最大并以球面对称的贝塞尔函数、高斯或 其它平滑或分段平滑曲线向周边逐渐衰减，如图4所示。
[0019] 再次参见图1，CT扫描器14采集扫描数据，该扫描数据由CT重建处理器或处理过 程56重建成3D CT体素图像并存储在CT图像存储器58内。CT数据优选重建成由3D笛 卡尔坐标系定义的矩形体素。典型地，CT体素远远小于PET角顶点间距，例如，1 一 3mm对 3 - 6mm。插值处理器或处理过程62将PET团块图像直接插值到CT体素空间内，S卩，具有 与CT图像相同的3D笛卡尔网格和体素尺寸的图像。一旦PET和CT图像处于同一空间，它 们就易于相加或者合并。所得到的合成图像存储在合成图像存储器64内。视频处理器66 处理接收到的合成PET - CT图像数据以在一个或多个监视器68上显示。
[0020] 继续参见图1，在扫描之前的某个时间，通过采用包含CT核标记72的模型70确定 成像系统10的CT和PET部分之间的变换矩阵。示范性CT核标记72包括盘状或球状致密 CT成像材料和位于其中心的放射性同位素标记。该标记可以由玻璃填充的特氟隆或其它合 适的材料构建而成，这些材料具有与模型明显不同的CT数，从而该标记在CT扫描期间易于 识别。典型地，放射性同位素标记是装有同位素的小瓶，该同位素具有相对较长半衰期，优 选大于100天，以防止频繁更换。对于PET - CT联合装置，优选的同位素为具有2. 6年半衰 期的Na-22。然而，也可同样采用具有较短半衰期的同位素，如Ge-68。对于PET - CT联合装 置，用于放射性标记的同位素选自半衰期大于1〇〇天且强能量峰值在50keV和600keV之间 的同位素。合适的SPECT可成像同位素实例为Go-57，Gd-153，Cd-109，Ce-139，Am-241， Cs-137 和 Ba-133。
[0021] 带有模型70的床32由床移动装置42移入到检查区域40内以由PET扫描器28 生成3D图像。PET重建处理器或处理过程50将电子数据重建成PET图像。确定每个放射 性同位素标记的质心坐标。
[0022] 接下来，床移动装置42移动床32以将模型70定位到CT扫描器检查区域24内， 在此进行CT成像。CT重建处理器56将电子数据重建成CT标记质心的3D CT图像。利用 PET和CT扫描器的已知几何形状和机械对准将为PET图像计算的放射性同位素标记位置坐 标平移到CT图像空间内。确定相应标记的坐标。接下来确定变换参数或矩阵，例如线性移 位量、比例、旋转和可选择的非线性平移，以使PET图像和CT图像准确对准。变换参数存储 在校准存储器80内并由对准处理器或装置82用于使PET和CT图像在后续扫描的插值处 理过程62中彼此对准。例如，对准装置82可向PET团块图像的角顶点和中心点施加确定 的变换以使该团块图像与CT图像对准。
[0023] 继续参见图1并进一步参见图5,带有物体的床32由床移动装置42移动到检查区 域40内以由PET扫描器28生成3D图像。PET重建处理器50将电子数据重建到团块空间 内。
[0024] 通常，3D扫描物体的连续分布f(x，y，z)可近似为基函数Φ (X)的伸缩和移位叠加 的副本的总和，基函数Φ (X)的中心设置在具有网格点P(xn，yn，zn)的网格上：
【权利要求】
1. 一种诊断成像系统，包括： 用于获得对象感兴趣区域的第一组电子数据的第一扫描器； 用于获得对象感兴趣区域的第二组电子数据的第二扫描器； 用于将第一组电子数据重建成非体素图像空间内的非体素图像空间第一扫描器图像 表不的第一重建装置； 用于将第二组电子数据重建成第二图像空间内的第二图像表示的第二重建装置； 用于将第一组电子数据的非体素图像空间第一扫描器图像表示转换到第二图像空间 的转换装置；以及 用于将第一组电子数据重建成团块的处理装置。
2. 根据权利要求1所述的诊断成像系统，其中所述团块布置在体心立方网格中并且 第二图像空间是笛卡尔体素图像空间。
3. 根据权利要求2所述的诊断成像系统，其中所述团块的尺寸设置成彼此重叠。
4. 根据权利要求1所述的诊断成像系统，其中第二图像空间是体素图像空间，并且该 转换装置包括： 用于确定体素图像空间的每个体素的中心点的第一确定装置；和 用于确定非体素图像空间图像表示内对应于每个所确定的体素中心点的点以定义体 素图像空间第一扫描器图像的第二确定装置。
5. 根据权利要求4所述的诊断成像系统，其中该非体素图像空间包括布置在非笛卡 尔网格上的团块，所述团块部分重叠，并且该转换装置还包括： 用于计算在非体素图像空间图像表示内重叠每个所确定的点的全部团块的总和以重 建合成体素图像f (X，y，z)的求和装置：
其中f(x，y，z)是合成3D图像表示，^^1(*是非体素图像空间图像表示，Φν° Μ?是体素 图像空间第一扫描器图像表示，η是非体素图像空间图像内的采样点数，m是体素图像内的 采样点数，Λ2是团块尺寸或采样间隔，Λ 1是体素图像空间内的体素尺寸或采样间隔，(^是 在每个采样点η处的图像系数，tm是在每个采样点m处的图像系数，采样点m取决于覆盖体 素图像空间的体素中心的团块数。
6. 根据权利要求4所述的诊断成像系统，还包括： 对准装置，用于向以下的一个施加事先确定的系统变换矩阵以使非体素图像空间图像 表示与体素图像空间图像表示彼此配准：（a)在确定非体素图像空间图像表示内对应于体 素图像空间的体素中心的点之前的体素图像空间第一扫描器图像表示，（b)体素图像空间， 和（c)非体素图像空间第一扫描器图像表不。
7. 根据权利要求6所述的诊断成像系统，还包括： 用于使体素图像空间第一扫描器图像表示与第二扫描器图像表示合并的合并装置。
8. 根据权利要求1所述的诊断成像系统，其中第一和第二扫描器包括至少以下中的 两个：PET、SPECT、MRI、超声、荧光透视、CT和数字X射线。
9. 一种诊断成像系统，包括： 用于获得对象感兴趣区域的第一组电子数据的第一扫描器； 用于获得对象感兴趣区域的第二组电子数据的第二扫描器； 用于将第一组电子数据重建成非体素图像空间内的非体素图像空间第一扫描器图像 表不的第一重建装置； 用于将第二组电子数据重建成第二图像空间内的第二图像表示的第二重建装置；以及 用于将第一组电子数据的非体素图像空间第一扫描器图像表示转换到第二图像空间 的转换装置，其中该非体素图像空间是非笛卡尔图像空间并且该第二图像空间是笛卡尔体 素图像空间。
10. -种诊断成像系统，包括： 用于获得对象感兴趣区域的第一组电子数据的第一扫描器； 用于获得对象感兴趣区域的第二组电子数据的第二扫描器； 用于将第一组电子数据重建成非体素图像空间内的非体素图像空间第一扫描器图像 表示的第一重建装置，该非体素图像空间包括多个团块，每个团块的中心在体心立方网格 的体心点上； 用于将第二组电子数据重建成第二图像空间内的第二图像表示的第二重建装置，该第 二图像空间是基于体素的笛卡尔网格图像空间；以及 用于将非体素图像空间第一扫描器图像表示直接转换到第二图像空间的转换装置。
11. 根据权利要求10所述的诊断成像系统，其中所述团块由Kaiser-Bessel函数b (r) 描述：
其中r是距团块中心的径向距离，Ιπ表示第一种级数m的修正贝塞尔函数，a是团块维 持的半径，α是控制团块形状的非负实数，以及m是控制团块在其边界r=a处的平滑度的 非负整数。
12. 根据权利要求11所述的诊断成像系统，其中所述处理装置包括： 用于选择最佳参数m、a和α以控制团块形状的参数选择装置。
13. -种诊断成像方法，包括： 从第一扫描器接收对象感兴趣区域的第一组电子数据； 从第二扫描器接收对象感兴趣区域的第二组电子数据； 将第一组电子数据重建成非体素图像空间内的非体素图像空间第一扫描器图像表 示； 将第二组电子数据重建成基于体素的笛卡尔图像空间内的第二图像表示； 确定基于体素的笛卡尔图像空间内的每个体素的中心点； 确定非体素图像空间第一扫描器图像表示内对应于每个所确定的体素中心点的点； 确定对应于非体素图像空间图像表示内的每个所确定的中心点的数据值，从而生成体 素图像空间第一扫描器图像表示；以及 在显示装置上显示下述的至少其一： 基于体素的笛卡尔图像空间第二图像表示； 体素图像空间第一扫描器图像表示；以及 基于体素的笛卡尔图像空间图像表示和体素图像空间第一扫描器图像表示的组合。
14. 根据权利要求13所述的方法，其中重建第一组电子数据包括将第一组电子数 据重建成团块，所述团块布置在体心立方网格上并且至少重叠邻近团块，并且该方法还包 括： 计算团块图像表示内重叠对应于体素中心的每个所确定的点的全部团块的总和；以及 重建合成体素图像f(x，y，z):
其中f(x，y，z)是合成3D图像表示，ΦωΛ是团块图像表示，φ?μ1是体素图像空间第 一扫描器图像表示，η是团块图像内的采样点数，m是体素图像内的采样点数，Λ2是团块尺 寸或采样间隔，Λ 1是体素图像空间内的体素尺寸或采样间隔，(^是在每个采样点η处的图 像系数，1是在每个采样点m处的图像系数，采样点m取决于覆盖体素中心的团块数。
15. 根据权利要求13所述的方法，还包括： 使体素图像空间第一扫描器图像表示与第二扫描器图像表示合并。
16. 根据权利要求13所述的方法，其中第一和第二扫描器包括至少以下中的两个： PET、SPECT、MRI、超声、荧光透视、CT和数字X射线。
17. -种用于执行权利要求13的方法的诊断成像系统。
18. -种多模态成像系统，包括： 第一扫描器， 第二扫描器， 编程为执行权利要求13的方法的数据处理器。
【文档编号】A61B8/00GK104107065SQ201410319237
【公开日】2014年10月22日 申请日期:2005年11月4日 优先权日:2004年11月19日
【发明者】D.加弄, W.王 申请人:皇家飞利浦电子股份有限公司