estimationyahoo.cn

yahoo.cn  时间:2021-05-21  阅读:()
SVDBasedKalmanParticleFilterforRobustVisualTrackingXiaoqinZhang1,WeimingHu1,ZixiangZhao2,Yan-guoWang1,XiLi1,QingdiWei11NationalLaboratoryofPatternRecognition,InstituteofAutomation,Beijing,China{xqzhang,wmhu,ygwang,lixi,qdwei}@nlpr.
ia.
ac.
cn2BeijingUniversityofAeronautics&Astronauticszhaozx531@yahoo.
cnAbstractObjecttrackingisoneofthemostimportanttasksincomputervision.
Theunscentedparticlelteralgorithmhasbeenextensivelyusedtotacklethisproblemandachievedagreatsuccess,becauseitusestheUKF(un-scentedKalmanlter)togenerateasophisticatedpro-posaldistributionswhichincorporatesthenewestob-servationsintothestatetransitiondistributionandthusovercomesthesampleimpoverishmentproblemsufferedbytheparticlelter.
However,UKFoftenencounterstheill-conditionedproblemwhensolvingthesquarerootofthecovariancematrixinpractice.
Inthispaper,weproposeanovelKalmanparticlelterbasedonSVD(singularvaluedecomposition),andapplyitforvisualtracking.
Experimentalresultsdemonstratethat,com-paredwiththeparticlelterandtheunscentedparticlelter,theproposedalgorithmismorerobustintrackingperformance.
1.
IntroductionObjecttrackinghasreceivedsignicantattentionduetoitscrucialvalueinvisualapplicationsinclud-ingsurveillance,human-computerinteraction,intelli-genttransportation,augmentedrealityandvideocom-pression.
Theparticlelter[1,2]hasbeenextensivelystud-iedinthetrackingliteratureduetoitseffectivenessandexibility.
FromaBayesianview,particlelterisessen-tiallyasequentialMonteCarloapproachtosolvethere-cursiveBayesianlteringproblem,whichcombinestheMonteCarlosamplingtechniqueswithBayesianinfer-ence.
ItrelaxesthelinearityandGaussianityconstraintsoftheKalmanlterandprovidesatractablesolutiontonon-linearandnon-Gaussiansystems.
Thebasicideaofparticlelteristouseanumberofindependentran-domvariablescalledparticles,sampleddirectlyfromaproposaldistribution,torepresenttheposteriorprob-ability,andupdatetheposteriorbyinvolvingthenewobservations.
Althoughithasachievedaconsiderablesuccessinthetrackingliterature,itisfacedwithafatalproblem-sampleimpoverishmentduetoits'suboptimalsampling'mechanism.
Fortheconventionalparticlel-ter,theparticlesaredirectlysampledfromstatetransi-tiondistribution.
However,itisnotthe'optimal'pro-posalsamplingdistribution.
Whenthestatetransitiondistributionliesinthetailoftheobservationlikelihooddistribution,theweightsofmostparticlesarelow,lead-ingtothepoorperformanceinpractice.
Muchefforthasbeenexpendedtoovercomethisproblemandimprovetheperformanceofparticlel-terinrecentyears[3,4,5,6,7,8].
Amongthem,theunscentedparticlelter[4]isthesuccessfulone.
Intheunscentedparticlelter,theUKFbasedproposaldistributionisintroducedasfollows.
Firstly,asetthesigmasamplesaregeneratedbyUT(unscentedtrans-formation)withcorrespondingweights,andthenarepropagatedthroughthestatetransitionmodel,nallytheweightedmeanandcovariancearefurthercalcu-latedtoformabetterproposaldistribution.
ComparedwiththeEKF(extendedKalmanlter)whichapproxi-matestotherst-orderaccuracyfornon-Gaussiandata,theestimationaccuracyofUKFisimprovedtoatleastsecond-order.
However,UKFoftenencounterstheill-conditionedproblemwhensolvingthesquarerootofthecovariancematrixinpractice.
Toovercomethisproblem,weproposeaSVDbasedKalmanparticlel-ter,wherethesigmasamplesaregeneratedbySVDoftheeigen-covariancematrix.
Whilemaintainingthesamecomputationalcomplexity,theproposedtrackingalgorithmperformsquiterobustlyintrackingperfor-mance.
Thefollowingpaperisarrangedasfollows.
Section2presentstheunscentedparticlelteringframework.
ThedetailoftheproposedSVDbasedKalmanlterisdescribedinSection3.
Section4introducestheincre-mentalsubspaceleaningbasedappearancemodel.
Ex-perimentalresultsareshowninSection5,andSection6isdevotedtoconclusion.
2.
UnscentedParticleFilteringFrameworkTomakethispaperself-contained,werstbrieyreviewtheparticlelteranditsmajorlimitation,andthenpresenttheunscentedparticlelterindetail.
2.
1.
ParticleFilterParticlelter[2]isanonlineBayesianinferencepro-cessforestimatingtheunknownstatextattimetfromasequentialobservationsy1:tperturbedbynoises.
Adynamicstate-spaceformemployedintheBayesianin-ferenceframeworkisshownasfollows,xt=f(xt1,t)p(xt|xt1)(1)yt=h(xt,νt)p(yt|xt)(2)wherext,ytrepresentsystemstateandobservation,t,νtarethesystemnoiseandobservationnoise.
f(.
,.
)andh(.
,.
)arethestatetransitionandobservationmod-els,whichcharacterizethestatetransitiondistributionp(xt|xt1)andtheobservationdistributionp(yt|xt)re-spectively.
Thekeyideaofparticlelteristoapprox-imatetheposteriorprobabilitydistributionp(xt|y1:t)byasetofweightedsamples{xit,wit}Ni=1,whicharesampledfromaproposaldistributionq(·),i.
e.
xitq(xt|xit1,y1:t),(i=1,N),andtheneachparti-cle'sweightissettowit∝p(yt|xit)p(xit|xit1)q(xt|xit1,y1:t)(3)Finally,theposteriorprobabilitydistributionisapprox-imatedasp(xt|y1:t)=Ni=1witδ(xtxit),whereδ(·)istheDiracfunction.
Doucetetal.
[9]provethatthe'optimal'proposaldistributionisp(xt|xit1,yt)inthesenseofminimizingthevarianceoftheimportanceweights.
Sotheques-tionis,howtoincorporatethecurrentobservationytintothetransitionmodelp(xt|xt1)toformaneffec-tiveproposaldistribution.
2.
2.
UnscentedParticleFilterInordertoutilizethecurrentobservations,Freitasetal.
[4]proposeahigh-performanceunscentedparti-clelter(UPF)byusingUKFtogeneratetheproposaldistribution.
Intheimplementation,thestatespaceisexpandedas:xat1=[xTt1Tt1νTt1],whosedimensionandco-variancematrixareNa=Nx+N+NνandPat1respectively.
Considerthenonlineartrackingproblemmodeledbythestate-spaceequations(1)and(2),thepseudo-codeoftheunscentedKalmanlterispresentedasfollows.
1.
Calculate2Nasigmapointsasin[4]X(i)a0,t1=x(i)at1X(i)aj,t1=[x(i)at1x(i)at1±(na+λ)Pat1j]W(m)0=λNa+λ,W(c)0=λNa+λ+(1α2+β)W(m)j=W(c)j=12(Na+λ),λ=α2(Na+κ)Naj=1,2Na2.
Timeupdate:X(i)xj,t|t1=f(X(i)xj,t1,X(i)j,t1),x(i)t|t1=2Naj=0W(m)jX(i)xj,t|t1P(i)t|t1=2Naj=0W(c)j[X(i)xj,t|t1x(i)t|t1][X(i)xj,t|t1x(i)t|t1]TY(i)j,t|t1=h(X(i)xj,t1,X(i)νj,t1),y(i)t|t1=2Naj=0W(m)jY(i)j,t|t13.
Measurementupdate:Pyt,yt=2Naj=0W(c)j[Y(i)j,t|t1y(i)t|t1][Y(i)j,t|t1y(i)t|t1]TPxt,yt=2Naj=0W(c)j[X(i)xj,t|t1x(i)t|t1][Y(i)j,t|t1y(i)t|t1]TKt=Pxt,ytP1yt,yt,x(i)t=x(i)t|t1+Kt(yty(i)t|t1)P(i)t=P(i)t|t1KtPyt,ytKTtAsaresult,theproposaldistributionisobtainedasq(xit|xit1,y1:t)=N(x(i)t,P(i)t),andtheunscentedpar-ticlelterisanaturalcombinationoftheUKFproposaldistributionandtraditionalparticlelteraspresentedinSection2.
1.
3.
SVDBasedKalmanFilterHowever,UKFoftenencounterstheill-conditionedproblemwhensolvingthesquarerootofthecovariancematrixPat1inpractice.
Therefore,weproposeanSVDbasedKalmanltertoovercomethisproblem.
Togiveaclearview,theowchartoftheSVDbasedKalmanlterframeworkisschematicallyshowninFig.
1.
TheSVDbasedKFsharesaclosespirittoUKF,rstly,themeanstateandeigen-covariancematrixofthesigmasamplesattimet1arecalculated,andweapplytheSVDtotheeigen-covariancematrixtoobtain/86:22:=95-52;+82;35/0,.
=>82;35.
5B2:=95Figure1.
OverviewoftheSVDbasedKalmanlteritseigenvectors.
Thentheobtainedmeanandeigen-vectorsarecombinedtogeneratenewsigmasamples.
Finally,thenewsamplesarelteredbythestandardKalmanlter.
ThethedetailSVD-basedKalmanl-terprocessispresentedasfollows.
1.
ComputetheSVDoftheeigen-pointcovariancema-trixPat1=Ut1St1VTt12.
Calculatenewsigmasamples:X(i)a0,t1=x(i)at1X(i)aj,t1=[x(i)at1x(i)at1±ρUj,t1√sj,t1]whereUj,t1,sj,t1arethejtheigenvectorandeigen-valuerespectively,andρisthescaleparameter.
3.
ThefollowingstepisthesameasthestandardKal-manltering.
ThebasicmotivationbehindSVD-KFisthatthecovariancematrixcanbecharacterizedbyitseigen-vectors,andSVDismorenumericallyrobustthanCholeskyfactorizationintheunscentedtransformation.
4.
IncrementalSubspaceLeaningBasedAppearanceModelInourpaper,weadoptasubspacebasedappearancemodel[10]forobservationevaluation,whichmodelstheappearanceofanobjectbyincrementallylearningalow-ordereigenspacerepresentation.
ObservationLikelihood:Asshownin[10],giventhelearnedthesubspaceUandthenewobservationyt,theobservationlikelihoodisbasedonthereconstructioner-roroftheobservationyiintheobjectsubspace,whichisdenedasfollows.
RE=||ytUUTyt||2(4)Asaresult,theobservationlikelihoodisnaturallyformedasp(yt|xt)=exp(RE)(5)IncrementallySubspaceLearning:GiventheSVDofthepreviousappearancedataA={I1,It},i.
e.
A=UΣVT,whereeachcolumnIiistheobservationoftheobjectintheithframe.
Aftertrackingkframes,wehaveobtainedknewestobservationsoftheobjectE={It+1,It+k},theR-SVDalgorithm[11]ef-cientlycomputestheSVDofthematrixA=(A|E)=UΣVTbasedontheSVDofAasfollows:1.
ApplyQRdecompositiontoandgetorthonormalba-sisEofE,andU=(U|E).
2.
LetV=V00IkwhereIkisak*kidentitymatrix.
Itfollowsthen,Σ=UTAV=UTE(A|E)V00Ik=UTAVUTEETAVETE=ΣUTE0ETE3.
ComputetheSVDofΣ=UΣVTandtheSVDofAisA=U(UΣVT)VT=(UU)Σ(VTVT)Inthisway,theR-SVDalgorithmcomputestheneweigenbasisefciently.
5.
ExperimentalResultsInourexperiment,theobjectisinitializedman-uallyandafnetransformationsisconsideredonly.
Specically,themotionischaracterizedbys=(tx,ty,a1,a2,a3,a4)where{tx,ty}denotethe2-Dtranslationparametersand{a1,a2,a3,a4}aredeforma-tionparameters.
Eachcandidateimageisrectiedtoa20*20patch,andthefeatureisa400-dimensionvectorwithzero-mean-unit-variancenormalization.
Inordertodemonstratetheeffectivenessofourap-proach,weconductacomparisonexperimentamongtheSVDbasedKPF(Kalmanparticlelter),astandardPF(particlelter)1andUPF[4]onavideowithmanu-allylabeledgroundtruth.
TheDavidsequence2issampledalternatelytoformarapidmotiontestingsequence.
Inourimplementa-tion,theparametersaresetto{N=200,var()=1Here,aGaussiantransitiondistributionxtN(xt1,Σ)istakenastheproposaldistribution2WeacknowledgetotheauthorofthesourcedataavailableattheURL:http://www.
cs.
toronto.
edu/dross/ivt/TrackingMethodFramesTrackedMSE(bypixels)PF16/6126.
9481UPF61/617.
1875SVDbasedKPF61/613.
9868Table1.
QuantitativeresultsofSVDbasedKPFtrackeranditscomparisonwithPFtrackerandUPFtracker[52,52,0.
012,0.
022,0.
0022,0.
0012]}correspondingtothenumberofparticlesandthecovariancematrixofthetransitiondistributionrespectively.
AsshownintherstcolumnofFig.
2,theparticlelterbasedtrackerfailstotracktheobjectatframe31,becausetheparticlesaresampledfromthetransitiondistributiontocatchtheob-jectmotion.
Whentheobjecthasrapidandarbitrarymotion,theparticlesdrawnfromthisdistributiondonotcoverasignicantregionofthelikelihood,andthustheweightsofmostparticlesarelow,leadingtothetrackingfailure.
Moreparticlesandanenlargementforthediag-onalelementsofthecovariancematrixwouldimproveitsperformance,butthisstrategyinvolvesmorenoisesandaheavycomputationalload.
ThesecondcolumnofFig.
2showsthetrackingperformanceoftheunscentedparticlelter,fromwhichwenoticethatthetrackerfollowstheobjectthroughoutthesequence.
However,Choleskyfactorizationisnotnumericallyrobustandof-tenencounterstheill-conditionedproblem,therebyre-sultingtotheinaccuratelocalizationandsize.
Incom-parison,ourmethodachievesthemoreaccurateresults,becausethecovariancematrixisfullycharacterizedbyitseigenvectors,andSVDismorenumericallyrobustthanCholeskyfactorization.
Meanwhile,wehavecon-ductedaquantitativeevaluationofthesealgorithms,andhaveacomparisoninthefollowingaspects:framesofsuccessfultracking,MSE(meansquareerror)betweentheestimatedpositionandthelabeledgroundtruth.
Intable1,itisclearthatthePFtrackerfailsatframe31whiletheUPFandSVDbasedKPFtrackerssucceedintrackingthroughoutthesequence.
Additionally,theSVDbasedKPFtrackeroutperformstheUPFtrackerintermofaccuracy.
6.
ConclusionThispaperpresentsanSVDbasedKalmanparticlelterforvisualtracking.
Inouralgorithm,asetofsigmasamplesaregeneratedbySVDofthecovariancema-trix,andthenthesesigmapointsarepropagatedbythestandardKalmanltertogenerateasophisticatedpro-posaldistribution.
Theobtainedproposaldistributionisincorporatedintotheparticleltertoformarobusttrackingalgorithm.
Experimentalresultsdemonstratetheeffectivenessandpromisingofourapproach.
Figure2.
Thetrackingresults(rstcol-umn:PF,secondcolumn:UPF,thirdcol-umn:SVDbasedKPF)7.
AcknowledgmentThisworkispartlysupportedbyNSFC(GrantNo.
60672040,60705003)andtheNational863High-TechR&DProgramofChina(GrantNo.
2006AA01Z453).
References[1]M.
IsardandA.
Blake,"Condensation:conditionaldensitypropagationforvisualtracking",IJCV.
,29(1):5-28,1998.
[2]M.
Arulampalam,S.
Maskell,N.
Gordon,andT.
Clapp,"Atutorialonparticlesltersforonlinenonlinear/non-gaussianbayesiantracking",IEEETrans.
onSignalProcessing,50(2):174-188,2002.
[3]J.
Deutscher,A.
Blake,I.
Reid,"Articulatedbodymotioncap-turebyannealedparticleltering",Proc.
CVPR'00,pp.
126-133,2000.
[4]R.
Merwe,A.
Doucet,N.
Freitas,andE.
Wan,"Theunscentedparticlelter",TechnicalReportCUED/F-INFENG/TR380,CambridgeUniversityEngineeringDepartment,2000.
[5]Y.
RuiandY.
Chen,"Betterproposaldistributions:objecttrackingusingunscentedparticlelter",Proc.
CVPR'01,pp.
786-793,2001.
[6]Y.
Wu,T.
Huang,"RobustVisualTrackingbyIntegratingMul-tipleCuesBasedonCo-InferenceLearning",IJCV.
,58(1):55-71,2004.
[7]B.
Han,Y.
Zhu,D.
ComaniciuandL.
Davis,"Kernel-BasedBayesianFilteringforObjectTracking",Proc.
CVPR'05,pp.
227-234,2005.
[8]X.
Zhang,W.
Hu,G.
Luo,andS.
Maybank,"Kernel-BayesianFrameworkforObjectTracking",Proc.
ACCV'07,pp.
821-831,2007.
[9]A.
Doucet,S.
Godsill,andC.
Andrieu,"OnsequentialMonteCarlosamplingmethodsforBayesianltering",StatisticsandComputing,10(3):197-208,2000.
[10]J.
Lim,D.
Ross,R.
S.
Lin,andM.
H.
Yang,"Incrementallearn-ingforvisualtracking",InAdvancesinNeuralInformationProcessingSystems,pp.
793-800,2004,TheMITPress.
[11]A.
LevyandM.
Lindenbaum,"SequentialKarhunen-Loeveba-sisextractionanditsapplicationtoimages",IEEETrans.
onImageProcessing,9(8):1371-1374,2000.

数脉科技香港自营,10Mbps CN2物理机420元/月

数脉科技怎么样?数脉科技品牌创办于2019,由一家从2012年开始从事idc行业的商家创办,目前主营产品是香港服务器,线路有阿里云线路和自营CN2线路,均为中国大陆直连带宽,适合建站及运行各种负载较高的项目,同时支持人民币、台币、美元等结算,提供支付宝、微信、PayPal付款方式。本次数脉科技给发来了新的7月促销活动,CN2+BGP线路的香港服务器,带宽10m起,配置E3-16G-30M-3IP,...

舍利云30元/月起;美国CERA云服务器,原生ip,低至28元/月起

目前舍利云服务器的主要特色是适合seo和建站,性价比方面非常不错,舍利云的产品以BGP线路速度优质稳定而著称,对于产品的线路和带宽有着极其严格的讲究,这主要表现在其对母鸡的超售有严格的管控,与此同时舍利云也尽心尽力为用户提供完美服务。目前,香港cn2云服务器,5M/10M带宽,价格低至30元/月,可试用1天;;美国cera云服务器,原生ip,低至28元/月起。一、香港CN2云服务器香港CN2精品线...

炭云188元/年,上海CN2 VPS/2核/384MB内存/8GB空间/800GB流量/77Mbps端口/共享IP

炭云怎么样?炭云(之前的碳云),国人商家,正规公司(哈尔滨桓林信息技术有限公司),主机之家测评介绍过多次。现在上海CN2共享IP的VPS有一款特价,上海cn2 vps,2核/384MB内存/8GB空间/800GB流量/77Mbps端口/共享IP/Hyper-v,188元/年,特别适合电信网络。有需要的可以关注一下。点击进入:炭云官方网站地址炭云vps套餐:套餐cpu内存硬盘流量/带宽ip价格购买上...

yahoo.cn为你推荐
参数winrar5neaseset支持ipad重庆网通重庆联通宽带ipadwifiipad wifi信号差怎么办ipadwifiIPAD连上了WIFI,但是无法上网,急!!x-router设置路由器是我的上网设置是x怎么弄联通版iphone4s苹果4s怎么分移动版联通版电信版?苹果5.1完美越狱iOS5.1.1完美越狱教程fastreport2.5pm10和pm2.5的区别,pm10和pm2.5哪个危害大
vps vps安全设置 80vps 42u标准机柜尺寸 mysql主机 日本空间 qq数据库 福建天翼加速 qingyun 阿里云浏览器 工信部icp备案号 200g硬盘 老左来了 股票老左 北京双线 徐正曦 亚马逊香港官网 1g空间 国外免费asp空间 免费phpmysql空间 更多