摘要:
為了實現(xiàn)對骨關(guān)節(jié)空間角度的數(shù)字化測量與功能客觀評價,建立了骨關(guān)節(jié)角度數(shù)字化測量與功能評估系統(tǒng)。對該系統(tǒng)所采用的骨關(guān)節(jié)空間角度測量、功能評估算法進行研究。首先,根據(jù)Kinect數(shù)據(jù)采集原理及空間向量計算方法介紹了骨關(guān)節(jié)空間角度測量算法,即在Kinect深度數(shù)據(jù)流基礎上將骨關(guān)節(jié)點空間位置信息轉(zhuǎn)換為對應的三維空間坐標,根據(jù)空間向量夾角公式計算骨關(guān)節(jié)空間角度。然后,在分析比較分類評價方法的性能基礎上,說明了采用K-means聚類方法進行骨關(guān)節(jié)功能評估的算法。實驗結(jié)果表明,系統(tǒng)可以為用戶提供骨關(guān)節(jié)的數(shù)字化測量與功能評估結(jié)果,評估指標purity、RI與F-meaures等能穩(wěn)定在0.8及其以上,基本滿足骨關(guān)節(jié)功能評價的非接觸、客觀高效、適應能力強、準確度高、操作簡便和成本低等要求。
?
0 引言
目前有關(guān)骨關(guān)節(jié)疾病的診斷與功能評估大都基于醫(yī)生臨床查體、醫(yī)學影像成像技術(shù),通過分析MRI、CT、X線檢查結(jié)果,根據(jù)經(jīng)驗對病情做出大致判斷與分析,制定相應救治方案,并以此作為后期康復指導的依據(jù),總體來說以醫(yī)生的主觀判斷與評價為基礎,缺乏對病人客觀、精確的功能參數(shù)分析[1]。另外,現(xiàn)有醫(yī)學影像設備有輻射且價格相對昂貴,不適合長期康復治療的效果評估與藥物指導。
近年,3D實感攝像設備在國內(nèi)外發(fā)展很快,其中Kinect作為一款具有視覺和深度信息采集的圖像聲音傳感器,可以通過跟蹤定位人體20個主要關(guān)節(jié)點的位置深度信息而實現(xiàn)人體三維骨架的構(gòu)建[2-3],將它們應用在骨關(guān)節(jié)疾患的運動功能監(jiān)測中,可為醫(yī)務人員在診斷、治療方案確立、治療前后功能對比評價以及康復指導等過程提供一種更客觀有效的依據(jù)[4-5]。
本文將Kinect應用于骨關(guān)節(jié)空間角度的實時測量之中建立了骨關(guān)節(jié)空間角度數(shù)字化評估系統(tǒng),并討論了測量中的骨關(guān)節(jié)空間角度測量、功能分類評估等關(guān)鍵算法。實驗證明3D實感攝像技術(shù)應用于骨關(guān)節(jié)空間角度測量與功能評估領域,在保證數(shù)字化測量要求的同時,可使系統(tǒng)具有客觀功能評價、非接觸、適應能力強、快速高效、準確、操作簡便和成本低等特點。
1 系統(tǒng)構(gòu)成及其工作原理
圖1為骨關(guān)節(jié)空間角度數(shù)字化評估系統(tǒng)框圖,主要由Kinect骨關(guān)節(jié)數(shù)據(jù)采集、骨關(guān)節(jié)空間角度計算、骨關(guān)節(jié)功能分類評估等組成。系統(tǒng)基于Unity開發(fā)平臺,利用體感設備Kinect進行深度數(shù)據(jù)采集;實時計算關(guān)節(jié)點空間角度并以數(shù)字化方式在界面上顯示;同時通過數(shù)據(jù)分析與對比,完成骨關(guān)節(jié)空間角度數(shù)字化評估工作。
2 系統(tǒng)實現(xiàn)
2.1 Kinect骨關(guān)節(jié)數(shù)據(jù)采集
基于Kinect的骨關(guān)節(jié)數(shù)據(jù)采集過程包括人物控制與骨骼點的綁定、骨架系統(tǒng)生成及關(guān)節(jié)點坐標表示。
2.1.1 人物控制與骨骼點的綁定
為了控制系統(tǒng)場景中人物角色的移動,需要添加兩個人物控制器分別對應于人物場景模型和人物骨架模型。其中,人物場景模型由一系列分別代表頭部、肩部、手等人體部位的20個關(guān)節(jié)點組成。人物骨架模型對應Kinect控制的模型,控制場景中的角色的移動。
將骨架模型中控制動作的關(guān)鍵預制體拖放于相對應的腳本變量上,即可完成人物控制與骨骼點的綁定。人物控制器與骨骼點綁定具體實現(xiàn)流程如圖2所示。
2.1.2 骨架系統(tǒng)生成及關(guān)節(jié)點坐標表示
通過人物控制與骨骼點的綁定,可實現(xiàn)模型與人物的同步,也可以間接地獲取關(guān)節(jié)點的三維坐標。識別出20個關(guān)節(jié)點位置,可以生成相應的骨架系統(tǒng)。
關(guān)節(jié)點的三維坐標獲取方式采用場景物體獲取函數(shù)GameObject.Find(),將Kinect空間位置信息轉(zhuǎn)換為三維坐標信息。
2.2 骨關(guān)節(jié)空間角度計算
骨關(guān)節(jié)空間角度的計算是依據(jù)獲取的關(guān)節(jié)點三維坐標,根據(jù)空間向量夾角公式進行計算。如圖3所示,假設A、B、C分別代表3個相鄰關(guān)節(jié)點的空間坐標,以計算關(guān)節(jié)點A的空間角度為例,具體操作步驟如下:
(1)計算A與相鄰兩點B、C組成的空間向量:
2.3 骨關(guān)節(jié)功能分類評估
功能分類算法很多,常用的如K-medoids[6]、CLARANS[7]、K-means[8]等??紤]到骨關(guān)節(jié)功能龐大的數(shù)據(jù)量,系統(tǒng)采用更適合大數(shù)據(jù)集分析的K-means聚類算法,首先對正常的骨關(guān)節(jié)數(shù)據(jù)進行分析,得出標準的數(shù)據(jù)范圍[9],再以正常的數(shù)據(jù)范圍為分析比對的標準,對不同關(guān)節(jié)狀況的骨關(guān)節(jié)角度數(shù)據(jù)進行收集整理。
其中,K-means算法以距離為分類的標準,根據(jù)某個距離函數(shù)完成分類,本系統(tǒng)以歐式距離為距離準則函數(shù),完成整體數(shù)據(jù)的聚類。在此過程中,每次分類都將重新確定一次聚類中心,其新的聚類中心使用取樣本均值的方法進行計算,具體如式(6)所示:
式中,zj表示當前類別j的均值結(jié)果,n表示類別的樣本數(shù)量,xi表示類別j的某一樣本數(shù)值。
系統(tǒng)采用誤差平方和作為目標函數(shù),即誤差準則函數(shù),具體定義如式(7)所示:
式中,E表示對應類別的誤差平方和;k為聚類類別數(shù);n表示類別j的樣本數(shù)量;xi表示類別j的某一樣本數(shù)值;zj表示類別j的均值結(jié)果;d為數(shù)據(jù)xi和zj的偏差平方和,即(xi-zj)2值。最后將k個類別的誤差平方和相加即為總體誤差平方和E。
聚類效果評估采用purity、RI與F-meaures[10]等指標,其中purity是計算正確聚類的樣本數(shù)占總樣本數(shù)的比例,RI是指度量正確的百分比,而F-measure由Precision(查準率)和Recall(查全率)兩個指標組合而。假設Fi為每個分類的F-measure值,Pi為每個分類的查準率,Ri為每個分類的查全率,則Fi的計算公式如式(8)所示。
總的F值由每個分類Fi的加權(quán)平均得到。
3 系統(tǒng)實驗與結(jié)果
3.1 數(shù)據(jù)測試與分析
利用骨關(guān)節(jié)空間角度測量方法,實測66例健康骨關(guān)節(jié)角度數(shù)據(jù)。通過對骨關(guān)節(jié)活動狀態(tài)的分析,確定正常人群關(guān)節(jié)角度屈伸范圍,如表1所示。它們可以作為后期正常樣本的功能分類標準。
在66例健康骨關(guān)節(jié)屈伸數(shù)據(jù)分析基礎上,增加測試數(shù)據(jù)至152例,通過聚類程序分析這152例數(shù)據(jù)。圖4所示為以左膝關(guān)節(jié)為代表的聚類結(jié)果,其中橫軸為伸展角度,縱軸為屈曲角度,數(shù)據(jù)被分為三類,分別是正常的一類、屈伸角度不足的一類、伸展角度不足的一類。
添加需要歸類的新樣本數(shù)據(jù),會在分類基礎上進行歸類判定,如圖4中以“□”號表示參與歸類的新的測量數(shù)據(jù)。
以相同方式,利用上述K-means聚類算法可以完成其他下肢關(guān)節(jié)樣本數(shù)據(jù)的聚類分析。表2記錄了人體下肢各關(guān)節(jié)角度數(shù)據(jù)的具體分類情況,包括各關(guān)節(jié)類別,聚類過程中的迭代次數(shù)、狀態(tài)、聚類中心、單個誤差平方和以及整體誤差平方和。其中單個誤差平方和可以反映每一個聚類結(jié)果的差異性,整體誤差平方和可以判斷不同聚類結(jié)果的樣本差異。
3.2 功能分類算法驗證
對下肢每一個關(guān)節(jié)測試并記錄正常、欠屈和欠伸3種情況下各10組數(shù)據(jù),一共30組數(shù)據(jù),在聚類結(jié)果數(shù)據(jù)中輸入測試的數(shù)據(jù),完成對數(shù)據(jù)的歸類,驗證測試結(jié)果如表3所示。
表3中Ai~Ii分別為各個關(guān)節(jié)角度數(shù)據(jù)的測試結(jié)果情況,其中i={1,2,3,4,5,6}分別對應了左髖關(guān)節(jié)、右髖關(guān)節(jié)、左膝關(guān)節(jié)、右膝關(guān)節(jié)、左踝關(guān)節(jié)、右踝關(guān)節(jié)的數(shù)據(jù)測試結(jié)果,如A1~I1代表的是左髖關(guān)節(jié)的數(shù)據(jù)測試結(jié)果。結(jié)果表明10個伸展困難的樣本中,有7個被正確歸類,2個被錯誤地歸類到正常,1個被錯誤地歸類到屈曲困難;10個正常的測試樣本中,有9個被正確歸類,1個被錯誤地歸類到屈曲困難;10個屈曲困難的測試樣本完全正確歸類于屈曲困難類。
結(jié)合聚類評估原理之中涉及的計算方法,可以計算相應關(guān)節(jié)的聚類指標,各關(guān)節(jié)聚類評估指標結(jié)果如表4所示。
由表4可知,利用K-means聚類算法能較好地完成骨關(guān)節(jié)空間角度數(shù)據(jù)的聚類和樣本歸類,以左髖關(guān)節(jié)數(shù)據(jù)的聚類分析和指標計算為例,3個聚類指標的值分別為0.867、0.837、0.86,即該關(guān)節(jié)數(shù)據(jù)的聚類準確度穩(wěn)定在0.83以上。就整個系統(tǒng)來而言,所有關(guān)節(jié)聚類評估準確度能穩(wěn)定在0.8及其以上。
3.3 系統(tǒng)測試
系統(tǒng)最終以Unity為開發(fā)平臺,運用C#編程方式,連接Kinect,實時采集數(shù)據(jù)并進行骨關(guān)節(jié)功能評估,圖5所示為系統(tǒng)測試界面圖。
圖5所示的系統(tǒng)測試結(jié)果圖中包括Kinect圖像實時采集與顯示,以及下肢左右髖關(guān)節(jié)、膝關(guān)節(jié)、踝關(guān)節(jié)空間角度顯示與功能評估等部分。系統(tǒng)能夠?qū)崟r獲取關(guān)節(jié)點之間的角度并進行最大屈伸數(shù)據(jù)的更新記錄,適用于多種參數(shù)的系統(tǒng)聚類,且可視化強,能夠非接觸地完成骨關(guān)節(jié)空間角度功能的數(shù)字化評估。
4 結(jié)論
本文研究了一種非接觸的骨關(guān)節(jié)空間角度測量方法,并介紹了骨關(guān)節(jié)空間角度數(shù)字化評估系統(tǒng)的結(jié)構(gòu)和工作原理。實驗結(jié)果證明:系統(tǒng)以可視化方式為用戶呈現(xiàn)骨關(guān)節(jié)功能的數(shù)字化分析與功能評估結(jié)果;聚類評估指標purity、RI與F-meaures等能穩(wěn)定在0.8及其以上,基本滿足骨關(guān)節(jié)空間角度數(shù)字化評估要求,可以為醫(yī)務人員在診斷、治療方案確立、治療前后功能對比評價以及康復指導等過程提供一種更客觀有效的依據(jù)。
-
控制系統(tǒng)
+關(guān)注
關(guān)注
41文章
6604瀏覽量
110576 -
數(shù)據(jù)采集
+關(guān)注
關(guān)注
38文章
6053瀏覽量
113620 -
數(shù)字化
+關(guān)注
關(guān)注
8文章
8708瀏覽量
61726
原文標題:【學術(shù)論文】骨關(guān)節(jié)角度數(shù)字化測量與功能評估系統(tǒng)設計
文章出處:【微信號:ChinaAET,微信公眾號:電子技術(shù)應用ChinaAET】歡迎添加關(guān)注!文章轉(zhuǎn)載請注明出處。
發(fā)布評論請先 登錄
相關(guān)推薦
評論