ct重建實(shí)驗(yàn)指導(dǎo)手冊(cè)20xx-在線瀏覽
2025-01-06 20:39本頁(yè)面
【正文】 t(1:256,f3)。 圖 14 垂直線束投影的一維傅立葉圖像 葉變換。%對(duì)原始圖像進(jìn)行二維傅立葉變換。%取結(jié)果的模 f3=fftshift(fab2D)。 %顯示圖像 所得的結(jié)果如圖 15所示。 ：通過(guò)對(duì)以上程序的理解，如何能證明在 ? = 090 也就是 x線平行 x軸時(shí)傅立葉切片定理的正確性 ，并編寫(xiě) matlab程序加以證明。 3. 實(shí)驗(yàn)步驟 1) 在 Matlab程序的 Command window內(nèi)輸入如下命令： help radon 在 Matlab程序的 Command window內(nèi)輸入如下命令： help iradon 了解 radon與 iradon函數(shù)的用法 。 C = radon(P,0:179)。 figure， imshow（ C， []） 。 得出結(jié)果 如圖 21所示 。 （ 1） 反投影重建 shepplogan體模圖像，并得到它的投影值。 figure,imshow(I)。 THETA=0:179。t lose anything when we rotate. [iLength, iWidth] = size(IMG)。 LengthPad = ceil(iDiag iLength) + 2。 padIMG = zeros(iLength+LengthPad, iWidth+WidthPad)。 % loop over the number of angles, rotate 90theta (because we can easily sum % if we look at stuff from the top), and then add up. Don39。 PR = zeros(size(padIMG,2), n)。bilinear39。crop39。 PR(:,i) = (sum(tmpimg))39。 得到的圖像如圖 21所示 圖 22 a原始體模 b所得的投影圖 22b進(jìn)行反投影重建。 sideSize = n。 % set up the image m = length(THETA)。 % find the middle index of the projections midindex = (n+1)/2。 y = 1:sideSize。 xpr = X (sideSize+1)/2。 % loop over each projection for i = 1:m tic disp([39。, num2str(THETA(i))])。 BPIa = ze
醫(yī)療健康相關(guān)推薦
鄂ICP備17016276號(hào)-1