基于matlab環境下的圖像重建及吸收率分析

李彪，張淼

摘 要：針對重建未知介質的幾何形狀及吸收率，綜合運用濾波反投影算法進行圖像重建，首先先引入中心切片定理，采集不同視角下的投影并求解其傅里葉變換，匯集成圖像的傅里葉變換再經傅里葉反變換，其次利用matlab iradon函數重建得未知介質的幾何形狀，再次利用坐標變化和標定參數得出未知介質在正方形托盤中的位置，最后通過matlab rgb2gra函數得到重建圖形的灰度值矩陣，由灰度值與吸收率的關系待測未知介質吸收率。

關鍵詞：圖像重建；濾波反投影法；matlab灰度值；吸收率

中圖分類號：TP391.41 文獻標志碼：A 文章編號：2095-2945（2018）06-0193-02

Abstract： In view of the geometric shape and absorptivity of the unknown medium， the filter back-projection algorithm is used to reconstruct the image. Firstly， the central slice theorem is introduced， and the projection in different angle of view is collected and the Fourier transform is solved. The Fourier transform which gathers the image is then transformed by the inverse Fourier transform， and then the geometric shape of the unknown medium is reconstructed by using the matlab iradon function. The position of unknown media in square pallet is obtained by coordinate change and calibration parameters. Finally， the gray value matrix of reconstructed graph is obtained by Matlab rgb2gra function. According to the relationship between gray value and absorptivity， the absorptivity of unknown medium is determined.

Keywords： image reconstruction； filtering inverse projection method； MATLAB gray value； absorption rate

引言

CT可以在不破壞樣品的情況下，利用樣品對射線能量的吸收特性對生物組織和工程材料的樣品進行斷層成像，由此獲取樣品內部的結構信息[1]。一種典型的二維CT系統[2]如圖1所示，相互平行的X射線垂直射入512個等距排列的探測器內。然后X射線沿旋轉中心逆時針旋轉180次得到180組的接受信息數據。

CT系統的參數即為CT系統旋轉中心在正方形托盤中的位置、探測器單元之間的距離以及該CT系統使用的X射線的180個方向。我們在已知標定參數的條件下，利用上述CT系統得到的某未知介質的接收信息512*180組數據，就可得出未知介質的幾何形狀及其內部吸收率，從而了解內部結構信息，這對現實中我們探索未知介質內部信息擁有很重要的意義。

1 未知介質的幾何形狀求解

1.1 濾波反投影重建模型

待建圖像為f（x，y），它的傅里葉變換為F（？棕1，？棕2）=

■（？籽，？茲），由中心切片定理有，■（？籽，？茲）可通過f（x，y）在不同視角？準下的投影p？準（xr）的傅里葉變換求得[4]。

待建圖像：

將其寫成空域變量為xr的傅里葉反變換式得到：

1.2 模型的算法

Step 1：利用matlab的iradon函數重建出介質的幾何形狀。

Step 2：利用第一問標定參數旋轉中心在托盤中的位置和坐標軸變換得到待測介質在正方形托盤的位置。

Step 3：將重建出來的圖形轉換為灰度圖。

其中白色正方形線段為100*100（mm）的托盤的大小，里面的灰色部分則為重建出來的未知介質形狀的灰度圖，圓圈表示托盤的中心位置。

2 內部吸收率求解

2.1 灰度值求解原理及算法

由灰度值[5]概念可知，灰度值范圍一般從0到255，白色為255，黑色為0，其間部分則為黑白部分，由得出的灰度圖圖2可知，亮度越亮，灰度值越大，灰度值為0的部分則表示吸收率最小，即為黑色，亮度越暗，灰度值越小，灰度值為1的部分則表示吸收率最大，即為白色，于是我們得出結論：灰度值與吸收率成正比。

利用matlab中imread函數讀取圖2白色方形區域內（托盤內）的灰度圖，利用gb2gray函數將其轉換為相應的灰度值矩陣，再利用相應的imresize函數將灰度值矩陣轉換為符合要求的大小為256*256的灰度值矩陣，可得出矩陣大小為256*256中各個位置的灰度值，由于灰度值與吸收率成正比，再通過0-1規劃得出吸收率矩陣。

2.2 求取未知介質內部吸收率

在已知灰度值與吸收率的關系后，我們得到了未知介質內部的所有灰度值，接下來我們求取1個特殊位置的吸收率。

由濾波反投影重建所得到的灰度圖的坐標軸與所給位置坐標的坐標軸存在90°的夾角，所以需將所給坐標點坐標變換得到新的坐標。

因為所給點的坐標為0.5的倍數，于是我們將大小為100*100的正方形托盤區域分為以0.5為間距，大小為200*200的點集，用imresize函數將256*256的灰度值矩陣轉換為200*200的灰度值矩陣，再通過0-1規劃得出吸收率矩陣。最后通過表1新坐標從轉換后的200*200的灰度值矩陣讀取該點的吸收率并將轉換后的坐標顯示出來得出結果見表1

3 結束語

在僅已知未知介質的能量接收信息，從而來確定其幾何形狀、介質的具體位置以及其吸收率時，在并沒有明確指出未知介質是否勻質，不能簡單的將其當做均勻物質處理，需將其吸收強度看作位置平面的函數。X射線平行照射時，實質上是將待測介質向某一軸投影，通過逐漸的變化角度，求得在不同方向上的投影關系，這種由投影確定原來圖形的方法是典型的圖像重建問題，由此引入中心切片定理和濾波反投影進行圖像重建，利用重建后的灰度圖得出灰度值，從而建立灰度值與吸收率的關系。

參考文獻：

[1]張旭良，童隆正，湯韋華，等.X線機械斷層的仿真系統一用

MATLAB虛擬X線斷層（二）[J].醫療裝備，2000（12）：1-2.

[2]馬繼明，張建奇，宋顧周，等.全變分約束迭代濾波反投影CT重建[J].光學學報，2015，35（02）：392-398.

[3]伍偉文，全超，劉豐林.相對平行直線掃描CT濾波反投影圖像重建[J].光學學報，2016，36（09）：165-175.

[4]朱 ，王富東.利用MATLAB實現二維圖像傅立葉變換算法[J].計算機應用與軟件，2006（12）：141-142.

[5]張曉克，許建剛.超分辨率圖像重建技術研究[J].科技創新與應用，2014（32）：68-69.endprint