CT6圖像重建2課件_第1頁(yè)
CT6圖像重建2課件_第2頁(yè)
CT6圖像重建2課件_第3頁(yè)
CT6圖像重建2課件_第4頁(yè)
CT6圖像重建2課件_第5頁(yè)
已閱讀5頁(yè),還剩80頁(yè)未讀, 繼續(xù)免費(fèi)閱讀

下載本文檔

版權(quán)說明:本文檔由用戶提供并上傳,收益歸屬內(nèi)容提供方,若內(nèi)容存在侵權(quán),請(qǐng)進(jìn)行舉報(bào)或認(rèn)領(lǐng)

文檔簡(jiǎn)介

CT圖像重建2

(CTImageReconstruction)學(xué)習(xí)內(nèi)容(Learningobjects)迭代法(Iterativealgorithm)

傅立葉變換法(Fouriertransformation)

反投影法(Backprojection)濾波反投影法(Filterbackprojection)知識(shí)回顧(Review)Radon變換投影與弦圖

采樣幾何形狀Shepp-Logan體模

圖像重建投影P(t)f(x,y)y理想圖像投影與弦圖弦圖理想圖像投影投影采樣幾何形狀Shepp-Logan體模r的單位是角度,μ為負(fù)值時(shí)表示削弱原有橢圓的衰減系數(shù)

S-L體模是CT圖像重建領(lǐng)域用于仿真計(jì)算的經(jīng)典頭部模型,于1974年由L.A.Shepp和B.F.Logan首次提出,可生成2D或者3D的標(biāo)準(zhǔn)投影數(shù)據(jù)。S-L體模通過橢圓來表征不同的形狀,不同的灰度用來模擬不同組織的衰減系數(shù),例如最外層的橢圓模擬頭骨,內(nèi)部的兩個(gè)小橢圓模擬大腦內(nèi)部特征或者腫瘤。直接矩陣求解法2x2matrixA11A12Po(0)=3

2A21A22Po(1)=7

4P9o(0)P9o(1)=4=6A11+A12 =3 A21+ A22 =7A11+ A21 =4A12 +A22 =61100 A11 30011 A12 71010 A21 41001 A22 6

F*A=PA=PF^(-1)=xywxwysg(s,1,2,…,8)1D-FouriertransformF[g(s)]μ(x,y)中心切片定理(CentralSliceTheorem)學(xué)習(xí)內(nèi)容(Learningobjects)迭代法(Iterativealgorithm)

傅立葉變換法(Fouriertransformation)反投影法(Backprojection)

濾波反投影法(Filterbackprojection)實(shí)例:反投影算法(Backprojection)弦圖反投影圖像反投影濾波濾波后圖像弦圖實(shí)例:濾波反投影算法(Filteredbackprojection)濾波反投影圖像理想圖像反投影圖像Projection反投影法(Backprojection)P(t)f(x,y)y在一個(gè)“空白”圖像里,把投影值從各個(gè)角度重新加進(jìn)去反投影法(Backprojection)由投影重建圖像的算法很多,而反投影法是其中最簡(jiǎn)單、最粗略,也是最基本的算法原理:“斷層平面中某一點(diǎn)的密度值可看作這一平面內(nèi)所有經(jīng)過該點(diǎn)的射線投影之和(的平均值)”第一步(firstguess)=0=2

=1=31(0+1)5(2+3)1(0+1)5(2+3)把90角度的投影值加進(jìn)空白圖像實(shí)例第二步(secondguess)021303331(0+1)5(2+3)1(0+1)5(2+3)+1(0+1)8(5+3)4(1+3)8(5+3)第四步(fourthguess)021332133(2+1)10(2+8)8(4+4)12(4+8)+6(3+3)12(2+10)9(1+8)15(3+12)0/36/33/39/36-612-69-615-6063902136129150213所有反投影的和算法舉例根據(jù)反投影算法x1=p5=5 x6=p2+p3+p5=18 …平均化處理,除以投影線數(shù)目

xi=xi/6000005200100000056237181271108136250.8310.3300.51.16321.160.061.661.330.160.510.330.83反投影重建后原像素值再除以投影線數(shù),平均化斷層平面中某一點(diǎn)的密度值可看作這一平面內(nèi)所有經(jīng)過該點(diǎn)的射線投影之和的平均值123456算法舉例反投影重建后,原來為0的點(diǎn)不再為0,形成偽影00000520010000000.8310.330.51.16321.160.061.661.330.160.510.330.83原像素值除以投影線數(shù)的平均化星狀偽影考慮孤立點(diǎn)源反投影重建中心點(diǎn)A經(jīng)n條投影線投影后,投影值均為1:

p1=p2=...=pn=1因此重建后而其他點(diǎn)均為1/n

這類偽跡成為星狀偽影1/n1/n1/n1/n11/n1/n1/n1/n000010000反投影重建的本質(zhì)是把取自有限物體空間的射線投影均勻地回抹(反投影)到射線所及的無(wú)限空間的各點(diǎn)之上,包括原先像素值為零的點(diǎn)反投影法的缺陷圖像產(chǎn)生模糊反投影法的改進(jìn)后處理方法:用2D濾波函數(shù)對(duì)反投影法所得圖像進(jìn)行處理,可實(shí)現(xiàn)圖像質(zhì)量的提升,但不能改變換準(zhǔn)確重建反投影法的改進(jìn)前處理方法:在反投影之前先濾波,即先對(duì)1D投影進(jìn)行濾波,再進(jìn)行反投影重建(理論基礎(chǔ):在線性系統(tǒng)中濾波算子可以交換順序)濾波反投影法學(xué)習(xí)內(nèi)容(Learningobjects)迭代法(Iterativealgorithm)

傅立葉變換法(Fouriertransformation)反投影法(Backprojection)

濾波反投影法(Filterbackprojection)濾波反投影法(FBP)的思想反投影法和濾波反投影法的區(qū)別:

濾波運(yùn)算或卷積運(yùn)算的引入在時(shí)域中相卷積相當(dāng)于在頻域中相乘,故濾波反投影法亦被稱為卷積方法(convolutionmethod)FBP的圖形解釋FBP公式推導(dǎo)目標(biāo)函數(shù)

f(x,y)可由傅立葉函數(shù)F(u,v)的逆變換獲得,即雅可比行列式FBP公式推導(dǎo)頻域中的笛卡爾坐標(biāo)與極坐標(biāo)的關(guān)系為:

FBP公式推導(dǎo)Let:F(cos,sin)=P(,)

FBP公式推導(dǎo)FBP公式推導(dǎo)FBP公式推導(dǎo)P(,)表示對(duì)應(yīng)于角度的單位投影的傅立葉變換;里層的積分是P(,)||的逆傅立葉變換,記為g(t,),在空間域,它表示單位投影被一頻域響應(yīng)為||的函數(shù)做濾波運(yùn)算,故稱之為濾波反投影1DFouriertransforminverse1DFouriertransformbackprojectionforallanglesfilterFBP公式推導(dǎo)FBP的計(jì)算機(jī)實(shí)現(xiàn)濾波(Filter)或卷積(Convolution)影像濾波,即為與一個(gè)適當(dāng)?shù)暮耍╧ernel)作卷積(convolution)運(yùn)算–3x3averagingkernel…–7x7Gaussiankernel…卷積為一個(gè)積分運(yùn)算,通常用符號(hào)表示卷積對(duì)兩個(gè)函數(shù)進(jìn)行運(yùn)算獲得第三個(gè)函數(shù)的過程,其中一個(gè)輸入函數(shù)被稱之為卷積核卷積核通常為3x3的小矩陣卷積核(Kernel)卷積核(Kernel)卷積核(Kernel)卷積核(Kernel)卷積核(Kernel)以下依此類推…….卷積核(Kernel)濾過器的設(shè)計(jì)和選擇在空域中,與P(,)對(duì)應(yīng)的函數(shù)是測(cè)量得到的平行投影p(t,),濾波器||與對(duì)應(yīng)的沖擊響應(yīng)(t)之間是一個(gè)傅立葉變換對(duì):將t=0代入上式計(jì)算得到(0),即曲線||以下的面積。當(dāng)時(shí),(0),所以上式是無(wú)法直接計(jì)算的,必須另想它法,引入限帶函數(shù)(band-limitingfunction)濾波器是個(gè)無(wú)限頻帶的濾波函數(shù),由于即積分是發(fā)散的,根據(jù)佩利--維納準(zhǔn)則,這一理想濾波器是不可實(shí)現(xiàn)的。實(shí)際數(shù)值計(jì)算通常采用加窗的濾波函數(shù)。運(yùn)用不同的窗函數(shù)可以得到不同的濾波器濾過器的設(shè)計(jì)和選擇濾過器的設(shè)計(jì)和選擇Ram-Lak:usingtherectangularwindow矩形窗Shepp-Logan:usingasincwindow正弦窗Cosine:usingacosinewindow余弦窗Hamming:usingageneralizedHammingwindow通用Hamming窗濾過器的設(shè)計(jì)和選擇濾過器的設(shè)計(jì)和選擇-Ramp濾波頻率與幅度成正比關(guān)系(a)圖為理想濾波器(b)圖為修正后濾波器

亦理論上濾波器亦稱為Ramp濾波器,其高頻分量是無(wú)限延伸的,但實(shí)際實(shí)現(xiàn)時(shí)必須截?cái)嗵幚?,如圖(b)圖中虛線所示,相當(dāng)于在帶寬之外突然衰減為零,在重建圖像的邊緣時(shí)會(huì)出現(xiàn)環(huán)狀震蕩條紋,稱之為Gibbs現(xiàn)象。為有效地消除此現(xiàn)象,我們需對(duì)Ramp濾波器稍作平滑處理,如將之與作卷積,得到Shepp-Logan濾波器;

濾過器的設(shè)計(jì)和選擇-Ramp濾波濾過器的設(shè)計(jì)和選擇-Shepp-Logan濾波器平滑了圖像,損失了部分高頻信息濾過器的設(shè)計(jì)和選擇-Hamming濾波器降低了高頻噪聲,可得到Hamming濾波器和Hanning濾波器骨濾過器和軟組織濾過器GELightSpeedSoft軟組織Standard標(biāo)準(zhǔn)Detail細(xì)節(jié)根據(jù)診斷需求可選用不同的濾波函數(shù)平滑用于觀察軟組織銳利用于觀察高分辨力影像Lung肺Bone骨Edge邊緣FBP中的補(bǔ)0運(yùn)算原始濾波運(yùn)算包含一個(gè)非周期卷積運(yùn)算,變到頻域后就是周期卷積,直接計(jì)算將產(chǎn)竹干涉?zhèn)斡?,即所謂的warp-around效應(yīng)。因此必須在傅立葉變換和濾波操作之前給每一個(gè)投影補(bǔ)0,才能避免偽影產(chǎn)生由于補(bǔ)零前后數(shù)據(jù)長(zhǎng)度不一樣,它們的分辨率(分別為df1,df2)不一樣,在頻域中譜線所代表的頻率也不一樣,所以這兩個(gè)頻譜所描述的對(duì)象也不相同特殊情況:M=(2^n)N時(shí),補(bǔ)零后的頻譜相當(dāng)于在補(bǔ)零前的頻譜中插入(2^n)-1條譜線。與補(bǔ)零前的頻譜中相重合的譜線,它們的幅值和相位完全一致FBP中的補(bǔ)0運(yùn)算這里舉一個(gè)例子,fs=1000HZ,補(bǔ)零前后數(shù)據(jù)長(zhǎng)度N=500和M=800,對(duì)應(yīng)的df1=2,df2=1.25。補(bǔ)零前的頻譜是對(duì)應(yīng)于0,2,4,...,500HZ的頻譜,而補(bǔ)零后的頻譜是對(duì)應(yīng)于0,1.25,2.5,...,500HZ的頻譜,所以兩頻譜中對(duì)應(yīng)頻率不同,描述當(dāng)然不同測(cè)量每個(gè)視角的投影p(t,)對(duì)投影p(t,)補(bǔ)0,得到p(t,)傅立葉變換,得到P(,)乘以斜坡濾波器H(),得到G(,)傅立葉逆變換,得到濾波投影g(t,)。反投影g(t,),并加入圖像f(x,y)。循環(huán)結(jié)束?FBP實(shí)現(xiàn)步驟(平行束)扇形束重建(FanbeamReconstruction)學(xué)習(xí)內(nèi)容(Learningobjects)扇形束到平行束的轉(zhuǎn)換等角線束重建等距線束重建扇形束重建成像幾何扇形束重建扇束情況下的重建算法較為復(fù)雜,但實(shí)質(zhì)沒有改變??刹捎闷叫惺闆r下的算法實(shí)現(xiàn),只需加以適當(dāng)?shù)匦拚纯芍嘏潘惴?把一個(gè)視圖中采得的扇形數(shù)據(jù)重新組合成平行的

射線投影數(shù)據(jù),然后采用平行束重建算法重建直接重建算法:不必?cái)?shù)據(jù)重排,只需適當(dāng)加權(quán)即可運(yùn)用與平

行束類似的算法重建扇形束重建-從扇形束向平行束轉(zhuǎn)換橫軸表示X射線離開中心點(diǎn)的距離,縱軸表示X射線與x軸形成的夾角在平行投影重建中,射線由兩參數(shù)t和唯一決定,其中t是射線到等中心的距離,是投影角一條投影線映射到弦空間中就是一個(gè)點(diǎn)

一個(gè)投影集在弦空間就映射成一個(gè)均勻排列的網(wǎng)格網(wǎng)格中的某一行表示其中一個(gè)投影,如a圖中虛線長(zhǎng)方形所示。如果把一個(gè)扇形投影集映射到同一弦空間,那么每一扇形投影映射成傾斜的一行點(diǎn)陣,如b圖中虛線矩形所示。每個(gè)樣本的角度不同。扇形束投影射線與等中心點(diǎn)的距離也不是按照固定的規(guī)律增減的。扇形束重建-從扇形束向平行束轉(zhuǎn)換扇形束重建-EA等角扇形束重建當(dāng)同樣大小的探測(cè)器單元沿著中心為X射線焦點(diǎn)的弧排列時(shí),就形成等角采樣扇形束的每一條射線可由β和γ確定,其中是射線與中心射線(假想的通過X射線源和等中心的直線)的夾角,稱為探測(cè)器角;是中心射線與y軸的夾角,稱為投影角

扇形束重建-EA等角扇形束重建投影乘以探測(cè)器角的余弦,濾波后的樣本隨著到光源的距離的增長(zhǎng)而增長(zhǎng)重建公式可由用(t,θ)坐標(biāo)確定(γ,β)坐標(biāo)上的每個(gè)樣本來得到。扇形投影中的投影樣本q(,)就轉(zhuǎn)化為平行投影中的投影樣本p(t,)

平行線束重建vs扇形束重建與平行投影重建不同,扇形投影重建在濾波操作之前先乘以Dcos,它獨(dú)立于投影角,所以可在重建之前先進(jìn)行計(jì)算并保存

第二個(gè)區(qū)別是扇形投影重建中用到了加權(quán)反投,反投影沿扇形進(jìn)行,尺度因子L-2隨著像素的不同而不同

等角扇形線束重建流程圖平行束vs等角扇束對(duì)單個(gè)濾波投影做反投影實(shí)驗(yàn)圖示,扇形束反投影重建的圖像不但其形狀隨離開X射線源的距離而改變(圖中是6點(diǎn)位置),而且強(qiáng)度也在變化平行束扇形束

扇形束重建-ED等角扇形束重建扇形束重建-Matlab函數(shù)Matlab函數(shù):fanbeam()ifanbeam()CT圖像重建小結(jié)圖像重建算法比較傅立葉變換法反投影法濾波反投影法算法比較算法比較傅里葉變換重建方法:對(duì)于每次測(cè)得的投影數(shù)據(jù)先作一維傅里葉變換,根據(jù)中心切片定理,可將此變換結(jié)果看成二維頻域中同樣角度下過原點(diǎn)的直線上的值。在不同投影角度下所得的一維變換函數(shù)可在頻域中構(gòu)成完整的二維傅里葉變換函數(shù),將此二維變換函數(shù)進(jìn)行逆變換,就得到了所要求的空間域中的密度函數(shù)。傅里葉變換的方法重建圖像時(shí),投影函數(shù)的一維傅里葉變換在頻域中表現(xiàn)為極坐標(biāo)的形式,把極坐標(biāo)形式的數(shù)據(jù)通過插補(bǔ)運(yùn)算轉(zhuǎn)換為直角坐標(biāo)形式的數(shù)據(jù)時(shí),計(jì)算的工作量比較

溫馨提示

  • 1. 本站所有資源如無(wú)特殊說明,都需要本地電腦安裝OFFICE2007和PDF閱讀器。圖紙軟件為CAD,CAXA,PROE,UG,SolidWorks等.壓縮文件請(qǐng)下載最新的WinRAR軟件解壓。
  • 2. 本站的文檔不包含任何第三方提供的附件圖紙等,如果需要附件,請(qǐng)聯(lián)系上傳者。文件的所有權(quán)益歸上傳用戶所有。
  • 3. 本站RAR壓縮包中若帶圖紙,網(wǎng)頁(yè)內(nèi)容里面會(huì)有圖紙預(yù)覽,若沒有圖紙預(yù)覽就沒有圖紙。
  • 4. 未經(jīng)權(quán)益所有人同意不得將文件中的內(nèi)容挪作商業(yè)或盈利用途。
  • 5. 人人文庫(kù)網(wǎng)僅提供信息存儲(chǔ)空間,僅對(duì)用戶上傳內(nèi)容的表現(xiàn)方式做保護(hù)處理,對(duì)用戶上傳分享的文檔內(nèi)容本身不做任何修改或編輯,并不能對(duì)任何下載內(nèi)容負(fù)責(zé)。
  • 6. 下載文件中如有侵權(quán)或不適當(dāng)內(nèi)容,請(qǐng)與我們聯(lián)系,我們立即糾正。
  • 7. 本站不保證下載資源的準(zhǔn)確性、安全性和完整性, 同時(shí)也不承擔(dān)用戶因使用這些下載資源對(duì)自己和他人造成任何形式的傷害或損失。

評(píng)論

0/150

提交評(píng)論