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.
2021年恒创科技618活动香港美国服务器/云服务器/高防全场3折抢购,老客户续费送时长,每日限量秒杀。云服务器每款限量抢购,香港美国独服/高防每款限量5台/天,香港节点是CN2线路还不错。福利一:爆品秒杀 超低价秒杀,秒完即止;福利二:云服务器 火爆机型 3折疯抢;福利三:物理服务器 爆款直降 800元/月起;福利四:DDOS防护 超强防御仅 1750元/月。点击进入:2021年恒创科技618活...
易探云香港vps主机价格多少钱?香港vps主机租用费用大体上是由配置决定的,我们选择香港vps主机租用最大的优势是免备案vps。但是,每家服务商的机房、配置、定价也不同。我们以最基础配置为标准,综合比对各大香港vps主机供应商的价格,即可选到高性能、价格适中的香港vps主机。通常1核CPU、1G内存、2Mbps独享带宽,价格在30元-120元/月。不过,易探云香港vps主机推出四个机房的优惠活动,...
mineserver怎么样?mineserver是一家国人商家,主要提供香港CN2 KVM VPS、香港CMI KVM VPS、日本CN2 KVM VPS、洛杉矶cn2 gia端口转发等服务,云服务器网(yuntue.com)介绍过几次,最近比较活跃。现在新推出了3款特价KVM VPS,性价比高,香港CMI/洛杉矶GIA VPS,2核/2GB内存/20GB NVME/3.5TB流量/200Mbps...
小米3手机真实图片为你推荐
浏览器哪个好什么浏览器最好用?录音软件哪个好录音软件哪个好牡丹江教育云空间登录牡丹江教育云平台学生注册错了怎么办?willyunlee找一部关于摩托车的电影`dns服务器地址如何知道自己的IP地址和DNS服务器地址?360云网盘下载我有别人的360云盘里面的东西的链接,我要怎么下载他的这个东西?360云盘资源360云盘资源被和谐了怎么办360云安全中心为什么我的360云安全中心不能成功连接?360云盘下载速度慢怎么办360网盘下载速度好慢,我这几天下载怎么才30多KB。广东联通彩铃广东联通炫铃 怎么定制彩铃
免费虚拟主机申请 成都主机租用 拜登买域名批特朗普 香港bgp机房 vpsio paypal认证 创宇云 中国特价网 丹弗 七夕促销 怎么测试下载速度 太原网通测速平台 华为云服务登录 空间登入 lick 上海电信测速 服务器防火墙 wordpress中文主题 测试网速命令 后门 更多