FULLY3D'03TUAM2–31ExtendedParallelBackprojection(EPBP)forArbitraryConeAngleandArbitraryPitch3DandPhase–Correlated4DCTReconstructionMarcKachelrie,WilliKalenderAbstract—RecentdevelopmentsinmedicalCTaimatfasterrotationspeedsandahighernumberofsimultane-ouslyacquiredslices.
TheseeortsarepushedfurtherbycardiacCTwhichiscurrentlythemostprominentspecial–purposeapplicationinCT[1].
Today,16-slicescannersarestate–of–the–art.
ButCTmanufacturershavealreadyan-nouncedscannerswithfarmoreslicesandsomeevenpro-moteprototypeswithupto256slices.
MedicalCTmustsupportthecircularscantrajectory(se-quencescan)andthespiraltrajectory(spiralscan).
Arbi-trarypitchselectionisofhighimportance.
Inanycase,thecompleteareaofthedetectoristobeexposedandeachmeasuredrayshouldcontributetotheimage.
Onlythen,optimizeddoseusagecanbeachieved.
Theserequirementscannotbefullledbycurrentre-constructionapproaches.
Exactcone–beamreconstruction,thatiscapableofreconstructinglargecone–angledata,can-notcopewitharbitrarypitchandphase–correlateddataseg-ments.
Onlyapproximatereconstructionapproacheshavethepotentialtohandleallrequirements.
Currently,theonlyknownapproachthatcanhandlephase–correlatedtruecone–beamdataisanextensiontotheAdvancedSingle–SliceRebinning(ASSR)algorithm[2],[3].
However,thisgeneralizedapproachASSRCIislimitedto32slices.
WehavethereforedevelopedanewapproximateFeldkamp–typealgorithm,theextendedparallelbackpro-jection(EPBP)[4].
Itsmainfeaturesareaphase–weightedbackprojectionandavoxel–by–voxel180normalization.
Therstfeatureensures3Dand4Dcapabilitieswithoneandthesamealgorithm,thesecondensures100%detectorusage(eachraycounts!
).
Thealgorithmisevaluatedusingsimulateddataofathoraxphantomandacardiacmotionphantomforscannerswithupto256slices.
Thestandardreconstructions(EPBPStd)areofexcellentqualityevenforasmanyas256slices.
Thecardiacrecon-structions(EPBPCI)areofhighqualityaswellandshownosignicantdeteriorationofobjectsevenfarothecenterofrotation.
SinceEPBPCIusesthecardiointerpolation(CI)phaseweightingthetemporalresolutionisequivalenttothatofthewellknownsingle–sliceandmulti–slicecardiacapproaches180CI,180MCI,andASSRCI,respectively,andliesintheorderof50msto100msforrotationtimesbetween0.
4sand0.
5s.
I.
IntroductionMEDICALcomputedtomographyiscurrentlyevolv-ingfasterthanever.
Increasedspatialresolution,de-creasedscantime,increasedtemporalresolution,decreasedpatientdose,andincreasedvolumecoveragearesomeoftheimportanttrendstomention.
Aslittleasveyearsago,single–slicespiralCTwasthestate–of–the–art.
Then,therst4–slicescannersbecameavailable.
Alreadyin2001,16–slicescannersstartedtore-InstituteofMedicalPhysics,UniversityofErlangen–N¨urnberg,Krankenhausstr.
12,91054Erlangen.
Correspondingauthor:PDDr.
MarcKachelrie,E–mail:marc.
kachelriess@imp.
uni–erlangen.
deplacethe4–slicemachines.
Apparently,thenearfuturewillshiftthenumberofsimultaneouslyscannedslicesto32,64andevenmore.
Slicethickness,andtherebyspatialresolu-tion,willcontinuetodecreasetofurtherimprovediagnosticaccuracy.
Atthesametime,doseutilizationwillincreasetokeeptheeectivepatientdoseatanacceptablelevel:dosemodulationtechniques,automaticexposurecontrolandimproveddetectormaterialswillhelptodoso[1].
Besidesimprovedspatialresolution,improvedcontrastresolutionandlowpatientdoseoneisfurtherinterestedinhighesttemporalresolutiontoallowimagingtheheart.
Thisisdoneusingshortrotationtimescombinedwithde-dicatedphase–correlatedreconstructionalgorithmsastheyareavailablesince1996/1997[5].
Therelevanceforimagereconstructionalgorithmsliesintheexpectedincreaseinconeangleandinthedemandforphase–correlatedreconstruction.
Neglectingthecone–angleofthescannerasdoneinall4–sliceandinallcom-mercialcardiac16–slicereconstructionalgorithmswillyieldunacceptableimageartifactsforfuturescanners.
Therearefastandecientcone–beamalgorithmswithacceptableimagequalityavailableforcone–beamspiralCT.
However,theydonotensure100%doseusageand/ortheydonotworkforarbitraryspiralpitch.
Further,noneofthemiscapableofcombiningthesedemandswithphase–correlatedimagingatwiderconeangles.
Tollthemissinggap,weproposeEPBP,anewapprox-imatecone–beamreconstructionthatallowsforarbitrarypitchandforphase–correlatedreconstructionandthaten-sures100%detectorusage.
DetailsofcardiacCTscanningandtherestrictionsonthemaximumpitchasafunctionofthepatient'sheartratecanbefoundinreferences[6],[7],[8].
Thephaseselectionmechanismsdescribedthere,namelythecardiointerpolation(CI)andthecardiodelta(CD)approachareusedforEPBPaswell.
EspeciallyEPBPCIthatusesCI,themostpromisingmulti–phaseweightingcurrentlyavailable,isevaluatedhere.
Inthispaper,wewilloutlinetheEPBPStdandtheEPBPCIalgorithmandgivesomedescriptiveexamples.
II.
SimulationsToevaluateournewapproachwehavesimulatedspiralcone–beamdatacorrespondingtothein–planegeometryofatypicalmedicalCTscanner(1160projectionsperro-tation,672detectorchannelsperdetectorrow,andafanangleΦ=52)usingadedicatedx–raysimulationtool(ImpactSim,VAMPGmbH,M¨ohrendorf,Germany).
Two2FULLY3D'03TUAM2–3phantomshavebeensimulated:thethoraxphantomde-scribedinthephantomdatabasehttp://www.
imp.
uni–erlangen.
de/forbildandthecardiacmotionphantomde-scribedin[7].
Thesimulatedscanprotocoluses0.
42srotationtime(143rpm),M*0.
75mmcollimationwithM=2msimul-taneouslyscannedsliceswherem=4,8,andapitchofp=0.
375.
Thetableincrementcanbefoundasd=MSp.
III.
GeometryA.
ScanGeometryThescangeometryassumedhereisafan–beamgeome-trywithcylindricaldetectorsandaspiralfocustrajectory.
TheEPBPapproachisbasedonarebinningtoparallelgeometry.
Othergeometries,suchasatdetectors,canbeeasilyincorporatedbymodifyingthecorrespondingrebin-ningandtransformequations.
Notethatinthelimitofp=0,thespiralreducestoacirculartrajectory.
EPBPcopeswithsequencedataaswellasitdoeswithspiraldata.
Thesourcetrajectoryisparameterizedbytheviewangleαass(α)=RFsinαcosα0+d001α2π.
(1)RFdenotestheradiusofthefocalspottrajectoryandddenotesthetableincrementperrotation.
Thecoordinatevectorofthedetectorelement(α,β,b)isgivenasr(α,β,b)=s(α)+RFDsin(α+β)cos(α+β)0+b001;βandbdenotethetransaxialandlongitudinaldetectorcomponents,respectively.
Torebinthetransaxialcomponentsofthecone–beamdatatoparallelgeometryweparameterizearaybyitsdis-tanceξtotheaxisofrotationandbyitsanglewithre-specttothenegativey–axis.
Thenormalformoftheray'sx–y–componentsisgivenasxcos+ysinξ=0.
Thisdenitionwaschosentohavethecentralraysforfan–beam(β=0)andforparallelbeam(ξ=0)coincidingforα=.
Therelationbetweenarayinfan–beamcoordinates(α,β)andaparallel–beamray(,ξ)isthefamiliartransform=α+βξ=RFsinβandα=+arcsinξ/RFβ=arcsinξ/RF.
(2)B.
PointProjectionForbackprojection,weneedtoknowthedetectorcoor-dinates(ξ,b)thatresultfromprojectingthepoint(x,y,z)froms(α)ontothecylindricaldetector.
Theradialcoordi-nateisgivenasξ=xcos+ysin.
Thelongitudinaldetectorcoordinatecanbecomputedus-ingtheintersectiontheorem.
ThetransaxialdistanceoftherespectivevoxeltothesourceisgivenasD2=(RFsinαx)2+(RFcosα+y)2=R2F2RFrsin(α)+r2with(x,y)=(rcos,rsin)or,equivalently,D=RFcosβ+η=R2Fξ2+ηwithη=ycosxsin.
Now,ndbbyscalingtheaxialdistancezdα2πfromDtoRFD:b=RFDD(zdα2π).
(3a)Andwendanotherrepresentationofξ:ξ=RFsinβ=RFxcosα+ysinαD.
(3b)IV.
ReconstructionTheextendedparallelbackprojectionalgorithmconsistsofthefollowingvesteps:-azimuthalrebinning:p(α,β,b)→p(,β,b),-longitudinalrebinning:p(,β,b)→p(,β,l),-radialrebinning:p(,β,l)→p(,ξ,l),-convolution:p(,ξ,l)→p(,ξ,l),-weightingandbackprojection:p(,ξ,l)→f(x,y,z).
A.
AzimuthalRebinningTheoriginalprojectiondatap(α,β,b)areconvertedfromfan–beamtofan–parallelgeometryusing(2)asfollows:p(,β,b)=p(α,β,b)withα=β.
B.
LongitudinalRebinningConvolvingspiraldatainthedetectorrowdirection(con-stantb)yieldsseverecone–beamartifacts.
AsindicatedbyASSR[2],SMPR[9],andexactcone–beamreconstruc-tion[10]theoptimaldirectionofconvolutionisthetangentds(α)/dα.
Toalignthefan–paralleldetectorrowswiththeoptimaldirectionofconvolutionalongitudinalrebinningisrequired.
Therefore,weareinterestedintherelationshipbetweenbandξwhenmovingalongds.
Using(3)onendsdbdξ=d(bD)d(ξD)=RFDdzRF(dxcosα+dysinα)=RFDd2πR2F;inthelaststep(1)wasusedtoinsertthecomponentsofds.
Now,wedeneanewlongitudinalvariablelasb=l+λξwithλ=dbdξ=dRFD2πR2Fsuchthatdl/dξ=0inthedirectionofds.
Then,dothelongitudinalrebinningp(,β,l)=p(,β,b)withb=l+λξ=lλRFsinβtoswitchtolasthenewindependentvariable.
FULLY3D'03TUAM2–33Wheneverbexceedsthedetectorlimitsbminandbmaxforsomeβweextrapolatebyrepeatingtheoutermostde-tectorrow.
Thevaluesthataremadeupbyextrapolationarerequiredduringconvolutionwhichoperatesoncom-pletedetectorrows,always.
Duringbackprojectiontheseextrapolatedpointsarenotaccessed;backprojectionratherrespectsthephysicaldetectorarea!
C.
RadialRebinningTheradialrebinningconvertstoequidistantparallelco-ordinates.
Weuse(2)tondp(,ξ,l)=p(,β,l)withβ=arcsinξ/RF.
D.
ConvolutionNow,convolutionofthedetectorrowsisperformedusingastandardconvolutionkernelk(ξ),as,forexample,theShepp–Logankernel.
p(,ξ,l)=p(,ξ,l)k(ξ)yieldstheconvolveddatapneededforbackprojection.
E.
WeightingandBackprojectionInthisstep,weregardthebackprojectionofaxedvoxel,sayonelocatedatr=(x,y,z).
LetVdenotethesetofviewanglesunderwhichrismeasured.
AssumeatemporalwindowTthatcomprisesallthatcorrespondtoalloweddata.
Forthestandardreconstruc-tionEPBPStd,alldataacquiredarealloweddataandthereforeT=R.
Forthereconstructionofcardiacdata,TcanbedenedbyspecifyingacardiacmotionphasecRthatcountsrelativetosomesynchronizationpeaksandaphasewidth0EPBPCIchoosescassmallasallowedbythecompletenesscondition(seebelow).
Otherdenitionsmayincludeabsolutetiminginformationortherestrictiontoonlyonetemporalintervaloflengthπ(single–phasere-construction),ortwointervals(bi–phasereconstruction).
RegardlessofwhatconventionisusedtodeneT,theintersectionI=V∩T,thatcomprisesallviewstobeused,mustbe180–complete:k(I+kπ)=R.
Now,assumeaweightingfunctionw()whosesupportequalsI,i.
e.
w(R\I)={0},andkw(+kπ)=0.
ThelastconditioncaneasilybeachievedbyusingpositiveweightsonIonly.
ForEPBPCIweuseamulti–triangularweightfunction:trianglefunctionslocatedoneachofI'sdisjunctintervals.
Bynormalizingwasw()=w()kw(+kπ)weachievekw(+kπ)=1anddw()=π.
Backprojectionf(x,y,z)=dp(,ξ,l)w()withξ=ξ(x,y,)=xcos+ysinα=β=+arcsinξ/RFl=b(x,y,z,α)λξthenyieldsthedesiredvoxelvalueat(x,y,z).
V.
ResultsFigure1showsthatimagequalityofthethoraxphantomisexcellentwithEPBP,evenforasmanyas256slices.
Asindicatedbytheribs,ASSR(whichisinfactdesignedforuptoaboutonly60slices[2])cannotcopewiththislargecone–angle;thesameappliestothehighlyrelatedAMPRalgorithmdenedinreference[11].
Fig.
1.
Thorax,scannedwith256*0.
75mmcollimationandd=72mmtableincrement.
Aheartrateof120min1wassimulatedforEPBPCI.
(0/500)ConsideringthatEPBPCIusesonlyafractionofthedataavailable(here,roughly25%),dependingonthelocalheartrateandonthereconstructionposition,theimagesarealmostasgoodastheEPBPStdreconstructions,apartfromtheincreasedimagenoise.
Theonlyexceptionisaslightvariationinthereconstructeddensityclosetothevertebrae.
Figure2showsreconstructionsofthecardiacmotionphantomfora16–sliceanda256–slicescanner.
Sincethe4FULLY3D'03TUAM2–3(a)16*0.
75mmcollimation(b)256*0.
75mmcollimationFig.
2.
TransaxialslicesandMPRsofthecardiacmotionphantomreconstructedwithvariousalgorithmsusingastandard16–slicescannerandawidecone–angle256–sliceCT.
ThesagittalMPRs(bottom)showadditionalδ–objectsusedtomeasureresolution.
(0/500)eldofviewshowsonlythecentralpartsofthepatient,thein–planeimagesoftheASSRapproachareacceptableevenforthe256–slicescanner.
However,with256–sliceASSRthemultiplanarreformations(MPRs)tendtobeblurredinthez–directionandfullwidthathalfmaximumFWHMzoftheslicesensitivityproleisincreasedsignicantlywhereasthein–planeresolutionFWHMxyisthesameasforthe16–slicecase.
EPBP,incontrast,behavesverywellforallsimulatedscanners(16,32,64,128and256slices).
Spatialresolu-tionisslightlyhigherthanforthesingle–slicerebinningalgorithms.
ImagenoiseincreasesforEPBPCIduetothephase–weighting.
Thisobservationisvalidforalltheothersimulatedgeometriesandheartrates(wehavelookedintofH=40min1140min1).
EPBPgenerallybehavesequaltoorbetterthanASSR.
Finally,gure3givesanexampleofreconstructedpa-tientdata.
Thedatashownarecorrelatedtothepatientmotionfunction,theso–calledkymogram,whichcandi-rectlybederivedfromtheacquiredrawdata[12].
ThestandardreconstructionsofASSRStdandEPBPStdarecomparableduetothelownumberofslices;thephase–correlatedEPBPCIimagesareofhighimagequalityandcorrespondtothegold–standardincardiacCTimaging.
VI.
DiscussionTheextendedparallelbackprojectionappearstobeade-quateformedicalCTimagereconstructioninallrespects.
EPBPimagequalityisequivalenttoexisting4–or16–slicestandardandcardiacalgorithmsforawiderangeofFULLY3D'03TUAM2–35Fig.
3.
ASSRandEPBP,12*0.
75mmcollimation,3.
375mmtableincrement.
Thephase–correlatedimagesarereconstructedrelativetothekymogrampeaksat0%and50%ofK-K,respectively.
(0/500)simultaneouslyscannedslices.
EvendatawithM=256slicesyieldsexcellentimagequality.
Forstandardrecon-structionsthisisnotsurprisingsinceEPBPStdissimilartootherFeldkampalgorithms(aslongastheseperformconvolutionalongthetangentdirection).
Forwideconeanglecardiacdata,wherenootherphase–correlatedcone–beamalgorithmisavailableyet,EPBPCIperformsverywellevenforobjectsfarotheisocenter(ribsingure1).
Feldkamp–typealgorithmsaresuperiortoASSR[2],AMPR[11],orSMPR[9]forlargeM.
Wideconeanglecar-diacCTiscurrentlyonlypossiblewithEPBP.
Itsuniqueweightingstrategythatassignsindividualdatarangestoeachvoxel,ensures100%datausageandthusthemax-imumdoseutilizationpossible.
ThefutureofmedicalsequentialandspiralCTwillcertainlyincludetheideaofphase–correlated/phase–weighted3DbackprojectionofEPBP–type,ormodicationsthereof.
References[1]W.
A.
Kalender,ComputedTomography.
Wiley&Sons,2003.
[2]M.
Kachelrie,S.
Schaller,andW.
A.
Kalender,"Advancedsingle–slicerebinningincone–beamspiralCT,"Med.
Phys.
,vol.
27,pp.
754–772,Apr.
2000.
[3]M.
Kachelrie,T.
Fuchs,R.
Lapp,D.
-A.
Sennst,S.
Schaller,andW.
Kalender,"ImagetovolumeweightinggeneralizedASSRforarbitrarypitch3Dandphase–correlated4Dspiralcone–beamCTreconstruction,"Proc.
ofthe6thInt.
MeetingonFully3DImageReconstruction,pp.
179–182,Nov.
2001.
[4]M.
KachelrieandW.
A.
Kalender,"Extendedparallelbackpro-jectionforcardiaccone–beamCTforupto128slices,"Radiol-ogy,vol.
225(P),p.
310,Nov.
2002.
[5]M.
KachelrieandW.
A.
Kalender,"ECG–basedphase–orientedimagereconstructionfromsubsecondspiralCTscansoftheheart,"Radiology,vol.
205(P),p.
215,Nov.
1997.
[6]M.
KachelrieandW.
A.
Kalender,"Electrocardiogram–correlatedimagereconstructionfromsubsecondspiralcomputedtomographyscansoftheheart,"Med.
Phys.
,vol.
25,pp.
2417–2431,Dec.
1998.
[7]M.
Kachelrie,S.
Ulzheimer,andW.
A.
Kalender,"ECG–correlatedimagereconstructionfromsubsecondmulti–slicespi-ralCTscansoftheheart,"Med.
Phys.
,vol.
27,pp.
1881–1902,Aug.
2000.
[8]M.
Kachelrie,S.
Ulzheimer,andW.
A.
Kalender,"ECG–correlatedimagingoftheheartwithsubsecondmultisliceCT,"IEEETransactionsonMedicalImaging,vol.
19,pp.
888–901,Sept.
2000.
[9]K.
Stierstorfer,T.
Flohr,andH.
Bruder,"Segmentedmultipleplanereconstruction–anovelapproximatereconstructionformulti–slicespiralCT,"Phys.
Med.
Biol.
,vol.
47,pp.
2571–2851,2002.
[10]K.
Sourbelle,PerformanceEvaluationofExactandApprox-imateCone–BeamAlgorithmsinSpiralComputedTomogra-phy.
PhDThesis,Friedrich–Alexander–Universit¨atErlangen–N¨urnberg,2002.
[11]S.
Schaller,K.
Stierstorfer,H.
Bruder,M.
Kachelrie,andT.
Flohr,"Novelapproximateapproachforhigh-qualityimagereconstructioninhelicalconebeamCTatarbitrarypitch,"SPIEMedicalImagingConferenceProc.
,vol.
4322,pp.
113–127,2001.
[12]M.
Kachelrie,D.
-A.
Sennst,W.
Maxlmoser,andW.
A.
Kalen-der,"Kymogramdetectionandkymogram–correlatedimagere-constructionfromsubsecondspiralcomputedtomographyscansoftheheart,"Med.
Phys.
,vol.
29,pp.
1489–1503,July2002.
香港服务器多少钱一个月?香港服务器租用配置价格一个月多少,现在很多中小型企业在建站时都会租用香港服务器,租用香港服务器可以使网站访问更流畅、稳定性更好,安全性会更高等等。香港服务器的租用和其他地区的服务器租用配置元素都是一样的,那么为什么香港服务器那么受欢迎呢,香港云服务器最便宜价格多少钱一个月呢?阿里云轻量应用服务器最便宜的是1核1G峰值带宽30Mbps,24元/月,288元/年。不过我们一般选...
乐凝网络怎么样?乐凝网络是一家新兴的云服务器商家,目前主要提供香港CN2 GIA、美国CUVIP、美国CERA、日本东京CN2等云服务器及云挂机宝等服务。乐凝网络提供比同行更多的售后服务,让您在使用过程中更加省心,使用零云服务器,可免费享受超过50项运维服务,1分钟内极速响应,平均20分钟内解决运维问题,助您无忧上云。目前,香港HKBN/美国cera云服务器,低至9.88元/月起,支持24小时无理...
paypal贝宝可撸$10的代金券!这两天paypal出了活动,本次并没有其他的限制,只要注册国区的paypal,使用国内的手机号和62开头的银联卡,就可以获得10美元的代金券,这个代金券购买产品需要大于10.1美元,站长给大家推荐几个方式,可以白嫖一年的VPS,有需要的朋友可以看看比较简单。PayPal送10美元活动:点击直达活动sfz与绑定卡的号码可以重复用 注册的邮箱,手机号与绑的银联卡必须...
小米3手机真实图片为你推荐
桌面背景图片大全手机壁纸图片大全麒麟990和骁龙865哪个好骁龙八65和麒麟990谁强一点?小说软件哪个好用免费什么软件看小说全免费游戏盒子哪个好游戏盒子哪个好?核芯显卡与独立显卡哪个好独立显卡和核芯显卡有什么区别网校哪个好有什么网校比较好云盘哪个好云盘有哪些,哪个云盘好360云盘登陆怎样登陆360我的云盘360云u盘360云u盘无法连接怎么回事?360云盘下载速度慢怎么办360云盘为什么下载速度如此慢?
vps交流 什么是域名地址 国内免备案主机 优惠码 la域名 godaddy优惠码 mobaxterm lighttpd 一元域名 丹弗 绍兴高防 777te 促正网秒杀 中国电信测网速 tna官网 免费测手机号 ca187 drupal安装 免费外链相册 网站加速软件 更多