光学投影断层成像运动伪影校正方法技术

本发明专利技术公开了光学投影断层成像运动伪影校正方法，使用不基于特征点的估计样品运动轨迹的方法，其实现步骤如下，(1)获取投影数据；(2)计算投影数据的零阶矩；(3)计算投影数据的一阶矩；(4)计算投影数据的质心；(5)获取运动参数；(6)计算运动量；(7)运动伪影校正。本发明专利技术中利用数据一致性条件建立样品与投影数据之间的对应关系，采用多项式估计被扫描样品的运动，进而直接从投影数据中估计出样品的运动信息。本发明专利技术可应用于样品的投影断层重建，可以提高光学断层成像系统的空间分辨率，减少图像伪影。

【技术实现步骤摘要】

本专利技术属于医学影像处理领域，更进一步涉及一种对光学投影断层成像系统的运动伪影进行校正的方法。该方法可用于光学投影断层成像的图像处理。
技术介绍
光学投影断层成像(以下简称OPT)是一种新型的光学分子影像成像技术，其成像原理和X射线计算机断层成像的原理类似。OPT可获得扫描样品结构像，又能利用荧光染料或荧光蛋白进行分子特异性标记，实现分子特性成像，并且设备成本低，使用方便。在OPT 扫描过程中，被扫描样品的运动以及OPT系统的机械不稳定性都会导致重建图像中产生运动伪影。运动伪影会导致系统空间分辨率下降、产生图像伪影甚至使重建图像无法正常使用。为了抑制运动伪影，需要在图像重建过程中进行校正。主要方法通过加权以减少一致性最差的投影数据对图像的贡献，消除或减弱运动带来的影响。伊利克股份有限公司在其申请的专利“CT扫描中运动伪影的减少”(专利申请号 200780053157. 9，公开号CN101842807A)中公开了一种丢弃包含运动结构的投影图像以获得投影图像的一致子集的方法。该方法是通过自动检测进而标识出需要被丢弃的有角区域，该方法可以在处理偶发性运动的情况下进行图像伪影校正。该方法仍存在的不足是，这种方法只能抑制小幅度运动引起的伪影，由于丢弃了包含运动结构的投影图像，会降低图像信噪比。U. J. Birk 等人在 “Correction for specimen movement and rotation errors for in—vivo Optical Projection Tomography, Biomedical Optics Express, vol.1， no. l，pp. 87-96,2010. ”中提出了通过准确估计出样品的运动轨迹，从而在图像重建过程中进行补充以实现运动伪影的校正。主要方法是使用基于特征点的运动估计方法。该方法需要样品的投影数据中具有明显易辨识的特征点，然后利用手工、自动或半自动的方法提取这些特征点，利用特征点的运动轨迹估计扫描样品的运动信息。该方法的不足是，需要估计特征点的运动轨迹，不适用于不易辨识的特征点。
技术实现思路
本专利技术的目的在于克服上述现有技术的不足，提出了一种光学投影断层成像运动伪影校正的方法。采用多项式函数对被扫描样品的运动进行建模，利用数据一致性条件建立样品与投影数据之间的对应关系，进而直接从投影数据中估计出样品的运动信息，应用于样品的投影断层重建。为实现上述目的，本专利技术的具体步骤如下(1)获取投影数据la)照射源对固定在电控旋转台上的样品进行水平投影断层扫描；lb)使用探测器采集样品的投影数据。(2)计算机计算下式中投影数据的零阶距CQ{t,)= [ g{e{t,),l)dlJ-co其中，CtlUi)为投影数据的零阶距，￡为正负无穷范围内的一维积分，g( θ (ti)，l)为步骤(1)中采集的投影数据，θ (ti)为在ti时刻的样品的扫描角度，1为探测器各个像素与探测器中心像素的距离。(3)计算机计算下式中投影数据的一阶距C1(^i)= Γ g(e(t,),l)ldlJ-oO其中，C1Ui)为投影数据的一阶距，￡为正负无穷范围内的一维积分，g( θ (ti)，l)为步骤(1)中采集的投影数据，θ (ti)为在ti时刻的样品的扫描角度，1为探测器各个像素与探测器中心像素的距离。(4)计算机计算下式中投影数据的质心Q,=缴其中，Qi为投影数据的质心，C1Ui)为投影数据的一阶距，CtlUi)为投影数据的零阶距。(5)获取运动参数5a)计算机按照下式建立投影数据质心与扫描样品质心在探测器上投影的一致性条件方程Q1 = (/^+^/,+…+ /^//^(^^(。+ (朽,。+ ^/, +…+ 朽乂广)sin 0(0其中，&为投影数据的质心，plj0, P1j1, · · ·，P1,N，p2j0, P2j1, · · ·，P2,N是待确定的运动参数，、为扫描时刻，N为运动参数多项式的阶数，cos θ (ti)为、时刻扫描角度的余弦值，sin θ (ti)为、时刻扫描角度的正弦值；5b)当投影数据个数大于待定的运动参数的个数时，使用最小二乘法从步骤5a) 中得到的一致性条件方程中获取运动参数。(6)当样品的运动曲线是光滑连续的并且样品只发生平动不发生转动时，计算机按照下面的公式计算每个扫描时刻样品的运动量dAO = PiJi + .·. + PuAn< Ndy{t,)= P2Ji+ ■.. + P2^t,其中，C^ti)和(IyUi)分别为样品的水平运动分量和竖直运动分量的多项式表达式，、为扫描时刻，pia, ...，Pl,N，p2a, . . .，P2,N为步骤(5)中获取的运动参数，N为运动参数多项式的阶数。(7)运动伪影校正将样品的投影重建图像中减去步骤(6)中获取的样品运动量， 实现光学投影断层重建运动伪影校正。本专利技术与现有的技术相比具有以下优点第一，本专利技术采用估计样品运动轨迹的方法对投影图像进行处理，由于重建时使用了所有的投影数据，克服了现有技术降低图像信噪比的缺点，使得本专利技术重建图像的信噪比较高。第二，本专利技术估计样品运动轨迹时使用的是不基于特征点的运动估计方法，克服了现有技术需要明显易辨识特征点的缺点，能够直接从投影数据中提取样品的运动参数信息，不要求投影数据中具有明显易于辨识的特征点，因此本专利技术使用起来比较方便、快捷。附图说明图1为本专利技术的流程图；图2为本专利技术获取的样品运动量随时间变化的估计曲线；图3为未经运动伪影校正的重建结果；本文档来自技高网...

【技术保护点】

【技术特征摘要】
1. 一种光学投影断层成像运动伪影校正方法，包括如下步骤(1)获取投影数据la)照射源对固定在电控旋转台上的样品进行水平投影断层扫描； lb)使用探测器采集样品的投影数据；(2)计算机计算下式中投影数据的零阶距2.根据权利要求1所述的光学投影断层成像运动伪影校正方法，其特征在于照射源采用激光器。3.根据权利要...

【专利技术属性】
技术研发人员：朱守平，陈冬梅，梁继民，陈多芳，屈晓超，赵恒，李军，陈雪利，侯彦宾，田捷，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：