[go: up one dir, main page]

KR101886650B1 - 무릎인대 물성 도출방법 - Google Patents

무릎인대 물성 도출방법 Download PDF

Info

Publication number
KR101886650B1
KR101886650B1 KR1020170047941A KR20170047941A KR101886650B1 KR 101886650 B1 KR101886650 B1 KR 101886650B1 KR 1020170047941 A KR1020170047941 A KR 1020170047941A KR 20170047941 A KR20170047941 A KR 20170047941A KR 101886650 B1 KR101886650 B1 KR 101886650B1
Authority
KR
South Korea
Prior art keywords
knee
physical property
ligament
load
photographing
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.)
Active
Application number
KR1020170047941A
Other languages
English (en)
Inventor
고용곤
손주현
강경탁
Original Assignee
주식회사 티제이씨라이프
손주현
강경탁
Priority date (The priority date 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 date listed.)
Filing date
Publication date
Application filed by 주식회사 티제이씨라이프, 손주현, 강경탁 filed Critical 주식회사 티제이씨라이프
Priority to KR1020170047941A priority Critical patent/KR101886650B1/ko
Application granted granted Critical
Publication of KR101886650B1 publication Critical patent/KR101886650B1/ko
Active legal-status Critical Current
Anticipated expiration legal-status Critical

Links

Images

Classifications

    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B5/00Measuring for diagnostic purposes; Identification of persons
    • A61B5/45For evaluating or diagnosing the musculoskeletal system or teeth
    • A61B5/4533Ligaments
    • AHUMAN NECESSITIES
    • A61MEDICAL OR VETERINARY SCIENCE; HYGIENE
    • A61BDIAGNOSIS; SURGERY; IDENTIFICATION
    • A61B6/00Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
    • A61B6/02Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
    • A61B6/03Computed tomography [CT]
    • A61B6/032Transmission computed tomography [CT]

Landscapes

  • Health & Medical Sciences (AREA)
  • Life Sciences & Earth Sciences (AREA)
  • Engineering & Computer Science (AREA)
  • Medical Informatics (AREA)
  • Heart & Thoracic Surgery (AREA)
  • Surgery (AREA)
  • Physics & Mathematics (AREA)
  • Veterinary Medicine (AREA)
  • Biophysics (AREA)
  • Pathology (AREA)
  • Public Health (AREA)
  • Biomedical Technology (AREA)
  • General Health & Medical Sciences (AREA)
  • Animal Behavior & Ethology (AREA)
  • Molecular Biology (AREA)
  • Oral & Maxillofacial Surgery (AREA)
  • Rehabilitation Therapy (AREA)
  • Dentistry (AREA)
  • Rheumatology (AREA)
  • Orthopedic Medicine & Surgery (AREA)
  • Pulmonology (AREA)
  • Theoretical Computer Science (AREA)
  • High Energy & Nuclear Physics (AREA)
  • Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
  • Optics & Photonics (AREA)
  • Radiology & Medical Imaging (AREA)
  • Apparatus For Radiation Diagnosis (AREA)

Abstract

무릎인대 물성 도출방법이 개시된다. 본 발명의 실시예에 따른 무릎인대 물성 도출방법(S100)은, (a) 측정대상환자를 옆으로 눕힌 상태에서, 무릎을 특정 각도로 구부린 후, 무릎에 소정 크기의 하중을 특정 방향으로 가하는 하중인가단계(S110); (b) 특정 각도로 구부러진 무릎에 대해 CT 촬영을 실시하여 촬영데이터를 획득하는 촬영단계(S120); (c) 촬영단계로부터 획득한 데이터를 바탕으로 무릎에 대한 유한요소해석을 위한 3차원 모델을 생성시키고, 무릎인대의 물성값은 기공지된 평균 물성값으로 설정하는 무릎모델 생성단계(S130); (d) 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하는 민감부위 선별단계(S140); 및 (e) 하중인가단계에서 인가된 가압력 값과 촬영단계로부터 획득한 촬영데이터를 바탕으로 민감부위 선별단계로 통해 선별된 인대의 물성값을 최적화하는 물성 최적화단계(S150);를 포함하는 것을 구성의 요지로 한다.
본 발명에 따르면, 사체 실험을 통하지 않고 CT 장비와 TELOS 장비를 연동하여 특정하는 방법과 X-RAY-TELOS 측정을 하여 하중을 받지 않는 CT 장비와 정합함으로써, 측정대상환자의 무릎인대에 대한 물성을 정확히 도출할 수 있는 무릎인대 물성 도출방법을 제공을 제공할 수 있다.

Description

무릎인대 물성 도출방법{Material Properties Evaluaion Method for Ligaments in the Tibiofemoral}
본 발명은 무릎인대 물성 도출방법에 관한 것으로서, 보다 상세하게는 사체 실험을 통하지 않고 CT 장비와 TELOS 장비를 연동하여 특정하는 방법과 X-RAY-TELOS 측정을 하여 하중을 받지 않는 CT 장비와 정합하는 무릎인대 물성 도출방법에 관한 것이다.
도 1에는 종래 기술에 따른 웨이트 베어링 MRI(Weight-bearing MRI) 장비를 이용하여 측정대상환자를 촬영하고 있는 모습을 나타내는 사진이 도시되어 있다. 또한, 도 2에는 측정대상환자의 무릎에 대한 MRI촬영 사진으로서, 좌측사진은 일반 MRI 촬영 사진이고, 우측사진은 웨이트 베어링 MRI 촬영 사진이 도시되어 있다.
한편, 무릎환자에게 사용되는 인공 무릎인대는 기성복과 같이 평균화된 물성으로 제작된다.
경우에 따라서, 무릎환자의 정상인 무릎 부분에 대한 무릎인대 물성을 측정하여, 이와 동일한 물성을 가지는 인공무릎인대를 제작하여 무릎환자에게 시술하기도 한다. 이때, 무릎환자의 정상인 무릎에 대한 무릎인대 물성예측은 여러 가지 방법이 활용된다.
종래 기술에 따른 신체의 무릎인대 물성예측은, 사체(cadaver)를 이용하거나 생체외(in vitro) 방법이 사용된다.
생체외(in vitro) 방법의 경우, 예측된 물성은 생체내(in vivo) 물성값과 차이를 가진다.
생체내(in vivo) 방법으로 무릎인대에 대한 물성을 정확히 예측하려면, 측정대상 환자의 체중(weight)을 무릎에 부여한 상태에서 MRI촬영을 실시해야 한다. 이때 사용되는 장비가, 도 1에 도시된 웨이트 베어링 MRI(Weight bearing MRI) 장비이다.
도 1에 도시된 웨이트 베어링 MRI(Weight bearing MRI) 장비를 이용하여 도 2에 도시된 MRI촬영 사진을 획득할 수 있다.
도 2에 도시된 MRI촬영 사진을 통해 무릎인대에 대한 물성을 간접적으로 측정할 수 있다.
그러나, 도 1에 도시된 웨이트 베어링 MRI(Weight bearing MRI) 장비는 국내에 널비 보급되지 않은 고가의 장비로서, 전국에 분포된 여러 측정대상 환자에 대한 효율적인 촬영작업을 수행할 수 없다.
더욱이, 도 3에 도시된 바와 같이, 웨이트 베어링 MRI를 수행하기 위해서는 무릎에 인위적인 하중을 가하는 장비를 착용해야 하며, 장비 착용에 있어 많은 시간과 인력이 소요된다. 또한 도 3에 도시된 장비를 착용한 후 수행되는 웨이트 베어링 MRI 촬영 역시 일반 MRI 촬영 대비 많은 시간이 필요하며 숙련된 장비 운용 능력이 필요하다.
한편, 무릎인대를 포함한 무릎연조직은 24개의 연조직(meniscus 14개, cartilage 2개, ligament 8개)으로 구성되는 바, 종래 기술에 따른 무릎인대 및 연조직 물성 측정 방법에서는, 무릎인대를 포함한 무릎연조직을 구성하고 있는 모든 연조직에 대해 최적화 작업이 필요하여, 측정작업에 긴 시간이 소요된다는 문제점을 가지고 있다.
결과적으로, 종래 기술로는 무릎환자에게 최적화된 물성을 가지는 인공연골조직을 용이하게 제작할 수 없다.
따라서, 종래 무릎인대 물성 측정방법에 따른 문제점을 해결할 수 있는 기술이 필요한 실정이다.
한국공개특허 제10-2007-0058372호 (2007년 06월 08일 공개)
본 발명의 목적은, 사체 실험을 통하지 않고 CT 장비와 TELOS 장비를 연동하여 특정하는 방법과 X-RAY-TELOS 측정을 하여 하중을 받지 않는 CT 장비와 정합함으로써, 측정대상환자의 무릎인대에 대한 물성을 정확히 도출할 수 있는 무릎인대 물성 도출방법을 제공하는 것이다.
이러한 목적을 달성하기 위한 본 발명의 일 측면에 따른 무릎인대 물성 도출방법은, 측정대상환자의 무릎인대에 대한 물성을 도출하는 방법으로서, (a) 측정대상환자를 옆으로 눕힌 상태에서, 무릎을 특정 각도로 구부린 후, 무릎에 소정 크기의 하중을 특정 방향으로 가하는 하중인가단계; (b) 특정 각도로 구부러진 무릎에 대해 CT 촬영을 실시하여 촬영데이터를 획득하는 촬영단계; (c) 촬영단계로부터 획득한 데이터를 바탕으로 무릎에 대한 유한요소해석을 위한 3차원 모델을 생성시키고, 무릎인대의 물성값은 기공지된 평균 물성값으로 설정하는 무릎모델 생성단계; (d) 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하는 민감부위 선별단계; 및 (e) 하중인가단계에서 인가된 가압력 값과 촬영단계로부터 획득한 촬영데이터를 바탕으로 민감부위 선별단계를 통해 선별된 인대의 물성값을 최적화하는 물성 최적화단계;를 포함하는 구성일 수 있다.
본 발명의 일 실시예에 있어서, 상기 하중인가단계에서, 무릎에 가해지는 하중은, 무릎의 전방, 후방, 좌측방, 우측방 방향 가해질 수 있다.
이때, 상기 하중인가단계에서, 무릎에 하중을 가하는 도구는, 두 개의 지지부재와 하나의 가압부재 및 가압부재에 장착된 가압력 측정부재를 포함하는 구성일 수 있다.
경우에 따라서, 상기 무릎에 하중을 가하는 도구는, 텔로스 스트레스 디바이스(Telos Stress Device)일 수 있다.
본 발명의 일 실시예에 있어서, 상기 하중인가단계에서, 무릎에 가해지는 하중의 크기는 100 내지 200 N일 수 있다.
또한, 상기 하중인가단계에서, 무릎을 구부리는 각도는 30도, 60도 및 90도로 구성되는 군에서 하나 이상 선택되는 것일 수 있다.
본 발명의 일 실시예에 있어서, 상기 무릎인대 물성 도출방법은, 촬영단계로부터 획득한 촬영데이터를 바탕으로, 가압력에 따른 경골과 대퇴골의 이격 거리, 및 경골과 대퇴골 각각의 위치변경 거리를 측정하는 단계를 더 포함하는 구성일 수 있다.
본 발명의 일 실시예에 있어서, 상기 무릎모델 생성단계는, 인터 옵저버 리프로듀서빌러티(Inter-observer reproducibility) 방법을 이용하여 생성된 3차원 모델을 검증하는 3차원모델 복원검증단계를 포함하는 구성일 수 있다.
본 발명의 일 실시예에 있어서, 상기 민감부위 선별단계는, 몬테카를로(Monte Carlo) 기법을 이용하여 민감도해석을 수행할 수 있다.
본 발명의 일 실시예에 있어서, 상기 물성 최적화단계는, 3차원 모델의 선별된 인대의 유한요소해석 결과값이 CT 촬영 데이터와 동일한 거동을 보이도록 인대에 대한 물성값을 도출할 수 있다.
이상에서 설명한 바와 같이, 본 발명의 무릎인대 물성 도출방법에 따르면, 하중인가단계, 촬영단계, 무릎모델 생성단계, 민감부위 선별단계 및 물성 최적화단계를 포함하고 있어, 사체 실험을 통하지 않고 CT 장비와 TELOS 장비를 연동하여 특정하는 방법과 X-RAY-TELOS 측정을 하여 하중을 받지 않는 CT 장비와 정합함으로써, 측정대상환자의 무릎인대에 대한 물성을 정확히 도출할 수 있는 무릎인대 물성 도출방법을 제공할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 소정각도만큼 무릎을 구부린 상태(flexion)에서 특정방향으로 특정크기의 하중을 인가한 후, CT 촬영을 실하여 특정데이터를 획득하여, 이를 바탕으로 무릎인대의 물성을 손쉽게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, CT 촬영 데이터를 이용하여 생성된 3차원 유한요소 모델을 이용하여 무릎인대의 물성을 정확하게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하여 이에 대한 물성값을 최적화함으로써, 종래 기술 대비 짧은 시간 내에 무릎인대에 대한 물성을 신속하고 정확하게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 측정대상환자의 무릎인대의 물성을 손쉽게 측정할 수 있고, 결과적으로 측정대상환자의 무릎인대 물성과 정확히 일치하는 물성을 가지는 인공인대를 손쉽게 제작할 수 있다.
도 1은 종래 기술에 따른 웨이트 베어링 MRI(Weight-bearing MRI) 장비를 이용하여 측정대상환자를 촬영하고 있는 모습을 나타내는 사진이다.
도 2는 측정대상환자의 무릎에 대한 MRI촬영 사진으로서, 좌측사진은 일반 MRI 촬영 사진이고, 우측사진은 웨이트 베어링 MRI 촬영 사진이다.
도 3은 측정대상환자의 무릎에 대한 웨이트 베어링 MRI(Weight-bearing MRI) 촬영을 위해 무릎에 인위적인 하중을 가하는 장비를 착용한 상태를 나타내는 모식도이다.
도 4는 본 발명의 일 실시예에 따른 무릎인대 물성 도출방법을 나타내는 흐름도이다.
도 5는 도 4에 도시된 촬영단계를 통해 획득한 CT촬영 데이터를 나타내는 그림이다.
도 6은 도 4에 도시된 무릎모델 생성단계를 통해 생성된 3차원 모델을 나타내는 그림이다.
도 7은 도 4에 도시된 민감부위 선별단계를 통해 무릎인대를 구성하는 각 부분에 대한 민감도를 측정한 후 결과값을 나타낸 표이다.
도 8은 전방경골변형 상태(ATT)에서 무릎은 30도 구부린 상태에서 측정된 민감도 결과값을 나타낸 표이다.
도 9는 전방경골변형 상태(ATT)에서 무릎은 90도 구부린 상태에서 측정된 민감도 결과값을 나타낸 표이다.
도 10은 후방경골변형 상태(PTT)에서 무릎은 30도 구부린 상태에서 측정된 민감도 결과값을 나타낸 표이다.
도 11은 후방경골변형 상태(PTT)에서 무릎은 90도 구부린 상태에서 측정된 민감도 결과값을 나타낸 표이다.
이하 도면을 참조하여 본 발명의 바람직한 실시예를 상세히 설명하기로 한다. 이에 앞서, 본 명세서 및 청구범위에 사용된 용어나 단어는 통상적이거나 사전적인 의미로 한정하여 해석되어서는 아니되며, 본 발명의 기술적 사상에 부합하는 의미와 개념으로 해석되어야 한다.
본 명세서 전체에서, 어떤 부재가 다른 부재 "상에" 위치하고 있다고 할 때, 이는 어떤 부재가 다른 부재에 접해 있는 경우뿐 아니라 두 부재 사이에 또 다른 부재가 존재하는 경우도 포함한다. 본 명세서 전체에서, 어떤 부분이 어떤 구성요소를 "포함" 한다고 할 때, 이는 특별히 반대되는 기재가 없는 한 다른 구성요소를 제외하는 것이 아니라 다른 구성 요소를 더 포함할 수 있는 것을 의미한다.
도 4에는 본 발명의 일 실시예에 따른 무릎인대 물성 도출방법을 나타내는 흐름도가 도시되어 있고, 도 5에는 도 4에 도시된 촬영단계를 통해 획득한 CT촬영 데이터를 나타내는 그림이 도시되어 있으며, 도 6에는 도 4에 도시된 무릎모델 생성단계를 통해 생성된 3차원 모델을 나타내는 그림이 도시되어 있다.
도 5에 도시된 3D 모델은 랙서티 테스트(laxity test)를 위해 생성한 모델이다. 종래 기술에 따르면 반월판(meniscus)의 물성을 예측하기 위해 반월판의 움직임을 측정하고 그 수치에 맞도록 물성 최적화 작업을 수행하였다. 그러나, 본원발명에 따르면 인대(ligament)의 물성을 예측하기 위해 랙서티 테스트(laxity test)를 통해 ATT, PTT 수치에 맞도록 최적화 작업을 수행한다.
이들 도면을 도 7 내지 도 11과 함께 참조하면, 본 실시예에 따른 무릎인대 물성 도출방법(S100)은, 하중인가단계(S110), 촬영단계(S120), 무릎모델 생성단계(S130), 민감부위 선별단계(S140) 및 물성 최적화단계(S150)를 포함하는 구성일 수 있다.
구체적으로, 하중인가단계(S110)는, 측정대상환자를 옆으로 눕힌 상태에서, 무릎을 특정 각도로 구부린 후, 무릎에 소정 크기의 하중을 특정 방향으로 가하는 단계이다.
이때, 무릎에 가해지는 하중의 방향은 무릎의 전방, 후방, 좌측방, 우측방일 수 있다. 무릎에 가해지는 하중은 무릎인대의 변형 거리를 측정하기 위한 것으로서, 정확한 방향과 하중량을 제어할 수 있도록 특정 도구를 활용함이 바람직하다. 또한, 무릎에 가해지는 하중의 크기는, 무릎 인대의 변형 거리를 손쉽게 측정할 수 있는 범위 내의 하중 값이라면 특별히 제한되는 것은 아니며, 예를 들어 100 내지 200 N의 범위 내의 값일 수 있다.
구체적으로, 하중인가단계(S110)에서, 무릎에 하중을 가하는 도구(14)는, 두 개의 지지부재와 하나의 가압부재 및 가압부재에 장착된 가압력 측정부재를 포함하는 구조일 수 있다. 이에 대해 대표적인 예로서, 텔로스 스트레스 디바이스(Telos Stress Device)를 들 수 있으며, 이에 한정되지 않음은 물론이다.
종래 기술에 따르면, 수직하중이 가해졌을 때, 반월판(meniscus)가 눌리면서 발생하는 움직임을 통하여 반월판의 물성을 예측하였다. 그러나, 본원발명에 따르면, 인대 손상 여부를 판단하는 랙서티 테스트(laxity test)를 통하여 굴곡(flexion) 시 획득할 수 있는 ATT 값, PTT 값을 통하여 인대(ligament) 물성을 예측할 수 있다.
또한, 하중인가단계(S110)에서 무릎을 구부리는 각도는, 30도, 60도 및 90도로 구성되는 군에서 하나 이상 선택되는 것일 수 있다. 무릎을 구부리는 각도에 따라 무릎 인대의 변형 거리를 측정함으로써 무릎 인대의 물성을 더욱 정확히 도출해 낼 수 있다.
한편, 무릎모델 생성단계(S130)는, 촬영단계로부터 획득한 데이터를 바탕으로 무릎에 대한 유한요소해석을 위한 3차원 모델을 생성시키고, 무릎인대의 물성값은 기공지된 평균 물성값으로 설정하는 단계이다.
이때, 사용되는 프로그램은 촬영 데이터를 이용하여 3차원 모델을 생성시킬 수 있는 프로그램이라면 특별히 제한되는 것은 아니다. 예를 들어, Materialise사의 Mimics 프로그램을 사용할 수 있으나, 이에 한정되는 것은 아니다.
종래 기술에 따르면 랙서티 테스트(laxity test)는 텔로스장비(Telos 장비)를 이용하여 X-Ray 촬영을 통하여 수행된다. 그러나, 본원발명에 따르면, 굴곡(flexion) 시 발생하는 무릎관절의 내외측 회전을 고려하여, 보다 정확하게 ATT값과 PTT값을 측정하기 위하여 3차원 모델을 구축한다. 따라서, 이를 위해 CT 장비를 이용하여 랙서티 테스트(laxity test)를 수행한다. 경우에 따라서, 추가로 X-Ray 촬영 후 CT 촬영 데이터와 정합하여 3차원 모델을 구축할 수 있다.
또한, 무릎모델 생성단계(S130)에서, 3차원모델 복원검증단계가 수행될 수 있다. 구체적으로, 3차원모델 복원검증단계는 인터 옵저버 리프로듀서빌러티(Inter-observer reproducibility) 방법을 이용하여 생성된 3차원 모델을 검증하는 단계이다.
인터 옵저버 리프로듀서빌러티(Inter-observer reproducibility) 방법이란, 관찰자간 재현성에 관한 방법으로서, 둘 이상의 독립적인 관찰자들이 동일한 대상을 관찰함에 있어 그 측정결과가 얼마나 일치하는가를 나타내는 정도를 이용하는 방법이다. 관찰자간 일치도를 측정하는 전통적인 방법은 관찰자 점수들 간의 상관관계나 일치도 계수(예, Kappa coefficient)를 계산하는 방법이다. 점수들 간의 상관관계가 높거나 일치도 계수가 높게 추정되면 관찰자간 일치도가 높은 것을 가정한다. 경우에 따라서, 관찰자에 의해 측정된 결과값의 표준편차(CV, coefficient of Variation)를 구하여 관찰자간 일치도를 판단할 수 있다. 이때, 표준편차(CV) 값이 작을수록 관찰자간 일치도가 높다고 판단한다. 본 발명에서는 표준편차(CV) 값을 이용하여 관찰자간 일치도를 판단하였다.
만약 상관관계가 낮거나 일치도 계수가 낮게 추정되면 일치도가 낮은 것으로 해석하는데, 이 때는 여러 가지 원인이 있을 수 있다. 만약 관찰시 사용한 측정도구의 이 모호한 경우라면 그 기준을 보다 분명히 수정하여야 할 것이다. 측정도구상의 문제가 아니라면 관찰자 훈련을 강화하여 측정도구에 대한 이해를 높임으로써 측정과정상의 오차를 줄일 필요가 있다. 관찰 상황 하에서의 일치도 및 신뢰도 문제분석은 단순지수의 계산보다는 (generalizability theory)의 활용을 통한 보다 면밀한 분석을 요한다
인터 옵저버 리프로듀서빌러티 방법을 이용하여 생성된 3차원 모델을 검증하여, 상관관계가 높을 경우 생성된 3차원 모델은 올바르게 생성된 것으로 판단하고 다음 단계를 수행한다.
또한, 본 실시예에 따른 무릎인대 물성 도출방법(S100)은, 촬영단계(S120)로부터 획득한 촬영데이터를 바탕으로, 가압력에 따른 경골과 대퇴골의 이격 거리, 및 경골과 대퇴골 각각의 위치변경 거리를 측정하는 단계를 더 수행할 수 있다. 이때 측정된 데이터를 바탕으로 무릎인대에 대한 민감도해석을 수행할 수 있다.
다음 단계로서, 민감부위 선별단계(S140)가 수행될 수 있다. 구체적으로, 민감부위 선별단계(S140)는, 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하는 단계이다.
이때, 민감부위 선별단계(S140)는 몬테카를로(Monte Carlo) 기법을 이용하여 민감도해석을 수행할 수 있다. 민감도해석을 수행하여 민감도가 다른 인자보다 현저히 높은 인자를 선별할 수 있다.
이때 상기 언급한 몬테카를로(Monte Carlo) 기법이란, 확률론적 시뮬레이션 기법에 관한 것으로서, 불확실한 상황하에서의 의사결정을 목적으로 확률적시스템의 모의 실험에 이용되는 절차를 말한다. 몬테카를로 기법의 핵심은 모형의 확률요소들에 대한 실험인데, 이는 확률적 또는 우연결과를 발생시켜 주는 도구를 이용하여 수행된다. 이 도구는 모형에서 가정한 확률분포에 따라 무작위표본추출에 의해서 우연결과를 발생시켜 주는데 이용된다. 따라서 몬테카를로 기법을 모의적 표본 추출법(simulated sampling technique)이라고도 한다. 우연결과 또는 확률적 결과를 발생시켜 주는데 이용되는 도구로는 주사위나 룰렛바퀴(roulett wheel), 복권추첨에 이용되는 숫자공 등이 있으며, 가장 일반적으로 쓰이는 것은 난수(random number)이며 컴퓨터에 의해 발생되는 의사난수(pseudo-random number) 등이 있다.
몬테카를로 기법의 시행을 위한 절차는, (1) 확률변수의 확률분포를 얻는 단계, (2) 누적확률분포를 설정하는 단계, (3) 확률변수의 값이나 값의 범위를 나타내기 위해서 적절한 난수의 집합/난수의 구간을 할당하는 단계, (4) 무작위표본추출을 이용하여 시뮬레이션 실험을 실시하는 단계, (5) 행동방안을 설계/시행하고 통제하는 단계로 구성될 수 있다.
몬테카를로 기법을 이용하여, 인체의 무릎 연조직을 구성하고 있는 수많은 연조직 중 민감도가 현저히 높은 특정 인자를 선별할 수 있다.
인체의 무릎 연조직을 구성하고 있는 각 연조직에 대한 물성값은, 인종별, 대상환자별, 성별별로 큰 차이를 가지고 있으며, 결과적으로 상당히 넓은 표준편차를 가지고 있다. 이 때문에 모든 인자에 대해 물성값 최적화 단계를 수행할 경우, 많은 시간과 인력이 소요되는 문제점이 발생한다.
따라서, 함수로 분포되어 있는 무수한 인자들의 민감도를 분석하기 위해 몬테카를로 기법을 사용하고, 무수히 많은 인자들 중 민감도가 현저히 높은 특정 인자를 손쉽게 선별할 수 있다.
한편, 민감부위 선별단계(S140)는, 무릎의 완전이완(Full extension) 상태에서 무릎연조직에 대한 민감도해석이 수행됨이 바람직하다.
도 7 내지 도 11에 도시된 표와 같이, 민감부위 선별단계(S140)를 통해, 무릎인대를 구성하는 각 부분에 대한 민감도를 측정한 후 결과값을 바탕으로 소정 수치값 이상인 민감부위를 선별할 수 있다.
다음 단계로서, 물성 최적화단계(S150)가 수행될 수 있다. 구체적으로, 물성 최적화단계(S150)는, 3차원 모델의 무릎을 소정각도만큼 구부리며(flexion) 민감부위 선별단계를 통해 선별된 무릎인대의 물성을 최적화하는 단계이다.
구체적으로, 물성 최적화단계(S150)는 3차원 모델의 선별된 무릎인대의 유한요소해석 결과값이 CT 촬영데이터와 동일한 거동을 보이도록 무릎인대에 대한 물성값을 도출하는 단계이다.
무릎모델 생성단계(S130)로부터 생성된 3차원 모델의 무릎인대에는 기공지된 평균 물성값이 설정되어 있다. 이러한 3차원 모델에 대한 유한요소 해석과정을 실시한 후, 해설 결과값이 CT 촬영데이터와 일치함을 판단한다. 만약 일치하지 않을 경우 선별된 무릎인대에 대한 물성값을 재설정한 후 재차 유한요소 해석과정이 실시된다. 재차실시되는 유한요소 해석과정과 선별된 무릎인대에 대한 물성값 재설정 과정에 의해, 무릎인대에 대한 물성값을 최적화 될 수 있다.
본 실시예에 따른 물성 최적화 단계는, 앞서 수행된 민감부위 선별단계(S140)를 통해 선별된 무릎인대에 대한 최적화 작업을 수행한다. 따라서, 무릎인대를 구성하고 있는 모든 조직에 대한 최적화 작업을 수행하지 않아도 되므로, 모든 조직에 대해 최적화 작업을 수행하는 종래 기술에 따른 방법 대비 빠른 시간 내에 무릎인대에 대한 물성을 측정할 수 있다.
이상에서 설명한 바와 같이, 본 발명의 무릎인대 물성 도출방법에 따르면, 하중인가단계, 촬영단계, 무릎모델 생성단계, 민감부위 선별단계 및 물성 최적화단계를 포함하고 있어, 사체 실험을 통하지 않고 CT 장비와 TELOS 장비를 연동하여 특정하는 방법과 X-RAY-TELOS 측정을 하여 하중을 받지 않는 CT 장비와 정합함으로써, 측정대상환자의 무릎인대에 대한 물성을 정확히 도출할 수 있는 무릎인대 물성 도출방법을 제공할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 소정각도만큼 무릎을 구부린 상태(flexion)에서 특정방향으로 특정크기의 하중을 인가한 후, CT 촬영을 실하여 특정데이터를 획득하여, 이를 바탕으로 무릎인대의 물성을 손쉽게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, CT 촬영 데이터를 이용하여 생성된 3차원 유한요소 모델을 이용하여 무릎인대의 물성을 정확하게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하여 이에 대한 물성값을 최적화함으로써, 종래 기술 대비 짧은 시간 내에 무릎인대에 대한 물성을 신속하고 정확하게 도출할 수 있다.
또한, 본 발명의 무릎인대 물성 도출방법에 따르면, 측정대상환자의 무릎인대의 물성을 손쉽게 측정할 수 있고, 결과적으로 측정대상환자의 무릎인대 물성과 정확히 일치하는 물성을 가지는 인공인대를 손쉽게 제작할 수 있다.
이상의 본 발명의 상세한 설명에서는 그에 따른 특별한 실시예에 대해서만 기술하였다. 하지만 본 발명은 상세한 설명에서 언급되는 특별한 형태로 한정되는 것이 아닌 것으로 이해되어야 하며, 오히려 첨부된 청구범위에 의해 정의되는 본 발명의 정신과 범위 내에 있는 모든 변형물과 균등물 및 대체물을 포함하는 것으로 이해되어야 한다.
즉, 본 발명은 상술한 특정의 실시예 및 설명에 한정되지 아니하며, 청구범위에서 청구하는 본 발명의 요지를 벗어남이 없이 본 발명이 속하는 기술 분야에서 통상의 지식을 가진 자라면 누구든지 다양한 변형 실시가 가능하며, 그와 같은 변형은 본 발명의 보호 범위 내에 있게 된다.
10: 측정대상환자의 무릎
11: 연골(cartilage)
12: 인대(ligament)
13: 반월판(meniscus)
14: 무릎에 인위적인 하중을 가하는 장비
S100: 무릎인대 물성 도출방법
S110: 하중인가단계
S120: 촬영단계
S130: 무릎모델 생성단계
S140: 민감부위 선별단계
S150: 물성 최적화단계

Claims (10)

  1. 측정대상환자의 무릎인대에 대한 물성을 도출하는 방법으로서,
    (a) 측정대상환자를 옆으로 눕힌 상태에서, 무릎을 특정 각도로 구부린 후, 무릎에 소정 크기의 하중을 특정 방향으로 가하는 하중인가단계(S110);
    (b) 특정 각도로 구부러진 무릎에 대해 CT 촬영을 실시하여 촬영데이터를 획득하는 촬영단계(S120);
    (c) 촬영단계로부터 획득한 데이터를 바탕으로 무릎에 대한 유한요소해석을 위한 3차원 모델을 생성시키고, 무릎인대의 물성값은 기공지된 평균 물성값으로 설정하는 무릎모델 생성단계(S130);
    (d) 무릎인대에 대한 민감도해석을 통해 민감도가 소정 수치값 이상인 인대를 선별하는 민감부위 선별단계(S140); 및
    (e) 하중인가단계에서 인가된 가압력 값과 촬영단계로부터 획득한 촬영데이터를 바탕으로 민감부위 선별단계로 통해 선별된 인대의 물성값을 최적화하는 물성 최적화단계(S150);
    를 포함하는 것을 특징으로 하는 무릎인대 물성 도출방법(S100).
  2. 제 1 항에 있어서,
    상기 하중인가단계(S110)에서, 무릎에 가해지는 하중은, 무릎의 전방, 후방, 좌측방 및 우측방 방향으로 가해지는 것을 특징으로 하는 무릎인대 물성 도출방법.
  3. 제 1 항에 있어서,
    상기 하중인가단계(S110)에서, 무릎에 하중을 가하는 도구(14)는, 두 개의 지지부재와 하나의 가압부재 및 가압부재에 장착된 가압력 측정부재를 포함하는 것을 특징으로 하는 무릎인대 물성 도출방법.
  4. 제 3 항에 있어서,
    상기 무릎에 하중을 가하는 도구(14)는, 텔로스 스트레스 디바이스(Telos Stress Device)인 것을 특징으로 하는 무릎인대 물성 도출방법.
  5. 제 1 항에 있어서,
    상기 하중인가단계(S110)에서, 무릎에 가해지는 하중의 크기는 100 내지 200 N인 것을 특징으로 하는 무릎인대 물성 도출방법.
  6. 제 1 항에 있어서,
    상기 하중인가단계(S110)에서, 무릎을 구부리는 각도는 30도, 60도 및 90도로 구성되는 군에서 하나 이상 선택되는 것을 특징으로 하는 무릎인대 물성 도출방법.
  7. 제 1 항에 있어서,
    상기 촬영단계(S120)는, 획득한 촬영데이터를 바탕으로, 가압력에 따른 경골과 대퇴골의 이격 거리, 및 경골과 대퇴골 각각의 위치변경 거리를 측정하는 것을 특징으로 하는 무릎인대 물성 도출방법.
  8. 제 1 항에 있어서,
    상기 무릎모델 생성단계(S130)는, 인터 옵저버 리프로듀서빌러티(Inter-observer reproducibility) 방법을 이용하여 생성된 3차원 모델을 검증하는 3차원모델 복원검증단계를 포함하는 것을 특징으로 하는 무릎인대 물성 도출방법.
  9. 제 1 항에 있어서,
    상기 민감부위 선별단계(S140)는, 몬테카를로(Monte Carlo) 기법을 이용하여 민감도해석을 수행하는 것을 특징으로 하는 무릎인대 물성 도출방법.
  10. 제 1 항에 있어서,
    상기 물성 최적화단계(S150)는, 3차원 모델의 선별된 인대의 유한요소해석 결과값이 CT 촬영 데이터와 동일한 거동을 보이도록 인대에 대한 물성값을 도출하는 것을 특징으로 하는 무릎인대 물성 도출방법.
KR1020170047941A 2017-04-13 2017-04-13 무릎인대 물성 도출방법 Active KR101886650B1 (ko)

Priority Applications (1)

Application Number Priority Date Filing Date Title
KR1020170047941A KR101886650B1 (ko) 2017-04-13 2017-04-13 무릎인대 물성 도출방법

Applications Claiming Priority (1)

Application Number Priority Date Filing Date Title
KR1020170047941A KR101886650B1 (ko) 2017-04-13 2017-04-13 무릎인대 물성 도출방법

Publications (1)

Publication Number Publication Date
KR101886650B1 true KR101886650B1 (ko) 2018-08-10

Family

ID=63229548

Family Applications (1)

Application Number Title Priority Date Filing Date
KR1020170047941A Active KR101886650B1 (ko) 2017-04-13 2017-04-13 무릎인대 물성 도출방법

Country Status (1)

Country Link
KR (1) KR101886650B1 (ko)

Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6419645B1 (en) * 2000-11-02 2002-07-16 Arie M. Rijke Device and method for evaluating injuries to ligaments
KR20070058372A (ko) 2004-02-04 2007-06-08 엘디알 메디칼 추간 연골 인공 기관
JP2016096889A (ja) * 2014-11-19 2016-05-30 株式会社東芝 画像解析装置、画像解析方法、およびプログラム
WO2017020069A1 (en) * 2015-07-31 2017-02-09 360 Knee Systems Pty Ltd Post-operative prediction
KR101707646B1 (ko) * 2016-06-08 2017-02-20 주식회사 티제이씨라이프 무릎연조직 물성 도출방법

Patent Citations (5)

* Cited by examiner, † Cited by third party
Publication number Priority date Publication date Assignee Title
US6419645B1 (en) * 2000-11-02 2002-07-16 Arie M. Rijke Device and method for evaluating injuries to ligaments
KR20070058372A (ko) 2004-02-04 2007-06-08 엘디알 메디칼 추간 연골 인공 기관
JP2016096889A (ja) * 2014-11-19 2016-05-30 株式会社東芝 画像解析装置、画像解析方法、およびプログラム
WO2017020069A1 (en) * 2015-07-31 2017-02-09 360 Knee Systems Pty Ltd Post-operative prediction
KR101707646B1 (ko) * 2016-06-08 2017-02-20 주식회사 티제이씨라이프 무릎연조직 물성 도출방법

Similar Documents

Publication Publication Date Title
Kang et al. Probabilistic evaluation of the material properties of the in vivo subject‐specific articular surface using a computational model
CN107847275B (zh) 用于术前规划的用以定位软组织的系统和方法
EP2805277B1 (en) Method and system for human joint treatment plan and personalized surgery planning using 3-d kinematics, fusion imaging and simulation
Szczerbik et al. The influence of knee marker placement error on evaluation of gait kinematic parameters
Lopes et al. Reliability and validity of frontal plane kinematics of the trunk and lower extremity measured with 2-dimensional cameras during athletic tasks: a systematic review with meta-analysis
KR102323492B1 (ko) 영상 분석을 통한 근골격계 균형 진단 및 이를 이용한 운동처방 방법
Rachmat et al. Generating finite element models of the knee: How accurately can we determine ligament attachment sites from MRI scans?
US20160331465A1 (en) System and method for simulating reconstructive surgery of anterior cruciate ligament using medical images
Belyea et al. Validity and intrarater reliability of a 2-dimensional motion analysis using a handheld tablet compared with traditional 3-dimensional motion analysis
CN114431957A (zh) 基于深度学习的全膝关节置换术后翻修术前规划方法
Relph et al. Criterion-related validity of knee joint-position-sense measurement using image capture and isokinetic dynamometry
Kang et al. In vivo evaluation of the subject-specific finite element model for knee joint cartilage contact area
Alahmari et al. Concurrent validity of two-dimensional video analysis of lower-extremity frontal plane of movement during multidirectional single-leg landing
King et al. Reliability of using a handheld tablet and application to measure lower-extremity alignment angles
KR101707646B1 (ko) 무릎연조직 물성 도출방법
Coyne et al. Marker location and knee joint constraint affect the reporting of overhead squat kinematics in elite youth football players
Price et al. A model‐based motion capture marker location refinement approach using inverse kinematics from dynamic trials
KR101118372B1 (ko) 가상 마커 방식을 이용한 비만/체형을 평가하는 시스템 및 방법
KR101886650B1 (ko) 무릎인대 물성 도출방법
Elyasi et al. Biomechanical lower limb model to predict patellar position alteration after medial open wedge high tibial osteotomy
CN116158774B (zh) 用于测量胫骨平台后倾角的方法、设备和计算机存储介质
Lukas et al. Automated measurement of joint space width in small joints of patients with rheumatoid arthritis
KR102021942B1 (ko) 골강도 및 예상골절조건 산출 방법 및 이를 적용한 예상골절조건 산출 시스템
Kahmann et al. A combined experimental and finite element analysis of the human elbow under loads of daily living
Lafayette et al. Hybrid solution for motion capture with kinect v2 to different biotypes recognition

Legal Events

Date Code Title Description
PA0109 Patent application

Patent event code: PA01091R01D

Comment text: Patent Application

Patent event date: 20170413

PA0201 Request for examination
E902 Notification of reason for refusal
PE0902 Notice of grounds for rejection

Comment text: Notification of reason for refusal

Patent event date: 20180627

Patent event code: PE09021S01D

E701 Decision to grant or registration of patent right
PE0701 Decision of registration

Patent event code: PE07011S01D

Comment text: Decision to Grant Registration

Patent event date: 20180716

GRNT Written decision to grant
PR0701 Registration of establishment

Comment text: Registration of Establishment

Patent event date: 20180802

Patent event code: PR07011E01D

PR1002 Payment of registration fee

Payment date: 20180803

End annual number: 3

Start annual number: 1

PG1601 Publication of registration
PC1903 Unpaid annual fee
PR0401 Registration of restoration

Patent event code: PR04011E01D

Patent event date: 20220310

Comment text: Registration of Restoration

PR1001 Payment of annual fee

Payment date: 20220311

Start annual number: 4

End annual number: 4

PR1001 Payment of annual fee

Payment date: 20220729

Start annual number: 5

End annual number: 5