基于改進(jìn)動態(tài)時間規(guī)整算法的奶牛步態(tài)分割方法

蘇力德 張 永 王 健 尹 玉 宗哲英 鞏彩麗

(1.內(nèi)蒙古農(nóng)業(yè)大學(xué)機(jī)電工程學(xué)院， 呼和浩特 010018； 2.內(nèi)蒙古大學(xué)電子信息工程學(xué)院， 呼和浩特 010021)

0 引言

內(nèi)蒙古自治區(qū)是我國重要的奶牛養(yǎng)殖優(yōu)勢區(qū)域，其奶牛存欄數(shù)、牛奶產(chǎn)量，以及乳制品企業(yè)的規(guī)模、效益和影響力均位于全國前列。據(jù)我國奶業(yè)年鑒統(tǒng)計(jì)，截止到2017年，奶牛存欄數(shù)已達(dá)到216萬頭，牛奶產(chǎn)量達(dá)到730萬噸，乳制品產(chǎn)量達(dá)到337萬噸。奶牛疾病阻礙了奶業(yè)的健康發(fā)展，據(jù)統(tǒng)計(jì)，我國大中型牧場臨床跛行率和嚴(yán)重跛行率分別為(31±12)%和(10±6)%[1]。跛行不僅導(dǎo)致奶牛身體狀況下降、牛奶產(chǎn)量降低，還將使奶牛過早地被淘汰。因此，準(zhǔn)確分割奶牛步態(tài)、識別奶牛跛行對減少養(yǎng)殖場的經(jīng)濟(jì)損失、提高福利化養(yǎng)殖水平具有重要意義。

步態(tài)由周期性的步幅組成，如何從步態(tài)序列中分割單個步幅已成為運(yùn)動分析領(lǐng)域的研究熱點(diǎn)。國內(nèi)外學(xué)者主要采用峰值檢測法[2-3]、視頻分割法[4]、模板匹配方法分割步態(tài)。峰值檢測法必須滿足一定的閾值條件才可以進(jìn)行分割；視頻分割方法工作量大、自動化程度不高，限制了在實(shí)際中的應(yīng)用；模板匹配方法通過計(jì)算兩個時間序列之間相似程度對步態(tài)進(jìn)行分割。國內(nèi)外基于模板匹配方法的步態(tài)分割主要集中在人體步態(tài)研究，對奶牛步態(tài)分割的研究尚未見報(bào)道。文獻(xiàn)[5-6]利用基于模板匹配的互相關(guān)方法對人體步態(tài)進(jìn)行分割，但步幅模板長度和形式固定，不適用于不同步幅長度序列的分割。

動態(tài)時間規(guī)整算法(Dynamic time warping,DTW)是一種非線性匹配算法，該算法來源于動態(tài)規(guī)整原理，解決時間長度不一致時的模板匹配問題。文獻(xiàn)[7]利用動態(tài)時間規(guī)整算法對人體步態(tài)進(jìn)行分割，通過研究步態(tài)周期數(shù)據(jù)間的差異性，證明該方法在臨床疾病檢測上優(yōu)于其他方法。文獻(xiàn)[8]利用傳感器技術(shù)采集步態(tài)信息，以動態(tài)時間規(guī)整算法分割步態(tài)，通過k近鄰分類器對偏癱患者進(jìn)行識別，得出曲線下面積為0.94。文獻(xiàn)[9]利用動態(tài)時間規(guī)整算法對帕金森患者步態(tài)進(jìn)行分割，獲得100%的分割準(zhǔn)確率。文獻(xiàn)[10]利用動態(tài)時間規(guī)整算法對25名老年人的步態(tài)進(jìn)行分割，分割準(zhǔn)確率達(dá)到97%。

常規(guī)的動態(tài)時間規(guī)整算法樣本序列是一維數(shù)據(jù)，而奶牛步態(tài)加速度信號的時間和幅度是周期性變化的。本文利用三維加速度傳感器采集奶牛趾蹄加速度信號，將時間和幅度組成的二維數(shù)據(jù)作為步態(tài)序列，提出一種改進(jìn)的動態(tài)時間規(guī)整算法，對奶牛步態(tài)進(jìn)行分割，并對每個步幅提取步態(tài)特征，利用Wilcoxon秩和檢驗(yàn)分析健康奶牛和跛行奶牛特征值之間的關(guān)系，選取與跛行相關(guān)的6個特征，利用邏輯回歸方法建立跛行識別模型，以準(zhǔn)確識別奶牛的跛行。

1 材料和方法

1.1 測量設(shè)備

實(shí)驗(yàn)場地為內(nèi)蒙古旗幟牧業(yè)有限公司，選擇牧坤電子科技有限公司生產(chǎn)的M63X系列三維加速度傳感器，測量范圍±8g，采樣頻率50 Hz。將傳感器封裝在密閉防水盒，固定在制作好的綁帶上，形成可穿戴式設(shè)備。數(shù)據(jù)通過ZigBee協(xié)議無線傳輸?shù)交荆?shí)時上傳到計(jì)算機(jī)。奶牛進(jìn)入長、寬分別為23、0.8 m的測試通道，射頻識別閱讀器識別身份信息，使用PCO.dimax系列高速攝像機(jī)拍攝奶牛行走視頻，如圖1所示。由于奶牛85%的趾蹄疾病發(fā)生在后趾外側(cè)爪上[11]，因此將加速度傳感器綁在后肢外側(cè)，為避免奶牛應(yīng)激反應(yīng)產(chǎn)生的測量誤差，每頭牛綁上加速度傳感器行走10 min后開始采集，以此適應(yīng)測量設(shè)備。

圖1 測試環(huán)境及設(shè)備Fig.1 Test environment and equipment1.加速度測量裝置 2.捆綁式加速度傳感器 3.無線接收基站

1.2 數(shù)據(jù)采集與預(yù)處理

數(shù)據(jù)采集使用實(shí)驗(yàn)室虛擬儀器工程平臺(Laboratory virtual instrumentation engineering workbench，LabVIEW)實(shí)時顯示在上位機(jī)，LabVIEW與Access數(shù)據(jù)庫連接實(shí)現(xiàn)對數(shù)據(jù)的管理，數(shù)據(jù)處理及算法實(shí)現(xiàn)通過Matlab軟件編程，并與LabVIEW的DLL(Dynamic-link library)及ActiveX等技術(shù)實(shí)現(xiàn)混合編程。選取30頭荷斯坦奶牛，其中健康奶牛21頭、跛行奶牛9頭，在測量通道內(nèi)連續(xù)采集6 d，每天上午和下午擠完奶各采集一次，每次重復(fù)采集一次。為確保數(shù)據(jù)的準(zhǔn)確性高速攝像機(jī)與加速度傳感器系統(tǒng)實(shí)現(xiàn)同步采集，其中高速攝像機(jī)以940 f/s速率拍攝奶牛行走視頻，并以該視頻數(shù)據(jù)作為黃金標(biāo)準(zhǔn)。3名獸醫(yī)通過觀察每頭奶牛行走視頻，以5分制[12]評分規(guī)則對其進(jìn)行評分，取平均值作為跛行評分值。在三維加速度傳感器中，定義X軸為垂直方向、Y軸為前進(jìn)方向、Z軸為水平方向。

由于噪聲的存在，加速度數(shù)據(jù)不能直接應(yīng)用于步態(tài)分割中，需對其進(jìn)行預(yù)處理。小波變換去噪方法能有效區(qū)分信號突變部分和噪聲[13]，本文通過小波變換方法對加速度數(shù)據(jù)進(jìn)行預(yù)處理。利用Matlab軟件在Daubechies小波基框架內(nèi)，選擇DB(Daubechies)6小波進(jìn)行3層分解，通過信號重構(gòu)去除干擾信號。

2 改進(jìn)動態(tài)時間規(guī)整算法步態(tài)分割實(shí)現(xiàn)

將連續(xù)的步態(tài)序列分割為單個步幅稱為步態(tài)分割。對于周期性步態(tài)的分割相對簡單實(shí)現(xiàn)，然而在現(xiàn)實(shí)中步幅并不完全是周期性的，不同跛行評分、不同體質(zhì)量和泌乳的奶牛步幅行走模式均有所不同。這些異構(gòu)的步態(tài)序列是步態(tài)分割的主要難點(diǎn)之一，而動態(tài)時間規(guī)整算法主要解決序列長度不一致時模板匹配的問題，適用于不同行走模式下的步態(tài)分割。奶牛步幅定義為同一肢體連續(xù)2次接觸地面之間的間隔，對應(yīng)的時間差為一個步態(tài)周期。步態(tài)周期分為支撐階段和擺動階段，其中，在支撐階段，趾蹄與地面接觸開始，并在推開后終止。在擺動階段，趾蹄需在地面以上抬起。一個步態(tài)周期中，擺動時間占31.6%，支撐階段占68.4%[14]，平均步態(tài)周期為(1.22±0.05) s[15]。通過前期實(shí)驗(yàn)分析可知，利用垂直方向加速度信號分割步態(tài)容易存在誤差，而水平加速度很大程度上受到乳房的影響，整個步行過程中變化并不明顯。前進(jìn)方向加速度信號極大值點(diǎn)突出容易區(qū)分，本文采用該方向加速度信號對步態(tài)進(jìn)行分割。

利用動態(tài)時間規(guī)整算法分割步態(tài)是通過求模板樣本序列映射到測試樣本序列的成本距離，利用遞推公式計(jì)算滿足最小累積成本的規(guī)整路徑，尋找與模板樣本相似的子序列，最終從連續(xù)的步態(tài)序列中分割單個步幅。本文算法的步態(tài)分割整個流程如圖2所示。

圖2 步態(tài)分割算法流程圖Fig.2 Flow chart of gait segmentation algorithm

2.1 步態(tài)序列特征點(diǎn)提取

奶牛一個步態(tài)周期中趾蹄觸地和足底離地時刻前進(jìn)方向加速度信號均有極值出現(xiàn)，利用小波變換過零點(diǎn)方法提取信號極值點(diǎn)[16-17]，并將該點(diǎn)作為樣本特征點(diǎn)構(gòu)成樣本序列。信號突變點(diǎn)檢測是對原始信號進(jìn)行“磨光”的過程，磨光可以理解為對信號進(jìn)行平滑處理[18-19]。磨光函數(shù)計(jì)算公式為

(1)

式中θ(t)——低通濾波器

假設(shè)θ(t)二次可導(dǎo)，對其求導(dǎo)表示為

(2)

函數(shù)ψ(1)(t)和ψ(2)(t)滿足小波容許條件，可作為小波基函數(shù)，以此定義的卷積型小波變換為

(3)

式中 *——卷積運(yùn)算符號s——尺度因子

θs(t)——θ(t)在s下的伸縮

f(t)——原始信號

2.2 模板樣本建立

將連續(xù)步態(tài)序列定義為測試樣本序列，向量表示為Q=(q(1),q(2),…,q(m))，其中m為序列長度。模板樣本向量表示為P=(p(1),p(2),…,p(n))，其中n為序列長度。由于不同跛行程度的奶牛行走模式不同，將步態(tài)樣本分為健康、輕度跛行、中度跛行和嚴(yán)重跛行4個組。根據(jù)高速攝像機(jī)拍攝的視頻手動分割步態(tài)，將每個步幅起點(diǎn)和終點(diǎn)手動標(biāo)記為模板和黃金標(biāo)準(zhǔn)。為建立能代表大范圍步幅的模板序列，每組手動標(biāo)記的步幅通過Matlab的interp1函數(shù)線性內(nèi)插至各組樣本中，并取平均值生成模板序列，待分割的測試樣本及生成的模板樣本如圖3所示。

圖3 模板樣本與測試樣本Fig.3 Template sample and test sample

2.3 累積成本矩陣計(jì)算

為計(jì)算2個不同樣本序列的最小規(guī)整路徑，構(gòu)造了一個N×M的最小距離矩陣D。該矩陣中的每行表示模板樣本的一個序列到測試樣本的一個序列之間的距離。如果2個樣本序列相似，則局部距離d較小，規(guī)整成本較低。相反，序列不相似，則d較大，規(guī)整成本較大。本文使用歐氏距離計(jì)算局部距離，計(jì)算公式為

(4)

式中p(n)——模板樣本序列元素

q(m)——測試樣本序列元素

為了尋找模板樣本與測試樣本間最小規(guī)整路徑，定義了一個N×M的累積成本矩陣C。該矩陣表示模板樣本映射到測試樣本各個部分的累積成本，是通過距離矩陣建立累積成本矩陣來搜索最佳規(guī)整路徑的。將距離矩陣D的第1列元素值依次求和得到累積成本矩陣C的第1列元素。而矩陣C的第一行元素由矩陣D的第1行構(gòu)成，定義該行為起始行，用于計(jì)算模板樣本映射到測試樣本的累積成本。累積成本矩陣其余元素計(jì)算公式為

C(n,m)=d(n,m)+
min(C(n-1,m-1)，C(n-1,m)，C(n,m-1))

(5)

由式(5)可知，其可以通過將累積成本矩陣C3個方向的最小成本元素與距離矩陣D的各個元素相加來計(jì)算。因此，累積成本矩陣中元素表示模板樣本映射到測試樣本的最小累積成本。

2.4 算法改進(jìn)

由于奶牛前進(jìn)方向加速度信號極值是隨時間周期性變化的，即每個步態(tài)周期特征點(diǎn)是由時間序列和幅度序列構(gòu)成的二維序列。本文將該二維序列距離作為改進(jìn)動態(tài)時間規(guī)整算法的局部距離。設(shè)模板樣本的時間為pt，幅值為pa，測試樣本的時間為qt，幅值為qa，因此，樣本序列分別為Ps={(pt(1),pa(1)),(pt(2),pa(2)),…,(pt(n),pa(n))},Qs={(qt(1),qa(1)),(qt(2),qa(2)),…,(qt(m),qa(m))}，將式(4)中的p(n)和q(m)用該二維序列來代替，二維歐氏距離計(jì)算公式為

(6)

在將式(6)代入式(5)中，計(jì)算模板樣本映射到測試樣本的累積成本矩陣。

2.5 最小規(guī)整路徑確定

最小規(guī)整路徑反映2個樣本相似性的同時可以確定步幅的分割起點(diǎn)及終點(diǎn)。將規(guī)整路徑定義為W，該路徑由成本距離矩陣的每一行元素局部最小值構(gòu)成，計(jì)算公式為

W=(w1,w2,…,wk，…，wK)

(7)

式中K——路徑長度

wk——路徑中的點(diǎn)元素

k——測試樣本序列長度

假設(shè)w0為規(guī)整路徑第1行元素，計(jì)算式為

w0=min {C(1,1),C(1，2)，…,C(1,m)}

(8)

動態(tài)時間規(guī)整算法對規(guī)整路徑W做了一定限制，只有滿足下面3個條件才可成立。

(1)邊界約束

(9)

任何2個序列都有可能變化，但其各部分的先后次序不會改變，因此所選路徑必定是從網(wǎng)格左下角出發(fā)，在右上角結(jié)束。

(2)連續(xù)性約束

對路徑中wk-1=(n′,m′)和wk=(n,m) 2點(diǎn)，需滿足連續(xù)性約束條件

(10)

式(10)表明在計(jì)算路徑時不可能跨過某個點(diǎn)去匹配，只能和相鄰點(diǎn)對齊。這可以保證測試樣本和模板樣本中每個坐標(biāo)都出現(xiàn)在W中。

(3)單調(diào)性約束

對路徑中wk-1=(n′,m′)和wk=(n,m) 2點(diǎn)，要求

(11)

式(11)限制了W上面的點(diǎn)必須要隨時間單調(diào)進(jìn)行。由式(10)、(11)可知，對于任意一個網(wǎng)格點(diǎn)(n,m)，下一個通過的網(wǎng)格點(diǎn)只可能是(n,m-1)、(n-1,m-1)和(n-1,m)。最小規(guī)整成本路徑求解過程如圖4所示。

圖4 最小規(guī)整路徑Fig.4 Minimum warping path

通過計(jì)算最小規(guī)整路徑從連續(xù)的步態(tài)序列中分割單個步幅。路徑W的起點(diǎn)和終點(diǎn)分別代表步幅的起點(diǎn)和終點(diǎn)，每個分割點(diǎn)的累積成本必須小于分割閾值，因此，選取合適的分割閾值對步態(tài)進(jìn)行分割。

3 跛行識別模型建立

3.1 步態(tài)特征提取

對于每個分割的步幅分別提取統(tǒng)計(jì)學(xué)和運(yùn)動學(xué)步態(tài)特征參數(shù)，包括均值A(chǔ)[20]、方差[20]、合加速度Ml[20]、運(yùn)動變化量Mv[20]、信號幅度面積Sma[21]、平均強(qiáng)度Al[21]、偏度Sk[22]、峰度Fk[22]、步態(tài)周期Ts、步幅數(shù)Fs、支持時間St和步幅長度Ls。步幅數(shù)為連續(xù)序列中分割的步幅數(shù)量，支撐時間通過前進(jìn)加速度數(shù)據(jù)變化率計(jì)算，步態(tài)周期為步幅分割點(diǎn)之間的時間差，而步幅長度由步行總距離除去步幅數(shù)得出。

由于健康奶牛和跛行奶牛樣本數(shù)量不一致，且可能存在非正態(tài)分布，因此采用Wilcoxon秩和檢驗(yàn)對提取的特征進(jìn)行分析，確定跛行奶牛與健康奶牛特征參數(shù)間關(guān)系。

3.2 跛行識別模型

目前，國內(nèi)跛行識別研究主要集中在基于計(jì)算機(jī)視覺和壓力板測量方法上[23-25]，國外主要通過基于各種傳感器及壓敏墊方法采集數(shù)據(jù)進(jìn)行跛行識別的研究[26-30]。本文利用邏輯回歸法建立奶牛跛行識別模型[31]，根據(jù)統(tǒng)計(jì)檢驗(yàn)結(jié)果選取與跛行相關(guān)的6個步態(tài)特征參數(shù)，建立單因素邏輯回歸模型。將輕度、中度和嚴(yán)重度跛行合并為一個類別，稱之為“跛行”。以二分類變量為模型輸出，“1”表示跛行，“0”表示健康。

對于二分類問題可設(shè)置一個閾值，如果G(x)≥0.5，則預(yù)測y=1，即該奶牛患有跛行。如果G(x)≤0.5，則預(yù)測y=0，即該奶牛為健康。為預(yù)測輸出類似于概率一樣介于0～1之間，將sigmoid函數(shù)作為映射函數(shù)得到邏輯回歸模型,計(jì)算公式為

(12)

式中xi——步態(tài)特征值

α——回歸截距β——回歸系數(shù)

G——跛行發(fā)生概率

yi——觀測值標(biāo)簽i——樣本序號

邏輯回歸是非線性模型，利用最大似然估計(jì)法對模型系數(shù)進(jìn)行估計(jì)。設(shè)Gi=G(yi=1|xi)為給定xi條件下得到的跛行yi=1概率，而yi=0的條件概率為G(yi=0|xi)=1-Gi。因此，邏輯回歸模型的似然函數(shù)計(jì)算公式為

(13)

對式(13)兩邊取對數(shù)，得到邏輯回歸對數(shù)似然函數(shù)，并利用極大似然估計(jì)法求最佳回歸系數(shù)。

4 結(jié)果與分析

4.1 步態(tài)分割結(jié)果對比

本文共采集30頭奶牛976個樣本數(shù)據(jù)，其中健康奶牛(得分1)21頭共613個樣本、輕度跛行(得分2)奶牛2頭共102個樣本、中度跛行(得分3)奶牛4頭共137個樣本、嚴(yán)重跛行(得分大于等于4)奶牛3頭共124個樣本。從每組生成模板樣本，計(jì)算每組模板樣本與對應(yīng)測試樣本間累積成本矩陣，根據(jù)動態(tài)規(guī)整原則尋找最小規(guī)整路徑，確定分割起點(diǎn)和終點(diǎn)，最終從連續(xù)的步態(tài)序列中分割單個步幅。

圖5為利用前進(jìn)方向加速度信號分割的單個步幅，步態(tài)周期開始于趾蹄接觸地面終止于趾蹄下一次接觸地面，圖5中紅色虛線之間為一個步態(tài)周期。在趾蹄接觸地面時刻前進(jìn)方向加速度值反向最大，之后有一個正方向的加速度信號后進(jìn)入支撐階段，在支撐階段趾蹄相對靜止，因此前進(jìn)方向加速度值變化不大，圖中黑色虛線之間表示支撐階段。支撐階段后進(jìn)入擺動階段，擺動階段開始于趾蹄抬起地面，終止于趾蹄接觸地面，趾蹄抬起時刻前進(jìn)方向加速度值正向最大，圖中黑色虛線和紅色虛線之間為擺動階段。通過分析可知，前進(jìn)方向加速度信號下峰值均出現(xiàn)在每個步態(tài)周期開始時刻，因此該信號下峰值對應(yīng)步幅的起點(diǎn)和終點(diǎn)。

圖5 分割的步幅示意圖Fig.5 Schematic of segmented stride

步態(tài)分割主要以減少步幅丟失和最小化錯誤分割為目的。本文利用混淆矩陣來對分割結(jié)果進(jìn)行分析，精確度U、靈敏度V、準(zhǔn)確率E計(jì)算公式為

(14)

(15)

(16)

式中Tp——正確分割為步幅的樣本數(shù)量

Fp——錯誤分割為步幅的樣本數(shù)量

FN——無法識別為步幅的樣本數(shù)量

TN——無效步幅的樣本數(shù)量

對30頭奶牛改進(jìn)前后的步態(tài)分割結(jié)果對比如表1所示。將每個分割的步幅與黃金標(biāo)準(zhǔn)進(jìn)行比較來驗(yàn)證分割的準(zhǔn)確性。

表1 算法改進(jìn)前分割結(jié)果對比Tab.1 Comparison of segmentation results %

由表1結(jié)果可知，本文算法的精確度、靈敏度、準(zhǔn)確率平均值可提高到89.53%、95.51%、87.49%，比改進(jìn)前分別提高了5.31、4.48、8.43個百分點(diǎn)。其中，輕度跛行奶牛步態(tài)分割靈敏度、準(zhǔn)確率提高最為明顯，較改進(jìn)前分別提高7.22、12.75個百分點(diǎn)。而中度跛行奶牛步態(tài)分割精確度提高到86.44%，較改進(jìn)前提高了10.81個百分點(diǎn)。

為進(jìn)一步驗(yàn)證該算法性能，利用峰值檢測法、自相關(guān)函數(shù)法和本文算法對30頭奶牛進(jìn)行步態(tài)分割試驗(yàn)，結(jié)果表明，峰值檢測法總準(zhǔn)確率87.44%，自相關(guān)函數(shù)法總準(zhǔn)確率88.82%，改進(jìn)算法在步態(tài)分割準(zhǔn)確率方面有所提高，總體準(zhǔn)確率提高到90.57%，相較自相關(guān)函數(shù)法和峰值檢測法分別提高了1.75、3.13個百分點(diǎn)。

4.2 步態(tài)特征提取結(jié)果

為研究健康奶牛和跛行奶牛特征參數(shù)間關(guān)系，從分割的每個步幅提取步態(tài)特征參數(shù)，并取均值。利用Wilcoxon秩和檢驗(yàn)分析各特征值，結(jié)果如表2所示，表中特征值用差異性高到低依次進(jìn)行排列。分析可知，信號幅度面積、平均強(qiáng)度、步態(tài)周期、步幅長度、支撐時間、前進(jìn)方向加速度均值、運(yùn)動變化量和方差在健康奶牛和跛行奶牛之間均表現(xiàn)顯著性差異(P<0.05)。而峰度和偏度在2個方向上均不存在顯著性差異。

表2 Wilcoxon秩和檢驗(yàn)分析結(jié)果Tab.2 Analysis results of Wilcoxon rank sum test

4.3 跛行預(yù)測模型的測試

根據(jù)統(tǒng)計(jì)學(xué)分析結(jié)果，選取支撐時間、步幅長度、平均強(qiáng)度、信號幅度面積、前進(jìn)方向加速度均值和運(yùn)動變化量等6個特征，利用邏輯回歸法建立單因素回歸模型。模型輸入為特征值，輸出為二分類結(jié)果，其中“1”表示跛行，“0”表示健康，以0.5作為模型分類閾值。將樣本數(shù)據(jù)的70%作為訓(xùn)練集來獲得最佳回歸系數(shù)，30%作為測試集進(jìn)行識別模型驗(yàn)證，通過極大似然估計(jì)法分別求出各模型回歸系數(shù)，確定邏輯回歸方程。

利用接受者操作特征曲線(Receiver operating characteristic curve，ROC)下的面積[32]評估各預(yù)測模型，該面積范圍在0～1之間，值越大表明識別率越高，各預(yù)測模型識別率與ROC曲線下的面積見表3。結(jié)果可知，單因素邏輯回歸模型中以前進(jìn)方向加速度均值為自變量的模型獲得較高的識別準(zhǔn)確率，為89.45%，ROC曲線下的面積為0.902。其他5個單因素回歸模型識別率均不小于81.72%，ROC曲線下的面積均不小于0.804。

表3 跛行預(yù)測模型結(jié)果Tab.3 Results of lameness prediction model

5 結(jié)論

(1)利用三維加速度傳感器和無線數(shù)據(jù)采集設(shè)備采集奶牛后趾蹄加速度信號，并對其進(jìn)行預(yù)處理，通過分析前進(jìn)方向加速度信號可以為奶牛步態(tài)的分割提供理論依據(jù)。

(2)提出采用改進(jìn)的動態(tài)時間規(guī)整算法對奶牛步態(tài)進(jìn)行分割，步態(tài)分割精確度、靈敏度、準(zhǔn)確率的均值較改進(jìn)前分別提高了5.31、4.48、8.43個百分點(diǎn)。總體準(zhǔn)確率為90.57%，比自相關(guān)函數(shù)法和峰值檢測法總準(zhǔn)確率分別提高了1.75、3.13個百分點(diǎn)。因此，本文算法準(zhǔn)確率高，能獲得較好的分割效果。

(3)依據(jù)提取的步態(tài)特征，建立跛行識別模型，獲得了較高的識別率，其中以前進(jìn)方向加速度均值、信號幅度面積、平均強(qiáng)度、運(yùn)動變化量、支撐時間和步幅長度建立的模型識別率分別為89.45%、86.81%、86.15%、85.71%、83.44%、81.72%。為后期奶牛跛行識別研究奠定了良好的基礎(chǔ)。