工業 CT 系統幾何參數校正

陳志祥 須穎 石錦洋

摘要：工業計算機斷層成像技術（computer tomography，CT）系統幾何參數不準確會影響重建圖像的質量和尺寸測量的精度。為校正工業 CT 幾何位姿失配參數，提出基于標準球板模體的幾何校正方法。首先，利用球板在兩個不同位置的投影圖求解出射線源到探測器的距離以及射線源到運動平臺的距離。再在360°范圍內間隔一定角度均勻采集球板投影信息，擬合球心投影軌跡，解析求解探測器偏角η、φ和探測器中心點坐標。最后根據重建出的體素模型的球心距測量結果，利用牛頓迭代法，以工業 CT 尺寸測量最大允許誤差為約束條件迭代求解探測器偏角θ。實驗結果顯示：該文幾何校正方法可以求解出所有工業 CT 幾何位姿失配參數。參數校正后，提升圖像質量，消除重建圖像偽影，且工業 CT 球心距測量最大誤差由44.8μm 降低到16.3μm，表明該幾何校正方法可以滿足實際工業 CT 系統校正的要求。

關鍵詞：工業 CT;幾何校正;校正模體;圖像重建

中圖分類號： TP391.4文獻標志碼： A文章編號：1674–5124（2021）12–0023–06

Geometric parameter calibration of industrial CT system

CHEN Zhixiang，XU Ying，SHI Jinyang

（School of Mechanical and Electrical Engineering， Guangdong University of Technology， Guangzhou 510006， China）

Abstract： Incorrect geometric parameters of the industrial CT （computer tomography） system will affect the qualityof thereconstructedimageandtheaccuracyof dimensionalmeasurement. Inordertocalibrate geometric mismatching parameters of industrial CT， a geometric calibration method based on the standard spherical plate phantom was proposed. Firstly， the distance from the ray source to the detector and the distance from the ray source to the rotating platform were calculated by using the projection of the spherical plate at two different positions. Then， the projection information of the spherical plate was uniformly collected at a certain angle interval within 360° range， fitted the sphere center projection trajectory， and the detector deflection angle η ， φ and the detector center point coordinates were calculated analytically. Finally， based on the spherical- center distance measurement results of the voxel model， the Newton iteration method was used to solved the detectordeflectionangle θ bytakingtheMaximumPermissibleErrorof theindustrialCTdimensional measurementastheconstraintcondition. Theexperimentalresultsshowedthatthegeometriccorrection method could solved all the geometric position mismatch parameters of industrial CT. After the parameters were calibrated， the image quality was improved， and the reconstructed image artifacts were eliminated. The maximum error of industrial CTspherical distance measurement was reduced from 44.8μm to 16.3μm，indicated that the geometric calibration method could met the requirements of actual industrial CT system calibration.

Keywords： industrial computer tomography; geometry calibration; calibration phantom; image reconstruction

0引言

工業計算機斷層成像技術（ computer tomography，CT）是一門集光學、機械、控制和計算機科學于一體的數字化 X 射線成像技術，其利用物體對 X 射線吸收率的不同來獲取物體的多角度投影圖，并根據這些投影圖重建出物體的三維圖像[1]。目前這項技術正在向尺寸測量領域發展，但因為工業 CT 通常采用 FDK 算法對投影數據進行重建，而 FDK（ Feldkamp-Davis-Kress）算法要求工業 CT系統幾何位姿滿足理想條件[2]，實際情況下，工業 CT 系統難于滿足理想條件，系統幾何參數會存在偏差。在圖像重建過程中，如果不進行校正，會影響工業 CT 尺寸測量的精度。

針對這個問題，Noo[3]利用其設計的兩點模體，通過獲取模體多角度的投影數據，擬合橢圓運動軌跡，根據橢圓的相關參數求解系統幾何參數，在忽略探測器偏角θ對其他參數的影響的情況下，精確求解其余參數，但實際工業 CT 系統中偏角θ是一定存在的，且也會影響圖像的重建質量。Jacobson[4]設計了一種金屬絲校準模型，利用兩條不共面的直線在不同角度的投影平面內的交點位置的不同，根據這些交點的軌跡來計算出系統幾何參數，但該方法需要準確知道校正模體中每一根金屬絲準確的空間位置，在實際過程中，很難精確定位空間中金屬絲的位置，且對空間中金屬絲線位置測量的難度明顯高于球板球心位置測量的難度。Zhao[5]利用含有九個標準球的正方形模體，采集模體單個角度的投影數據，根據標準球的投影位置關系迭代求出系統幾何參數，雖能通過迭代求解出系統所有的失配參數，但該方法僅根據球板一個位置的投影關系求解，其求解精度難以保證。目前，工業 CT 系統參數校正的方法[6-14]都是利用特定校正模體，采集模體不同角度的投影數據，根據投影數據建立起與工業CT 系統幾何位姿參數的映射關系，再求解出相應的參數。

本文采用高精度球板模體，采集球板模體在不同角度，不同位置下的二維投影圖，求出不同情況下球板各球體球心在探測器上的投影坐標，擬合運動軌跡，再根據投影關系，利用解析法精確求解除θ外的失配參數。之后，用帶參數的 FDK 重建算法軟件重建出的球板三維圖像，將其導入可視化軟件 VG Studio MAX 3.0中并提取邊緣，測量球體球心距，分析工業 CT 測量誤差，以工業 CT 尺寸測量最大允許誤差為約束條件，利用牛頓迭代法求解系統最后一個失配參數。相較于其他校正方法，本方法可求解出工業 CT 系統所有的失配參數，且用工業 CT 尺寸測量誤差來對比校正前后的效果，更為直觀。

1工業 CT 系統和幾何參數

工業 CT 系統主要由射線源、旋轉平臺和探測器構成，根據 FDK 算法的要求，理想的工業 CT 系統幾何結構應滿足：射線源焦點和探測器中心點的連線經過旋轉軸，且垂直于探測器平面;旋轉軸與探測器中心列平行。工業 CT 系統理想位姿關系如圖1所示。其中f 為射線源到旋轉軸的距離;d 為射線源到探測器中心的距離;（u0 ， v0）為射線源焦點在探測器上的垂直投影點。

為更好地研究工業 CT 系統的幾何位姿關系，建立三維坐標系O? XYZ，如圖1所示。X 軸沿射線方向;Y 軸沿水平方向垂直于射線方向且與探測器中心行平行;Z 軸由右手坐標系原則確定。同樣在探測器所在平面上建立二維坐標系O?uv，u 軸沿探測器中心行方向;v 軸沿探測器中心列方向。工業 CT 系統在安裝調試過程中，系統位姿關系一定會存在偏差，工業 CT 的校正即校正這些參數。經國內外眾多研究人員對工業 CT 校正的研究，需要校正的參數可概括為：探測器繞 X、Y、Z 軸的轉動角度η、θ、φ;射線源焦點在探測器上的垂直投影點（u0 ， v0）和在 u、v 方向上的偏移Δu、Δv;源到旋轉中心的距離以及源到探測器中心的距離。其中和探測器有關的5個參數如圖2所示。

2參數校正

2.1校正模體

本文使用的校正模體是球板標準器，如圖3所示，由12個高精度的紅寶石球體和碳纖維平板組成，部分球體之間的球心距已知。

2.2計算d、f

在保持管電壓、管電流和曝光時間等掃描參數不變的情況下，僅改變探測器到運動平臺的距離，采集兩組不同位置時的球板二維投影圖。擬合球板各球體在二維投影圖上的圓心位置，再計算圓心距 L，L 為像素數量與像素大小的乘積，即

式中：L——被測量，mm;

n——像素數量;

M——放大比;

D——像元尺寸，mm。

根據式（1），利用上述采集的兩組投影圖數據，利用球板各球體間球心距，即可建立多個關于 d 和 f 的二元一次方程，在精確知道運動平臺沿射線源方向上的運動距離的情況下，可求解出多個 d、f值，最后取 d、f 的平均值作為最終校正結果。

2.3計算η、u0、v0

采集 N（N 為偶數）張1號球和12號球每隔一定角度的投影圖，計算出1號球和12號球在不同角度投影圖下的球心位置（ui ， vi）（i =1， 2，3， ··· ， N）。并根據這些球心坐標擬合運動軌跡，可以得到兩個橢圓（橢圓1和橢圓2），擬合出的橢圓可用橢圓方程式表示如下：

式中：（ ， ）——橢圓中心點坐標;

a、b、c——橢圓方程參數。

由式（3）分別計算出1號球和12號球球心連線在每隔180°投影圖下的中心點位置（1 ， 1）、（2 ， 2），如圖4所示。

其中j = i+ N/2;i =1， 2， ··· ， N/2。

探測器面內扭轉角η即可由式（4）計算得出。探測器平面內射線源焦點的垂直投影點（u0 ， v0）可由式（5）、（6）計算出。

2.4計算φ

Sun[15]等利用四球模體校正系統幾何參數，提出平行于轉軸的兩球投影圓心距之間滿足以下關系：

式中：——平行于轉軸的兩球投影圓心距，

——實際圓心距，mm。

根據這一原理，將 1 號球到 10 號球的球心距和 3 號球到 12 球的球心距 SD3-12 代入式 （7），可計算出。

2.5 計算

根據上述計算出的探測器失配參數對在理想坐標系下的 1 號球球心坐標進行坐標變換得，則：

其中：

經過坐標變換后得到：

射線源焦點坐標S =（0， 0， 0）經坐標變換后，射線源焦點在理想坐標系下的坐標應為S＼=（xs ，ys ，zs）;將 S＼與p 相連交探測器于點p = v ＼），p ＼即為1號球球心經 X 射線投射后在探測器上的投影點。其中：

同理可計算出其他球體球心經 X 射線投射后在探測器上的投影點坐標。

此時使用帶參數的 FDK 重建算法軟件，利用上面所計算出的6個系統失配參數，重建出此時球板模體的三維圖像，再利用 VG Studio MAX 3.0圖像可視化軟件測量各球體間的球心距，找出此時球心距誤差最大的兩球體。以1號球和12號球為例，此時兩球球心在探測器上的投影點間距離為：

式（12）中，除θ外的其他參數均已知，故S D1?12= f （θ）是一個關于θ的函數，解析求解過于復雜，可用牛頓迭代法求解，以工業 CT 尺寸測量最大允許誤差為迭代約束條件，迭代求解出θ。

3實驗驗證

3.1實驗設備及參數

為驗證上述方法的可行性，在實際工業 CT 系統中進行了實驗，實驗所使用的設備為天津三英精密儀器有限公司生產的 nanoVoxel-4000，如圖5所示。射線源為開管發射式，最大管電壓為220 kV，最大管電流為1.0 mA;探測器成像面積為427 mm ×427 mm，像素矩陣為3072×3072;像素細節分辨力為500 nm，空間分辨率為2?m。

3.2校正過程

首先驗證 d 和f，在 d 固定的情況下，改變f其他參數如表1所示，采集兩組投影圖。表1 初始參數

通過計算可得，射線源到探測器的距離誤差為–5.2029 mm，射線源到旋轉軸的距離誤差為–4.3907 mm。

在調整工業 CT 系統 d 和f后，再次掃描校正模體，保持 d=1500 mm，f=500 mm，在其他參數不變的情況下360°范圍內采集球板模體的投影圖，然后根據本文方法求解圓心，擬合橢圓軌跡如圖6所示，橢圓參數計算出工業 CT 其余幾何參數，最后利用迭代法經4次迭代后計算出工業 CT 需校正的最后一個參數。

至此，所有參數校正完畢，最終校正結果如表2所示。

3.3校正前后效果對比

實際重建過程中，利用帶參數的 FDK 重建算法軟件重建球板模體，校正前后的重建圖像對比如圖7所示，可以看出，校正后球板模體的重建圖像質量有了明顯提升，消除了圖像偽影，提升了圖像清晰度。

利用 VG Studio MAX 3.0圖像可視化軟件分別對校正前后重建的球板模體三維圖像進行尺寸測量，兩次測量結果如圖8所示?？芍Ｕ蠊I CT 尺寸測量的最大誤差由44.8?m 降低到16.3?m，效果明顯，本文工業 CT 系統參數校正的方法降低了工業 CT 尺寸測量的誤差。

4結束語

本文基于高精度球板模體，提出了一種工業 CT 幾何參數校正的方法，利用球板在不同角度下的投影圖像，獲取球心投影運動軌跡，擬合橢圓，再根據橢圓參數和工業 CT 尺寸測量的最大誤差，結合解析法與迭代法，可求解出工業 CT 系統所有幾何失配參數。實驗表明，校正后重建圖像的質量有了明顯的提升，減少了圖像偽影，校正后工業 CT 尺寸測量最大誤差由44.8?m 降低到16.3?m，能滿足實際工業 CT 系統參數的校正要求。但是，本文幾何校正方法需要高精度的校正模板，需要準確知道球板模體各球心距，也受限于運動平臺回轉精度和沿射線源方向的運動精度，否則會影響校正參數的準確性。

參考文獻

[1]張朝宗.工業 CT 技術和原理[M].北京：科學出版社， 2009.[2] FELDKAMP L A， DAVIS L C， KRESS J W. Practical cone-beam algorithm[J]. Journal of the Optical Society of America A， 1984， 1（6）：612-619.

[3] NOO F， CLACKDOYLE R， MENNESSIER C， et al. Analyticmethodbasedonidentificationofellipseparametersfor scanner calibration in cone-beam tomography[J]. Physics in Medicine and Biology， 2000， 45（11）：3489-3508.

[4] JACOBSON M W， KETCHA M D， CAPOSTAGNO S， et al.A line fiducial method for geometric calibration of cone-beam CTsystemswithdiversescantrajectories[J]. Physicsin Medicine and Biology， 2018， 63（2）：025030.

[5] ZHAO J T， HU X D， ZOU J，etal. Geometric parametersestimationandcalibrationincone-beammicro-CT[J]. Sensors， 2015， 15（9）：22811-22825.

[6]周凌宏， 李翰威， 徐圓， 等.錐束 CT 圓軌道掃描的幾何校正[J].光學精密工程， 2014， 22（10）：2847-2854.

[7] SMEKALLV，KACHELRIEM，STEPINAE，etal.Geometricmisalignmentandcalibrationincone-beam tomography[J]. Medical Physics， 2004， 31（12）：3242.

[8] CHO Y， MOSELEY D J，SIEWERDSEN J H， et al. Accuratetechniqueforcompletegeometriccalibrationof cone-beam computedtomographysystems[J]. MedicalPhysics， 2005，32（4）：968-983.

[9] XIANG Q， WANG J， CAI Y F. A geometric calibrationmethodforconebeamCTsystem[C]//EighthInternational ConferenceonDigitalImageProcessing （ICDIP 2016）. International Society for Optics and Photonics， 2016.

[10] PATEL V， CHITYALA R N， HOFFMANN K R， et al. Self- calibrationofacone-beammicro-CTsystem[J]. MedicalPhysics， 2008， 36（1）：48-58.

[11]張才鑫， 黃魁東， 陳華.基于標定量具的快速錐束 CT 幾何校準方法[J].儀器儀表學報， 2019， 40（1）：175-182.

[12] ZHOU K， HUANG Y X， MENG X X， et al. A new method for cone-beamcomputedtomographygeometricparameters estimation[J]. JournalofComputerAssistedTomography， 2016， 40（4）：639-648.

[13] LI G， LUO S H， YOU C Y， et al. A novel calibration method incorporatingnon-linearoptimizationandball-bearing markers for cone-beam CT with a parameterized trajectory[J]. Medical Physics， 2019， 46（1）：152-164.

[14] LI H C， HONG W， LIU Y， et al. A novel method of micro- tomographygeometricanglecalibrationwithrandom phantom[J]. Journal of X-ray science and technology， 2017，25（4）：641-652.

[15] SUN Y， HOU Y， ZHAO F Y， et al. A calibration method for misalignedscannercogeometryinne-beamcomputed tomography[J]. NDT and E International， 2006， 39（6）：499-513.

（編輯：莫婕）