CN111415394B - 骨硬化伪影校正方法、装置、计算机设备和可读存储介质 - Google Patents

Abstract

本申请涉及一种骨硬化伪影校正方法、装置、计算机设备和可读存储介质，包括获取原始数据以及与原始数据对应的重建图像；根据重建图像获取水投影数据；根据原始数据以及水投影数据，计算得到误差投影数据；根据误差投影数据和原始数据，得到校正图像。上述方法通过在投影域上对原始数据进行校正，可以省掉一个反向投影和正向投影的过程，减小了计算量。另外，由于x射线可能会同时经过水和骨，则水的存在也会导致x射线能谱的改变，进而会影响骨硬化导致的误差，因此骨硬化误差是由射线穿过水和骨的厚度共同决定的。上述方法通过同时使用水投影数据和原始数据可以更为精确地计算出误差投影数据，提高了骨硬化伪影的校正精度。

Description

骨硬化伪影校正方法、装置、计算机设备和可读存储介质

技术领域

本发明涉及医学影像处理技术，特别是涉及一种骨硬化伪影校正方法、装置、计算机设备和可读存储介质。

背景技术

计算机断层成像(Computed Tomography，简称CT)是最常用的医学诊断方法之一。众所周知，由于球管产生的x射线并不是单色的，当x射线穿过物体时，会发生射束硬化这种现象，在CT成像时会引入硬化伪影。因此为了得到能够用于临床诊断的CT图像，在图像重建过程中需要进行射束硬化校正。对人体图像而言，主要是水硬化校正和骨硬化校正。尤其对于头部，临床上对图像的要求更高，在此种情况下，骨硬化校正必不可少。倘若不进行骨硬化校正，骨硬化伪影会导致在软组织和骨组织边界处的图像的CT值增高，产生模糊的边界，还会导致致密物体之间和骨的延长线方向的暗色条带，这在颅底更为明显。

目前常用的骨硬化校正方法是一种图像后处理方法，主要包括以下步骤：首先对重建好的图像进行分割，设置合适的阈值将图像划分为水图和骨图；其次对水图和骨图分别进行正向投影，得到投影数据；然后对投影数据进行多项式计算得到骨硬化伪影的投影数据；最后对上步得到的投影数据进行反向投影重建出骨硬化伪影，并将原始图像减去骨硬化伪影图像得到校正后的图像。

但是上述方法需要单独地对每张图像进行校正。对于探测器排数很多的CT，比如320排，由于准直宽度较大，对于边缘的排，X射线经过物体的路径有很大倾斜，会同时穿过多张图像，而不在单张图像的平面上，由于几何上的不匹配，骨硬化校正的效果不理想。此外，上述方法需要两个正投过程，需要分别对分割得到的水图和骨图进行正向投影，计算量较大。

发明内容

本申请提供一种骨硬化伪影校正方法、装置、计算机设备和可读存储介质，以至少解决骨硬化伪影校正精度低和计算量大的问题。

第一方面，本申请实施例提供了一种骨硬化伪影校正方法，包括：

获取原始数据以及与所述原始数据对应的重建图像；

根据所述重建图像获取水投影数据；

根据所述原始数据以及水投影数据，计算得到误差投影数据；

根据所述误差投影数据和所述原始数据，得到校正图像。

在其中一些实施例中，所述获取所述原始数据对应的重建图像包括：

对所述原始数据进行降采样处理；

对降采样处理后得到的数据进行重建，得到所述重建图像。

在其中一些实施例中，所述根据所述重建图像获取水投影数据包括：

根据所述重建图像获取骨头图像；

对所述骨头图像进行正投影得到骨头投影数据；

将所述原始数据减去所述骨头投影数据得到所述水投影数据。

在其中一些实施例中，所述对所述骨头图像进行正投影得到骨头投影数据包括：

对所述骨头图像进行3D平行束正向投影或3D锥束正向投影，得到所述骨头投影数据。

在其中一些实施例中，所述根据所述重建图像获取水投影数据包括：

根据所述重建图像获取水图图像；

对所述水图图像进行正投影得到所述水投影数据。

在其中一些实施例中，所述根据所述原始数据以及水投影数据，计算得到误差投影数据包括：

根据所述原始数据以及水投影数据建立误差模型；

查找所述水投影数据所对应的校正系数；

将所述校正系数以及所述原始数据代入所述误差模型，得到所述误差投影数据。

在其中一些实施例中，所述查找所述水投影数据所对应的校正系数包括：

根据所述水投影数据得到水的厚度值；

根据所述水的厚度值，得到所述校正系数。

在其中一些实施例中，所述根据所述误差投影数据和所述原始数据，得到校正图像包括：

将所述原始数据减去所述误差投影数据得到校正数据；

对所述校正数据重建得到所述校正图像。

在其中一些实施例中，所述根据所述误差投影数据和所述原始数据，得到校正图像包括：

对所述原始数据进行重建得到重建实际图像；

对所述误差投影数据进行重建得到重建伪影图像；

将所述重建实际图像减去所述重建伪影图像得到所述校正图像。

第二方面，本申请实施例提供了一种骨硬化伪影校正装置，所述装置包括：

第一获取模块，用于获取原始数据以及与所述原始数据对应的重建图像；

第二获取模块，用于根据所述重建图像获取水投影数据；

计算模块，用于根据所述原始数据以及水投影数据，计算得到误差投影数据；

校正模块，用于根据所述误差投影数据和所述原始数据得到校正图像。

第三方面，本申请实施例提供了一种计算机设备，包括存储器、处理器以及存储在所述存储器上并可在所述处理器上运行的计算机程序，所述处理器执行所述计算机程序时实现如上述第一方面所述的骨硬化伪影校正方法。

第四方面，本申请实施例提供了一种计算机可读存储介质，其上存储有计算机程序，该程序被处理器执行时实现如上述第一方面所述的骨硬化伪影校正方法。

相比于相关技术，本申请实施例提供的骨硬化伪影校正方法，包括获取原始数据以及与所述原始数据对应的重建图像；根据所述重建图像获取水投影数据；根据所述原始数据以及水投影数据，计算得到误差投影数据；根据所述误差投影数据和所述原始数据，得到校正图像，上述方法通过在投影域上对原始生数据进行校正，解决了骨硬化伪影校正精度低和计算量大的问题。

本申请的一个或多个实施例的细节在以下附图和描述中提出，以使本申请的其他特征、目的和优点更加简明易懂。

附图说明

此处所说明的附图用来提供对本申请的进一步理解，构成本申请的一部分，本申请的示意性实施例及其说明用于解释本申请，并不构成对本申请的不当限定。在附图中：

图1为一实施例提供的CT扫描系统的结构示意图；

图2为一实施例提供的骨硬化伪影校正方法的流程图；

图3为另一实施例提供的骨硬化伪影校正方法的流程图；

图4为再一实施例提供的骨硬化伪影校正方法的流程图；

图5为一个实施例中骨硬化伪影校正装置的结构框图；

图6为一个实施例中计算机设备的内部结构图。

具体实施方式

为了使本申请的目的、技术方案及优点更加清楚明白，以下结合附图及实施例，对本申请进行描述和说明。应当理解，此处所描述的具体实施例仅仅用以解释本申请，并不用于限定本申请。基于本申请提供的实施例，本领域普通技术人员在没有作出创造性劳动的前提下所获得的所有其他实施例，都属于本申请保护的范围。

显而易见地，下面描述中的附图仅仅是本申请的一些示例或实施例，对于本领域的普通技术人员而言，在不付出创造性劳动的前提下，还可以根据这些附图将本申请应用于其他类似情景。此外，还可以理解的是，虽然这种开发过程中所作出的努力可能是复杂并且冗长的，然而对于与本申请公开的内容相关的本领域的普通技术人员而言，在本申请揭露的技术内容的基础上进行的一些设计，制造或者生产等变更只是常规的技术手段，不应当理解为本申请公开的内容不充分。

在本申请中提及“实施例”意味着，结合实施例描述的特定特征、结构或特性可以包含在本申请的至少一个实施例中。在说明书中的各个位置出现该短语并不一定均是指相同的实施例，也不是与其它实施例互斥的独立的或备选的实施例。本领域普通技术人员显式地和隐式地理解的是，本申请所描述的实施例在不冲突的情况下，可以与其它实施例相结合。

除非另作定义，本申请所涉及的技术术语或者科学术语应当为本申请所属技术领域内具有一般技能的人士所理解的通常意义。本申请所涉及的“一”、“一个”、“一种”、“该”等类似词语并不表示数量限制，可表示单数或复数。本申请所涉及的术语“包括”、“包含”、“具有”以及它们任何变形，意图在于覆盖不排他的包含；例如包含了一系列步骤或模块(单元)的过程、方法、系统、产品或设备没有限定于已列出的步骤或单元，而是可以还包括没有列出的步骤或单元，或可以还包括对于这些过程、方法、产品或设备固有的其它步骤或单元。本申请所涉及的“连接”、“相连”、“耦接”等类似的词语并非限定于物理的或者机械的连接，而是可以包括电气的连接，不管是直接的还是间接的。本申请所涉及的“多个”是指两个或两个以上。“和/或”描述关联对象的关联关系，表示可以存在三种关系，例如，“A和/或B”可以表示：单独存在A，同时存在A和B，单独存在B这三种情况。字符“/”一般表示前后关联对象是一种“或”的关系。本申请所涉及的术语“第一”、“第二”、“第三”等仅仅是区别类似的对象，不代表针对对象的特定排序。

本申请中描述的各种技术可以应用于CT扫描系统中，如图1所示，CT扫描系统包括数据采集设备101、扫描床102、控制设备103和重建设备104。其中，控制设备103即通常所称的主控台，可以接收CT扫描参数，例如，控制设备103可以是一个装载有CT扫描控制软件的计算机。医师在对病人进行扫描时，可以在该控制软件的操作界面上设置本次扫描的CT扫描参数，包括扫描方式、螺旋螺距、图像厚度等。医师设置的CT扫描参数可以转换为扫描控制命令，用以控制数据采集设备101对病人进行扫描。数据采集设备101可以是一个内部包含有球管、探测器等装置的机架，即通常所称的机架设备，用于数据采集和传输，扫描床102可以移动至该机架内以对病人进行CT扫描。并且，对病人进行的扫描是依据控制设备103获取的CT扫描参数，按照该参数指定的扫描方式、螺旋螺距等进行扫描。数据采集设备101获取到的采样数据，比如对病人扫描时得到的X射线穿过病人的衰减信息，可以传输至重建设备104，重建设备104即通常所称的建像机，由重建设备104将该采样数据存储在硬盘上并进行CT扫描的图像重建，以得到病人的扫描图像。此外，控制设备103可以将接收到的CT扫描参数的全部或部分传输至重建设备104，至少包括重建设备104在图像重建时需要使用到的一些参数，例如，图像间隔、图像厚度等建像参数，以使得重建设备104根据该CT扫描参数和由数据采集设备101得到的采样数据进行建像。

图2为一实施例提供的骨硬化伪影校正方法的流程图，如图2所示，骨硬化伪影校正方法包括步骤210至步骤240，其中：

步骤210，获取原始数据以及与原始数据对应的重建图像。

根据不同的扫描要求，可以对被扫描物体的不同部位设置不同的CT扫描参数。根据CT扫描参数，在对被扫描物体的不同部位扫描时，使用对应该部位的CT扫描参数进行扫描，获取原始数据。例如，控制设备在接收到设置的CT扫描参数后，可以根据该参数控制数据采集设备进行数据采集，得到原始数据。

数据采集设备获取到的原始数据，比如对被扫描物体扫描时得到的X射线穿过被扫描物体的衰减信息，可以传输至重建设备，重建设备可以将该原始数据存储在硬盘上并进行CT扫描的图像重建，以得到被扫描物体的重建图像。

在其中一些实施例中，在得到原始数据后，对原始数据进行预处理。预处理过程具体包括对原始数据进行校正和重排。因为实际的CT扫描系统并不是理想的，会造成采集到的信号带有缺陷或几何上有偏差，因此在重建之前需要校正成所需要的理想信号。而数据重排与图像重建使用的反投方法相关，一般重建使用的反投方法是基于平行束射线的，需要将原始数据重排成等效的平行束，如果是基于锥束，则可以不用重排。

步骤220，根据重建图像获取水投影数据。

人软组织的衰减性质与水十分相近，近似认为人体由水和骨组成。因此对人体扫描得到的重建图像包括水图和骨图，通过选择合适的阈值对重建图像进行划分，得到水投影数据。

在其中一些实施例中，根据重建图像获取水投影数据包括：

根据重建图像获取骨头图像；

对骨头图像进行正投影得到骨头投影数据；

将原始数据减去骨头投影数据得到水投影数据。

具体地，将重建好的图像通过阈值进行划分提取出骨图。然后将骨图进行正投，得到骨的投影数据，将原始数据减去骨的投影数据得到水投影数据。本实施例中，对骨头图像进行3D平行束正向投影或3D锥束正向投影，得到骨头投影数据。使得投影更符合真实的x射线路径，对于准直宽度更大的多排CT，校正效果更为理想。

在其中一些实施例中，根据重建图像获取水投影数据包括：

根据重建图像获取水图图像；

对水图图像进行正投影得到水投影数据。

需要说明的是，可以直接利用水图图像的投影数据进行骨硬化校正，但用于提取水图图像的FOV应比较大，要包含x射线从球管到探测器经过的路径上穿过的所有物体，比如头托或床板。

步骤230，根据原始数据以及水投影数据，计算得到误差投影数据。

可以利用水投影数据和原始数据进多项式计算得到骨的误差投影数据。

在其中一些实施例中，根据原始数据以及水投影数据，计算得到误差投影数据包括：

根据原始数据以及水投影数据建立误差模型；

查找水投影数据所对应的校正系数；

将校正系数以及原始数据代入误差模型，得到误差投影数据。

在其中一些实施例中，查找水投影数据所对应的校正系数包括：

根据水投影数据得到水的厚度值；

根据水的厚度值，得到校正系数。

骨硬化伪影校正的目的是将多色x射线穿过物体的效果校正成x射线是单能时的效果。对已知能谱的x射线，骨硬化导致的误差Err是实际投影值P_real与理想投影值P_ideal的差，这由x射线穿过水和骨的厚度决定，因为P_real和P_bone包含了水和骨厚度的信息，因此误差模型Err为P_real和P_bone的函数，即：

Err＝P_real-P_ideal＝Err(P_real,P_bone)

其中，P_ideal是假设单能的x射线穿过物体时的理想投影值，P_real是球管实际发出的x射线穿过物体时的投影值，P_bone是骨图通过正投得到的投影值。对于水的厚度一定，则误差模型仅是P_real的函数，可以将之写为以下多项式的形式，多项式的阶数不限，可以为3至5阶，误差模型可以写为下式：

由该式可以看出，误差模型式与原始数据和水投影数据有关。其中，C(m；L_water)为校正系数。校正系数可以是通过理论计算或用标准模体定标得到，对应不同的探测器，校正系数的数值也是不同的。在做骨硬化校正时，将P_real减去P_bone得到P_water，再将P_water除以一个因子λ近似得到水的厚度，已知水的厚度挑选出相应的校正系数C(:,L_water)，利用上式计算误差投影数据Err。或者也可将误差模型Err写为下式形式，

上式C的数值与水的厚度无关，水的信息隐含在P和P_bone中。以上两种方式的效果是一样的。

步骤240，根据误差投影数据和原始数据，得到校正图像。

本实施提供的骨硬化伪影校正方法，通过获取原始数据以及与原始数据对应的重建图像；根据重建图像获取水投影数据；根据原始数据以及水投影数据，计算得到误差投影数据；根据误差投影数据和原始数据，得到校正图像。上述方法通过在投影域上对原始数据进行校正，可以省掉一个反向投影和正向投影的过程，减小了计算量。另外，由于x射线可能会同时经过水和骨，则水的存在也会导致x射线能谱的改变，进而会影响骨硬化导致的误差，因此骨硬化误差是由射线穿过水和骨的厚度共同决定的。如果x射线同时经过了水和骨，则水的存在会也导致x射线能谱的改变，也会影响骨硬化导致的误差，也即上面提到的骨硬化误差Err是由射线穿过水和骨的厚度共同决定的，单独骨的投影值并不能完全决定Err，只能作为一种近似。因此，本申请同时使用骨投影数据和原始数据(或水的投影数据)可以更为精确地计算出误差投影数据Err，校正效果更好。

在其中一些实施例中，获取原始数据对应的重建图像包括：

对原始数据进行降采样处理；

对降采样处理后得到的数据进行重建，得到重建图像。

本实施例中，首先对原始数据进行降采样处理，可以同时降低采样的角度数和通道数。将降采样处理后的数据进行图像重建，图像大小应包括完整的头部，重建的图像可以是低分辨图像。例如若采集到的原始数据进行降采样处理，降采样因子为2，则可以将重建图像的橡素数量减小为原来的1/2。例如原始重建图像矩阵大小为[512，512]，则降采样处理后重建的图像矩阵大小为[256，256]。由于骨硬化伪影主要校正的是低频部分，因此，首先对原始数据进行降采样处理，可以在不影响校正质量的同时，可以减小计算量，从而提高校正效率。

需要说明的是，根据降采样处理后的原始数据得到的重建图像计算得到误差投影数据后，误差投影数据需要通过插值进行过采样处理，以使误差投影数据的个数与原始生数据相同，然后将原始数据减去误差投影数据，从而得到校正后的数据。

在其中一些实施例中，如图3所示，根据误差投影数据和原始数据，得到校正图像包括：

将原始数据减去误差投影数据得到校正数据；

对校正数据重建得到校正图像。

本实施例提供的骨硬化伪影方法的具体实现如下：获取原始数据并对原始数据进行预处理(重排)，对预处理后的原始数据进行降采样处理，然后对降采样处理后的原始数据进行重建得到低分辨率的重建图像。通过阈值划分，提取出骨图并对骨图进行正投影得到骨投影数据，利用骨投影数据和原始数据进行多项式计算得出误差投影数据。对误差投影数据进行过采样处理后，将原始数据减去过采样处理后的原始数据得到校正数据，对校正数据进行重建得到校正图像。

本实施例中，骨硬化伪影校正方法分为两个分支，一个可以称为主分支，另一个称为骨硬化分支。其中，骨硬化分支用于获取误差投影数据，主分支用于对原始数据和误差投影数据进行处理得到校正图像。具体地，主分支在得到骨硬化分支获取到的误差投影数据后，将原始数据减去误差投影数据得到校正数据，再对校正数据进行重建得到校正图像。

本实施例通过在投影域上直接对原始数据进行校正，校正效果更好；其次可以省略重建骨硬化伪影图像这个反投的过程，有利于算法性能的提升。

在其中一些实施例中，如图4所示，根据误差投影数据和原始数据，得到校正图像包括：

对原始数据进行重建得到重建实际图像；

对误差投影数据进行重建得到重建伪影图像；

将重建实际图像减去重建伪影图像得到校正图像。

本实施例提供的骨硬化伪影方法的具体实现如下：获取原始数据并对原始数据进行预处理(重排)，对预处理后的原始数据进行降采样处理，然后对降采样处理后的原始数据进行重建得到低分辨率的重建图像。通过阈值划分，提取出骨图并对骨图进行正投影得到骨投影数据，对骨投影数据进行重建得到伪影图，同时对原始数据进行重建得到原始数据的重建图像，将原始数据的重建图像减去伪影图，得到校正图像。

本实施例中，将骨硬化分支和主分支并行，如果骨硬化分支进行得足够快，则主分支不会因为要等待骨硬化校正过程的完成而停留，相当于骨硬化分支不会给整个流程增加运行时间，提高了校正效率。此外，此方法也可适用于生数据没有经过重排的情况，相应地，对骨图进行投影时应改用3D锥束正投。

应该理解的是，虽然图2至图4的流程图中的各个步骤按照箭头的指示依次显示，但是这些步骤并不是必然按照箭头指示的顺序依次执行。除非本文中有明确的说明，这些步骤的执行并没有严格的顺序限制，这些步骤可以以其它的顺序执行。而且，图2至图4中的至少一部分步骤可以包括多个子步骤或者多个阶段，这些子步骤或者阶段并不必然是在同一时刻执行完成，而是可以在不同的时刻执行，这些子步骤或者阶段的执行顺序也不必然是依次进行，而是可以与其它步骤或者其它步骤的子步骤或者阶段的至少一部分轮流或者交替地执行。

在一个实施例中，如图5所示，提供了一种骨硬化伪影校正装置，包括：第一获取模块510、第二获取模块520、计算模块530和校正模块540，其中：

第一获取模块510，用于获取原始数据以及与原始数据对应的重建图像；

第二获取模块520，用于根据重建图像获取水投影数据；

计算模块530，用于根据原始数据以及水投影数据，计算得到误差投影数据；

校正模块540，用于根据误差投影数据和原始数据得到校正图像。

本申请涉及一种骨硬化伪影校正装置，包括第一获取模块510、第二获取模块520、计算模块530和校正模块540，通过第一获取模块510

获取原始数据以及与原始数据对应的重建图像；第二获取模块520根据重建图像获取水投影数据；计算模块530根据原始数据以及水投影数据，计算得到误差投影数据；校正模块540根据误差投影数据和原始数据，得到校正图像。上述装置通过在投影域上对原始数据进行校正，可以省掉一个反向投影和正向投影的过程，减小了计算量。另外，由于x射线可能会同时经过水和骨，则水的存在也会导致x射线能谱的改变，进而会影响骨硬化导致的误差，因此骨硬化误差是由射线穿过水和骨的厚度共同决定的。上述方装置通过同时使用水投影数据和原始数据可以更为精确地计算出误差投影数据，提高了骨硬化伪影的校正精度。

在其中一些实施例中，第一获取模块510还用于对原始数据进行降采样处理；

对降采样处理后得到的数据进行重建，得到重建图像。

在其中一些实施例中，第二获取模块520还用于根据重建图像获取骨头图像；

对骨头图像进行正投影得到骨头投影数据；

将原始数据减去骨头投影数据得到水投影数据。

在其中一些实施例中，第二获取模块520还用于对骨头图像进行3D平行束正向投影或3D锥束正向投影，得到骨头投影数据。

在其中一些实施例中，第二获取模块520还用于根据重建图像获取水图图像；

对水图图像进行正投影得到水投影数据。

在其中一些实施例中，计算模块530还用于根据原始数据以及水投影数据建立误差模型；

查找水投影数据所对应的校正系数；

将校正系数以及原始数据代入误差模型，得到误差投影数据。

在其中一些实施例中，计算模块530还用于根据水投影数据得到水的厚度值；

根据水的厚度值，得到校正系数。

在其中一些实施例中，校正模块540还用于将原始数据减去误差投影数据得到校正数据；

对校正数据重建得到校正图像。

在其中一些实施例中，校正模块540还用于对原始数据进行重建得到重建实际图像；

对误差投影数据进行重建得到重建伪影图像；

将重建实际图像减去重建伪影图像得到校正图像。

关于骨硬化伪影校正装置的具体限定可以参见上文中对于骨硬化伪影校正方法的限定，在此不再赘述。上述骨硬化伪影校正装置中的各个模块可全部或部分通过软件、硬件及其组合来实现。上述各模块可以硬件形式内嵌于或独立于计算机设备中的处理器中，也可以以软件形式存储于计算机设备中的存储器中，以便于处理器调用执行以上各个模块对应的操作。

在一个实施例中，提供了一种计算机设备，该计算机设备可以是终端，其内部结构图可以如图6所示。该计算机设备包括通过系统总线连接的处理器、存储器、网络接口、显示屏和输入装置。其中，该计算机设备的处理器用于提供计算和控制能力。该计算机设备的存储器包括非易失性存储介质、内存储器。该非易失性存储介质存储有操作系统和计算机程序。该内存储器为非易失性存储介质中的操作系统和计算机程序的运行提供环境。该计算机设备的网络接口用于与外部的终端通过网络连接通信。该计算机程序被处理器执行时以实现一种骨硬化伪影校正方法。该计算机设备的显示屏可以是液晶显示屏或者电子墨水显示屏，该计算机设备的输入装置可以是显示屏上覆盖的触摸层，也可以是计算机设备外壳上设置的按键、轨迹球或触控板，还可以是外接的键盘、触控板或鼠标等。

本领域技术人员可以理解，图6中示出的结构，仅仅是与本申请方案相关的部分结构的框图，并不构成对本申请方案所应用于其上的计算机设备的限定，具体的计算机设备可以包括比图中所示更多或更少的部件，或者组合某些部件，或者具有不同的部件布置。

在一个实施例中，提供了一种计算机设备，包括存储器和处理器，存储器中存储有计算机程序，该处理器执行计算机程序时实现以下步骤：

获取原始数据以及与所述原始数据对应的重建图像；

根据所述重建图像获取水投影数据；

根据所述原始数据以及水投影数据，计算得到误差投影数据；

根据所述误差投影数据和所述原始数据，得到校正图像。

在一个实施例中，提供了一种计算机可读存储介质，其上存储有计算机程序，计算机程序被处理器执行时实现以下步骤：

获取原始数据以及与所述原始数据对应的重建图像；

根据所述重建图像获取水投影数据；

根据所述原始数据以及水投影数据，计算得到误差投影数据；

根据所述误差投影数据和所述原始数据，得到校正图像。

本领域普通技术人员可以理解实现上述实施例方法中的全部或部分流程，是可以通过计算机程序来指令相关的硬件来完成，所述的计算机程序可存储于一非易失性计算机可读取存储介质中，该计算机程序在执行时，可包括如上述各方法的实施例的流程。其中，本申请所提供的各实施例中所使用的对存储器、存储、数据库或其它介质的任何引用，均可包括非易失性和/或易失性存储器。非易失性存储器可包括只读存储器(ROM)、可编程ROM(PROM)、电可编程ROM(EPROM)、电可擦除可编程ROM(EEPROM)或闪存。易失性存储器可包括随机存取存储器(RAM)或者外部高速缓冲存储器。作为说明而非局限，RAM以多种形式可得，诸如静态RAM(SRAM)、动态RAM(DRAM)、同步DRAM(SDRAM)、双数据率SDRAM(DDRSDRAM)、增强型SDRAM(ESDRAM)、同步链路(Synchlink)DRAM(SLDRAM)、存储器总线(Rambus)直接RAM(RDRAM)、直接存储器总线动态RAM(DRDRAM)、以及存储器总线动态RAM(RDRAM)等。

以上所述实施例的各技术特征可以进行任意的组合，为使描述简洁，未对上述实施例中的各个技术特征所有可能的组合都进行描述，然而，只要这些技术特征的组合不存在矛盾，都应当认为是本说明书记载的范围。

以上所述实施例仅表达了本申请的几种实施方式，其描述较为具体和详细，但并不能因此而理解为对申请专利范围的限制。应当指出的是，对于本领域的普通技术人员来说，在不脱离本申请构思的前提下，还可以做出若干变形和改进，这些都属于本申请的保护范围。因此，本申请专利的保护范围应以所附权利要求为准。

Claims (11)

1.一种骨硬化伪影校正方法，其特征在于，所述方法包括：

获取原始数据以及与所述原始数据对应的重建图像；

根据所述重建图像获取水投影数据；

根据所述原始数据以及水投影数据，计算得到误差投影数据；

根据所述误差投影数据和所述原始数据，得到校正图像；

所述根据所述原始数据以及水投影数据，计算得到误差投影数据包括：根据所述原始数据以及水投影数据建立误差模型；查找所述水投影数据所对应的校正系数；将所述校正系数以及所述原始数据代入所述误差模型，得到所述误差投影数据。

2.根据权利要求1所述的方法，其特征在于，所述获取所述原始数据对应的重建图像包括：

对所述原始数据进行降采样处理；

对降采样处理后得到的数据进行重建，得到所述重建图像。

3.根据权利要求1所述的方法，其特征在于，所述根据所述重建图像获取水投影数据包括：

根据所述重建图像获取骨头图像；

对所述骨头图像进行正投影得到骨头投影数据；

将所述原始数据减去所述骨头投影数据得到所述水投影数据。

4.根据权利要求3所述的方法，其特征在于，所述对所述骨头图像进行正投影得到骨头投影数据包括：

对所述骨头图像进行3D平行束正向投影或3D锥束正向投影，得到所述骨头投影数据。

5.根据权利要求1所述的方法，其特征在于，所述根据所述重建图像获取水投影数据包括：

根据所述重建图像获取水图图像；

对所述水图图像进行正投影得到所述水投影数据。

6.根据权利要求1所述的方法，其特征在于，所述查找所述水投影数据所对应的校正系数包括：

根据所述水投影数据得到水的厚度值；

根据所述水的厚度值，得到所述校正系数。

7.根据权利要求1所述的方法，其特征在于，所述根据所述误差投影数据和所述原始数据，得到校正图像包括：

将所述原始数据减去所述误差投影数据得到校正数据；

对所述校正数据重建得到所述校正图像。

8.根据权利要求1所述的方法，其特征在于，所述根据所述误差投影数据和所述原始数据，得到校正图像包括：

对所述原始数据进行重建得到重建实际图像；

对所述误差投影数据进行重建得到重建伪影图像；

将所述重建实际图像减去所述重建伪影图像得到所述校正图像。

9.一种骨硬化伪影校正装置，其特征在于，所述装置包括：

第一获取模块，用于获取原始数据以及与所述原始数据对应的重建图像；

第二获取模块，用于根据所述重建图像获取水投影数据；

计算模块，用于根据所述原始数据以及水投影数据，计算得到误差投影数据；

校正模块，用于根据所述误差投影数据和所述原始数据得到校正图像；

所述计算模块，还用于根据所述原始数据以及水投影数据建立误差模型；查找所述水投影数据所对应的校正系数；将所述校正系数以及所述原始数据代入所述误差模型，得到所述误差投影数据。

10.一种计算机设备，包括存储器和处理器，所述存储器存储有计算机程序，其特征在于，所述处理器执行所述计算机程序时实现权利要求1至8中任一项所述方法的步骤。

11.一种计算机可读存储介质，其上存储有计算机程序，其特征在于，所述计算机程序被处理器执行时实现权利要求1至8中任一项所述的方法的步骤。