CT系統(tǒng)參數(shù)標(biāo)定

1、CT系統(tǒng)參數(shù)標(biāo)定問題研究摘要本文通過對(duì)體模系統(tǒng)和投影點(diǎn)坐標(biāo)的分析，提出了較為充分合理的假設(shè)，并對(duì)數(shù)據(jù)進(jìn)行了擬合處理。先后建立了三個(gè)模型（理想狀態(tài)下的模型，只考慮載物臺(tái)傾斜的模型，考慮載物臺(tái)傾斜與探測(cè)器的俯仰角的模型）來求解CT系統(tǒng)參數(shù)標(biāo)定的問題。在對(duì)問題一的分析中，只考慮在理想情況下：即載物臺(tái)水平和探測(cè)器無偏轉(zhuǎn)，并不考慮機(jī)械誤差的情況下各個(gè)參數(shù)之間的關(guān)系，以X光放射源為原點(diǎn)建立坐標(biāo)系，沿水平方向射向探測(cè)器的水平直線為X軸,豎直方向?yàn)閅軸,建立模型一,繪出側(cè)視圖如一所示。運(yùn)用幾何知識(shí)作圖，通過建立相似圖形的等比例方程來研究確定幾個(gè)系統(tǒng)參數(shù)之間的關(guān)系（X射線源的位置參數(shù)、載物中心線的位置參數(shù)和探測(cè)

2、器的位置參數(shù)）。同時(shí)求出載物臺(tái)半徑等參數(shù)，最后討論了引入像素、載物臺(tái)偏角、觀測(cè)器俯仰角和數(shù)據(jù)不確定度對(duì)誤差的影響。在分析問題二時(shí)，考慮到所給橢圓點(diǎn)集可能有偏差，于是首先運(yùn)用擬合的方法對(duì)數(shù)據(jù)行處理，得到更為確信的數(shù)據(jù)并由此求出修正后的“橢圓”方程（觀察均方根誤差,得知此方程能符合條件，可以使用）。由于考慮多個(gè)誤差因素將使模型計(jì)算變得極為復(fù)雜，在此不同時(shí)考慮。假設(shè)探測(cè)器無偏轉(zhuǎn)，只有載物臺(tái)傾斜，建立模型二。易知傾斜的載物臺(tái)與安裝沒有偏差的水平觀測(cè)器可以等效為水平的載物臺(tái)與只縱向傾斜的觀測(cè)器。然后根據(jù)擬合后的數(shù)據(jù)求出兩類橢圓長(zhǎng)軸斜率，從而求出載物臺(tái)傾斜角為0.1146°,并依據(jù)此進(jìn)一步修正數(shù)

3、據(jù)。在模型二修正數(shù)據(jù)的基礎(chǔ)上，考慮探測(cè)器的俯角因素（易知0不影響模型二的求解）建立模型三。由修正后的圖知兩“橢圓”長(zhǎng)軸長(zhǎng)等。由上方長(zhǎng)軸長(zhǎng)大于下方長(zhǎng)軸長(zhǎng)知探測(cè)器向前倒。經(jīng)過一定數(shù)據(jù)近似后，再根據(jù)等比例關(guān)系可求得探測(cè)器俯角q為1.623。近一步推導(dǎo)可得坐標(biāo)及光源到旋轉(zhuǎn)臺(tái)的旋轉(zhuǎn)軸的距離為263.97mm,光源在成像平面上的水平直射點(diǎn)，光源到探測(cè)器水平直射點(diǎn)的間距為688.57mm,。最后在問題三，對(duì)模型三的機(jī)械誤差進(jìn)行了分析計(jì)算。機(jī)械誤差主要來源為像素、載物臺(tái)偏角，探測(cè)器的俯仰角，圓心距離的誤差。并進(jìn)行了靈敏度分析。關(guān)鍵詞：坐標(biāo)修正載物臺(tái)偏角探測(cè)器俯仰角機(jī)械誤差一、問題的重述掃描物件圍繞某一固定轉(zhuǎn)軸

4、旋轉(zhuǎn)，每隔一定角度采集一張圖像，然后根據(jù)采集的圖像采用3D圖像重建算法即可將原始3D物件重建出來。在實(shí)際中，由X射線源發(fā)出X射線，經(jīng)過掃描物體衰減后照射在探測(cè)器上，探測(cè)器根據(jù)接收到的光子數(shù)的計(jì)數(shù)實(shí)現(xiàn)光電轉(zhuǎn)換,從而形成灰度圖像。一般平板探測(cè)器大小為3000X2000像素，每個(gè)像素為0.127mm。在CT系統(tǒng)安裝過程中往往存在機(jī)械誤差，而這些誤差對(duì)于物件重建的準(zhǔn)確性往往是致關(guān)重要的，實(shí)際中就需要對(duì)安裝好的CT系統(tǒng)進(jìn)行參數(shù)標(biāo)定，系統(tǒng)參數(shù)主要包括X射線源的位置、載物中心線的位置和探測(cè)器的位置參數(shù)等。對(duì)于CT系統(tǒng)參數(shù)標(biāo)定精度要求是不夠的（通常CT重建過程中系統(tǒng)參數(shù)的機(jī)械誤差不允許超過一個(gè)像素）。通過掃描

5、已知參數(shù)的體模，分析投影數(shù)據(jù)，估計(jì)系統(tǒng)參數(shù)的方法。針對(duì)工業(yè)CT要求解決如下問題：（1）建立合適的坐標(biāo)系，正確描述CT系統(tǒng)的各種參數(shù)和機(jī)械誤差，并建模分析這些參數(shù)的關(guān)系和可能的機(jī)械誤差。（2）通常采用軸承鋼球作為CT參數(shù)標(biāo)定體模。因?yàn)樗哂懈飨蛲队耙恢?，邊緣清晰等?yōu)點(diǎn)，在里采用兩個(gè)鋼球?qū)崿F(xiàn)CT的參數(shù)標(biāo)定實(shí)驗(yàn).兩個(gè)鋼球置于有機(jī)玻璃管中（X射線容易透過有機(jī)玻璃，容易后期球心獲取的圖像處理），鋼球直徑為8mm±0.0008mm,兩球心距離約為100mm±lmmo將有機(jī)玻璃管固定在旋轉(zhuǎn)載物臺(tái)上，載物臺(tái)攜帶有機(jī)玻璃管以均勻速度旋轉(zhuǎn)，與此同時(shí)X射線源發(fā)出的射線，探測(cè)器采集數(shù)據(jù)（旋轉(zhuǎn)一周，

6、等間距采集180張）。兩個(gè)鋼球球心投影的軌跡是兩個(gè)橢圓。試根據(jù)實(shí)驗(yàn)數(shù)據(jù)估計(jì)該CT系統(tǒng)的參數(shù)值，給出CT系統(tǒng)的標(biāo)定。（3）在參數(shù)標(biāo)定過程中可能存在多種可能的機(jī)械誤差，對(duì)參數(shù)標(biāo)定可能的誤差進(jìn)行分析。二、問題的分析問題一：本題誤差主要來源于工業(yè)CT系統(tǒng)安裝過程。除了裝置本身的由于在裝系統(tǒng)時(shí)，可能會(huì)導(dǎo)致探測(cè)器扭轉(zhuǎn)、前后俯仰，載物臺(tái)的傾斜等，都會(huì)導(dǎo)致機(jī)械誤差的存在。我們?cè)谟?jì)算時(shí)會(huì)考慮這些因素,并建立相應(yīng)的模型。所以，在標(biāo)定時(shí)應(yīng)該給出其機(jī)械誤差大小?？梢愿鶕?jù)體模的投影數(shù)據(jù)，利用幾何知識(shí)來求解主要包括X射線源的位置、載物中心線的位置和探測(cè)器的位置并求出他們的關(guān)系。對(duì)于載物臺(tái)來說，由于載物臺(tái)要進(jìn)行旋轉(zhuǎn)，所以其

7、傾斜方式只有一種。問題二：根據(jù)分析，載物臺(tái)的傾斜可以等效為探測(cè)器投影繞一個(gè)點(diǎn)在豎直方向旋轉(zhuǎn)（圖三所示）。由附錄一中圖分析可知，兩個(gè)橢圓的長(zhǎng)軸不是水平的，這是由于載物臺(tái)的傾斜造成的：上橢圓與下橢圓的長(zhǎng)軸不一樣長(zhǎng)，這是由于探測(cè)器的前后俯仰造成的:“橢圓”圖形左右不對(duì)稱，是由于探測(cè)器的扭轉(zhuǎn)造成的。我們分別建立模型將數(shù)據(jù)一步一步進(jìn)行處理，得到機(jī)械誤差來源，再進(jìn)行修正，從而給出準(zhǔn)確的CT參數(shù)標(biāo)定。問題三：本題中考慮的機(jī)械誤差來源有兩圓心距離誤差，載物臺(tái)傾斜，探測(cè)器前俯后仰，利用問題二中所得數(shù)據(jù)，可求得總機(jī)械誤差。三、模型的假設(shè)與符號(hào)說明1 .模型的假設(shè)（1）考慮到此處光子的波動(dòng)不明顯，故忽略光的衍射現(xiàn)象

8、，以光的直線傳播進(jìn)行處理；（2）忽略X光輻射源的體積和形狀，將其視為質(zhì)點(diǎn)：（3） 一些題目中未說明但實(shí)際易得的參數(shù)，如載物臺(tái)的半徑等，假設(shè)可以測(cè)量得到，此處不考慮具體方法及相關(guān)誤差；（4）數(shù)據(jù)能客觀反映實(shí)際情況：（5）載物臺(tái)與有機(jī)玻璃管始終垂直；（6）兩小球圓心視為質(zhì)點(diǎn)；（7）在數(shù)據(jù)采集過程中，放射源點(diǎn)，載物臺(tái)，探測(cè)器參數(shù)不再變化。2 .符號(hào)說明X光射線源到載物臺(tái)中心垂線的距離；12：X光射線源到探測(cè)器的距離；13：探測(cè)器俯仰偏轉(zhuǎn)后下方橢圓最下沿點(diǎn)到Y(jié)軸距離：14：探測(cè)器俯仰偏轉(zhuǎn)后上方橢圓最上沿點(diǎn)到Y(jié)軸距離；r：有機(jī)玻璃管和小球的軌跡圓半徑；h：兩個(gè)小球重心之間的距離；%：橢圓間最近兩點(diǎn)的距離

9、：h,：橢圓間最遠(yuǎn)兩點(diǎn)的距離：n1：上投影橢圓短軸長(zhǎng)度；n2：上投影橢圓短軸長(zhǎng)度；m：投影的豎直方向最大長(zhǎng)度；載物臺(tái)的傾斜角；探測(cè)器偏轉(zhuǎn)的俯仰角；tl：大橢圓的長(zhǎng)軸長(zhǎng)度；t2：小橢圓的長(zhǎng)軸長(zhǎng)度；口2：兩橢圓中心之間的距離；的：當(dāng)偏差為一個(gè)像素時(shí)候，探測(cè)器的俯仰角度：At1：當(dāng)偏差為一個(gè)像素時(shí)候片的變化量；At3：當(dāng)偏差為一個(gè)像素時(shí)候t?的變化量；At12：當(dāng)偏差為一個(gè)像素時(shí)候G的變化量四、模型的建立和求解問題一：1、模型一的建立首先考慮在理想情況下，即載物臺(tái)水平和探測(cè)器無偏轉(zhuǎn)和不考慮機(jī)械誤差的情況下各個(gè)參數(shù)之間的關(guān)系。建立坐標(biāo)系時(shí)選定X光輻射源為坐標(biāo)原點(diǎn)，沿水平方向射向探測(cè)器的水平直線為X軸，

10、豎直方向?yàn)閅軸，繪出側(cè)視圖如一所示。圖二理想模型下的俯視圖由平面幾何中相似圖形成比例可以得出下列方程：hL+rm-(n1+n2)12(2)若考慮引入機(jī)械誤差，則可能存在以下幾種機(jī)械誤差：（1）載物臺(tái)與地面有一定傾斜角，此角會(huì)使投影偏轉(zhuǎn)。（2）探測(cè)器地面間存在一定俯仰角，此偏角的存在會(huì)使原本應(yīng)是橢圓的投影產(chǎn)生一定偏差：（3）x放射點(diǎn)與載物臺(tái)軸線成面與探測(cè)器成角不成90度，會(huì)使原本應(yīng)是橢圓的投影產(chǎn)生一定偏差：（4）像素不夠大引起的誤差，一般平板探測(cè)器大小為3000x2000,每個(gè)像素為0.127nin,若兩個(gè)投影點(diǎn)之間的水平或豎直距離小于0.127in】，在探測(cè)器上就會(huì)區(qū)分不出，例如以下2組數(shù)據(jù)：

11、上方球心坐標(biāo)（X,Y）：（1415,1639）,（1375,1639）；（5）題目中給出的數(shù)據(jù)（小球半徑、球心距離等）不精確產(chǎn)生的誤差：?jiǎn)栴}一.1、誦二的建立暫只考慮載物臺(tái)的傾斜角。的情況下，建立此模型。傾斜的載物臺(tái)與安裝沒有偏差的觀測(cè)器可以等效為水平的載物臺(tái)與只縱向傾斜的觀測(cè)器，如圖二由于數(shù)據(jù)可能有較大的誤差，我們用matlab對(duì)數(shù)據(jù)進(jìn)行擬合。2000180016001400120010008006003000圖四：修正前的圖形2000修正后1800160014001200100080060040020050010003000150020002500圖五：擬合修正后的圖形再由圖中所示方法求出

12、長(zhǎng)軸短軸大小及對(duì)應(yīng)各個(gè)切點(diǎn)?！捌叫星芯€之間的距肉最短為短軸圖六：橢圓的切線與橢圓軸的關(guān)系再根據(jù)對(duì)應(yīng)切點(diǎn)確定長(zhǎng)軸斜率tanx=(lxj-rxj)/(lx-rx);%下方橢圓長(zhǎng)軸斜率tans=(Isj-rsj)/(Is-rs);專上方橢圓長(zhǎng)軸斜率tan=(tans+tanx)/2;得q)2=arctaiiO.OO2根據(jù)相應(yīng)切點(diǎn)求出圓心軌跡所成圓之圓心，水平射于探測(cè)器的點(diǎn)橫縱坐標(biāo)xsO=(ls+rs)/2;ysO=(lSj4-rsj)/2;%上方圓心軌跡圓圓心水平射于探測(cè)器的點(diǎn)橫縱坐標(biāo)xxO=(lx+rx)/2;yxO=(lxj+rxj)/2;%卜方圓心軌跡圓圓心水平射于探測(cè)器的點(diǎn)橫縱坐標(biāo)bili=

13、(gs-ds)/(gs-ds)+(gx-dx);根據(jù)附件一中所給數(shù)據(jù)，計(jì)算出X放射點(diǎn)水平射于探測(cè)器的橫縱坐標(biāo)yO=ds-(gx-ds)*乂放射點(diǎn)水平射于探測(cè)器的縱坐標(biāo)xO=dsi-(gxi-dsi)放射點(diǎn)水平射于探測(cè)修正后坐標(biāo)軸不聲前坐標(biāo)軸、(Xo.Vo>1/cos6=y2圖七修正的原理圖再根據(jù)對(duì)應(yīng)切點(diǎn)確定長(zhǎng)軸斜率h=abs(sqrt(tanA2+l)-(tan*(x(izj)-Is)-y(irj)+lsj)/sqrt(tanA2+l);%上方橢圓上任意一點(diǎn)至上方橢圓長(zhǎng)軸的距離y(irj)=yO+(ysO-yO)*sqrt(tanA2+l)-(tan*(x(izj)-Is)-y(irj)

14、+lsj)/sqrt(tanA2+l);%核正后的縱坐標(biāo)x(i,j)=xO+fuhao*hh;宅修正后的縱橫坐標(biāo)修正后圖形如圖五。詳細(xì)數(shù)據(jù)見附錄一。2、模型三的建立我們可以建立在模型二的基礎(chǔ)上考慮探測(cè)器的俯仰角度W，如圖八所示，由于載物臺(tái)的偏轉(zhuǎn)角度c與探測(cè)器俯仰角度W不影響，所以可以把他們結(jié)合起來計(jì)算?？紤]他們偏轉(zhuǎn)角的誤差。由圖八和圖可知，探測(cè)器向前偏轉(zhuǎn)了0,導(dǎo)致上橢圓和下橢圓大小不一樣，根據(jù)幾何知識(shí)，可以求得方程(4)和方程(5)：于是，我們就可以求出在探測(cè)器前后俯仰時(shí)候的偏轉(zhuǎn)角q與體模投影點(diǎn)坐標(biāo)之間的關(guān)系。找出偏轉(zhuǎn)角q與其他參數(shù)之間的關(guān)系，進(jìn)行修正。得到準(zhǔn)確的載物中心線的位置和探測(cè)器的位置

15、參數(shù)，較精確對(duì)CT進(jìn)行標(biāo)定。而“與t2的值非常相近，即q非常小。相似三角形對(duì)應(yīng)邊要小于對(duì)應(yīng)1cos/,但由與W非常小，相差不大，二者可約等于。則由圖三和圖四可以得到下列方程如下：21=皿叫)hjcos弘h2cos.(tjt)lnll(4)(6)w非常小。則可以將L，h，。的長(zhǎng)度看成一樣。即式子可改為:Acos（flh2cos.（tt»）L%5m.二J-；ll于是，我們就可以求出在探測(cè)器前后俯仰角例=1.62%由模型三的修正，我們可以得到準(zhǔn)確的數(shù)據(jù)，就可以求得：h=787m=1905=347n2=215%=2708t3=2685=1625.5帶入模型一的方程,由模型一，X射線源的位置在

16、扇點(diǎn)，即坐標(biāo)為（0：0）,由射區(qū)源到載物中心線的距離為1=2078.55（像素）=263.97皿，即位置坐標(biāo)為（0,263.97）,由射線源到探測(cè)器的距離為12=5421.85（像素）=688.5711111】，即位置坐標(biāo)為（0,688.57）,我們還求出在載物臺(tái)上物體繞載物臺(tái)中心線旋轉(zhuǎn)的半徑,由于小球在載物臺(tái)邊緣，則旋轉(zhuǎn)半徑約等于載物臺(tái)半徑r戶468.56（像素）=59.5mm問題三：根據(jù)前面的描述，誤差的來源主要有三種：像素、偏轉(zhuǎn)角和球心距。對(duì)于CT系統(tǒng)參數(shù)標(biāo)定精度要求是其機(jī)械誤差不超過一個(gè)像素，而由式（1）（2）（5）知球心距只與1?有直接聯(lián)系，與像素?zé)o直接聯(lián)系，故此處主要研究偏轉(zhuǎn)角與投

17、影像素之間的關(guān)系。有次可進(jìn)一步推出它們的關(guān)系如下：-tjTAt9;（8）cos一%=倏；（9）.（tj+Atj1->At、）1,tj2sin%=-=二-G+M；（10）計(jì)算（8）得到=2.009。；傷一%=0.386。：H算（10）得藥價(jià)=1.764°：的一%=0.1410；比較兩次結(jié)果，取最小值0.141。，即要保證投影在探測(cè)板上的點(diǎn)位誤差不超過一個(gè)相位，探測(cè)板的俯仰角浮動(dòng)范圍最大不能超過0.141。五、模型的推廣與改進(jìn)方向考慮三維直角坐標(biāo)系，能更精確的求出機(jī)械誤差，具有更好的實(shí)際應(yīng)用?？紤]探測(cè)器的扭轉(zhuǎn)角對(duì)體模投影圖像的影響?？紤]多種誤差同時(shí)影響。六、模型的優(yōu)缺點(diǎn)優(yōu)點(diǎn)：1、建

18、立多個(gè)模型，模型的建立由易到難，比較有層次，能更清晰的分析問題，找到它們之間的關(guān)系：2、對(duì)體模的投影點(diǎn)坐標(biāo)進(jìn)行數(shù)據(jù)的修正，計(jì)算的結(jié)果更加準(zhǔn)確：3、考慮到的誤差比較合理，有較大的實(shí)際意義；缺點(diǎn)：1、對(duì)探測(cè)器的扭轉(zhuǎn)角考慮不全面，僅在理論上進(jìn)行分析。2、由于像素造成的機(jī)械誤差比較大，像素精度不高，此方法解決這個(gè)問題效果欠佳；七、參考文獻(xiàn)1姜啟源.數(shù)學(xué)模型Ml.北京：高等教育出版社，1999.2韓中庚.數(shù)學(xué)建模方法及其應(yīng)用M.北京：高等教育出版社，2009.八、附錄附錄一：投影點(diǎn)坐標(biāo)的修正zhen=zeros(2000,3000);fori=l:180forj=l:2zhen(y(i,.i),x(i,

19、i)=l:endend%矩陣之生成rs=1000;lx=1000;rx=1000;fori=l:3000forj=l:2000ifzhenfj,i)=1ifi<lsifj>1400ls=i;lsj=j;endendifi<lxifj<500lx=i;lxj=j;endendifi>rsifj>1400rs=i;rsj=j；endendifi>rxifj<500rx=i;rxj=j;endendendendendgs=1800;ds=1800;fori=l:3000forj=l:2000ifzhenij,i)=1ifj>1000ifj<

20、dsds=j;dsi=i;endendendifzhen；j,i)=lifj>gsss=j;gsi=i；endendendendgx=200,dx=200；fori=l:3000forj=1:2000ifzhen;j,i)=1ifj<dxdx=j;dxi=i;endendifzhen'.j,i)=1ifj<1000ifj>gxgx=j；gxi=i;endendendendend%下方橢圓長(zhǎng)軸斜率%上方橢圓長(zhǎng)軸斜率%上方圓心軌跡圓圓心水平射于探測(cè)器的點(diǎn)橫縱坐%下方圓心軌跡圓圓心水平射于探測(cè)器的點(diǎn)橫縱tanx=(lxj-rxj)/(lx-rx);tans=(lsj-rsj)/(ls-rs);tan=(tans+tanx)/2;xs0=(ls+rs)/2;ys0=(lsj+rsj)/2;林xxO=：1x+rx)/2;yxO=(lxj+rxj)/2;坐標(biāo)bili=(gs-ds)/(gs-ds)+(gx-dx);yO=ds-(gx-ds)*bili;%x放射點(diǎn)水平射于探測(cè)器的縱坐標(biāo)forj=l:2xO=dsi-(gxi-dsi)*bi1i;%x放射點(diǎn)水平射手探測(cè)器的橫坐標(biāo)fori=l:180if(y(i,j)>1500)h=(tan*(x(i,j)-xsO)