一种CT投影中心的自动校正方法

本发明提供一种CT投影中心的自动校正方法，包括如下步骤：1)使用CT对被测物进行全周扫描，得到在不同角度下被测物在探测器阵列的各探测器单元的投影原始数据；2)对所述投影原始数据进行数据处理，找出在所述被测物边缘的投影点在所述探测器阵列中的位置；3)根据所述被测物边缘的投影点在所述探测器阵列中的位置计算出CT投影中心并对所述投影原始数据进行校正，得到最终的二维断层数据。本发明无需采用专用模体，无需移动被测物，即可对CT的中心投影位置进行校正，而且本发明的各校正步骤可以编入重建前处理程序自动运行，节约了测量时间，延长了射线源及探测器的有效工作寿命。

1.一种CT投影中心的自动校正方法，包括如下步骤：1)使用CT对被测物进行全周扫描，得到在不同角度下被测物在探测器阵列的各探测器单元的投影原始数据；2)对所述投影原始数据进行数据处理，找出在所述被测物边缘的投影点在所述探测器阵列中的位置；3)根据所述被测物边缘的投影点在所述探测器阵列中的位置计算出CT投影中心并对所述投影原始数据进行校正，得到最终的二维断层数据；所述步骤2)包括如下子步骤：21)将投影原始数据以矩阵P的形式存取，所述矩阵P中的元素以p_ab表示；所述矩阵P中行向量a为同一角度下不同探测器单元的数据，列向量b为同一探测器单元不同角度下的数据；找出所有投影原始数据中的最大值p_max；22)根据该最大值p_max得出未发生截断扫描的衰减阈值p_valve；设置未发生截断扫描的阈值系数k，所述衰减阈值p_valve＝kp_max；23)当探测器单元的投影原始数据小于衰减阈值p_valve时，判定该探测器单元探测到所述被测物；遍历矩阵P，求得符合条件p_ab＜p_valve时，b的极小值b_min和极大值b_max；所述步骤3)包括如下子步骤：31)计算投影中心p_center＝(b_max+b_min)/2；32)计算N_offset＝p_center-d/2，d为探测器单元的总数，N_offset为需对原始数据进行插值的列数；33)对投影原始数据矩阵P进行插值，插值方法是：在矩阵P中插入一N_offset×b的插值矩阵，该插值矩阵中所有元素的值等于p_max；N_offset为正时，表示在矩阵P的右端插值，N_offset为负时，则表示在矩阵P的左端插值；34)对插值后的矩阵进行对数运算，并进行数值重建，得到最终的二维断层数据。
2.根据权利要求1所述的CT投影中心的自动校正方法，其特征在于，所述阈值系数k的取值范围为：0.5＜k＜1。