ct重建實驗指導手冊20xx-文庫吧資料
2024-11-11 20:39本頁面
【正文】 PIa。 newfiltIndex = filtIndex(spota)。 spota = find((filtIndex 0) amp。 filtIndex = round(midindex + xpr*sin(th(i)) ypr*cos(th(i)))。On angle 39。 ypr = Y (sideSize+1)/2。 [X,Y] = meshgrid(x,y)。 % set up x and y matrices x = 1:sideSize。 BPI = zeros(sideSize,sideSize)。 % convert THETA to radians th = (pi/180)*THETA。 % back projection n = size(PR,1)。 THETA(i) toc end figure,imshow(PR,[])。)。, 39。 for i = 1:n tic tmpimg = imrotate(padIMG, 90THETA(i), 39。t perform any % interpolation on the rotating. n = length(THETA)。 padIMG(ceil(LengthPad/2):(ceil(LengthPad/2)+iLength1), ... ceil(WidthPad/2):(ceil(WidthPad/2)+iWidth1)) = IMG。 WidthPad = ceil(iDiag iWidth) + 2。 iDiag = sqrt(iLength^2 + iWidth^2)。 % This MATLAB code takes an image matrix and vector of angles and then % finds the 1D projection (Radon transform) at each of the angles. It returns % a matrix whose columns are the projections at each angle. % pad the image with zeros so we don39。 IMG=double(I)。 %form the image I=phantom(256)。 圖 21 a iradon重建圖像 b 投影圖 2) 首先， 我們 分別對前面所創(chuàng)建的 shepplogan進行 平行束 反投影重建以及濾波反投影重建。 figure， i
醫(yī)療健康相關推薦
鄂ICP備17016276號-1