An Uncertainty-Driven Hybrid of Intensity-Based and Feature-Based Registration with Applica

  • 格式:pdf
  • 大小:296.17 KB
  • 文档页数:8

AnUncertainty-DrivenHybridof

Intensity-BasedandFeature-BasedRegistration

withApplicationtoRetinalandLungCT

Images󰀂

CharlesV.Stewart,Ying-LinLee,andChia-LingTsai

RensselaerPolytechnicInstituteTroy,NY12180-3590USA

Abstract.Anewhybridoffeature-basedandintensity-basedregistra-

tionispresented.Thealgorithmreflectsanewunderstandingoftheroleofalignmenterrorinthegenerationofregistrationconstraints.Thisleadstoaniterativeprocesswheredistinctiveimagelocationsfromthemovingimagearematchedagainsttheintensitystructureofthefixedimage.Thesearchrangeofthismatchingprocessiscontrolledbyboththeuncertaintyinthecurrenttransformationestimateandtheproper-tiesoftheimagelocationstobematched.TheresultinghybridalgorithmisappliedtoretinalimageregistrationbyincorporatingitasthemainestimationenginewithinourrecentlypublishedDual-BootstrapICPal-gorithm.Thehybridalgorithmisusedtoalignserialand4dCTimagesofthelungusingaB-splinebaseddeformationmodel.

1Introduction

Feature-basedandintensity-basedregistrationalgorithmsdifferconsiderablyin

theimage-basedfeaturesthatdrivethealignmentprocess[4].Intensity-basedtechniquesuseallimagepixelsanddonotrequireexplicitfeatureextraction.

Theytendtobemorestablearoundminimaoftheobjectivefunctionbecause

theydon’trelyonuncertainfeaturelocationsoroncorrespondenceswhichmay

fluctuatewithslightchangesinthetransformation.Ontheotherhand,feature-

basedtechniquestendtobefaster,haveawidercapturerange,andallowalign-

menttobefocusedononlyselectedsubsetsoftheimagedata.Thesestrengths

andweaknessesarestudiedexperimentallyin[9].

Agrowingsetofpapershasbeguntoaddresstheissueofcombiningfeature-

basedandintensity-basedregistration.Asexamples,in[5]registrationisdriven

bothbyintensitysimilarityerrorandbyerrorsinthepositionsofmatched

features.InFeldmaretal.[7]intensityistreatedasa4thdimensionforICP

matching,whileSharpetal.[13]combinesICPwithinvariantfeatures.Inamuch

󰀁ThesupportoftheGEGlobalResearchCenterandoftheNationalScienceFoun-dationthroughtheCenterforSubsurfaceSensingandImagingSystemsisgratefullyacknowledged.

C.Barillot,D.R.Haynor,andP.Hellier(Eds.):MICCAI2004,LNCS3216,pp.870–877,2004.c󰀂Springer-VerlagBerlinHeidelberg2004AnUncertainty-DrivenHybrid871

differentapproach,Aylwardetal.[1]locatetubularstructuresinoneimageand

thenalignimagesusinggradient-descentofa“tubularness”measureevaluated

atthetransformedlocationsofthesestructuresintheotherimage.Ourselin,et

al.[10]useblock-matchingofintensitiestogetherwithrobustregression.Shum

andSzeliski[15]useblockmatchingofintensitiestorefinevideomosaics.The

PASHAalgorithm[4]combinescorrespondenceanddeformationprocessesina

globalobjectivefunction.IntheHAMMERalgorithm[14]registrationisdriven

bythealignmentoffeaturepositionsbasedonahierarchicaldescriptionofthe

surroundingintensitydistribution.

Thispaperpresentsanewhybridofintensity-basedandfeature-basedregis-

trationandappliestheresultingalgorithmintwocontexts:aligninglow-overlap

retinaimagesandspline-basedalignmentoflungCTvolumes.Featurepoints

foundinthemovingimagearematchedagainsttheintensitystructureofthe

fixedimage.Importantly,thesearchrangeforthematchisdictatedbyboth

thepropertiesofthefeatureandtheuncertaintyinthecurrenttransforma-

tionestimate.Thisdiffersfromcommonintensity-basedtechniqueswherecon-

straintsaredrivenbylocalchangesinthesimilaritymeasure.Italsodiffersfrom

correspondence-basedmethodslikeICP[2]wherematchingispurelyanearest-

pointsearch.

Thiscorealgorithmisbuiltintotwodifferentoverallregistrationalgorithms.

ThefirstisanextensionofourrecentDual-BootstrapICP[17]algorithmfor

aligningretinalimages.Replacingthefeature-to-featureICPalgorithmwith

thenewfeature-to-intensitysimilaritymatchingincreasestheeffectivenessfor

extremelydifferentcases.Thesecondoverallalgorithmisanewtechniquefor

non-rigid,B-splinealignmentoflungCTvolumes.Smallregionswithsufficient

intensityvariationinthemovingimagearematchedagainstthefixedimage.

ThesematchescontroltheestimateofhierarchicalB-splinedeformations.The

algorithmisappliedtoeffectivelyalignserialCTvolumes.

2TheRoleofUncertaintyinRegistration

Wemotivatethenewhybridalgorithmbyconsideringtheimportanceofun-

certaintyinregistration.LetS(p,q)measuretheregion-basedsimilarityerror

betweenmovingimageImatlocationpandfixed(target)imageIfatlocation

q.LetT(p;θ)bethetransformationmappingpfromImontoIfbasedonthe

(tobeestimated)parametervectorθ.Finally,letPbeasetoflocationsinImwheretransformationconstraintsareapplied.

Intensity-basedalgorithmssearchforthetransformationminimizingtheag-

gregatesimilarityerror:

E(θ)=󰀂

pi∈PS(pi,T(pi,θ)),(1)

Regularizingconstraintsmaybeplacedonθ,butwewillignorethesefornow.E(θ)ismostfrequentlyminimizedthroughagradient-descenttechnique[8].

Theestimate,ˆθ,remainsuncertainthroughouttheminimizationprocess.This