0 引言

    目前有關骨關節(jié)疾病的診斷與功能評估大都基于醫(yī)生臨床查體、醫(yī)學影像成像技術,，通過分析MRI,、CT、X線檢查結果,，根據(jù)經(jīng)驗對病情做出大致判斷與分析,，制定相應救治方案，并以此作為后期康復指導的依據(jù),，總體來說以醫(yī)生的主觀判斷與評價為基礎,，缺乏對病人客觀、精確的功能參數(shù)分析^[1],。另外,，現(xiàn)有醫(yī)學影像設備有輻射且價格相對昂貴，不適合長期康復治療的效果評估與藥物指導,。

    近年,，3D實感攝像設備在國內(nèi)外發(fā)展很快，其中Kinect作為一款具有視覺和深度信息采集的圖像聲音傳感器,，可以通過跟蹤定位人體20個主要關節(jié)點的位置深度信息而實現(xiàn)人體三維骨架的構建^[2-3],，將它們應用在骨關節(jié)疾患的運動功能監(jiān)測中，可為醫(yī)務人員在診斷,、治療方案確立,、治療前后功能對比評價以及康復指導等過程提供一種更客觀有效的依據(jù)^[4-5]。

    本文將Kinect應用于骨關節(jié)空間角度的實時測量之中建立了骨關節(jié)空間角度數(shù)字化評估系統(tǒng),，并討論了測量中的骨關節(jié)空間角度測量,、功能分類評估等關鍵算法。實驗證明3D實感攝像技術應用于骨關節(jié)空間角度測量與功能評估領域,，在保證數(shù)字化測量要求的同時,，可使系統(tǒng)具有客觀功能評價,、非接觸、適應能力強,、快速高效,、準確、操作簡便和成本低等特點,。

1 系統(tǒng)構成及其工作原理

    圖1為骨關節(jié)空間角度數(shù)字化評估系統(tǒng)框圖,，主要由Kinect骨關節(jié)數(shù)據(jù)采集、骨關節(jié)空間角度計算,、骨關節(jié)功能分類評估等組成,。系統(tǒng)基于Unity開發(fā)平臺，利用體感設備Kinect進行深度數(shù)據(jù)采集,；實時計算關節(jié)點空間角度并以數(shù)字化方式在界面上顯示,；同時通過數(shù)據(jù)分析與對比，完成骨關節(jié)空間角度數(shù)字化評估工作,。

2 系統(tǒng)實現(xiàn)

2.1 Kinect骨關節(jié)數(shù)據(jù)采集

    基于Kinect的骨關節(jié)數(shù)據(jù)采集過程包括人物控制與骨骼點的綁定、骨架系統(tǒng)生成及關節(jié)點坐標表示,。

2.1.1 人物控制與骨骼點的綁定

    為了控制系統(tǒng)場景中人物角色的移動,，需要添加兩個人物控制器分別對應于人物場景模型和人物骨架模型。其中,，人物場景模型由一系列分別代表頭部,、肩部、手等人體部位的20個關節(jié)點組成,。人物骨架模型對應Kinect控制的模型,，控制場景中的角色的移動。

    將骨架模型中控制動作的關鍵預制體拖放于相對應的腳本變量上,，即可完成人物控制與骨骼點的綁定,。人物控制器與骨骼點綁定具體實現(xiàn)流程如圖2所示。

2.1.2 骨架系統(tǒng)生成及關節(jié)點坐標表示

    通過人物控制與骨骼點的綁定,，可實現(xiàn)模型與人物的同步,，也可以間接地獲取關節(jié)點的三維坐標。識別出20個關節(jié)點位置,，可以生成相應的骨架系統(tǒng),。

    關節(jié)點的三維坐標獲取方式采用場景物體獲取函數(shù)GameObject.Find()，將Kinect空間位置信息轉(zhuǎn)換為三維坐標信息,。

2.2 骨關節(jié)空間角度計算

    骨關節(jié)空間角度的計算是依據(jù)獲取的關節(jié)點三維坐標,，根據(jù)空間向量夾角公式進行計算。如圖3所示,，假設A,、B,、C分別代表3個相鄰關節(jié)點的空間坐標，以計算關節(jié)點A的空間角度為例,，具體操作步驟如下：

    (1)計算A與相鄰兩點B,、C組成的空間向量：

2.3 骨關節(jié)功能分類評估

    功能分類算法很多，常用的如K-medoids^[6],、CLARANS^[7],、K-means^[8]等?？紤]到骨關節(jié)功能龐大的數(shù)據(jù)量,，系統(tǒng)采用更適合大數(shù)據(jù)集分析的K-means聚類算法，首先對正常的骨關節(jié)數(shù)據(jù)進行分析,，得出標準的數(shù)據(jù)范圍^[9],，再以正常的數(shù)據(jù)范圍為分析比對的標準，對不同關節(jié)狀況的骨關節(jié)角度數(shù)據(jù)進行收集整理,。

    其中,，K-means算法以距離為分類的標準，根據(jù)某個距離函數(shù)完成分類,，本系統(tǒng)以歐式距離為距離準則函數(shù),，完成整體數(shù)據(jù)的聚類。在此過程中,，每次分類都將重新確定一次聚類中心,，其新的聚類中心使用取樣本均值的方法進行計算，具體如式(6)所示：

式中,，z_j表示當前類別j的均值結果,，n表示類別的樣本數(shù)量，x_i表示類別j的某一樣本數(shù)值,。

    系統(tǒng)采用誤差平方和作為目標函數(shù),，即誤差準則函數(shù)，具體定義如式(7)所示：

式中,，E表示對應類別的誤差平方和,；k為聚類類別數(shù)；n表示類別j的樣本數(shù)量,；x_i表示類別j的某一樣本數(shù)值,；z_j表示類別j的均值結果；d為數(shù)據(jù)x_i和z_j的偏差平方和,，即(x_i-z_j)²值,。最后將k個類別的誤差平方和相加即為總體誤差平方和E。

    聚類效果評估采用purity,、RI與F-meaures^[10]等指標,，其中purity是計算正確聚類的樣本數(shù)占總樣本數(shù)的比例,，RI是指度量正確的百分比，而F-measure由Precision（查準率）和Recall（查全率）兩個指標組合而,。假設F_i為每個分類的F-measure值,，P_i為每個分類的查準率，R_i為每個分類的查全率,，則F_i的計算公式如式(8)所示,。

    總的F值由每個分類F_i的加權平均得到。

3 系統(tǒng)實驗與結果

3.1 數(shù)據(jù)測試與分析

    利用骨關節(jié)空間角度測量方法,，實測66例健康骨關節(jié)角度數(shù)據(jù),。通過對骨關節(jié)活動狀態(tài)的分析，確定正常人群關節(jié)角度屈伸范圍,，如表1所示,。它們可以作為后期正常樣本的功能分類標準。

    在66例健康骨關節(jié)屈伸數(shù)據(jù)分析基礎上,，增加測試數(shù)據(jù)至152例,，通過聚類程序分析這152例數(shù)據(jù)。圖4所示為以左膝關節(jié)為代表的聚類結果,，其中橫軸為伸展角度,，縱軸為屈曲角度，數(shù)據(jù)被分為三類,，分別是正常的一類、屈伸角度不足的一類,、伸展角度不足的一類,。

    添加需要歸類的新樣本數(shù)據(jù)，會在分類基礎上進行歸類判定,，如圖4中以“□”號表示參與歸類的新的測量數(shù)據(jù),。

    以相同方式，利用上述K-means聚類算法可以完成其他下肢關節(jié)樣本數(shù)據(jù)的聚類分析,。表2記錄了人體下肢各關節(jié)角度數(shù)據(jù)的具體分類情況,，包括各關節(jié)類別，聚類過程中的迭代次數(shù),、狀態(tài),、聚類中心、單個誤差平方和以及整體誤差平方和,。其中單個誤差平方和可以反映每一個聚類結果的差異性,，整體誤差平方和可以判斷不同聚類結果的樣本差異。

3.2 功能分類算法驗證

    對下肢每一個關節(jié)測試并記錄正常,、欠屈和欠伸3種情況下各10組數(shù)據(jù),，一共30組數(shù)據(jù),，在聚類結果數(shù)據(jù)中輸入測試的數(shù)據(jù)，完成對數(shù)據(jù)的歸類,，驗證測試結果如表3所示,。

    表3中A_i～I_i分別為各個關節(jié)角度數(shù)據(jù)的測試結果情況，其中i={1,，2,，3，4,，5,，6}分別對應了左髖關節(jié)、右髖關節(jié),、左膝關節(jié),、右膝關節(jié)、左踝關節(jié),、右踝關節(jié)的數(shù)據(jù)測試結果,，如A₁～I₁代表的是左髖關節(jié)的數(shù)據(jù)測試結果。結果表明10個伸展困難的樣本中,，有7個被正確歸類,，2個被錯誤地歸類到正常，1個被錯誤地歸類到屈曲困難,；10個正常的測試樣本中,，有9個被正確歸類，1個被錯誤地歸類到屈曲困難,；10個屈曲困難的測試樣本完全正確歸類于屈曲困難類,。

    結合聚類評估原理之中涉及的計算方法，可以計算相應關節(jié)的聚類指標,，各關節(jié)聚類評估指標結果如表4所示,。

    由表4可知，利用K-means聚類算法能較好地完成骨關節(jié)空間角度數(shù)據(jù)的聚類和樣本歸類,，以左髖關節(jié)數(shù)據(jù)的聚類分析和指標計算為例,，3個聚類指標的值分別為0.867、0.837,、0.86,，即該關節(jié)數(shù)據(jù)的聚類準確度穩(wěn)定在0.83以上。就整個系統(tǒng)來而言,，所有關節(jié)聚類評估準確度能穩(wěn)定在0.8及其以上,。

3.3 系統(tǒng)測試

    系統(tǒng)最終以Unity為開發(fā)平臺，運用C#編程方式,，連接Kinect,，實時采集數(shù)據(jù)并進行骨關節(jié)功能評估,，圖5所示為系統(tǒng)測試界面圖。

    圖5所示的系統(tǒng)測試結果圖中包括Kinect圖像實時采集與顯示,，以及下肢左右髖關節(jié),、膝關節(jié)、踝關節(jié)空間角度顯示與功能評估等部分,。系統(tǒng)能夠?qū)崟r獲取關節(jié)點之間的角度并進行最大屈伸數(shù)據(jù)的更新記錄,，適用于多種參數(shù)的系統(tǒng)聚類，且可視化強,，能夠非接觸地完成骨關節(jié)空間角度功能的數(shù)字化評估,。

4 結論

    本文研究了一種非接觸的骨關節(jié)空間角度測量方法，并介紹了骨關節(jié)空間角度數(shù)字化評估系統(tǒng)的結構和工作原理,。實驗結果證明：系統(tǒng)以可視化方式為用戶呈現(xiàn)骨關節(jié)功能的數(shù)字化分析與功能評估結果,；聚類評估指標purity、RI與F-meaures等能穩(wěn)定在0.8及其以上,，基本滿足骨關節(jié)空間角度數(shù)字化評估要求,，可以為醫(yī)務人員在診斷、治療方案確立,、治療前后功能對比評價以及康復指導等過程提供一種更客觀有效的依據(jù),。

參考文獻

[1] CARLAK H F，GENCER N G,，BESIKCI C.Theoretical sssessment of electro-thermal imaging：a new technique for medical diagnosis[J].Infrared Physics and Technology,，2016，76：227-234.

[2] Ye Minxiang,，Yang Cheng,，VLADIMIR S.A depth camera motion analysis framework for tele-rehabilitation：motion capture and person-centric Kinematics analysis[J].IEEE Journal of Selected Topics in Signal Processing，2016,，10(5)：877-887.

[3] CHOI J S，KANG D W,，SEO J W.The development and evaluation of a program for leg-strengthening exercises and balance assessment using Kinect[J].Journal of Physical Therapy Science,，2016，28(1)：33-37.

[4] DOLATABADI E,，TAATI B,，MIHAILIDIS A.Concurrent validity of the microsoft kinect for Windows v2 for measuring spatiotemporal gait parameters[J].Medical Engineering and Physics，2016,，38(9)：952-958.

[5] 張克華,，王書平，尹曉紅,，等.復雜室內(nèi)環(huán)境下的單一Kinect導盲系統(tǒng)[J].光學精密工程,，2015,，23(8)：2419-2427.

[6] RAI A，UPADHYAY S H.Bearing performance degradation assessment based on a combination of empirical mode decomposition and k-medoids clustering[J].Mechanical Systems and Signal Processing,，2017,，93：16-29.

[7] 何童.不確定性目標的CLARANS聚類算法[J].計算機工程，2012,，38(11)：56-58.

[8] 李曉瑜,，俞麗穎，雷航,，等.一種K-means改進算法的并行化實現(xiàn)與應用[J].電子科技大學學報,，2017，46(1)：61-68.

[9] 葉燦華,，陳峰,，錢文偉，等.成人型髖關節(jié)發(fā)育不良的分型與治療[J].中華骨與關節(jié)外科雜志,，2017,，10(1)：70-75.

[10] 向培素.兩種聚類有效性評價指標的MATLAB實現(xiàn)[J].西南民族大學學報(自然科學版)，2013,，39(6)：1002-1005.

作者信息:

方艷紅,，楊雪梅，張紅英,，王學淵

（西南科技大學 信息工程學院,，四川 綿陽621010）