[发明专利]光学投影断层成像运动伪影校正方法

权利要求书

1.一种光学投影断层成像运动伪影校正方法，包括如下步骤：

(1)获取投影数据

1a)照射源对固定在电控旋转台上的样品进行水平投影断层扫描；

1b)使用探测器采集样品的投影数据；

(2)计算机计算下式中投影数据的零阶距：

C0(ti)=∫-∞∞g(θ(ti),l)dl]]>

其中，C₀(t_i)为投影数据的零阶距，为正负无穷范围内的一维积分，g(θ(t_i)，l) 为步骤(1)中采集的投影数据，θ(t_i)为在t_i时刻的样品的扫描角度，l为探测器各个 像素与探测器中心像素的距离；

(3)计算机计算下式中投影数据的一阶距：

C1(ti)=∫-∞∞g(θ(ti),l)ldl]]>

其中，C₁(t_i)为投影数据的一阶距，为正负无穷范围内的一维积分，g(θ(t_i)，l) 为步骤(1)中采集的投影数据，θ(t_i)为在t_i时刻的样品的扫描角度，l为探测器各个 像素与探测器中心像素的距离；

(4)计算机计算下式中投影数据的质心：

Qi=C1(ti)C0(ti)]]>

其中，Q_i为投影数据的质心，C₁(t_i)为投影数据的一阶距，C₀(t_i)为投影数据的零 阶距；

(5)获取运动参数

5a)计算机按照下式建立投影数据质心与扫描样品质心在探测器上投影的一致性 条件方程：

Qi=(p1,0+p1,1ti+...+p1,NtiN)cosθ(ti)+(p2,0+p2,1ti+...+p2,NtiN)sinθ(ti)]]>

其中，Q_i为投影数据的质心，p_1，0，p_1，1，...，p_1，N，p_2，0，p_2，1，...，p_2，N是待确定的运动参 数，t_i为扫描时刻，N为运动参数多项式的阶数，cosθ(t_i)为t_i时刻扫描角度的余弦值， sinθ(t_i)为t_i时刻扫描角度的正弦值；

5b)当投影数据个数大于待定的运动参数的个数时，使用最小二乘法从步骤5a) 中得到的一致性条件方程中获取运动参数；

(6)当样品的运动曲线是光滑连续的并且样品只发生平动不发生转动时，计算 机按照下面的公式计算每个扫描时刻样品的运动量：

dx(ti)=p1,1ti+...+p1,NtiNdy(ti)=p2,1ti+...+p2,NtiN]]>

其中，d_x(t_i)和d_y(t_i)分别为样品的水平运动分量和竖直运动分量的多项式表达 式，t_i为扫描时刻，p_1，1，...，p_1，N，p_2，1，...，p_2，N为步骤(5)中获取的运动参数，N为运动 参数多项式的阶数；

(7)运动伪影校正：将样品的投影重建图像中减去步骤(6)中获取的样品运动 量，实现光学投影断层重建运动伪影校正。