摘要

中子成像技術在類似考古（古董、古生物）、冶金金屬領域獲得了廣泛關注，但其受到同步輻射光源限制，更有效的調(diào)配好使用光源的時間是急需解決的問題。N-Cam中子相機系統(tǒng)相對傳統(tǒng)的中子相機，能在較短的曝光時間獲得詳細清晰的圖像，是一種新型的針對冷中子、熱中子的成像系統(tǒng)。N-Cam將20微米厚的Gd2O2S:Tb（Gadox）閃爍體直接應用于像增強器的輸入窗口。N-Cam系統(tǒng)在英國盧瑟福.阿普爾頓國家實驗室散裂中子源設備進行了實驗，在75mm靶面視場、5s曝光情況下獲得了良好的對比度效果圖。此篇分析文章通過相機模式傳遞函數(shù)法（MTF）在多個位置對空間分辨率進行了方向性依賴的分析。邊緣效應在X、Y方向都可以觀察到，但是Y方向更為明顯。同時此次實驗在中子通量為2*107n/cm2/s的情況下，600s內(nèi)測量了幾個毫米空間尺度下350個灰階的變化。最后基于平場圖像分析得出系統(tǒng)的量子效率約為16%。

前言介紹

中子成像由于是一種無損檢測技術（NDT/NDE），被廣泛應用于汽車、航天航空以及學術界領域【1.2.3】。隨著中子技術的不斷成熟發(fā)展進步，中子成像正逐漸成為一種標準檢測觀察技術。

中子成像常和X射線成像比較，這是因為這兩種技術非常相似，區(qū)別只是在于電離輻射的產(chǎn)生和衰減的方式 【3，4】。X射線與材料的電子云相互作用過程中密度較大的材料或較重的元素衰減。中子由于不帶有電荷因此與原子核相互作用、同時與較輕的元素作用衰減得更快，故中子可以拍攝到更微觀的結(jié)構(gòu)。中子成像技術包括頻閃觀測、能量選擇成像和三維層析成像等。

中子成像作為一種無損檢測技術的主要局限性之一是高通量、高準直性的中子源很難獲得同時成本高。在過去的十多年里，人們一直在努力使工業(yè)界和學術界更容易使用中子成像技術【3，4，5，6】。其中包括加強相關研究項目、改進中子研究反應堆及其為相關領域提供更多資金，以及將這項技術引入新產(chǎn)業(yè)等。增強型Gadox-sCMOS中子相機（N-Cam）的研制增強了中子成像設備的可行性（高靈敏度、良好空間分辨率、更快速的積分時間）。同時它還為低通量反應堆提供了實驗可能，其可以保證低通量條件下的中子成像。中子散裂源（ISIS）以及和μ光源的成形譜儀（IMAT）裝置在此次實驗中給樣品提供了相對穩(wěn)定的中子束流【7】。本篇應用報告通過此套設備表征N-Cam在不同工作模式下的性能。

實驗與分析方法

探測器部分

N-Cam探測器部分包含一個厚度為20 µm、密度為7.3mg/cm2的Gadox閃爍體同時具有75mm靶面的有效探測區(qū)域。一個厚度為1mm的發(fā)黑鋁材料（防止低能量X射線干擾，類似X射線相機防止可見光的鈹窗）被放置在Gadox閃爍體前端2.7mm處，這樣就使待測對象與Gadox閃爍體之間的最小距離為3.7 mm。Gadox閃爍體直接應用于Photek 的MCP175像增強器的光纖輸入窗口前端。光纖輸入的數(shù)值孔徑（NA）為1.0，保證了閃爍體光的高收集效率。典型的中子和閃爍體材料互相作用會產(chǎn)生30kev的內(nèi)部電子沉積【8，9，10】，從而在閃爍體的平均電子路徑上會形成1200個光子的計數(shù)。大約32%的光信號會通過光纖窗口傳輸?shù)焦怆婈帢O。在此次實驗中閃爍體由于沒有反射覆蓋層以最大化優(yōu)化空間分辨率，造成了一部分的光損失。同時在光纖窗口的內(nèi)表面沉積的S20光陰極在閃爍體545nm發(fā)射峰位置提供了10%的量子效率。通過上述的計算得出，每吸收一個中子，在像增強器端會產(chǎn)生38個光電子；這些產(chǎn)生的光電子會在光陰極和微通道板之間的300um空間，通過一個200v的電壓加速。微通道板MCP是一個有數(shù)百萬個孔的薄玻璃板，可以用來放大電子信號【11】。N-Cam中子相機系統(tǒng)配置有一級MCP,該MCP具有直徑為12um間距為15um的陣列孔，孔徑比L:D是60：1。MCP的增益是由電壓控制的，對于大多數(shù)測試設置為+900v。通過MCP放大后的光電子通過5.5kv的加速電場轟擊P46熒光屏并輸出成像。在上述的電壓增益過程中，增強器的總光增益可達104。

Figure 1: 中子像增強器示意圖

高壓模塊HVPS以及門控控制器GCU集成在相機系統(tǒng)內(nèi)部，用雙層3.2mm厚柔性硼材料包裹以防止中子對電路的干擾。同時Photek使用GCU去控制50ns門控所需要的高壓脈沖電源。光電陰極的電壓開關是由光陰極與MCP之間的-50v偏置電壓控制，以防止無效光電子被放大。像增強器模塊可以在50ns尺度進行通斷操作。門控器GCU通過USB連接到數(shù)據(jù)采集器以及電腦上，這樣可以通過軟件控制MCP增益、光電陰極開關、門控持續(xù)時間以及延遲時序控制。利用中子源成形譜儀IMAT的信號和GCU模塊進行同步從而來實現(xiàn)不同的采集模式。

探測器記錄

N-Cam探測器的成像部分集成的是制冷型sCMOS相機，這臺相機具有2048x2048陣列格式，同時單像素尺寸為6.5μm。 此臺相機滿畫幅情況下最快成像速度是35fps，通過選擇區(qū)域ROI可以獲得更高的幀頻。在此次實驗中成像速度介于0.1fps和30fps之間。相機的滿阱容量為30000e-，讀出噪聲＜2e-，16bit數(shù)位動態(tài)范圍。同時暗噪聲以及其他噪聲可以通過-15℃的制冷進行大幅度削弱。相機通訊采用USB3.0,并有專業(yè)的控制軟件。

使用一個45度的反射鏡將像增強器的輸出光反射到相機口成像，90度的探測設計可以有效減少中子束流對相機的干擾。為了進一步防止粒子干擾以及輻射損傷，在相機外還放置了雙層3.2mm厚的柔性硼材料作為中子屏蔽防護層。75mm的探測成像區(qū)域通過Schneider-Kreuznach Componon f2.8/40mm 的鏡頭成像在上述的sCMOS相機上。同時相機還被安裝在一個可以調(diào)整聚焦位置的平移臺上。在實驗中我們還通過鏡頭延長環(huán)前移相機靶面獲得了33mm的視場，有效的分辨率在16.3um。

Figure 2: Diagram of N-Cam with labelled components

圖像處理

使用ImageJ軟件對圖像記錄并進行處理，對每一個數(shù)據(jù)組、平均暗場信號、平場校正信號以及樣品成像進行編譯。首先從一組參考的背景圖中獲取出平均背景信號；為了估算出平場信號，我們從每一個平場信號中減去平均背景，平場信號減去背景而后進行歸一化和平均處理這樣就可以得到平均平場校正信號；為了得到有效的目標信號，對所有的目標圖像進行背景扣除、歸一化處理而后除以平均平場信號。調(diào)整對象堆棧的大小可以改變相機曝光時間。通過將平均目標圖像除以平場信號，這樣可以使增益帶來的變化以及固定模式噪聲最小化【7，12】。最后處理得到的圖像再進一步分析。

空間分辨率

為了標定中子相機系統(tǒng)的空間分辨率，我們從the Paul Scherrer Institute (PSI)實驗中心取得一組分辨率標定圖。PSI提供了三組不斷減小的線對組圖(200-1000um;25-150um;5-25um)以及用于模式傳輸函數(shù)MTF分析的3°的傾角盒。通過以上工具對75mm以及30mm視場進行了極限分辨率分析。

對于75mm的視場，IMAT中子源的針孔設置為40mm，產(chǎn)生了5x106n/cm2/s的中子束流，同時中子源輸出孔徑比為245：1。實驗中在0.2fps情況下采集了60組平場信號以及60組背景信號，單大組耗時300s。分辨率標定圖緊貼像增強器的前方中央位置，也記錄了60組。為了評估X、Y軸的依賴，分辨率圖也被放置在了像增強器邊緣位置的3、6、9以及12點方向，每個位置被記錄了60幀。對于30mm視場實驗，記錄了30張平場以及分辨率圖，其中分辨率圖只是放在像增強器中心進行了測試。

處理圖像后，在第一和第二組線對（25-1000um）成像中重現(xiàn)出多條線陣的輪廓。同時對于不同線對的尺寸用方程二進行對比度計算。計算圖像的MTF模式傳遞函數(shù)邊緣效應，使用了ImageJ軟件插件，針對邊緣ROI區(qū)域進行定義分析。

靈敏度

基于閃爍體的中子成像系統(tǒng)的靈敏度很難評估。根據(jù)最初原理 ，閃爍屏的外部量子效率如下：

其中nc項代表是中子捕獲效率，ice項表示捕獲產(chǎn)生電子的效率。本次實驗使用的Gadox閃爍體，理論上有85%的中子可以被捕獲，同時內(nèi)部電子轉(zhuǎn)換效率85%，故得出理論效率為47%。其中捕獲概率是通過NIST數(shù)據(jù)中心捕獲效率以及散射計算得出；IMAT中子成形光譜的單幀白光能量譜圖參考了Kochelmann等人的相關文獻【7】。

在實際過程中，由于探測器的非均勻性、光收集效率、每一次內(nèi)部中子沉積引起電子不確定性以及閃爍體內(nèi)部光產(chǎn)額，所以DQE(探測器量子效率)明顯小于QE(理論效率)。其中中子作用閃爍體的能量沉積就可以在接近0keV（這個能量有大概率從閃爍體表面產(chǎn)生和逃逸）到大于100keV(在閃爍體內(nèi)移動)。此部分，我們嘗試通過分析暴露在全白光束的信號來估算量子效率DQE，同時也考慮到系統(tǒng)內(nèi)部的光路系統(tǒng)會引起相鄰像素之間重復計算。對于一個均勻的中子通量源，每個像素點的平均強度標記為ηFF，每個像素點的平均偏差σFF ,在一定區(qū)域的探測器中對捕獲的中子產(chǎn)生的內(nèi)轉(zhuǎn)換電強度進行泊松分布函數(shù)表征如下方程：

其中N(x)表征的是在x大小區(qū)域內(nèi)產(chǎn)生信號的中子數(shù)量，F(xiàn)N是中子通量，T是積分時間，r是相關的標度長度，可以認為是中子源產(chǎn)生的光的半徑。通過對得到的全白光平場數(shù)據(jù)針對不同Binning進行整合處理，而后用方程4擬合得到DQE和r值。在標準的像增強器設置情況下，1s積分時間，獲得了IMAT光源（80mm，中子通量2x107n/cm2/s）的平場數(shù)據(jù)。在連續(xù)成像模式下，采集了150幀平場的參考圖像，同時在之前又采集了5幀信號圖像。數(shù)據(jù)的矯正處理通過上述提到的方程1進行處理.

對比度波動分析CNR是衡量成像系統(tǒng)灰度分辨率的指標。為了測量這一個對比度指標，使用了各10個不同厚度的C101銅以及082-T651級別的鋁材料（從5mm-14mm，每1mm遞增），每種材料的截面都是5mm*5mm。成像使用IMAT光源針孔束結(jié)構(gòu)，其中針孔為80mm，孔徑比125：1，產(chǎn)生中子通量為2x107n/cm2/s。相機設置積分時間為1s，達到了50%的像素滿阱強度。像增強模塊MCP電壓設置為900v，熒光屏電壓為5500v。在模型拍照前后拍攝兩組平場圖像，每組積分1s。同時在IMAT光源關閉且衰減器放置情況下獲得同樣曝光時間的背景圖像。同時獲得模型材料的圖像。

所有的平場圖像、背景圖像以及模型圖像都對像素靶面39872個像素區(qū)域進行了歸一化處理，以校正曝光期間的光束強度和增益變化。從得到的圖像進行分析，20個不同對比度材料的截面的邊緣附近不均勻，假設是由于邊緣效應，包括光束發(fā)散、中子散射和沿著每個截面之間的尖銳邊界的衍射而引起的。對于上述提到的5 mm正方形部分的模型材料，在其中心附近選擇32 x 64像素進行ROI分析，并計算像素值的平均值和標準偏差。相鄰兩段（A和B）之間的CNR定義為

2結(jié)果分析

空間分辨率

根據(jù)圖三得到的是300s積分情況下，75mm視場與30mm視場的中心成像區(qū)域的處理后的圖像，利用ImageJ測得的線輪廓，我們估算了30mm和75mm視場下，N-cam系統(tǒng)的極限分辨率（5%）分別是12.5lp/mm以及10lp/mm 。對于有效像素尺寸為38um的75mm視場，同時進行采樣處理（Nyquist），得到的極限分辨率為13lp/mm，和MCP175的16lp/mm分辨率結(jié)合分析，10lp/mm屬于合理的數(shù)據(jù)。對于有效像素尺寸16um的30mm視場，采樣處理后極限是30lp/mm，與16lp/mm像增強器結(jié)合得到12.5lp/mm與理論相差較大，當然這個數(shù)據(jù)還有其他我們沒有考慮到的影響因素，比如我們需要更詳細的測量來確定熒光屏對分辨率的影響。

(a) 中心分辨率圖像 (75mm FOV)       (b) 中心分辨率 (75mm 視場)線寬25-150µm

Figure 3: 中子成像得到的PSI分辨率圖樣

Figure 4: 不同尺寸線對的平均測量對比圖

如圖5所示，將曝光時間從300秒減少到30秒不會影響測量的對比度或計算的MTF曲線。將曝光時間減少到30秒以下，測量對比度和計算的MTF曲線的趨勢相同，但引入了更多的噪聲，如圖6所示。然而，MTF曲線小于5%的點保持在9.5–10.5 lp/mm之間 。

Figure 5: 不同尺寸線寬 、不同積分時間75mm視場下對比度     Figure 6: 75mm不同分辨率條件下75mm視場 MTF曲線（5、10、30、300s曝光時間）

DQE探測器量子效率

對成像區(qū)域的24mm*24mm進行平場成像區(qū)域進行校正, 曝光時間T=5s,中子流量FN = 2x107n/cm2/s , 并對此區(qū)域內(nèi)使用泊松分布平場進行了矯正。 通過方程4進行擬合，而后估算得到的DQE=16.1%,中子產(chǎn)生的光半徑r=81um。低于理論結(jié)果的47%，這是由于閃爍體表面的低能量損耗、探測器不均勻性以及儀器其他噪聲的綜合影響。r=81um表示距離像素邊緣一定距離的中子互相作用將對相機信號產(chǎn)生影響。考慮到基于熒光屏的光散射，這個數(shù)值并不奇怪，并且和上述的MTF測量值一致，在2.5lp/mm情況下，50%對比度，200um線寬，大致為2r。未來將研究Gadox閃爍體的替代材料，以確定它們對r的影響和系統(tǒng)的極限分辨率。

R2 = 0.999625 Figure 9: Fit to standard deviation as a function of pixel size using Equation 4, parameterized with QEeff and r.

對比度波動

圖10給出了600s積分情況下對比度分析成像組的平場對比度圖像。其中如像圖10b所示，有四個長方形的模塊以及一個鋁制支架、兩個鈦材料螺絲組成。每一個長方形的尺寸為5mm寬、25mm長，銅材料（左側(cè)10a）和鋁材料（右側(cè)10c）兩大組中分別是以1mm遞增的5-9mm以及10-14mm的厚度。

10a) Copper step-wedge from 5 to 14mm in 1 mm steps    10b) Contrast Phantom    10c) Aluminum Step-wedge from 5 to14 mm in 1 mm steps

Photek Limited

26 Castleham Road, St Leonards on Sea, East Sussex, TN38 9NS, United Kingdom

Revision 1 Last Updated: 03/04/202

圖11通過方程5分析測量了不同厚度情況下的透過噪聲曲線。 通過圖11可以得出1.0cm與1.1cm厚度鋁材料之間得CNR值是2.6。相同厚度的銅材料CNR是11.4。假定CNR=1的情況是最小可以檢測的分辨率對比度，我們發(fā)現(xiàn)對于鋁的可探測的最小對比度是0.26%，銅是0.3%，表明圖像中的灰階值應該是350。這個值同時也取決于中子束特性、樣品的成分以及幾何結(jié)構(gòu)。

Figure 11a. Transmission through the aluminum step-wedge with error bars representing 1 σ noise.

Figure 11a. Transmission through the copper step-wedge with error bars representing 1 σ noise.

結(jié)論

N-cam在75mm視場下分辨率為10lp/mm（300s曝光測試時間）。曝光時間減小到30s以上對分辨率幾乎沒有影響。將曝光時間減小到30s以下時，分辨率也在10lp/mm，但會影響到分辨率。因此，采集300s積分情況下的平場、暗場圖像，采集360度旋轉(zhuǎn)的成像，一次斷層掃描成像最快需要通過5次拍攝歷時40分鐘或者30次拍攝耗時3個小時。N-cam在測試過程中顯示了x,y方向的依賴性，其中y方向?qū)D像的影響大于x方向，但是其中的原因需要進一步的研究確定。同時系統(tǒng)的量子效率得出是16.1%，閃爍體的厚度會影響量子效率，同時也會影響小視場情況下的分辨率，這一點也有待研究。然而，對于75mm的大視場的影響可以忽略不計,與理論值近似，這主要是因為大視場情況下的有效像素點比較大所以影響較小。
References