CN112617816A - 一种基于运动学参数定位功能性肘关节旋转中心方法 - Google Patents
一种基于运动学参数定位功能性肘关节旋转中心方法 Download PDFInfo
- Publication number
- CN112617816A CN112617816A CN202011495353.5A CN202011495353A CN112617816A CN 112617816 A CN112617816 A CN 112617816A CN 202011495353 A CN202011495353 A CN 202011495353A CN 112617816 A CN112617816 A CN 112617816A
- Authority
- CN
- China
- Prior art keywords
- point
- reference axis
- upper arm
- vector
- mark point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Granted
Links
- 238000000034 method Methods 0.000 title claims abstract description 108
- 210000002310 elbow joint Anatomy 0.000 title claims abstract description 78
- 210000000245 forearm Anatomy 0.000 claims abstract description 74
- 238000007781 pre-processing Methods 0.000 claims abstract description 6
- 238000012360 testing method Methods 0.000 claims description 37
- 210000002758 humerus Anatomy 0.000 claims description 22
- 239000003550 marker Substances 0.000 claims description 16
- 238000005457 optimization Methods 0.000 claims description 15
- 238000013519 translation Methods 0.000 claims description 12
- 238000012937 correction Methods 0.000 claims description 5
- 238000012545 processing Methods 0.000 claims description 5
- 210000002659 acromion Anatomy 0.000 claims description 4
- 210000000707 wrist Anatomy 0.000 claims description 4
- 238000004364 calculation method Methods 0.000 abstract description 3
- 230000036544 posture Effects 0.000 description 33
- 210000000323 shoulder joint Anatomy 0.000 description 21
- 230000000877 morphologic effect Effects 0.000 description 11
- 238000002474 experimental method Methods 0.000 description 7
- 210000002683 foot Anatomy 0.000 description 7
- 210000000038 chest Anatomy 0.000 description 6
- 210000003414 extremity Anatomy 0.000 description 6
- 230000005021 gait Effects 0.000 description 6
- 238000007796 conventional method Methods 0.000 description 5
- 210000001503 joint Anatomy 0.000 description 4
- 210000002414 leg Anatomy 0.000 description 4
- 210000000689 upper leg Anatomy 0.000 description 4
- 210000004247 hand Anatomy 0.000 description 3
- 238000005259 measurement Methods 0.000 description 3
- 238000012986 modification Methods 0.000 description 3
- 230000004048 modification Effects 0.000 description 3
- 210000001991 scapula Anatomy 0.000 description 3
- 230000002194 synthesizing effect Effects 0.000 description 3
- 101100155017 Arabidopsis thaliana RCE2 gene Proteins 0.000 description 2
- 102100022359 CAAX prenyl protease 2 Human genes 0.000 description 2
- 101000824514 Homo sapiens CAAX prenyl protease 2 Proteins 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- PCHJSUWPFVWCPO-UHFFFAOYSA-N gold Chemical compound [Au] PCHJSUWPFVWCPO-UHFFFAOYSA-N 0.000 description 2
- 238000003384 imaging method Methods 0.000 description 2
- 238000002595 magnetic resonance imaging Methods 0.000 description 2
- 239000000203 mixture Substances 0.000 description 2
- 238000002559 palpation Methods 0.000 description 2
- 238000001356 surgical procedure Methods 0.000 description 2
- 230000009897 systematic effect Effects 0.000 description 2
- 210000001364 upper extremity Anatomy 0.000 description 2
- 208000031872 Body Remains Diseases 0.000 description 1
- 241001653121 Glenoides Species 0.000 description 1
- 238000012307 MRI technique Methods 0.000 description 1
- 206010037660 Pyrexia Diseases 0.000 description 1
- 241000746998 Tragus Species 0.000 description 1
- 210000003423 ankle Anatomy 0.000 description 1
- 230000015572 biosynthetic process Effects 0.000 description 1
- 230000007423 decrease Effects 0.000 description 1
- 230000004064 dysfunction Effects 0.000 description 1
- 230000000694 effects Effects 0.000 description 1
- 210000000527 greater trochanter Anatomy 0.000 description 1
- 238000009499 grossing Methods 0.000 description 1
- 210000004095 humeral head Anatomy 0.000 description 1
- 210000000629 knee joint Anatomy 0.000 description 1
- 210000003205 muscle Anatomy 0.000 description 1
- 210000004197 pelvis Anatomy 0.000 description 1
- 238000002360 preparation method Methods 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000005070 sampling Methods 0.000 description 1
- 238000012216 screening Methods 0.000 description 1
- 238000003786 synthesis reaction Methods 0.000 description 1
- 238000012956 testing procedure Methods 0.000 description 1
- 210000000115 thoracic cavity Anatomy 0.000 description 1
- 238000010200 validation analysis Methods 0.000 description 1
- 238000010792 warming Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1118—Determining activity level
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1116—Determining posture transitions
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/112—Gait analysis
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1121—Determining geometric values, e.g. centre of rotation or angular range of movement
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/103—Detecting, measuring or recording devices for testing the shape, pattern, colour, size or movement of the body or parts thereof, for diagnostic purposes
- A61B5/11—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb
- A61B5/1126—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique
- A61B5/1128—Measuring movement of the entire body or parts thereof, e.g. head or hand tremor, mobility of a limb using a particular sensing technique using image analysis
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06V—IMAGE OR VIDEO RECOGNITION OR UNDERSTANDING
- G06V40/00—Recognition of biometric, human-related or animal-related patterns in image or video data
- G06V40/20—Movements or behaviour, e.g. gesture recognition
-
- G—PHYSICS
- G16—INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR SPECIFIC APPLICATION FIELDS
- G16H—HEALTHCARE INFORMATICS, i.e. INFORMATION AND COMMUNICATION TECHNOLOGY [ICT] SPECIALLY ADAPTED FOR THE HANDLING OR PROCESSING OF MEDICAL OR HEALTHCARE DATA
- G16H20/00—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance
- G16H20/30—ICT specially adapted for therapies or health-improving plans, e.g. for handling prescriptions, for steering therapy or for monitoring patient compliance relating to physical therapies or activities, e.g. physiotherapy, acupressure or exercising
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- General Health & Medical Sciences (AREA)
- Physics & Mathematics (AREA)
- Public Health (AREA)
- Medical Informatics (AREA)
- Biophysics (AREA)
- Veterinary Medicine (AREA)
- Molecular Biology (AREA)
- Physiology (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Animal Behavior & Ethology (AREA)
- Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Heart & Thoracic Surgery (AREA)
- Computer Vision & Pattern Recognition (AREA)
- General Physics & Mathematics (AREA)
- Social Psychology (AREA)
- Human Computer Interaction (AREA)
- Multimedia (AREA)
- Psychiatry (AREA)
- Theoretical Computer Science (AREA)
- Geometry (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Radiology & Medical Imaging (AREA)
- Physical Education & Sports Medicine (AREA)
- Epidemiology (AREA)
- Primary Health Care (AREA)
- Measurement Of The Respiration, Hearing Ability, Form, And Blood Characteristics Of Living Organisms (AREA)
Abstract
本发明提供一种基于运动学参数定位功能性肘关节旋转中心方法,包括获取被试者的身体参数和数据预处理,还包括以下步骤:修正参考轴线方向;参考轴向量高度补偿;定位功能性肘关节旋转中心。本发明提出基于运动学参数定位功能性肘关节旋转中心方法,定位功能性肘关节旋转中心更加符合肘关节解剖学知识,在上臂与前臂两个方面精度更高,并且定位过程使上臂与前臂身体段落模型统一,使得模型中这两段身体段落首尾相连利于后续计算。
Description
技术领域
本发明涉及数字图像处理的技术领域,特别是一种基于运动学参数定位功能性肘关节旋转中心方法。
背景技术
依据身体段落运动学参数估计关节旋转中心位置的功能性方法是一种被广泛应用的分析手段,其在临床检测、步态分析、姿态采集等过程中非常实用。人体内部关节无法通过触诊方法被准确估计,即使检测员具有丰富***的解剖学知识,也会存在较大误差,尤其是在面对关节置换术术后患者时准确性更加难以保证。同时,金标准核磁成像手段虽然准确性最高,但成本高,操作复杂,每次扫描姿态单一。而功能性方法不需要***的解剖学知识,准确性也被学者Stokdijk等人通过测量尸体肩关节位置的方法被证明。功能性方法由于其实用性,被广泛应用于各种领域。
肘关节由于其远不如肩关节的旋转范围与灵活度,受到了很少的关注,甚至有学者认为功能性肩关节旋转中心定位方法可以直接推广到全身关节而将肘关节一带而过。但肘关节连接上臂与小臂,在上臂参考坐标系,功能性肘关节旋转中心是前臂运动轨迹中心,在前臂参考坐标系,功能性肘关节旋转中心是上臂运动轨迹旋转中心。并且,参考功能性肩关节旋转中心定位方法,两种参考系下求出的结果在全局坐标系下不是重合的。因此,尽管肘关节不如肩关节灵活,但通过移植功能性肩关节旋转中心定位方法来估计功能性肘关节旋转中心位置的方法是不合理的。
文献[C.A.C,L.D.G,A.J.A,and E.B.C,MRI development and validation oftwo new predictive methods of glenohumeral joint centre locationidentification and comparison with established techniques,Journal ofbiomechanic,vol.42,no.10,2009.],该文章提出了一种新的回归模型和简单的三维偏移量方法来估计肩肱关节的位置,采用易于管理的措施,并将估计值与磁共振成像(MRI)测量的已知GHJ位置进行了比较。该方法的缺点是当被试者躯干是运动的,如行走过程,MRI技术很难检测每一时刻的肩关节位置,功能性旋转中心也不能用一个点来代替。
文献[H.E.J.Veeger,"The position of the rotation center of theglenohumeral joint,"Journal of Biomechanics,vol.33,no.12,2000.]该文章公开了为了验证关节旋转中心可以根据关节的几何形状来描述的假设,比较了两种计算关节旋转中心的方法。这是基于瞬时螺旋轴计算的运动学估计,以及基于通过关节盂表面的球面拟合的几何估计。四个新鲜的尸体手臂固定在肩胛骨上,并安装了电磁传感器。每只手臂向不同方向移动,同时记录肱骨的方位。随后,对每个标本进行解剖,并对其肩胛盂和肱骨头表面进行数字化。结果表明,两种方法之间没有差异。该方法使用尸体验证了功能性方法的重复性与准确性,但尸体被动运动与人体主动运动时肢体末端可达域是不同的。
文献[K.Aoki,K.Kawachi,M.Kouchi,and M.Mochimaru,Functional jointrotation centers for whole body digital manikin,SAE Technical Papers vol.1,2008,doi:10.4271/2008-01-1859.]该文章提出了一种估计整个人体功能关节旋转中心的方法,利用这种方法,只要人体各节段的长度保持不变,就可以估计出运动中关节的位置,当数字模型的全身运动时,最远端的轨迹与实际人体的轨迹更接近。该方法的缺点是全身22个关节都采用肩关节功能性旋转中心定位方法进行定位,从而忽略了肘关节的重要性。
发明内容
为了解决上述的技术问题,本发明提出基于运动学参数定位功能性肘关节旋转中心方法,定位功能性肘关节旋转中心更加符合肘关节解剖学知识,在上臂与前臂两个方面精度更高,并且定位过程使上臂与前臂身体段落模型统一,使得模型中这两段身体段落首尾相连利于后续计算。
本发明提供一种基于运动学参数定位功能性肘关节旋转中心方法,获取被试者的身体参数和数据预处理,还包括以下步骤:
步骤1:修正参考轴线方向;
步骤2:参考轴向量高度补偿;
步骤3:定位功能性肘关节旋转中心。
优选的是,所述被试者的身体参数包括形态学参数和/或步态下相关运动学参数。
在上述任一方案中优选的是,所述形态学参数包括围度参数、宽度参数、高度参数和人体四肢的长度与高度。
在上述任一方案中优选的是,所述形态学参数的获取方法包括:
1)采用自然站立姿势测量所述围度参数和所述宽度参数,所述自然站立姿势为自然站立姿态下人体双脚打开与肩同宽,双手打开远离躯干;
2)采用立正姿势测量所述高度参数,所述立正姿势为人体双臂紧贴躯干,双腿并拢,双脚脚跟贴紧;
3)采用坐姿测量所述人体四肢的长度与高度,所述坐姿为通过调节升降椅的高矮,坐姿下的人体大腿与地面平行,右手手臂水平向前伸直,左手上臂垂直向下,小臂水平向前。
在上述任一方案中优选的是,所述步态下的相关运动学参数的获取方法为贴好标记点的受试者以自然站立姿势站在测试区域内部边缘处,约T初后开始行走,沿直线行走至测试区域另一边后缓慢转身,再沿原路返回行走,行走五次后站立在测试区域内部边缘,直至T末时结束,测试重复NT次结束,N个标记点的三维坐标T帧的数据被采集。
在上述任一方案中优选的是,所述数据预处理的方法为筛选一次实验中连续Ti帧数据,所述Ti帧数据连续且完整,对筛选出来的所述Ti帧数据进行平滑处理。
在上述任一方案中优选的是,所述步骤1包括以下子步骤:
步骤11:分别确定上臂中轴线和前臂圆柱体中轴线,以肱骨外上髁点为起点,肩峰点为终点的向量为上臂圆柱参考轴向量,为了方便后续的优化处理,将参考轴向量平移至起点与肱骨内外上髁点中点重合,记参考轴向量为以手腕标记点中点为起点,肱骨内外上髁点中点为终点的向量为前臂圆柱参考轴向量记为
在上述任一方案中优选的是,在任意测试时刻ta,系数C有以下性质:
其中,表示ta时刻上臂圆柱的参考轴向量,和分别表示ta时刻上臂第一标记点高度、第二标记点高度和第三标记点高度占参考轴向量总高度比例,和分别表示上臂第一标记点、第二标记点和第三标记点到参考轴向量的垂向量,和分别表示上臂第一标记点、第二标记点和第三标记点的位置坐标,表示肱骨内外上髁点的中点位置坐标,和分别表示ta时刻前臂第一标记点高度、第二标记点高度和第三标记点高度占参考轴向量总高度比例,前臂参考轴向量与系数CD同理。
在上述任一方案中优选的是,在整个测试时间段(ts-tp)内,同一标记点不同时刻的系数或不全相同,选取受试者站立姿态下双臂垂直向下的时刻ti时的系数为标准系数,对每一时刻参考轴向量加入方向补偿向量改变向量方向,通过修正参考轴向量方向的方式改变系数C,使得补偿后系数C保持恒定,并且与标准系数相等,
在上述任一方案中优选的是,选取受试者站立姿态下双臂垂直向下的时刻ti时上臂标记点到参考轴向量的距离为标准距离和对参考周向量的起点与终点同时修正,其中,MU表示肩峰点位置坐标,MD表示肱骨内外上髁点中点位置坐标。
在上述任一方案中优选的是,所述步骤2包括对前臂参考轴向量进行伸缩补偿,补偿过程如下式
在上述任一方案中优选的是,所述步骤3包括采用上臂圆柱参考轴线的起点和前臂圆柱参考轴线的终点来计算功能性肘关节旋转中心,在测试时间内任意时刻ta,上臂圆柱参考轴线起点为 前臂参考轴线终点为 其中,其中,Melb表示......,表示......,表示......,表示......,MD表示肱骨内外上髁点中点,Mwrist表示桡骨茎突与尺骨茎突中点。
在上述任一方案中优选的是,所述步骤3还包括对Elb1参考肘关节1进行优化,通过对上臂参考轴向量进行伸缩与平移减小误差的同时,使参考肘关节1向Elb2参考肘关节2靠近。
在上述任一方案中优选的是,在所述优化的过程中,高度改变对于上臂三标记点是同向的,而距离改变则是不完全同向的,将上臂参考轴向量最终平移量由上臂三标记点垂线的合成来表示。在测试时间内任意时刻ta,上臂参考轴平移量可表示为:
在上述任一方案中优选的是,所述优化的过程如下:
其中,表示ta时刻上臂第一标记点在上臂参考轴向量中高度,表示ta时刻高度变化,m1、m2和m3分别表示参考肘关节Elb1到上臂第一标记点、上臂第二标记点和上臂第三标记点距离的平均值,表示ta时刻修正后前臂参考轴线终点。
本发明提供了一种基于运动学参数定位功能性肘关节旋转中心方法,使得定位结果更符合解剖学知识,也保证了算法的重复性。
附图说明
图1为按照本发明的基于运动学参数定位功能性肘关节旋转中心方法的一优选实施例的流程图。
图2为按照本发明的基于运动学参数定位功能性肘关节旋转中心方法的标记点距离变化标准差的一实施例的示意图。
图3为按照本发明的基于运动学参数定位功能性肘关节旋转中心方法与传统方法精度对比的一优选实施例的示意图。
具体实施方式
下面结合附图和具体的实施例对本发明做进一步的阐述。
实施例一
如图1所示,执行步骤100,获取被试者的身体参数和数据预处理,被试者的身体参数包括形态学参数和/或步态下相关运动学参数。形态学参数包括围度参数、宽度参数、高度参数和人体四肢的长度与高度。形态学参数的获取方法包括:1)采用自然站立姿势测量所述围度参数和所述宽度参数,所述自然站立姿势为自然站立姿态下人体双脚打开与肩同宽,双手打开远离躯干;2)采用立正姿势测量所述高度参数,所述立正姿势为人体双臂紧贴躯干,双腿并拢,双脚脚跟贴紧;3)采用坐姿测量所述人体四肢的长度与高度,所述坐姿为通过调节升降椅的高矮,坐姿下的人体大腿与地面平行,右手手臂水平向前伸直,左手上臂垂直向下,小臂水平向前。步态下的相关运动学参数的获取方法为贴好标记点的受试者以自然站立姿势站在测试区域内部边缘处,约T初后开始行走,沿直线行走至测试区域另一边后缓慢转身,再沿原路返回行走,行走五次后站立在测试区域内部边缘,直至T末时结束,测试重复NT次结束,N个标记点的三维坐标T帧的数据被采集。
数据预处理的方法为筛选一次实验中连续Ti帧数据,所述Ti帧数据连续且完整,对筛选出来的所述Ti帧数据进行平滑处理。
执行步骤110,修正参考轴线方向。分别确定上臂中轴线和前臂圆柱体中轴线,以肱骨外上髁点为起点,肩峰点为终点的向量为上臂圆柱参考轴向量,为了方便后续的优化处理,将参考轴向量平移至起点与肱骨内外上髁点中点重合,记参考轴向量为以手腕标记点中点为起点,肱骨内外上髁点中点为终点的向量为前臂圆柱参考轴向量记为引入系数C对所述上臂中轴线和前臂圆柱体中轴线进行修正,所述系数C包括上臂标记点到参考轴向量的垂足在向量中高度占比系数CU和前臂标记点到参考轴向量的垂足在向量中高度占比系数CD。在任意测试时刻ta,系数C有以下性质:
其中,表示ta时刻上臂圆柱的参考轴向量,和分别表示ta时刻上臂第一标记点高度、第二标记点高度和第三标记点高度占参考轴向量总高度比例,和分别表示上臂第一标记点、第二标记点和第三标记点到参考轴向量的垂向量,和分别表示上臂第一标记点、第二标记点和第三标记点的位置坐标,表示肱骨内外上髁点的中点位置坐标,和分别表示ta时刻前臂第一标记点高度、第二标记点高度和第三标记点高度占参考轴向量总高度比例,前臂参考轴向量与系数CD同理。
在整个测试时间段(ts-tp)内,同一标记点不同时刻的系数或不全相同,选取受试者站立姿态下双臂垂直向下的时刻ti时的系数为标准系数,对每一时刻参考轴向量加入方向补偿向量改变向量方向,通过修正参考轴向量方向的方式改变系数C,使得补偿后系数C保持恒定,并且与标准系数相等,
选取受试者站立姿态下双臂垂直向下的时刻ti时上臂标记点到参考轴向量的距离为标准距离和对参考周向量的起点与终点同时修正,其中,MU表示肩峰点位置坐标,MD表示肱骨内外上髁点中点位置坐标。任意时刻ta下,上臂第一标记点到参考轴向量的距离为
执行步骤120,参考轴向量高度补偿。对前臂参考轴向量进行伸缩补偿,补偿过程如下式
执行步骤130,定位功能性肘关节旋转中心。采用上臂圆柱参考轴线的起点和前臂圆柱参考轴线的终点来计算功能性肘关节旋转中心,在测试时间内任意时刻ta,上臂圆柱参考轴线起点为 前臂参考轴线终点为 其中,Melb表示......,表示......,表示......,表示......,MD表示肱骨内外上髁点中点,Mwrist表示桡骨茎突与尺骨茎突中点。
对Elb1参考肘关节1进行优化,通过对上臂参考轴向量进行伸缩与平移减小误差的同时,使参考肘关节1向Elb2参考肘关节2靠近。在所述优化的过程中,高度改变对于上臂三标记点是同向的,而距离改变则是不完全同向的,将上臂参考轴向量最终平移量由上臂三标记点垂线的合成来表示。在测试时间内任意时刻ta,上臂参考轴平移量可表示为:
所述优化的过程如下:
其中,表示ta时刻上臂第一标记点在上臂参考轴向量中高度,表示ta时刻高度变化,m1、m2和m3分别表示参考肘关节Elb1到上臂第一标记点、上臂第二标记点和上臂第三标记点距离的平均值,表示ta时刻修正后前臂参考轴线终点。
实施例二
依据身体段落运动学参数估计关节旋转中心位置的功能性方法是一种被广泛应用的分析手段,其在临床检测、步态分析、姿态采集等过程中非常实用。人体内部关节无法通过触诊方法被准确估计,即使检测员具有丰富***的解剖学知识,也会存在较大误差,尤其是在面对关节置换术术后患者时准确性更加难以保证。同时,金标准核磁成像手段虽然准确性最高,但成本高,操作复杂,每次扫描姿态单一。而功能性方法不需要***的解剖学知识,准确性也被学者Stokdijk等人通过测量尸体肩关节位置的方法被证明。功能性方法由于其实用性,被广泛应用于各种领域。
肘关节由于其远不如肩关节的旋转范围与灵活度,受到了很少的关注,甚至有学者认为功能性肩关节旋转中心定位方法可以直接推广到全身关节而将肘关节一带而过。但肘关节连接上臂与小臂,在上臂参考坐标系,功能性肘关节旋转中心是前臂运动轨迹中心,在前臂参考坐标系,功能性肘关节旋转中心是上臂运动轨迹旋转中心。并且,参考功能性肩关节旋转中心定位方法,两种参考系下求出的结果在全局坐标系下不是重合的。因此,尽管肘关节不如肩关节灵活,但通过移植功能性肩关节旋转中心定位方法来估计功能性肘关节旋转中心位置的方法是不合理的。基于以上原因,改进上臂刚性圆柱体模型,设计了估计上臂与前臂运动轨迹共同中心的功能性肘关节旋转中心(RCE,RotationCenterofElbow)定位方法。
选择28名健康成年男性(18-55岁,形态学参数如表1所示)受试者加入了本次测试。
表1
在实际测量过程中,使用了AnthroscanBodyscan三维人体扫描仪(德国Humansolutions公司)与Qualisys三维运动采集与分析***(由运动捕捉摄像机、分析软件、获取单元、校准设备、标记球和设备固定装置组成,瑞典Qualisys公司)。实验中,测量人体形态学参数的三维扫描仪测量误差为0.5mm,测量运动中人体运动学参数的运动捕捉***包括17个摄像头:4个视频摄像头,13个捕捉摄像头均匀分布在实验场地四周。测量摄像头采样频率100Hz,每次实验的校准精度均保持在0.7mm以下。
1、准备阶段
形态学参数进行测量前,检查受试者上肢关节情况,排除上肢功能障碍等情况。形态学参数测量结束后安排受试者热身,尤其激活肩部肌肉,准备活动强度与量以全身发热而不感到疲惫为宜。热身结束后受试者穿着紧身衣裤,粘贴共69个标记点,其中45个骨形标志点,24个非骨形标志点。
骨形标志点包括:头顶点,眉间点,耳屏点,颈点,肩峰点,肱骨内上髁点,肱骨外上髁点,桡骨茎突点,尺骨茎突点,手部二五指中点,胸上点,胸下点,胸下点对应高度的胸椎点,胸下点对应高度胸宽点,胸下点到脐点中线最远两点,髂前上棘点,髂后上棘点,大转子点,膝关节,脚踝,脚长中线一五趾中点,右肩胛骨。
非骨形标志点为上臂、小臂、大腿和小腿共8处身体段落上的8组点,每组三个标记点。这类标记点的粘贴位置是根据文献规定的两大原则确定的:(1)三点不能呈一条直线;(2)三点之间的距离应尽可能的大。在本实验中,三点的粘贴位置不仅要符合以上两条原则,还要保证其在身体段落横截面的投影尽可能地将截面周长均分。
2测试步骤
首先对被试者进行三维扫描以获得其形态学参数。在扫描区域内,被试者先后以三种姿势接受扫描——自然站立姿势、立正姿势与坐姿。其中自然站立姿态下人体双脚打开与肩同宽,双手打开远离躯干。因为这种姿态下人体的躯干不与双臂紧挨,双腿也不相互干扰,所以主要的围度参数与宽度参数是在这种姿态下测量的;立正姿态下人体双臂紧贴躯干,双腿并拢,双脚脚跟贴紧,这类姿态下人体骨骼笔直向上适合测量高度参数;通过调节升降椅的高矮,坐姿下的人体大腿与地面平行,右手手臂水平向前伸直,左手上臂垂直向下,小臂水平向前,用于测量人体四肢的长度与高度。
其次,使用Qual isys***对被试者步态下相关运动学参数进行测量。测试开始时,贴好标记点的受试者以自然站立姿势站在测试区域内部边缘处,约2s后开始行走,沿直线行走至测试区域另一边后缓慢转身,再沿原路返回行走,行走五次后站立在测试区域内部边缘,直至30s测试时间结束,测试重复三次结束。69个标记点的三维坐标3000帧的数据被采集。
3、数据预处理
Qual i sys三维运动采集与分析***输出了实验过程中粘贴在运动人体表面标记点每一帧的三维坐标(100Hz,共3000帧,30秒)。每人三组数据,每组数据包含69个标记点3000帧的三维坐标。在这三组数据中筛选出一次实验中连续的500帧数据,要求这500帧数据连续且完整,避免标记点位置采集过程中出现漏点的情况影响实验结果。对筛选出的数据做平滑处理。
4、结果
在上臂刚性圆柱体模型中,上臂被抽象成圆柱体表面的点,而功能性肩关节旋转中心在圆柱体顶面圆心上。因此功能性肩关节旋转中心到标记点距离由圆柱体半径与标记点到圆柱顶面的高决定,其中,圆柱半径体现上臂维度,由上臂三个标记点决定,而标记点到圆柱顶面的高可以通过伸缩圆柱体高度来改变,所以采用伸缩圆柱体高度的方法来补偿上臂维度变化以保持标记点到圆柱顶面圆心距离固定,补偿后圆柱顶面圆心即功能性肩关节旋转中心。
考虑到上臂自旋对上臂皮肤的影响较大,功能性肩关节旋转中心需要在上臂圆柱的中轴线上以减小上臂自旋对定位结果的影响。而前臂自旋对前臂皮肤影响较小,且肘关节旋转范围小,因此功能性肘关节旋转中心并不是严格要求在上臂与前臂圆柱体中轴线上的。同时,由于标记点到对应身体段落中轴线距离与标记点到身体段落终点的高共同决定了功能性肘关节旋转中心到标记点的距离,可以对两者同时进行优化最终得到一个重合的点,它到上臂与前臂标记点距离均保持固定,即功能性肘关节旋转中心。
4.1修正参考轴线方向
首先分别确定上臂与前臂圆柱体中轴线,只有当中轴线正确,圆柱的半径才能用以描述上臂或前臂的维度。以肱骨外上髁点为起点,肩峰点为终点的向量为上臂圆柱参考轴向量,为了方便后续的优化处理,将参考轴向量平移至起点与肱骨内外上髁点中点重合,记参考轴向量为以手腕标记点中点为起点,肱骨内外上髁点中点为终点的向量为前臂圆柱参考轴向量,记为由于参考轴向量与真实轴向量之间的误差,上臂或前臂三标记点到对应身体段落参考轴向量的距离不固定,标记点到参考轴向量的垂足切割向量长度的比例也不同。因此需要对参考周向量进行修正使之与真实轴向量尽可能贴近。
引入上臂与前臂标记点到各自参考轴向量的垂足在向量中高度占比系数CU与CD,在任意测试时刻ta,系数C有以下性质:
其中,表示ta时刻上臂圆柱的参考轴向量,表示ta时刻上臂第一标记点高度占参考轴向量总高度比例,表示上臂第一标记点到参考轴向量的垂向量,表示第一标记点的位置坐标,表示肱骨内外上髁点的中点位置坐标。前臂参考轴向量与系数CD同理。
以最优解为上臂新参考轴向量,继续修正使得上臂标记点到参考轴的距离固定。同样地,选取受试者站立姿态下双臂垂直向下时(假设为时刻ti时)上臂标记点到参考轴向量的距离为标准距离 和由于这次修正只改变参考轴向量的位置而不改变方向,需要对参考周向量的起点与终点同时修正,
公式描述了任意时刻ta下,标记点1到参考轴向量的距离:
3.2参考轴向量高度补偿
在贴点过程中,为了保证肩关节估计精度,上臂标记点粘贴位置距离肘关节更近,因此仅对前臂参考轴向量进行伸缩补偿,补偿过程如式。
3.3定位RCE
考虑到肘关节位于上臂末端和前臂顶端附近,算法采用上臂圆柱参考轴线的起点和前臂圆柱参考轴线的终点来计算功能性肘关节旋转中心。在测试时间内任意时刻Ta,上臂圆柱参考轴线起点为前臂参考轴线终点为 其中MD表示肱骨内外上髁点中点,Mwrist表示桡骨茎突与尺骨茎突中点。
由于上臂参考轴没有做伸缩补偿处理,因此对Elb1参考肘关节1进行优化,通过对上臂参考轴向量进行伸缩与平移减小误差的同时,使参考肘关节1向Elb2参考肘关节2靠近。最终优化结果为功能性肘关节旋转中心估计结果。根据勾股定理,肘关节旋转中心到上臂标记点的距离可以表示为标记点到参考轴向量距离与标记点所在高度的平方和开根,因此只需满足不同时刻距离与高度的平方和固定即可。
在优化过程中,高度改变对于上臂三标记点是同向的,即要么共同增加要么共同减小。而距离改变则是不完全同向的,考虑到参考轴向量沿标记点垂线方向平移改变距离最快,将上臂参考轴向量最终平移量由上臂三标记点垂线的合成来表示。在测试时间内任意时刻ta,上臂参考轴平移量可表示为:
优化过程如下:
其中表示ta时刻上臂第一标记点在上臂参考轴向量中高度,表示ta时刻高度变化,m1表示参考肘关节Elb1到上臂第一标记点距离的平均值。上式描述了优化过程:优化结果到上臂标记点距离固定,优化结果与参考肘关节重合。
在5s测试时间内,二十八位被试者每人500帧功能性肘关节旋转中心通过算法被估计,功能性肘关节旋转中心到上臂与前臂标记点距离标准差的平均值为:到上臂第一标记点0.8666[mm],到上臂第二标记点0.3807[mm],到上臂第三标记点1.5473[mm],到前臂第一标记点2.0059[mm],到前臂第二标记点1.4804[mm],到前臂第三标记点1.767[mm]。图2展示了二十八位受试者功能性肘关节旋转中心估计结果到上臂与前臂标记点距离的标准差。
由于肘关节灵活度远小于肩关节,因此功能性肘关节旋转中心估计方法通常沿用了功能性肩关节旋转中心估计方法,但在几何上,肘关节不仅是前臂运动轨迹中心,也是上臂运动轨迹中心,功能性肩关节旋转中心估计方法往往不能兼顾上臂与前臂,将肘关节视为前臂运动轨迹中心与上臂运动轨迹中心往往会得到不同结果,并且这两个结果对前臂与上臂的精度也不同。算法研究了一种同时满足上臂与前臂运动轨迹中心的功能性肘关节旋转中心估计法。以一号被试者为例,图3展示了这种方法与传统方法精度对比,图3左展示了传统方法将功能性肘关节旋转中心视为前臂运动轨迹中心求解结果到上臂与前臂标记点距离变化,图3中展示了传统方法将功能性肘关节旋转中心视为上臂运动轨迹中心求解结果到上臂与前臂标记点距离变化,图3右展示了新方法估计的功能性肘关节旋转中心到上臂与前臂标记点的距离变化。定位功能性肘关节旋转中心RCE(Rotation Center ofElbow)新方法与传统方法估计结果精度对比。RCE1,传统方法中肘关节与前臂形成球杆模型,对球心求解结果;RCE2,传统方法中肘关节与上臂形成球杆模型,对球心求解结果;NewRCE,新方法中肘关节与前臂与上臂形成杆球杆模型,对球心求解结果。图上半部展示了不同方法求解的RCE到上臂三标记点距离标准差,下半部展示了到前臂三标记点距离标准差。由于原理不同,导致了RCE1到上臂精度高而对前臂精度低,RCE2到上臂精度低而对前臂精度高,新方法精度更平均。
可以看出,传统方法要么到前臂标记点距离标准差小而到上臂误差大,要么到前臂误差大而到上臂误差小。综合上臂与前臂来看,本申请中的算法误差最小。
为了更好地理解本发明,以上结合本发明的具体实施例做了详细描述,但并非是对本发明的限制。凡是依据本发明的技术实质对以上实施例所做的任何简单修改,均仍属于本发明技术方案的范围。本说明书中每个实施例重点说明的都是与其它实施例的不同之处,各个实施例之间相同或相似的部分相互参见即可。对于***实施例而言,由于其与方法实施例基本对应,所以描述的比较简单,相关之处参见方法实施例的部分说明即可。
可能以许多方式来实现本发明的方法、装置和***。例如,可通过软件、硬件、固件或者软件、硬件、固件的任何组合来实现本发明的方法和***。用于所述方法的步骤的上述顺序仅是为了进行说明,本发明的方法的步骤不限于以上具体描述的顺序,除非以其它方式特别说明。此外,在一些实施例中,还可将本发明实施为记录在记录介质中的程序,这些程序包括用于实现根据本发明的方法的机器可读指令。因而,本发明还覆盖存储用于执行根据本发明的方法的程序的记录介质。
本发明的描述是为了示例和描述起见而给出的,而并不是无遗漏的或者将本发明限于所公开的形式。很多修改和变化对于本领域的普通技术人员而言是显然的。选择和描述实施例是为了更好说明本发明的原理和实际应用,并且使本领域的普通技术人员能够理解本发明从而设计适于特定用途的带有各种修改的各种实施例。
Claims (10)
1.一种基于运动学参数定位功能性肘关节旋转中心方法,获取被试者的身体参数和数据预处理,其特征在于,还包括以下步骤:
步骤1:修正参考轴线方向;
步骤2:参考轴向量高度补偿;
步骤3:定位功能性肘关节旋转中心。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011495353.5A CN112617816B (zh) | 2020-12-17 | 2020-12-17 | 一种基于运动学参数定位功能性肘关节旋转中心方法 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN202011495353.5A CN112617816B (zh) | 2020-12-17 | 2020-12-17 | 一种基于运动学参数定位功能性肘关节旋转中心方法 |
Publications (2)
Publication Number | Publication Date |
---|---|
CN112617816A true CN112617816A (zh) | 2021-04-09 |
CN112617816B CN112617816B (zh) | 2021-09-21 |
Family
ID=75316370
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN202011495353.5A Active CN112617816B (zh) | 2020-12-17 | 2020-12-17 | 一种基于运动学参数定位功能性肘关节旋转中心方法 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN112617816B (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN115530813A (zh) * | 2022-10-20 | 2022-12-30 | 吉林大学 | 用于人体上身多关节三维运动测试分析的标记*** |
WO2024051735A1 (zh) * | 2022-09-10 | 2024-03-14 | 中南大学湘雅二医院 | 骨骼相对运动测量方法、关节旋转轴获取及定位方法 |
Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0345866A1 (en) * | 1988-06-07 | 1989-12-13 | Necchi Societa Per Azioni | Machine for supplying and sewing overlapped fabric layers |
CN1439352A (zh) * | 2003-03-11 | 2003-09-03 | 刘佰林 | 康复肢体运动功能用机械肢体组合系列机器人 |
EP2065771A2 (en) * | 2007-11-30 | 2009-06-03 | Fanuc Ltd | Numerical controller for controlling a five-axis machining apparatus |
DE102012016323A1 (de) * | 2012-08-10 | 2014-04-03 | Herbert Lutz | Meßvorrichtung |
CN104147770A (zh) * | 2014-07-24 | 2014-11-19 | 燕山大学 | 基于惯性传感器可穿戴式偏瘫康复设备及捷联姿态算法 |
CN105526848A (zh) * | 2016-01-11 | 2016-04-27 | 北京联合大学 | 一种姿态辅助测量装置及测量方法 |
CN108030496A (zh) * | 2017-12-02 | 2018-05-15 | 北京工业大学 | 一种人体上肢肩部盂肱关节旋转中心与上臂抬升角耦合关系测量方法 |
CN108836346A (zh) * | 2018-04-16 | 2018-11-20 | 大连理工大学 | 一种基于惯性传感器的人体步态分析方法和*** |
DE102018004119A1 (de) * | 2018-05-23 | 2019-11-28 | Jürgen Löschinger | Verfahren zur Erzeugung von Erinnerungsmeldungen oder Warnhinweisen aufgrund der Ähnlichkeit einer beliebigen Bewegung zu einer beliebigen, aufgezeichneten Referenzbewegung |
TWM604175U (zh) * | 2020-08-26 | 2020-11-21 | 王志堅 | 人工膝關節置換術之運動軸*** |
-
2020
- 2020-12-17 CN CN202011495353.5A patent/CN112617816B/zh active Active
Patent Citations (10)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
EP0345866A1 (en) * | 1988-06-07 | 1989-12-13 | Necchi Societa Per Azioni | Machine for supplying and sewing overlapped fabric layers |
CN1439352A (zh) * | 2003-03-11 | 2003-09-03 | 刘佰林 | 康复肢体运动功能用机械肢体组合系列机器人 |
EP2065771A2 (en) * | 2007-11-30 | 2009-06-03 | Fanuc Ltd | Numerical controller for controlling a five-axis machining apparatus |
DE102012016323A1 (de) * | 2012-08-10 | 2014-04-03 | Herbert Lutz | Meßvorrichtung |
CN104147770A (zh) * | 2014-07-24 | 2014-11-19 | 燕山大学 | 基于惯性传感器可穿戴式偏瘫康复设备及捷联姿态算法 |
CN105526848A (zh) * | 2016-01-11 | 2016-04-27 | 北京联合大学 | 一种姿态辅助测量装置及测量方法 |
CN108030496A (zh) * | 2017-12-02 | 2018-05-15 | 北京工业大学 | 一种人体上肢肩部盂肱关节旋转中心与上臂抬升角耦合关系测量方法 |
CN108836346A (zh) * | 2018-04-16 | 2018-11-20 | 大连理工大学 | 一种基于惯性传感器的人体步态分析方法和*** |
DE102018004119A1 (de) * | 2018-05-23 | 2019-11-28 | Jürgen Löschinger | Verfahren zur Erzeugung von Erinnerungsmeldungen oder Warnhinweisen aufgrund der Ähnlichkeit einer beliebigen Bewegung zu einer beliebigen, aufgezeichneten Referenzbewegung |
TWM604175U (zh) * | 2020-08-26 | 2020-11-21 | 王志堅 | 人工膝關節置換術之運動軸*** |
Non-Patent Citations (5)
Title |
---|
AIYANGAR, A 等: "Instantaneous centers of rotation for lumbar segmental extension in vivo", 《JOURNAL OF BIOMECHANICS》 * |
DE SMET, ARNE; VERREWAERE, DRIES; SYS, GWEN: "Enhancing rotational placement of reconstruction prostheses of the distal femur after sarcoma resection", 《MEDICAL ENGINEERING & PHYSICS 》 * |
DUAN X H , ALLEN R H , SUN J Q .: "A stiffness-varying model of human gait", 《MEDICAL ENGINEERING & PHYSICS》 * |
孙嘉琪 等: "捏力测量方法及影响因素研究进展", 《体育科研》 * |
曹天一 等: "可活动式外固定治疗肘关节结核", 《中国骨伤》 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
WO2024051735A1 (zh) * | 2022-09-10 | 2024-03-14 | 中南大学湘雅二医院 | 骨骼相对运动测量方法、关节旋转轴获取及定位方法 |
CN115530813A (zh) * | 2022-10-20 | 2022-12-30 | 吉林大学 | 用于人体上身多关节三维运动测试分析的标记*** |
CN115530813B (zh) * | 2022-10-20 | 2024-05-10 | 吉林大学 | 用于人体上身多关节三维运动测试分析的标记*** |
Also Published As
Publication number | Publication date |
---|---|
CN112617816B (zh) | 2021-09-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
Stagni et al. | Quantification of soft tissue artefact in motion analysis by combining 3D fluoroscopy and stereophotogrammetry: a study on two subjects | |
Scheys et al. | Calculating gait kinematics using MR-based kinematic models | |
CN112617816B (zh) | 一种基于运动学参数定位功能性肘关节旋转中心方法 | |
Hicks et al. | Clinical applicability of using spherical fitting to find hip joint centers | |
Clément et al. | Soft tissue artifact compensation in knee kinematics by multi-body optimization: Performance of subject-specific knee joint models | |
Kuo et al. | Influence of soft tissue artifacts on the calculated kinematics and kinetics of total knee replacements during sit-to-stand | |
AU2001290888B8 (en) | Assessing the condition of a joint and devising treatment | |
Tsai et al. | A novel dual fluoroscopic imaging method for determination of THA kinematics: in-vitro and in-vivo study | |
Solav et al. | A framework for measuring the time-varying shape and full-field deformation of residual limbs using 3-D digital image correlation | |
Monnet et al. | Comparison of the SCoRE and HA methods for locating in vivo the glenohumeral joint centre | |
Murakami et al. | In vivo kinematics of healthy male knees during squat and golf swing using image-matching techniques | |
Sangeux et al. | Quantification of the 3D relative movement of external marker sets vs. bones based on magnetic resonance imaging | |
Camomilla et al. | Soft tissue displacement over pelvic anatomical landmarks during 3-D hip movements | |
WO2015162158A1 (en) | Human motion tracking | |
US20200205900A1 (en) | Dynamic 3d motion capture for surgical implant orientation | |
JP2002532126A (ja) | 関節状態の評価及び損傷防止装置 | |
El Habachi et al. | A parallel mechanism of the shoulder—application to multi-body optimisation | |
Cereatti et al. | Scapular motion tracking using acromion skin marker cluster: in vitro accuracy assessment | |
WO2022116411A1 (zh) | 一种人体功能性关节旋转中心检测与定位分析方法 | |
Peters et al. | Determination of the optimal locations of surface-mounted markers on the tibial segment | |
CN112535475B (zh) | 基于刚性上臂模型的功能性肩关节旋转中心定位方法 | |
WO2024051735A1 (zh) | 骨骼相对运动测量方法、关节旋转轴获取及定位方法 | |
Donati et al. | Enhanced anatomical calibration in human movement analysis | |
Uzuner et al. | Dual fluoroscopic evaluation of human tibiofemoral joint kinematics during a prolonged standing: A pilot study | |
Lalone et al. | Effect of radial head implant shape on radiocapitellar joint congruency |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
GR01 | Patent grant | ||
GR01 | Patent grant |