尚志剛，馮平艷，劉新玉，牛曉可，萬 紅

(鄭州大學(xué)電氣工程學(xué)院，河南鄭州450001)

0 引言

朝向調(diào)諧特性是初級視覺皮層(Primary Visual Cortex，又稱V1區(qū))神經(jīng)元感受野的主要特性之一［1］.朝向調(diào)諧特性的傳統(tǒng)研究方法是通過使用神經(jīng)元spike發(fā)放率作為其朝向調(diào)諧指標(biāo).然而，由于在spike檢測和分類過程中經(jīng)常會發(fā)生假陽性和假陰性等問題［2］，嚴(yán)重影響了后續(xù)分析的結(jié)果.

局部場電位(Local Field Potential，LFP)是電極尖端附近神經(jīng)元響應(yīng)信號的總和，頻帶主要集中在1～200 Hz［3］.LFP低頻部分主要反映狀態(tài)信息，高頻部分尤其是γ頻帶(30～90 Hz)，主要反映了與視覺刺激相關(guān)的信息.Juergens等［4］在研究初級視覺皮層局部場電位信號和EEG信號的關(guān)系時發(fā)現(xiàn)γ頻帶的能量反映了視覺刺激的信息.Frien等［5］研究指出LFPγ頻帶的功率會隨著刺激朝向的改變而產(chǎn)生不同的調(diào)制作用;Siegel和K?nig［6］研究也發(fā)現(xiàn)γ頻帶的能量會隨著朝向的改變而改變;Henrie等［7］用運動光柵刺激研究發(fā)現(xiàn)局部場電位γ頻帶的能量在不同的朝向下有選擇性的加強.但文獻(xiàn)中沒有給出兩種信號對朝向調(diào)諧特性研究的定量比較.

以大鼠為研究對象，將微電極陣列植入大鼠V1區(qū)，以不同朝向的全屏光柵為視覺刺激，同時記錄LFP和spike信號，研究γ頻帶能量對光柵朝向的調(diào)制作用.結(jié)果表明:LFPγ頻帶能量和spike發(fā)放率獲得的神經(jīng)元偏愛朝向一致率達(dá)到了67.4%，朝向調(diào)諧一致率達(dá)到了87.2%.并且仿真研究結(jié)果發(fā)現(xiàn)，基于LFPγ頻帶能量的朝向調(diào)諧指數(shù)比spike發(fā)放率具有更強的抗噪能力.

1 原理與方法

1.1 朝向調(diào)諧指數(shù)的計算

不同朝向光柵刺激LFPγ頻帶能量朝向調(diào)諧指數(shù) ρLFP

［8］計算如下:

若ρLFP＞0.7，則該神經(jīng)元具有偏愛朝向.其中Pp為最優(yōu)朝向γ頻帶能量，Po為與最優(yōu)朝向正交的兩個朝向γ頻帶能量P1、P2的均值，即:Po=(P1+P2)/2.

不同朝向光柵刺激spike發(fā)放率朝向調(diào)諧指數(shù)ρs可類似計算:若ρs＞0.7，則該神經(jīng)元具有偏愛朝向.其中Pp為最優(yōu)朝向spike發(fā)放率;Po為與最優(yōu)朝向正交的兩個朝向 spike發(fā)放率 P3、P4的均值，即:Po=(P3+P4)/2.

偏愛朝向一致率與朝向調(diào)諧一致率定義為:設(shè)總神經(jīng)元個數(shù)為M，其中兩種信號都有偏愛朝向并且偏愛朝向一致的神經(jīng)元個數(shù)為N1，兩種都沒有偏愛朝向的神經(jīng)元個數(shù)為N2，偏愛朝向一致率為λ1=N1/M，朝向調(diào)諧一致率為λ2=(N1+N2)/M.

1.2 LFPγ頻帶的提取和功率譜估計

首先對原始數(shù)據(jù)使用通帶為0～250 Hz的二階Butterworth帶通濾波器對信號進(jìn)行濾波，然后對濾波后的信號重采樣，得到采樣頻率為2 kHz的局部場電位信號.

由于LFP為非平穩(wěn)非高斯信號，在檢測的過程中往往會出現(xiàn)一些短時瞬態(tài)脈沖，這些瞬態(tài)脈沖信號可能包含有刺激的信息.傳統(tǒng)的瞬態(tài)脈沖檢測方法中，匹配濾波［9］是最常用的一種，但是LFP非平穩(wěn)性極強，匹配模板需要不斷更新.因此，采用小波變換的方法，用逐漸精細(xì)的時間分辨率逼近信號的快變部分，用逐漸精細(xì)的頻率分辨率逼近信號的慢變部分.

小波基函數(shù)ψ(t)不具唯一性，不同的小波基函數(shù)分析得到的結(jié)果相差甚遠(yuǎn).筆者基于LFP信號的特點，選擇與LFP信號相似的‘db5’作為小波的基函數(shù).

小波分解公式為:

小波重構(gòu)公式為:

其中，小波基函數(shù)為:

式中:a0＞0;j，k∈Z;f(t)為原信號;a-j0為尺度因子，kτ0為平移因子.

LFP信號采樣頻率為2 kHz，為了精確獲得γ頻帶信號，按照式(3)的小波分解算法對原始信號進(jìn)行5層小波分解，經(jīng)過5層分解后，選擇其中62.5～125 Hz頻帶再分解后獲得62.5～93.75 Hz頻帶，利用此頻帶與第一次分解得到的31.25～62.5 Hz頻帶進(jìn)行重建，得到γ頻帶信號(30～90 Hz)，重構(gòu)出的γ頻帶信號如圖1所示.小波重構(gòu)出的γ頻帶較好地抑制了噪聲并且保留了響應(yīng)的瞬態(tài)特征.

筆者選用Welch功率譜估計對LFP進(jìn)行功率估計.Welch功率譜估計的關(guān)鍵是選取合適的窗和窗長以適應(yīng)不同分辨率的要求.因此筆者選擇分辨率和方差水平都較好的Hamming窗進(jìn)行譜估計.選用的Hamming窗分辨率為8 Hz，局部場電位頻帶信號帶寬約為60 Hz，采樣頻率2 kHz.對比研究發(fā)現(xiàn):分辨率8 Hz較好地抑制了噪聲，并且可以準(zhǔn)確地提取出頻帶的能量，其對應(yīng)的窗長N=500，其功率譜估計如圖2所示.由圖2可以看出，該分辨率下譜估計方差和噪聲水平都較好，頻帶能量主要集中在35～60 Hz頻帶內(nèi).

2 實驗分析

2.1 實驗數(shù)據(jù)采集

實驗采用年齡11周左右，體重200～250 g的成年雄性Long Evans(LE)大鼠.術(shù)前實驗動物用水合氯醛腹腔注射麻醉后，固定于Stoelting腦立體定位儀上，參考大鼠腦立體定位圖譜在V1區(qū)植入2×8鉑銥合金Microprobes微電極絲陣列.植入電極1.5 h后，對大鼠進(jìn)行12個朝向(0～330°，間隔30°)的全屏光柵視覺刺激，并用 CerebusTM數(shù)據(jù)采集系統(tǒng)記錄大鼠V1區(qū)初級視覺皮層信號.

2.2 LFPγ頻帶能量分析

為了研究γ頻帶能量是否能反應(yīng)刺激的信息，筆者首先分析了有光柵刺激和無光柵刺激γ頻帶的能量分布.圖3為其中一個通道的分析結(jié)果.由圖3可以看出:γ頻帶能量主要集中在35～60 Hz，并且有刺激時γ頻帶能量比無刺激時有較為明顯的變化，說明了γ頻帶的能量編碼了刺激的信息.

圖3 有刺激和無刺激γ頻帶能量圖Fig.3 Power of gamma band with simulation and without simulation

為了研究γ頻帶的能量對于光柵朝向的調(diào)制作用，筆者采用最佳參數(shù)光柵刺激(時間頻率4 Hz，空間頻率0.2 cpd)，隨機(jī)變換其朝向，其中一個通道偏愛朝向和正交朝向的分析結(jié)果如圖4所示.由圖4可知，不同的朝向引起的神經(jīng)元響應(yīng)不同，其中偏愛朝向刺激前后γ頻帶能量差變化最為顯著.

圖4 偏愛朝向和正交朝向刺激前后γ頻帶能量差Fig.4 Power increase of the preferred and orthogonal orientation

2.3 spike的提取和分類

對原始數(shù)據(jù)使用通帶0.25～3 kHz的二階帶通Butterworth濾波器濾波，然后進(jìn)行閾值檢測，閾值水平設(shè)定為信號標(biāo)準(zhǔn)差的5倍［10］.利用主成分分析(Principal Component Analysis，PCA)結(jié)合 kmeans聚類進(jìn)行spike的分類提?。?1］.

某通道第一和第二主成分利用k-means聚類后如圖5(b)所示.鋒電位被分為明顯的兩類，選用發(fā)放較多的一類(圖5(c))做研究，其朝向調(diào)諧曲線如圖5(d)所示.由圖可以看出，該神經(jīng)元對90°光柵表現(xiàn)出較為明顯的朝向選擇性.

2.4 LFPγ頻帶能量和spike的朝向調(diào)諧分析

2.4.1 LFP與spike朝向調(diào)諧對比

分別用兩種信號研究了172個神經(jīng)元的朝向調(diào)諧特性.其中兩種信號有150個神經(jīng)元(116個有偏愛朝向，34個無偏愛朝向)的分析結(jié)果保持一致，有偏愛朝向的一致率達(dá)到了67.4%，總體一致率達(dá)到了87.2%.這說明γ頻帶能量可以做為研究朝向選擇性的有效指標(biāo)，其中一個神經(jīng)元的對比結(jié)果如圖6所示.

圖6為某一通道兩種方法的朝向調(diào)諧曲線圖.圖6(a)為基于LFPγ頻帶能量的朝向調(diào)諧曲線，朝向調(diào)諧指數(shù)ρLFP=0.710 9;圖6(b)為基于spike發(fā)放率的朝向調(diào)諧曲線，朝向調(diào)諧指數(shù)ρs=0.704 1.可以看出，該神經(jīng)元LFP頻帶能量在研究朝向調(diào)諧特性時能達(dá)到和spike發(fā)放率一樣的效果.

2.4.2 不同信噪比下實驗數(shù)據(jù)分析

為研究利用LFP和spike兩種信號分析朝向調(diào)諧特性時的抗噪能力，筆者研究了不同信噪比下兩種信號的分析結(jié)果.實驗數(shù)據(jù)由實測信號加白噪聲組成.從原始信號中提取出spike信號和LFPγ頻帶信號，重構(gòu)成純凈的神經(jīng)放電信號，如圖7所示.

圖7 仿真數(shù)據(jù)的建立Fig.7 The simulation data generation

將重構(gòu)的原始純凈信號加入不同SNR的高斯白噪聲［12］，分析不同信噪比下兩種信號的朝向調(diào)諧特性.針對LFP頻帶能量和spike發(fā)放率，計算在不同信噪比時兩種信號的朝向調(diào)諧指數(shù)和偏愛朝向.兩種方法下偏愛朝向以及調(diào)諧指數(shù)的分析結(jié)果如表1所示.

表1 不同信噪比下LFPγ頻帶能量與spike發(fā)放率的比較Tab.1 The comparison between LFP Gamma band power and spike under different SNR

由表1可以看出:隨著信噪比的降低，兩種算法的朝向調(diào)諧指數(shù)都有所降低，其中基于spike發(fā)放率的朝向調(diào)諧指數(shù)降低更快.當(dāng)信噪比在0 dB以下時，基于spike發(fā)放率的偏愛朝向出現(xiàn)了改變.這說明，與spike發(fā)放率相比，LFPγ頻帶能量在計算神經(jīng)元感受野朝向調(diào)諧特性方面具有更強的抗噪能力.

3 結(jié)論

利用LFP信號γ頻帶能量分析了V1區(qū)神經(jīng)元感受野的朝向調(diào)諧特性.用小波變換結(jié)合Welch功率譜估計提取并計算了局部場電位γ頻帶在光柵朝向刺激下的能量分布，并對比分析了LFP和spike兩種信號下的研究結(jié)果.

研究表明，與傳統(tǒng)基于spike發(fā)放率的朝向調(diào)諧特性計算方法相比，基于LFPγ頻帶能量計算方法具有更強的抗噪能力，在低信噪比的情況下也能有效地估計出神經(jīng)元感受野的朝向調(diào)諧特性.因此，LFPγ頻帶能量在低信噪比的情況下，研究V1區(qū)神經(jīng)元感受野的朝向調(diào)諧特性方面比spike發(fā)放率更有優(yōu)勢，為研究大鼠V1區(qū)神經(jīng)元感受野的朝向調(diào)諧特性提供了另一種可行的方法.

［1］ MELIZA C D，DAN Y.Receptive field modification in rat visual cortex induced by paired visual stimulation and single cell spiking［J］.Neuron，2006，49:183-189.

［2］ HILL D N，MEHTA SB，KLEINFELD D.Quality metrics to accompany spike sorting of extracellular signals［J］.The Journal of Neuroscience，2011，31(24):8699-8705.

［3］ GAWNE T J.The local and nonlocal components of the local field potential in awake primate visual cortex［J］.Journal of Computational Neuroscience，2010，29:615-623.

［4］ JUERGENS E，GUETTLER A，ECKHORN R.Visual stimulation elicits locked and induced gamma oscillations in monkey intracortical and EEG potentials，but not in human EEG ［J］.Experimental Brain Research，1999，129:247-259.

［5］ FRIEN A，ECKHOMR，WOELBEMT，et al.Fast oscillations display sharper orientation tuning than slower components of the same recordings in striate cortex of the awake monkey［J］.European Journal of Neuroscience，2000，12:1453-1465.

［6］ SIEGE L M，KNIG P.Frequency distribution of stimulus dependent synchronization in area 18 of the awake cat［J］.The Journal of Neuroscience，2003，23:4251-4260.

［7］ HENNE J A，SHAPLEY R.LFP power spectra in V1 cortex:the graded effect of stimulus contrast［J］.

Journal of Neurophysiology，2005，94(1):479-490.［8］ CARANDINIA M，F(xiàn)ERSTE D.Membrane potential and

firing rate in cat primary visual cortex［J］.The Journal of Neuroscience，2000，20(1):470-484.

［9］ 萬紅，郜麗賽，牛曉可.基于匹配小波變換的初級視皮層神經(jīng)元鋒電位分類［J］.鄭州大學(xué)學(xué)報:工學(xué)版，2011，32(3):90-94.

［10］ VARGAS-IRWIN C，DONOGHUE J P.Automated spike sorting using density grid contour clustering and subtractive waveform decomposition［J］.Journal of Neuroscience Methods，2007，164(1):1-18.

［11］ LEWICKI MS.A review of methods for spike sorting:the detection and classification of neural action potentials［J］.Computation in Neural Systems，1998(9):53-78.

［12］ 王靜，封洲燕.多通道神經(jīng)元鋒電位檢測和分類的新方法［J］.生物化學(xué)與生物物理進(jìn)展，2009，36(5):641-647.