于海洋++金恩海++宋艾琪

摘要：簡單細(xì)胞對邊緣刺激具有選擇功能，大腦皮層中的感受野當(dāng)受到特定朝向或?qū)挾鹊拇碳r產(chǎn)生強(qiáng)烈的反應(yīng)，當(dāng)朝向逐漸偏離，反應(yīng)程度也呈現(xiàn)減小趨勢甚至完全消失。簡單細(xì)胞的感受野是如何模型化的？不同反應(yīng)程度的感受野是透過怎樣的機(jī)制整合得到精確的邊緣角度信息的？一直是研究中的難點(diǎn)。該論文以簡單細(xì)胞模型edge cell為例，針對現(xiàn)行模型在邊緣角度計算上存在的問題，給出兩種改良方案。利用高斯擬合大幅度降低其計算成本。利用高斯函數(shù)卷積疊加法模擬視覺系統(tǒng)的角度信息整合機(jī)制，解決了由少量的離散反應(yīng)值信息準(zhǔn)確的推算邊緣角度的問題。

關(guān)鍵詞：簡單細(xì)胞 邊緣檢測 高斯擬合 卷積 感受野

中圖分類號：TP394 文獻(xiàn)標(biāo)識碼：A 文章編號：1007-9416（2014）12-0099-02

1 緒論

輪廓線包含著物體的形狀信息，在圖像處理、模式識別領(lǐng)域，以構(gòu)成輪廓線的邊緣線段（Edge segments）為目標(biāo)的檢測具有重要意義。外界的光刺激投射到視網(wǎng)膜表面，通過視網(wǎng)膜轉(zhuǎn)換成的刺激信號經(jīng)由外側(cè)膝狀體中繼后傳導(dǎo)至大腦皮層V1野，再由V1野傳導(dǎo)至更高層視覺野進(jìn)行處理。V1野中的簡單細(xì)胞對邊緣刺激具有極強(qiáng)的選擇性，這種選擇包括對邊緣位置的選擇和對邊緣朝向的選擇。在過去的研究中作為簡單細(xì)胞的模型，廣泛認(rèn)可的有Gabor filters[1]， DOG filters [2]以及2G filters（laplacian-gaussian）[3]。這些模型中，采用兩個高斯函數(shù)差來描述外側(cè)膝狀體的數(shù)學(xué)模型DOG filters最貼合生理學(xué)實(shí)驗(yàn)數(shù)據(jù)，2G filters明確的定義了邊緣的位置，并通過唯一的參數(shù)σ明確了外側(cè)膝狀體細(xì)胞感受野與空間頻率的關(guān)聯(lián)。上述模型的邊緣檢測結(jié)果均為點(diǎn)排列，并不能呈現(xiàn)邊緣線段的結(jié)果。這是因?yàn)樯鲜瞿Ｐ途鶝]有對感受野的尺寸特別是長軸的尺寸做出規(guī)定，但感受野的大小和之間的距離對于輪廓線特征抽取是很重要的。近年提出的“Edge cell”簡單細(xì)胞模型彌補(bǔ)了這個缺點(diǎn)，該模型是由marr提出的簡單細(xì)胞模型改良而來。結(jié)合鏈?zhǔn)剿惴ǎ瑥倪吘壍囊稽c(diǎn)出發(fā)，一步一步跟蹤得到首尾相連的邊緣線段集合，其步長由空間周波數(shù)通過計算得到。

通過Edge cell模型所得到的邊緣線段，包含輪廓線局部的位置信息與朝向信息。其中位置信息基本可以忽略空間周波數(shù)變化帶來的影響，而朝向信息當(dāng)空間頻率較高時其精度受到影響。這種影響會延伸到后續(xù)利用朝向信息來檢測局部曲率的過程中。這是因?yàn)楝F(xiàn)行算法在對朝向進(jìn)行選擇時，只是簡單的通過比較細(xì)胞反應(yīng)值大小而造成的。本論文就朝向選擇問題提出兩種改良方法，并探討其各自的特點(diǎn)。

2 簡單細(xì)胞模型及角度推算方法

眾所周知簡單細(xì)胞具有方位選擇性，也就是說簡單細(xì)胞處于最匹配方位的時候反應(yīng)最為強(qiáng)烈，這個方位被稱為該簡單細(xì)胞的最適角度。Edge-cell的尺寸，如圖1所示， Edge-cell的高度對應(yīng)于2G*I負(fù)峰值和正峰值間的距離定義為2σ。Edge-cell的寬度定義為σ。這是出于這樣的一種考慮：外側(cè)膝狀體細(xì)胞的同心圓感受野最理想的排列方式是順序等間隔排列。像昆蟲復(fù)眼般的正六邊形排列是對平面覆蓋的一種有效方式。

圖1中所示，Pi，Qi（i=1，2，3）分別代表興奮中心型外側(cè)膝狀體細(xì)胞（on-center geniculate cell） 和抑制中心型外側(cè)膝狀體細(xì)胞（off-center geniculate cell）的中心坐標(biāo)。Cedge和θ表示edge cell感受野的中心和方位角。對于圖像I， edge cell的反應(yīng)函數(shù)定義如下：

RE（Cedge，θ）=2G（Q1）+2G（Q2）+2G（Q3）-2G（P1）-2G（P2）-2G（P3） （1）

此處2G（x）表示坐標(biāo)x位置的2G*I（laplacian-gaussian）值.

根據(jù)Marr提出的零交叉由負(fù)峰值、正峰值所決定的理論， RE（Cedge，θ）最大就意味著邊edge cell準(zhǔn)確的檢出了負(fù)峰值、正峰值，從而確定了位于負(fù)峰值和正峰值之間的零交叉線段的位置和朝向。從圖1中可以看到當(dāng)edge cell旋轉(zhuǎn)到和黑白境界線相同角度時對應(yīng)的RE（Cedge，θ）最大。從算式1中可以看出，edge cell的反應(yīng)函數(shù)值 RE（Cedge，θ） 是一個與位置和角度雙相關(guān)的值，在現(xiàn)行模型中最適角度的選擇是通過尋找RE（Cedge，θ）極大確定的，這里存在兩個問題。首先在RE（Cedge，θ）極大值附近的曲線變化均勻緩慢，從生理學(xué)角度出發(fā)，比較兩個近似值是難以實(shí)現(xiàn)的。其次，如果edge cell以1°為單位旋轉(zhuǎn)采樣，那就需要進(jìn)行360次算式1的計算，計算成本高昂。

3 高斯擬合算法與反應(yīng)值采樣

第一種解決方法，從減少計算成本角度考慮，利用少數(shù)離散的采樣值來推算最適方位。具體來說就是采用高斯曲線擬合法。edge cell的反應(yīng)值RE（Cedge，θ）曲線與高斯曲線非常相似（圖1），可以通過計算所擬合的高斯曲線來求取最適角度與最大值。高斯曲線方程如下：

其中A為高斯曲線的峰高，對應(yīng)RE（Cedge，θ）極大值。B對應(yīng)高斯曲線的中心位置，對應(yīng)最適角度。C為高斯曲線的幅寬。那么如果我們知道曲線上的3個點(diǎn)（θ1，R1）、（θ2，R2）、（θ3，R3）那么就可以通過求解聯(lián)立方程推算出最適角度B

具體做法是在0°～360°范圍內(nèi)等距離取采樣點(diǎn)n個，選取其中RE（Cedge，θ）值最大的3個采樣值來作為樣本，計算得出最適角度。此法可以大幅縮減計算量。

4 高斯卷積疊加法

高斯擬合算法利用等間隔離散的三個RE（Cedge，θ）最大樣本進(jìn)行計算。在視覺理論中，立方體模型中的朝向片不論反應(yīng)強(qiáng)弱均應(yīng)參與到方位角的計算中來，基于這種考慮RE（Cedge，θ）為正值的樣本都應(yīng)參加計算。通過采樣值集合fi[j]={RE1， RE2， … ， REn}中的正值樣本單獨(dú)與高斯函數(shù)（式（5））卷積，再對所有卷積結(jié)果求和得到高斯近似函數(shù)G（θ）的集合（式（6））。用來卷積的高斯函數(shù)區(qū)間定義在0°～180°，deg表示采樣間隔

用來實(shí)施卷積操作的高斯函數(shù)形狀由參數(shù)σ決定。當(dāng)選取適當(dāng)?shù)摩遥é?9.5）值卷積疊加后得到近似高斯曲線如圖3所示。如果σ（σ=5）值設(shè)定過小，卷積對象的高斯曲線半值幅過窄則不能得到預(yù)期的近似高斯曲線，也就無法準(zhǔn)確的推算最適角度。σ值的選取要參照RE（Cedge，θ）采樣間隔的設(shè)定，采樣間隔越大應(yīng)該選取更大的σ值來完成算法。

5 結(jié)語

本論文所提出的高斯擬合法與高斯卷積疊加法在邊緣檢測精度上均優(yōu)于現(xiàn)行模型?，F(xiàn)行模型的精度受制于采樣頻率，如果試圖通過提高采樣頻率來提高精度，還會帶來灰度差值的困難和計算成本的浪費(fèi)。高斯擬合法與高斯卷積疊加法對采樣間隔也有限制，當(dāng)采樣間隔超過45°后已不能保證至少三個樣本可以同時出現(xiàn)在高斯曲線上了，此兩種方法在采樣間隔相對較小的情況下結(jié)果更準(zhǔn)確。高斯擬合法在維持原模型檢測精度不變的基礎(chǔ)上將計算時間縮短至十分之一以下（表1），而高斯卷積疊加法利用離散的采樣值對邊緣角度進(jìn)行推算的方式也更貼近視覺處理的模式。

參考文獻(xiàn)

[1]Lee S.Image Representation Using 2D Gabor Wavelets [J]. IEEE Trans. Pattern Anal. Machine Intell.1996. 18， 959- 971.

[2]Daugmang D.Uncertainty relation for resolution in space， spatial frequency， and orientation optimized by two-dimensional visual cortical filters[J].Opt. Soc. Am.1985.2，1160-1169.

[3]Marr， D. and S. Ullman.， Directional selectivity and its use in early visual processing[J].C Proc. R. Soc. Lond. B，1981.211， 151-180.