基于Mumford-Shah模型的運動(dòng)目標檢測

摘要：使用Mumford-Shah模型進(jìn)行運行目標檢測以克服常規算法的缺點(diǎn)。利用改進(jìn)的水平集算法，使算法能夠快速收斂。為達到實(shí)時(shí)性的要求，利用多分辨方法進(jìn)一步提高算法的速度。使用改進(jìn)的區域生長(cháng)算法進(jìn)一步準確地檢測出運動(dòng)目標。關(guān)鍵詞：Mumford-Shah模型 運動(dòng)目標 水平集 區域生長(cháng) 魯棒性 運動(dòng)目標檢測在目標跟蹤、視頻監控和精確制導等領(lǐng)域有重要的應用。傳統的運動(dòng)目標檢測算法存在閾值確定困難、對噪聲較敏感等缺點(diǎn)?；顒?dòng)輪廓模型是解決靜止與運動(dòng)圖像分割和目標檢測問(wèn)題的一種有效方法，其主要缺點(diǎn)在于拓撲適應性較弱，即在演化過(guò)程中不能自適應地裂開(kāi)或合并。由Osher和Sethian提出的水平集方法解決了此問(wèn)題，將二維的閉合曲線(xiàn)嵌入一個(gè)三維的曲面，借助曲面的演化實(shí)現曲線(xiàn)的演化?；趲缀沃鲃?dòng)輪廓線(xiàn)模型的水平集算法僅利用圖像的邊緣信息，對邊緣模糊或存在離散狀邊緣的目標標難以得到理想的分割效果，而Chan-Vese提供的Mumford-Shah模型利用圖像的同質(zhì)區域的全局信息，可較好地分割出邊界模糊或離散邊界的目標。 本文對含有多個(gè)運動(dòng)目標的圖像序列相鄰幀圖像差建立Mumford-Shah模型，利用水平集算法求解此模型，利用改進(jìn)的偏微分方程和其數值解，使方程能夠快速收斂。為達到實(shí)時(shí)性的要求，利用多尺度方法提高算法的速度。采用改進(jìn)的區域生長(cháng)算法進(jìn)一步提高分割的準確性。本文圖像的背景相對固定，即攝像機相對固定。 1 水平集（Level set）方法 水平集方法主要是從界面傳播領(lǐng)域逐步發(fā)展起來(lái)的，是處理封閉運動(dòng)界面隨時(shí)間演化過(guò)程中幾何拓撲變化的有效計算工具。以水平集合函數φ所表達的曲線(xiàn)演化的最大特點(diǎn)是：即使隱含在φ中的零水平集曲線(xiàn)C發(fā)生了拓撲結構的變化，φ仍然保持一有效函數。 要使φ的演化與閉合曲線(xiàn)C的演化相關(guān)，φ的演化要滿(mǎn)足如下的Hamilton-Jacobi偏微分方程： аφ/аt=F|△φ| (1) 水平集方法中不需要顯式地求水平集函數φ，而是用圖像域初始閉合曲線(xiàn)C0生成的符號距離函數（signed distance function）代替，簡(jiǎn)記了SDF，即： φ(x,t=0)=%26;#177;d (2) x∈R2,d為圖像平面上的點(diǎn)x到曲線(xiàn)C的距離，若x在曲線(xiàn)的內部取正，則在曲線(xiàn)的外部取負。 2 運動(dòng)目標檢測算法 2．1 幀差法 定義灰度差分圖像D={d(x,y)}如下： d(x,y)=g(|I(x,y;t+1)-I(x,y;t)|) (3) I(x,y;t+1)和I(x,y;t)為圖像序列中相鄰幀圖像。函數g(x)定義為： 函數g的作用是調整圖像的對比度和均衡化圖像，ψ為調整系數。對分辨率較低場(chǎng)景中存在多個(gè)運動(dòng)目標，且其運動(dòng)速度和方向不相同的條件下，采用簡(jiǎn)單的閾值法分割效果較差。 2．2 基于Mumford-Shah模型的運動(dòng)區域檢測及水平集解 運動(dòng)區域檢測需要將灰度差分圖像D中的灰度一致域分離出來(lái)，即分割出運動(dòng)區域和靜止背景區域。傳統的幾何輪廓線(xiàn)模型侵害圖像的方法多采用活動(dòng)輪廓線(xiàn)所在位置的圖像局部信息，難于綜合圖像區域的全局信息，僅僅依靠進(jìn)化曲線(xiàn)C所在位置邊緣信息控制C的進(jìn)化。這種方法對邊緣模糊或離散邊緣的圖像分割效果不好。目前大多數視頻監控圖像的分辨率不高，采用傳統的幾何輪廓線(xiàn)的方法不能正確分割出圖像中的同質(zhì)區域。 Chan-Vese提出了一種簡(jiǎn)化Mumford-Shah的圖像分割模型，圖像I的定義域為Ω，C為Ω上的一閉合曲線(xiàn)，C將圖像I分割為目標和背景兩個(gè)同質(zhì)區域，定義如下的能量函數： 式（5）中c0、cb分別為圖像I在閉合曲線(xiàn)C內部和外部的灰度平均值。μ%26;#183;Length(C)、v%26;#183;Area(inside(C))為正則項，控制曲線(xiàn)的進(jìn)化。因此，最優(yōu)化圖像分割問(wèn)題轉化為求能量函數F（c0,cb,C）的最小值問(wèn)題?？梢钥闯?，只有C進(jìn)化到目標的邊界C0時(shí)，F（c1,c2,C）取最小值。 Chan-Vese以歐拉-拉格郎日法推導出水平集函數φ表達并滿(mǎn)足式（5）的偏微分方程： 式（6）的Ω為圖像函數和水平集函數的定義域。H（z）是Heaviside函數，δ(x)是Direc函數。此偏微分方程所涉及的圖像函數I的定義域為全圖數據，且方程中C0(φ)、cb(φ)也定義在全圖范圍內，Mumford-Shah的圖像分割模型的最大特點(diǎn)就是全局優(yōu)化。另外閉合曲線(xiàn)可以放置在圖像上任何位置。最后一個(gè)顯著(zhù)特點(diǎn)就是不依靠圖像的邊緣信息，即使圖像的邊緣模糊或離散，仍能得到較好的分割效果。 實(shí)驗表明，由于Dirac函數狹窄的定義范圍，限制了模型檢測圖像的全局性。Chan-Vese方法將此進(jìn)行了改進(jìn)，采用了正則化的Heaviside函數和Direc函數，定義為： 該函數保證了在圖像定義域范圍內，所有點(diǎn)的δε（z）值都是趨于零的正值。但當檢測曲線(xiàn)遠離檢測目標時(shí)，則δε（z）函數嚴重限制了對遠離進(jìn)化曲線(xiàn)C的目標的檢測，不能穩定地檢測出目標。 2．3 對Chan-Vese分割方法的改進(jìn) 為消除方程（6）中Dirac函數對檢測遠離進(jìn)化曲線(xiàn)C邊緣的抑制，將δ（φ）替換為|△φ|，使偏方程變?yōu)椋? 而co(φ)、cb(φ)按照方程（6）計算。Heaviside函數按照式（7）計算。由于|△φ|≈1，消除了Dirac函數對非零水平集的抑制，因此方程（8）比方程（6）有更好的全局優(yōu)化性能。 為保證式（8）解的穩定性，不采用Chan-Vese中的Jacobi方法，而采用有限差分方法： 改進(jìn)后的偏微分方程的求解過(guò)程如下： %26;#183;由φ0初始化φ0，n=0； %26;#183;由式（6）計算c0(φ)、cb(φ)； %26;#183;由式（10）計算曲率K； %26;#183;由式（9）計算φn+1； %26;#183;檢查解是否穩定，如不穩定，n=n+1，重復計算。 由于Mumford-Shah模型要計算整個(gè)圖像定義域的全局最優(yōu)解，不能使用常用的窄帶法求解水平集，需在整個(gè)定義域更新水平集函數，計算量較大，但由于方程（8）為全局優(yōu)化的偏微分方程，只需很少的幾次迭代，就可以得到理想的分割效果。 為減少算法的運行時(shí)間，采用多分辨方法，算法在保證檢測結果正確的前提下，可大大減少算法的計算時(shí)間。具體的計算時(shí)間比較見(jiàn)表1。表1 水平集算法改進(jìn)前后的計算代價(jià) 圖像序列等級（圖像大?。┍疚乃惴–han-Vese算法迭代次數收斂時(shí)間（s）迭代次數收斂時(shí)間（s）　0(320%26;#215;240)32.3100400129.3100圖像序列1(160%26;#215;120)30.6020706.8200　2（80%26;#215;60）30.1200100.43102．4 運動(dòng)變化分割基礎上的區域生長(cháng) 在圖像背景相對靜止的條件下，圖像變化區域C可表為： C（t,t+1）=O(t)∪O(t)∪O(t+1) (11) C(t-1,t)=O(t-1)∪O(t) 其中O（t）為t時(shí)刻屬于運動(dòng)目標的點(diǎn)集。則 C（t-1,t）∩C(t,t+1)=O(t)∪(O(t+1)∩O(t-1)) (12) 這意味著(zhù)兩個(gè)連續幀差圖像的交集更能代表運動(dòng)目標的準確位置，因此使用兩連續幀差圖像的交集作為運動(dòng)目標的初始位置，在此基礎上進(jìn)行區域生長(cháng)以更完整地檢測出運動(dòng)目標。 區域增長(cháng)法的基本思想是將具有相似性質(zhì)的像素集合起來(lái)構成區域。實(shí)驗中發(fā)現，采用檢測出的運動(dòng)區域為種子像素進(jìn)行區域生長(cháng)，難以確定相應的閾值，另外運動(dòng)區域的點(diǎn)可能包括背景像素點(diǎn)，這樣很易造成生長(cháng)錯誤。 本文采用一種基于圖像邊緣信息的區域生長(cháng)算法以克服上面的缺點(diǎn)，以?xún)蛇B續幀差圖像的交集作為運動(dòng)目標的種子像素點(diǎn)集，采用Sobel邊緣算子檢測種子像素周?chē)南袼攸c(diǎn)，如果為邊緣像素點(diǎn)，則標記為運動(dòng)目標點(diǎn)，將檢測到的邊緣像素點(diǎn)與運動(dòng)檢測結果進(jìn)行融合，能夠準確地檢測出運動(dòng)目標。 基于邊緣信息的區域生長(cháng)算法過(guò)程如下： （1）選擇水平集算法檢測出的運動(dòng)區域的運動(dòng)點(diǎn)為種子像素； （2）以該像素為中心檢測其8鄰域，采用Sobel算子進(jìn)行檢測，如為邊緣像素，則標記； （3）以新確定的像素為中心，反回步驟（2），檢查新像素的鄰域，直至遇到水平集算法確定的運動(dòng)點(diǎn)或超出沒(méi)定的區域范圍，返回（1），直到所有的種子像素都被檢查一遍，結束整個(gè)生長(cháng)過(guò)程。 在使用Sobel算子時(shí)采用各向同性的檢測模板，不采用固定閾值確定邊緣點(diǎn)，而是采用如下的判斷： G（x,y）/V(x,y)>Cof (13) G(x,y)為點(diǎn)（x,y）計算的邊緣強度，V（x,y）為點(diǎn)（xmy）圖像的灰度值。為檢測邊緣像素點(diǎn)較多，Cof取較小值。測試中Cof取0.1～0.3。圖13 實(shí)驗結果分析 采用實(shí)際的圖像序列進(jìn)行測試，檢測結果如圖1所示。利用方程（8）和方程（6）分別進(jìn)行計算，參數取值如下：λ0=λb=1,μ=0.2%26;#215;255,v=0,△t=0.1,Cof=0.2,μ取較大值以保證檢測較大目標。為加快SDF函數的初始化計算速度，初始化曲線(xiàn)選為圓，以圖像矩形中心為圓心，矩形的短邊的一半為半徑。編制程序列用不同等級分辨率的圖像進(jìn)行計算。計算機配置為：CPU為1.50GHz主頻率，內存為256MB，操作系統為Windows2000。 由表1和圖1的檢測結果可以看出，使用本文的算法對不同等級分辨率的圖像迭代三次均能正確地檢測圖像中的運動(dòng)車(chē)輛和行人。改進(jìn)后的算法與Chan-Vese算法相比大大減少了收斂的次數。由于改進(jìn)后的算法與Chan-Vese算法每次迭代的耗時(shí)基本相同，但改進(jìn)后的算法迭代數少，因此提高了算法的速度。改進(jìn)的算法對不同等級分辨率的圖像均能正確地檢測出運動(dòng)目標，這一點(diǎn)也驗證了算法的全局性能，圖像序列等級2的計算時(shí)間為0.1200s。Chan-Vese算法由于Dirac函數的影響，對初始化曲線(xiàn)的位置和曲線(xiàn)的長(cháng)度十分敏感，在實(shí)驗中也驗證了這一點(diǎn)。當初始化曲線(xiàn)遠離目標時(shí)，收斂次數大大增加，而改進(jìn)后的算法對初始化曲線(xiàn)的位置不敏感，算法可快速收斂。圖1（f）與（g）為水平集算法檢測結果，圖1（g）為兩者的交集，圖1（i）為采用Sobel邊緣檢測算子進(jìn)行區域生長(cháng)法的生長(cháng)結果，圖1（j）為融合兩者的分割結果，較好地檢測出運動(dòng)目標。對檢測出來(lái)的運動(dòng)目標的內部空洞可以直接填充以完整地檢測出運動(dòng)目標。 SDF函數的構造比較耗時(shí)，可利用當前幀水平集計算迭代收斂時(shí)的SDF函數，作為下一幀差圖像水平集計算時(shí)的SDF函數，避免了SDF函數的構造。由于連續幀間的目標運動(dòng)變化不大，因此在下一幀差圖像水平集計算的迭代次數更少，進(jìn)一步提高了算法的速度。在連續幀的測試中，迭代一到兩次均能檢測出運動(dòng)目標。 采用基于能量極小化的框架對含有多個(gè)運動(dòng)目標的圖像序列相鄰幀圖像差建立Mumford-Shah模型，為求解此模型，提出利用改進(jìn)的偏微分方程和水平集數值解法，使算法能夠快速收斂。為達到實(shí)時(shí)性的要求，利用多尺度方法進(jìn)一步提高算法的速度。使用改進(jìn)的區域生長(cháng)算法進(jìn)一步準確檢測出運動(dòng)目標。試驗結果表明，本文的算法能從較復雜的圖像序列中有效地檢測和提取出運動(dòng)目標并有較強的魯棒性。進(jìn)一步提高算法的魯棒性是今后的研究重點(diǎn)。