《MATLAB实现的医学图像FBP算法》由会员分享,可在线阅读,更多相关《MATLAB实现的医学图像FBP算法(15页珍藏版)》请在金锄头文库上搜索。
1、TheFilteredBackprojectionAlgorithminMATLAB,GregGallardo51:185,FilteredBackprojectionAlgorithm,MeasurethefamilyofprojectionsK=#ofprojections(angles),N=#ofrays(detectorsize)PerformFFTMultiplyFFTbyFFT.PerformIFFTproductfromstep3Performbackprojectionas:,1)Measurethefamilyofprojections,R,Xp=RADON(.)retur
2、nstwovariablesMatrixRcolumnsaretheRadontransformfortheangles.RowsaredetectorpositionVectorXp-containingtheradialcoordinatescorrespondingtoeachrowofR.,1)Measurethefamilyofprojections,Functionparameterssize():getsNandKnextpow2():getwidthforFFT,2)PerformFFT,MATLABstandsforMATrixLABoratoryMatrixoperatio
3、nsarefasterthanvisitingeachelementinaloop.,3)MultiplyFFTbyFFT,BackprojectionfiltersShepp-Logan:Ram-Lakmultipliedbysincfunction(seeiradonhelp)Cosine:Ram-LakmultipliedbycosineHamming:Ram-LakmultipliedbyHammingwindowHanning,Blackman,etc.,3)MultiplyFFTbyFFT,Ram-Lak,Shepp-LoganandCosinefiltersinfrequency
4、domain,3)MultiplyFFTbyFFT,FIRFilters.CommonlyusedWindows,Figurefrom“Discrete-TimeSignalProcessing”,Oppenheim&Schafer,Prentice-Hall,3)MultiplyFFTbyFFT,RectangularBartlett(triangular)Hanning(vonHann)HammingBlackman,from“Discrete-TimeSignalProcessing”,Oppenheim&Schafer,Prentice-Hall,3)MultiplyFFTbyFFT,3)MultiplyFFTbyFFT,3)MultiplyFFTbyFFT,Filterstep.Elementbyelementmultiplication,4)IFFTproductfromstep3,5)PerformBackprojectionas:,Results,Results,