发明内容
1.要解决的技术问题
基于现有的图像重建想要重建出达到临床诊断标准的动态图像序列,就必须对病人进行重复多次相同剂量CT扫描,这会给病人带来很大的X-射线辐射剂量,并存在潜在的引发被扫描部位辐射病变的危险的问题,本申请提供了一种图像重建方法及其应用。
2.技术方案
为了达到上述的目的,本申请提供了一种图像重建方法,所述方法包括如下步骤:
1)采集动态灌注图像序列和基线图像序列;
2)对基线图像投影数据序列进行降采样,在每一帧动态灌注图像序列中减去降采样的基线图像投影数据序列,得到增强投影数据时间序列;
3)对所述增强投影数据时间序列进行迭代预重建,得到预重建对比增强图像序列;
4)对每一个图像像素点均提取一个与之对应的像素特征向量;
5)基于所述像素特征向量构建时间序列信息先验矩阵;
6)将所述时间序列信息先验矩阵以“图像遮罩”的形式引入对比增强图像序列的重建迭代中,得到重建的对比增强图像;
7)将所述对比增强图像投影数据序列逐帧与所述基线图像投影数据相加,得到重建的动态灌注图像时间序列。
本申请提供的另一种实施方式为:所述步骤1)中基线图像序列为显像剂到达人体组织之前进行的扫描图像,所述动态灌注图像序列为显像剂到达人体组织之后进行的连续扫描图像。
本申请提供的另一种实施方式为:所述步骤1)中所述动态灌注图像采用稀疏的降采样扫描方案进行扫描得到;所述基线图像采用非稀疏的全采样扫描方案进行扫描得到。
本申请提供的另一种实施方式为:所述步骤3)中采用全变分正则化的最大似然期望最大化算法进行预重建。
本申请提供的另一种实施方式为:所述步骤4)中所述像素特征向量与显像剂时间演变相关。
本申请提供的另一种实施方式为:所述步骤4)中所述像素特征向量的提取方法为:将所述预重建对比增强图像序列中的每一帧图像的对应像素点像素值用于组成图像像素点的像素特征向量。
本申请提供的另一种实施方式为:所述步骤5)中构建时间序列信息先验矩阵包括:对于每一个像素特征向量在像素特征集合/>中确定k个与之建立联系的相近像素特征向量;对处在/>的k-最邻近簇中的k个像素特征向量与/>通过一个预先设计的高斯连接函数建立联系,不处在/>的k-最邻近簇中的其他像素特征向量与/>的连接函数值定义为0;将基于像素特征向量/>计算的所有连接函数值ψ((i,j),(p,q))重排为一个长度为512×512的列向量。
本申请提供的另一种实施方式为:相近像素特征向量的确定方法为:计算每一个像素特征向量与被考虑的像素特征向量/>之间的欧几里得距离,与/>的欧几里得距离最近的k个像素特征向量即构成/>的k-最邻近簇。
本申请提供的另一种实施方式为:所述步骤7)中基线图像为从全采样基线投影数据经过滤波反投影算法重建的正常剂量基线图像。
本申请还提供一种图像重建方法的应用,将所述的图像重建方法应用于动态心肌灌注扫描CT图像、脑部动态灌注CT成像与磁共振灌注成像。
3.有益效果
与现有技术相比,本申请提供的一种图像重建方法的有益效果在于:
本申请提供的图像重建方法,针对医学图像。
本申请提供的图像重建方法,为一种时间序列信息协助的动态灌注CT图像重建方法。
本申请提供的图像重建方法,为一种基于灌注增强图像信息与基线图像信息分离的时间序列协助下的稀疏低剂量心肌灌注动态CT图像(DMP-CT)重建方法。。
本申请提供的图像重建方法,用于精确重建通过稀疏低剂量扫描方案获取的动态灌注CT图像,在降低X射线扫描辐射剂量的条件下,重建图像达到临床冠状动脉疾病(CAD)诊断标准。
本申请提供的图像重建方法,解决从稀疏投影角度数据重建心肌灌注动态CT图像的技术难题,从重建算法层面支撑稀疏投影角度的心肌灌注动态CT扫描方案的临床应用。
本申请提供的图像重建方法,可以在减小心肌动态灌注CT扫描过程中的X射线扫描幅数的情况下获得达到临床诊断标准的重建图像,从而显著降低心肌动态灌注CT成像对病人的辐射剂量。
本申请提供的图像重建方法,在算法性能方面,本申请算法能够有效避免由于稀疏采样导致的重建图像伪影与错误信息问题,在稀疏角度CT的场景下,成像质量明显高于现有商用算法重建结果。
具体实施方式
在下文中,将参考附图对本申请的具体实施例进行详细地描述,依照这些详细的描述,所属领域技术人员能够清楚地理解本申请,并能够实施本申请。在不违背本申请原理的情况下,各个不同的实施例中的特征可以进行组合以获得新的实施方式,或者替代某些实施例中的某些特征,获得其它优选的实施方式。
使用稀疏角度对病人进行X射线是有效的降低CT扫描辐射剂量的方法,特别对于DMP-CT这种需要对病人进行多次扫描的医学影像检查。使用现有重建算法(滤波反投影、传统代数迭代重建)从稀疏采集的数据重建CT图像会由于数据的不完整性(incompleteness)和不一致性(inconsistency),导致重建图像中存在大量的伪影,严重影响诊断。
不失一般性地,CT成像过程在数学上可以用一个简单的矩阵方程描述:
y=Ax, (式1)
其中y为CT扫描得到的投影数据(数据校准与对数变换之后),A为与扫描系统对应的系统矩阵,表征从被扫描物体到投影数据的正投影过程,x为被扫描物(待重建图像)。CT成像的图像重建任务即是根据所获得的投影数据y,基于与成像系统对应的已知系统矩阵A,估算被扫描物体的X射线衰减系数值x。
对于心肌灌注动态CT成像(DMP-CT)而言,若对于每一帧图像的投影几何(投影角度,投影幅数,投影角度间隔,管电压与管电流等)是完全相同的,那么每一帧图像的成像过程均可用上述(式1)表示,则DMP-CT成像可用下式描述:
Y=AX, (式2)
其中,Y={y1,y2,…,yT}为投影数据时间序列,X={x1,x2,…,xT}为待重建图像序列。则DMP-CT成像任务可理解为基于成像系统矩阵A,根据投影数据时间序列Y,估计动态图像时间序列X。
参见图1~2,本申请提供一种图像重建方法,所述方法包括如下步骤:
1)采集动态灌注图像序列和基线图像序列;
2)对基线图像投影数据序列进行降采样,在每一帧动态灌注图像序列中减去降采样的基线图像投影数据序列,得到增强投影数据时间序列;
3)对所述增强投影数据时间序列进行迭代预重建,得到预重建对比增强图像序列;
4)对每一个图像像素点均提取一个与之对应的像素特征向量;
5)基于所述像素特征向量构建时间序列信息先验矩阵;
6)将所述时间序列信息先验矩阵以“图像遮罩”的形式引入对比增强图像序列的重建迭代中,得到重建的对比增强图像;
7)将所述对比增强图像投影数据序列逐帧与所述基线图像投影数据相加,得到重建的动态灌注图像时间序列。
进一步地,所述步骤1)中基线图像序列为显像剂到达人体组织之前进行的扫描图像,所述动态灌注图像序列为显像剂到达人体组织之后进行的连续扫描图像。
进一步地,所述步骤1)中所述动态灌注图像采用稀疏的降采样扫描方案进行扫描得到;所述基线图像采用非稀疏的全采样扫描方案进行扫描得到。
进一步地,所述步骤3)中采用全变分正则化的最大似然期望最大化算法进行预重建。
进一步地,所述步骤4)中所述像素特征向量与显像剂时间演变相关。
进一步地,所述步骤4)中所述像素特征向量的提取方法为:将所述预重建对比增强图像序列中的每一帧图像的对应像素点像素值用于组成图像像素点的像素特征向量。
进一步地,所述步骤5)中构建时间序列信息先验矩阵包括:对于每一个像素特征向量在像素特征集合/>中确定k个与之建立联系的相近像素特征向量;对处在/>k-最邻近簇中的k个像素特征向量与/>通过一个预先设计的高斯连接函数建立联系,不处在/>的k-最邻近簇中的其他像素特征向量与的连接函数值定义为0;将基于像素特征向量/>计算的所有连接函数值ψ((i,j),(p,q))重排为一个长度为512×512的列向量。
进一步地,所述相近像素特征向量的确定方法为:计算每一个像素特征向量与被考虑的像素特征向量/>之间的欧几里得距离,与/>的欧几里得距离最近的k个像素特征向量即构成/>的k-最邻近簇。
进一步地,所述步骤7)中基线图像为从全采样基线投影数据经过滤波反投影算法重建的正常剂量基线图像。
本申请还提供一种图像重建方法的应用,将所述的图像重建方法应用于动态心肌灌注扫描CT图像、脑部动态灌注CT成像与磁共振灌注成像。当然,本申请中只是列举,不限于其他类型图像的重建。
实施例
基于动态灌注增强CT图像的分离重建方法和动态图像序列的稀疏角度图像采集方案,对由显像剂引起的心脏组织对比增强图像借助动态图像序列中自然包含的显像剂时间演变信息进行重建,并将重建增强图像叠加到正常剂量采集方案得到的基线CT图像,获得最终的动态心肌灌注扫描CT图像序列。
图像数据采集:对心肌动态灌注CT成像的基线(baseline)图像采集,即显像剂到达人体心脏组织之前进行的CT扫描,和动态灌注图像序列采集,即显像剂到达人体心脏组织之后进行的连续CT扫描,采取不同剂量的CT图像采集方案。基线图像使用正常剂量,即非稀疏的全采样扫描方案进行扫描,动态灌注图像序列使用低剂量,即稀疏的降采样扫描方案进行扫描。该步骤可用下述式子描述:
基线投影数据分离:从所采集的稀疏的动态投影数据序列分理出对应于非显像剂引入的图像信息。具体实现方法为:将所采集基线投影数据y基线,全采样进行降采样,保证所得稀疏基线投影数据与动态灌注投影数据Y灌注,稀疏的每一帧都对应与相同的图像采集几何条件(即系统矩阵A)。在每一帧的动态灌注投影数据中减去降采样的稀疏基线投影数据,得到增强投影数据时间序列。该步骤可用下述式子描述:
对比增强投影数据的迭代预重建:对上述式(4)所得对比增强投影数据Y对比增强,稀疏使用全变分正则化的最大似然期望最大化算法(MLEM-TV)进行预重建,得到预重建的对比增强图像序列。该步骤可用下述式子描述:
显像剂时间演变特征提取:对于每一个CT图像像素点(i,j),均提取一个与之对应的与显像剂时间演变相关的像素特征向量像素特征向量的提取方法为:将上述步骤3)所得预重建对比增强图像序列中的每一帧图像的对应像素点像素值x1~T(i,j)用于组成像素点(i,j)的像素特征向量,具体可用下述式子描述:
其中,T表示动态对比增强图像的帧数。
时间序列信息先验矩阵的构建:基于所提取的关于每一个像素点像素特征向量构建时间序列信息先验矩阵。具体的构建方法为:
a.对于每一个像素特征向量在像素特征集合/>中确定k个与之建立联系的相近像素特征向量。相近像素特征向量的确定方法为:建立关于像素特征向量/>的k-最邻近簇,即计算每一个像素特征向量/>与被考虑的像素特征向量/>之间的欧几里得距离,与/>的欧几里得距离最近的k个像素特征向量即构成/>的k-最邻近簇。
b.处在的k-最邻近簇中的k个像素特征向量与/>通过一个预先设计的高斯连接函数建立联系,不处在/>的k-最邻近簇中的其他像素特征向量与/>的连接函数值定义为0。步骤b可用下述式描述:
c.将基于像素特征向量计算的所有连接函数值ψ((i,j),(p,q))重排为一个长度为512×512的列向量。因此,对于每一个像素特征向量/>都可以计算得一个由该像素特征向量与其他所有像素特征向量的连接函数值组成的列向量,所有的列向量即组成时间序列信息先验矩阵Ψ。该步骤可用下述式描述:
时间序列信息先验矩阵协助的增强图像重建:将上述步骤5)构建的时间序列信息先验矩阵以“图像遮罩”的形式引入对比增强图像序列的重建迭代中。将对比增强CT图像表示为一个图像遮罩与一个系数图像ζ的乘积,该图像遮罩即为前述步骤5)所构建的时间序列先验矩阵Ψ,即
则对应于时间帧t的对比增强图像成像任务可描述为:
将上述(式10)应用于极大似然期望值最大化迭代重建,则有关于系数图像ζ的迭代更新规则:
根据上式迭代计算所得系数图像ζ可由上述(式9)恢复为重建的对比增强图像
基线图像与对比增强图像的融合:将根据上述步骤重建得到的对比增强图像序列X对比增强,重建逐帧与从全采样基线投影数据经过滤波反投影算法(FBP)重建的正常剂量基线图像相加,即得到重建的动态灌注心肌CT图像时间序列X动态灌注,重建。该步骤可用下述时描述:
临床使用的心脏动态灌注CT检查:在显像剂到达人体组织达到对比增强效果之前对病人进行正常剂量的扫描,在显像剂到达人体组织形成对比增强后对病人实施低剂量的稀疏CT动态扫描。以此图像采集方案得到的原始图像数据即可使用本申请方法进行重建。
本申请方法经过猪心脏医学临床图像的仿真实验验证,效果明显:图2展示了使用本申请方法重建的心肌动态灌注CT图像重建结果。该实验在稀疏倍率为12的条件下进行(正常剂量全采样为984幅投影/360°,稀疏低剂量采样为82幅投影/360°)。由图可见本申请方法在稀疏角度采样心肌动态灌注CT图像重建中的优越性能。
本申请对于心肌动态灌注CT扫描的高低剂量配合扫描方案,即显像剂未到达之前使用正常剂量扫描方案与显像剂到达后的稀疏低剂量扫描方案切换。
本申请从预重建对比增强图像序列中提取像素特征向量,并基于像素特征向量构建时间序列先验矩阵的方法;
本申请将CT图像重建任务转化为图像遮蔽下的系数图像重建问题的方法,及最大似然期望最大化算法在系数图像重建中的应用;
本申请中灌注CT图像的分离重建思想,即将灌注图像的基线部分与对比增强部分分开重建,重建后再融合相加得到最终动态灌注CT图像。
尽管在上文中参考特定的实施例对本申请进行了描述,但是所属领域技术人员应当理解,在本申请公开的原理和范围内,可以针对本申请公开的配置和细节做出许多修改。本申请的保护范围由所附的权利要求来确定,并且权利要求意在涵盖权利要求中技术特征的等同物文字意义或范围所包含的全部修改。