【文章內(nèi)容簡(jiǎn)介】 檢測(cè)物體之間的若干相互作用過程的產(chǎn)物。例如光電吸收過程和康普頓散射過程。這些過程中的每一種都有它自己的發(fā)生幾率。幾率也是輻射能量的函數(shù) ,因?yàn)閺?X 線管產(chǎn)生的 X 線由全部能譜組成， 具有“線衰減系數(shù) u ”的某種組織的衰減性質(zhì)是一第二章 CT 圖像重建的基本原理 7 個(gè)復(fù)雜的函數(shù)，隨著輻射情況的變化，它可以表現(xiàn)出不同的數(shù)值。 對(duì)于能量為 E 的單能射線，假設(shè)其投射的物體是均勻且為同一介質(zhì)，被投射物體的線性衰減系數(shù)為 u ，厚度為 x ，則在入射線強(qiáng)度為 0I ，出射 X 射線強(qiáng)度為I 時(shí)有 uxeII ?? 0 即 )/ln( 0 IIux? （ ） 這 就是著名的 LambertBeer 定律。 實(shí)際被投射的物體是非均勻的，如果在 X 射線束通過的路徑 l 上，介質(zhì)是非均勻的，則可以將 X 射線穿過的介質(zhì)沿掃描路徑 l 人為劃分為大小相等的 n 個(gè)小方塊，每個(gè)小塊的厚 度為 d ，塊內(nèi)視為均勻介質(zhì)。這樣，每個(gè)小塊有一個(gè)均勻的線性衰減系數(shù)，如圖所示，設(shè)各個(gè)方塊的線性衰減系數(shù)分別為 nuuuu ,..., 321 。 X射線通過每一個(gè)方塊的衰減為 dueII 101 ?? ，射線通過第二個(gè)方塊的衰減為dueII 212 ?? ，通過第 n 個(gè)方塊的衰減為 dunn neIII ??? 。將上述各式代入 整理，有： 圖 X射線通過非均勻介質(zhì)時(shí)的體素劃分 duuuu neII )...(0 321 ?????? （ ） 即 dpuuuu n /...321 ????? （ ） 式中， )ln( 0IIp? 是通過射線強(qiáng)度檢測(cè)得到的檢測(cè)值，是一個(gè)已知量，稱為投影。這樣，上式中，等號(hào)右邊可視為一個(gè)常數(shù)，等待求解的是等式左邊的每一個(gè)方塊的衰減系數(shù)。由于未知量太多，這個(gè)方程沒有唯一解?？紤]到圍繞受檢體的一次掃描可以得到一個(gè)有關(guān)衰減系數(shù)的方程 ，如果將整個(gè)感興趣斷層都劃分為大小相等的小方塊，每個(gè)小方塊由于一個(gè)固定的衰減系數(shù)，這樣，上述單束射線環(huán)繞掃描得到的每個(gè)穿透強(qiáng)度值都可以形成一個(gè)關(guān)于衰減系數(shù)的方程。只要掃描線不重復(fù)且足夠多，這些衰減系數(shù)的確定值肯定可以通過解線性方程組的方法解第二章 CT 圖像重建的基本原理 8 出唯一實(shí)數(shù)解。把每個(gè)方塊的衰減系數(shù)值用灰度表示，就可以重建出以衰減系數(shù)為特征的斷層圖像。這種圖像重建方法稱為方程組法。 值得說明的是，當(dāng)上圖中 0?d 時(shí)，投影 p 就是 dlluIIp ??? )()ln ( 0 （ ） 式中，衰減系數(shù) )(lu 是隨路徑 l 連續(xù)變化的函數(shù)。上式說明，入射 X 射線強(qiáng)度與出射 X 射線強(qiáng)度之比經(jīng)對(duì)數(shù)運(yùn)算后，表示沿射線路徑上衰減系數(shù)的線積分，而投影 p 與射線穿越介質(zhì)的路徑長(zhǎng)度成正比。 然而，一般二維斷層圖像至少應(yīng)劃分出 160 160個(gè)像素，如要用方程組法重建這個(gè)斷層圖像，就要 25600個(gè)獨(dú)立的方程聯(lián)立求解 ，這是一個(gè)相當(dāng)費(fèi)時(shí)的任務(wù)。盡管 1967年世界首臺(tái) CT 試驗(yàn)機(jī)采用的就是這個(gè)圖像重建方法，但當(dāng)斷層被劃為更小更多的體素時(shí)，即使采用今天最先進(jìn)的計(jì)算機(jī)技術(shù)，解一組這樣的方程組也不是件容易的事。所以，這種方法在后來的實(shí)際應(yīng)用中已不再使用，而是采用下一章將要介紹的反投影解析法。 第三章 CT 圖像解析重建算法 9 第三章 CT 圖像解析重建算法 通過上一章已經(jīng)知道 CT 圖像重建的基本原理，但由于使用最簡(jiǎn)單直接的方程組法不易求得斷層的各個(gè)像素的衰減系數(shù)。這一章介紹現(xiàn)在 CT 最常用的圖像重建算法 —— 反投影解析法。主要介紹了由投影重建圖像的系統(tǒng)模型，引進(jìn)Radon 變換和傅立葉變換，它們對(duì)應(yīng)于 X 射線 CT 的二維重建問題；且給出了Radon 變換在圖像重建中的具體形式，即截面函數(shù)沿著直線的線積分等于其 Radon 變換，而圖像重建過程即是將截面函數(shù)沿許多不同角度下直線的線積分所產(chǎn)生的投影數(shù)據(jù)進(jìn)行逆 Radon 變換從而得到截面函數(shù)。另外由逆 Radon 變換的推導(dǎo)給出了兩種常用的圖像解析重建法：濾波反投影法和卷積反投影法。 反投影重建法（ BP） 反投影是一個(gè)求和的過程，它把所有角度上的數(shù)據(jù)都加在一起。它把投影域中的數(shù)據(jù)沿著原來的投影路徑“涂抹”回去，但不改變數(shù)據(jù)的值 ]3[ 。 直接反投影重建是指，斷層平面中某像素的線性衰減系數(shù)可由平面內(nèi)所有經(jīng)過該像素射線的投影和得出，而整幅圖像重建可以看做為所有方向上掃描投影的累加。通俗地講，所謂直接反投影，就是在對(duì)某體面一個(gè)方向的掃描完成后，以得到的投影為灰度值沿著掃描路徑經(jīng)過的體素回抹到體素對(duì)應(yīng)的像素上。改變方向后的多次掃描形成多次回抹，同一像素上多次回抹的灰度累加即完成圖像重建。 圖 斷層的像素和射線模型 第三章 CT 圖像解析重建算法 10 圖 反投影重建示例 如圖 所示，其中（ a）為原圖像像素，（ b）為反投影重建后圖像，（ c）為（ b）中像素值除以投影線數(shù)后得到的平均值。有上述例子可以看出原圖中像素值為零的像素，經(jīng)反投影重建后不再為零。即有偽跡。 直接反投影圖像產(chǎn)生星狀偽跡的原因在于，該方法是把取自有限物體空間的投影均勻地回抹（反投影）到了射線所及的無限空間的各個(gè)像素上，包括原來像素值為 0的點(diǎn)。 圖 反投影重建的等效系統(tǒng) 設(shè) ),( yxu 為處于 0,0 ?? yx 處的一個(gè)點(diǎn)源，即二維斷面 ),( yx 上的一個(gè)沖擊函數(shù) ),( yx? ，這時(shí)的輸出 ),( yxf 就是系統(tǒng)的沖擊響應(yīng) ),( yxh 。 已知 221),( yxyxh ?? ? （ ） 可見，反投影重建過程后得到的圖像不是原來的圖像 ),( yx? ，系統(tǒng)存在星狀偽跡。 根據(jù)信號(hào)與系統(tǒng)理論，對(duì)于圖 ： ),(**),(),( yxhyxuyxf ? （ ） 式中，兩個(gè)星號(hào)表示二維卷積 ]4[ 。 第三章 CT 圖像解析重建算法 11 如果想要去除星狀偽跡，可在系統(tǒng)的輸出端再加上一個(gè)濾波器，設(shè)其時(shí)域函數(shù)為 ),( yxq ，傳遞函數(shù)為 ),( 21 ??Q ，為了使該加了濾波器的系統(tǒng)輸出的圖像等于系統(tǒng)輸入的原圖 ),( yxu ，也就是要求 ),(),(**),( yxyxhyxq ?? （ ） 即 ),(1**),(22 yxyxyxq ?? ?? （ ） 對(duì)該式取二維傅立葉變換，得到： 11),(2221 ??? yxQ ??? （ ） 即 222121 ),( ????? ??Q （ ） 這是一個(gè)二維濾波器。即使加上某種近似，這種濾波器也不容易實(shí)現(xiàn)。如果21 ??和 的值取無窮，則是一個(gè)不可能實(shí)現(xiàn)的理想濾波器。為此提出了濾波反投影圖像重建法。 濾波反投影圖像重建法 （ FBP） 為了消除直接反投影法產(chǎn)生的星狀偽跡，提出了濾波反投影重建的方法，這種方法是在直接反投影重建方法的基礎(chǔ)上增加一個(gè)濾波器。增加濾波器的思路有兩種，如圖 所示。 第三章 CT 圖像解析重建算法 12 圖 兩種不同的去除星狀偽跡反投影圖像重建思路 所謂濾波反投影重建法是把獲得的投影函數(shù)先做卷積處理，即人為設(shè)計(jì)一種一維濾波器，用它對(duì)所得的投影函數(shù)進(jìn)行改造，然后把這些改造過的投影函數(shù)進(jìn)行反投影和累加等處理，就可以達(dá)到在一定程度上消除星狀偽跡的目的。 在介紹濾波反投影重建法之前先補(bǔ)充濾波反投影重建法要用到的數(shù)學(xué)知識(shí)和投影定理 —— 中心 切片定理。 Radon 變換與傅立葉變換 為了有效和快速地對(duì)圖像進(jìn)行處理，常需要將原定義在圖像空間的圖像以某種形式轉(zhuǎn)換到另外一些空間，并利用在這些空間的特有性質(zhì)方便地進(jìn)行一定的加工，最后再轉(zhuǎn)換回圖像空間以得到所需的效果 ]7[ 。通常 正變換：圖像空間到其他空間 反變換：其他空間到圖像空間 Radon 變換 Radon 變換 1917 年由奧地利數(shù)學(xué)家 Radon 提出，是 CT 圖像重建思想的數(shù)學(xué)表達(dá)。 Radon 從理論上證明了某物理量的二維分布函數(shù)可由該函數(shù)在 其定義域內(nèi)的所有方向上的線積分完全確定。 Radon 變換的意義在于，只要知道了一個(gè)未知二維分布函數(shù)的所有方向上的線積分，那么就能夠求得該二維分布函數(shù)。所謂CT 斷層的圖像重建，就是求取能夠反應(yīng)斷層內(nèi)部結(jié)構(gòu)和組成的某物理量的二維分布。斷層掃描獲得的投影實(shí)際上可視為被測(cè)物體斷層內(nèi)部結(jié)構(gòu)和組成的不同方向上的線積分，所以 Radon 變換的正變換和逆變換正好對(duì)應(yīng)了 CT 技術(shù)的射線投影和圖像二維分布函數(shù)的重建過程。 第三章 CT 圖像解析重建算法 13 對(duì) ),( yxf 的 Radon 變換 ),( ?pRf 定義為沿由 p 和 ? 定義的直線 l 的線積分。若設(shè) ),( yxf 表示物體的一個(gè)二維斷層分布，通過 Radon 變換的到 f 沿不同方向直線的線積分相當(dāng)于得到物體不同方向的投影。 圖 定義 Radon 變換的坐標(biāo)系統(tǒng) 上述線積分可寫為： ????? dlyxfpR f ),(),( ? （ ） 如果借助 Delta 函數(shù)，上述線積分還可寫為： ? ???? ??? ??? dx dyyxpyxfpR f )s i nc os(),(),( ???? （ ） 上式通常稱為 ),( yxf 的 Radon 變換。 由上述可知， ),( yxf 關(guān)于某直線的 Radon 變換就是 ),( yxf 沿該直線上的一維投影。 實(shí)際上當(dāng)將 ? 也視作可變參數(shù)時(shí)， ),( ?pRf 就是一個(gè)二元函數(shù)，此時(shí) Radon變換是空間域 ),( yx 到變換域 ),( ?p 的映射。域 ),( ?p 中的每個(gè)點(diǎn)對(duì)應(yīng)于空間域),( yx 中的一條直線，這里 ),( ?p 實(shí)際上就是《圖像分析》中直線的 Hough 變換 ]10[ 。 傅立葉變換的定義 ???? ???? dedeftf ii? ????????? ])([21)( （ ） 其中 第三章 CT 圖像解析重建算法 14 dtetfF ti?? ?????? )()( （實(shí)自變量的復(fù)值函數(shù)） （ ） 稱 為 )(tf 的 Fourier 變換，記為 )]([ tfF 。 ??? ? deF ti)(21 ???? 稱為 )(?F 的 Fourier 逆變換，記為 )]([1 ?fF? 。 中心切片定理 中心切片定理的含義是平行投影的一維傅立葉變換等同于原始物體的二維傅立葉變換的一個(gè)切片。將已知投影數(shù)據(jù)通過一個(gè)簡(jiǎn)單的二維傅立葉反變換可以得到物體截面的一個(gè)評(píng)估。 圖 中心切片定理示意圖 通過上圖可知：某斷層（或它對(duì)應(yīng)的圖像） ),( yxf 在視角為 ? 時(shí)得到的平行投影（函數(shù)）的一維傅立葉變換，等于 ),( yxf 二維傅立葉變換 ),( 21 ??F 過原點(diǎn)的一個(gè)垂直切片，且切片與軸 1? 相交成 ? 角 ]9[ 。 數(shù)學(xué)表達(dá)如下： 首先定義代表截面的函數(shù)的二維傅里葉變換： dx dyeyxfvuF vyuxj? ?????????? )(2),(),( ? （ ） 類似的定義某個(gè)角度 ? 下的一條投影 )(tP? 的傅里葉變換： 第三章 CT 圖像解析重建算法 15 dtetPS tj ???? ? 2)()( ?????? （ ） 舉個(gè)最簡(jiǎn)單的例子：沿著 ? =0 的直線方向，物體投影的傅里葉變換等于頻域里面 ? = 0的情形 dx dyeyxfuF uxj? ????????? ?2),()0,( （ ） 因?yàn)橄辔灰蜃硬辉僖蕾囉?y ，在此可以將積分分成兩部分， dxedyyxfuF uxj? ????????? ?2]),([)0,( （ ） 從平行投影的定義看，可以將上式括號(hào)中的部分看成是沿著常量 x 的積分 dyyxfxP ?????? ),()(0? （ ） 將其代入 ()得到 dxexPuF uxj??????? ?? 20 )()0,( （ ） 等式右邊表示的是投影 0??P 的一維傅里葉變換；因此可以得到垂直方向投影和投影函數(shù)的二維變換的關(guān)系式： )()0,( 0 uSuF ?? ? （ ） 這是傅里葉切片定理最簡(jiǎn)單的形式。很明顯物體與坐標(biāo)系統(tǒng)間的位置方向是獨(dú)立的。例如，在圖 中，如果 ),(st ,坐標(biāo)系統(tǒng)為 ),( yx ,坐標(biāo)系統(tǒng)旋轉(zhuǎn) ? ，式()定義的投影的傅里葉變換等于物體沿著一旋轉(zhuǎn)了 ? 角的直線的二維傅里葉變換。從而引出了成像系統(tǒng)中著名的中心（傅里葉）切片定理 ]2[ 。 即某圖像 ),( yxf 在視角 ? 時(shí)的平行投影的 傅里葉變換給出了圖像二維傅里葉變換 ),( vuF 在與 u 軸夾角為 ? 的一個(gè)切片，此切片通過原點(diǎn)。換句話說，的 )(tP?傅里葉變換在數(shù)值上對(duì)應(yīng)于 ),( vuF 沿著圖中所示直線 BB 的值。 通過將 ),( yx 坐標(biāo)系統(tǒng)旋轉(zhuǎn)得到 ),(st 坐標(biāo)系統(tǒng)后，傅里葉切片定理可以得到更廣泛的應(yīng)用。在 ),(st 系統(tǒng)中沿著常量為 t 的直線進(jìn)行