KR101179081B1 - Apparatus for designing custom-made prosthesis based on bone density measurement - Google Patents
Apparatus for designing custom-made prosthesis based on bone density measurement Download PDFInfo
- Publication number
- KR101179081B1 KR101179081B1 KR1020110012517A KR20110012517A KR101179081B1 KR 101179081 B1 KR101179081 B1 KR 101179081B1 KR 1020110012517 A KR1020110012517 A KR 1020110012517A KR 20110012517 A KR20110012517 A KR 20110012517A KR 101179081 B1 KR101179081 B1 KR 101179081B1
- Authority
- KR
- South Korea
- Prior art keywords
- bone
- image
- bone density
- prosthesis
- information
- Prior art date
Links
- 230000037182 bone density Effects 0.000 title claims abstract description 66
- 238000001739 density measurement Methods 0.000 title description 2
- 210000000988 bone and bone Anatomy 0.000 claims abstract description 70
- 238000000034 method Methods 0.000 claims abstract description 27
- 238000013461 design Methods 0.000 claims description 23
- 239000007943 implant Substances 0.000 claims description 13
- 238000004040 coloring Methods 0.000 claims description 2
- 238000002591 computed tomography Methods 0.000 description 30
- 238000012545 processing Methods 0.000 description 11
- 238000002595 magnetic resonance imaging Methods 0.000 description 6
- 208000029725 Metabolic bone disease Diseases 0.000 description 5
- 230000000877 morphologic effect Effects 0.000 description 4
- 210000004197 pelvis Anatomy 0.000 description 4
- 238000004088 simulation Methods 0.000 description 4
- 238000003325 tomography Methods 0.000 description 4
- 238000012800 visualization Methods 0.000 description 4
- 208000001132 Osteoporosis Diseases 0.000 description 3
- 241001270131 Agaricus moelleri Species 0.000 description 2
- 206010049088 Osteopenia Diseases 0.000 description 2
- 230000037118 bone strength Effects 0.000 description 2
- 238000001514 detection method Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 230000009977 dual effect Effects 0.000 description 2
- 210000001624 hip Anatomy 0.000 description 2
- 239000011159 matrix material Substances 0.000 description 2
- 238000005259 measurement Methods 0.000 description 2
- 238000011541 total hip replacement Methods 0.000 description 2
- 208000003076 Osteolysis Diseases 0.000 description 1
- 238000010521 absorption reaction Methods 0.000 description 1
- 210000003484 anatomy Anatomy 0.000 description 1
- 238000011882 arthroplasty Methods 0.000 description 1
- 230000005540 biological transmission Effects 0.000 description 1
- 230000008468 bone growth Effects 0.000 description 1
- 238000006243 chemical reaction Methods 0.000 description 1
- 238000003759 clinical diagnosis Methods 0.000 description 1
- 238000013170 computed tomography imaging Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000002059 diagnostic imaging Methods 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 238000010494 dissociation reaction Methods 0.000 description 1
- 230000005593 dissociations Effects 0.000 description 1
- 239000000284 extract Substances 0.000 description 1
- 238000000605 extraction Methods 0.000 description 1
- 230000004927 fusion Effects 0.000 description 1
- 238000011540 hip replacement Methods 0.000 description 1
- 238000003384 imaging method Methods 0.000 description 1
- 238000003780 insertion Methods 0.000 description 1
- 230000037431 insertion Effects 0.000 description 1
- 208000029791 lytic metastatic bone lesion Diseases 0.000 description 1
- 210000000276 neural tube Anatomy 0.000 description 1
- 230000005855 radiation Effects 0.000 description 1
- 238000011160 research Methods 0.000 description 1
- 238000001356 surgical procedure Methods 0.000 description 1
- 238000012546 transfer Methods 0.000 description 1
- 238000002604 ultrasonography Methods 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B5/00—Measuring for diagnostic purposes; Identification of persons
- A61B5/45—For evaluating or diagnosing the musculoskeletal system or teeth
- A61B5/4504—Bones
- A61B5/4509—Bone density determination
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T15/00—3D [Three Dimensional] image rendering
- G06T15/50—Lighting effects
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T19/00—Manipulating 3D models or images for computer graphics
- G06T19/20—Editing of 3D images, e.g. changing shapes or colours, aligning objects or positioning parts
Landscapes
- Engineering & Computer Science (AREA)
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Physics & Mathematics (AREA)
- Medical Informatics (AREA)
- Surgery (AREA)
- Public Health (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Veterinary Medicine (AREA)
- General Health & Medical Sciences (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- Computer Graphics (AREA)
- Molecular Biology (AREA)
- Heart & Thoracic Surgery (AREA)
- Pathology (AREA)
- Biomedical Technology (AREA)
- Software Systems (AREA)
- Optics & Photonics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Architecture (AREA)
- High Energy & Nuclear Physics (AREA)
- Computer Hardware Design (AREA)
- Radiology & Medical Imaging (AREA)
- General Engineering & Computer Science (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Orthopedic Medicine & Surgery (AREA)
- Rheumatology (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
본 발명에 따른 맞춤형 인공삽입물 설계 방법은 2D로 촬영된 뼈 이미지를 이용하여 3D(Three Dimensions) 영상을 생성하는 단계와; 상기 2D 뼈 이미지를 이용하여 골밀도 정보를 측정하는 단계와; 상기 골밀도 정보를 상기 3D 영상에 색채화하는 단계와; 상기 색채화된 3D 영상을 이용하여, 인공 삽입물을 설계하기 위한 인공삽입물의 크기, 위치, 골고정 나사 구멍의 위치 값을 산출하는 단계를 포함할 수 있다.According to the present invention, there is provided a method of designing a customized prosthesis, comprising: generating a 3D (Three Dimensions) image using a bone image photographed in 2D; Measuring bone density information using the 2D bone image; Colorizing the bone density information on the 3D image; Using the colored 3D image, the method may include calculating a size, a position, and a position value of the bone fixation screw hole for designing the prosthesis.
Description
본 발명은 인공삽입물의 설계에 관한 것이다.The present invention relates to the design of prostheses.
CT(Computed Tomography) 이미지는 x-ray를 여러 각도에서 인체에 투영하여 촬영하고 이를 컴퓨터로 재구성하여 인체 내부 단면의 이미지를 얻어내는 방법을 말한다. Computed Tomography (CT) image is a method of projecting x-rays to the human body from various angles and reconstructing them with a computer to obtain images of the internal cross section of the human body.
일반적으로 임상에서 진단용으로 사용되는 CT 이미지는 2~5㎜ 간격으로 촬영되며 연구 목적인 경우 1㎜ 이하의 간격으로 촬영되기도 한다. In general, CT images used for clinical diagnosis are taken at intervals of 2 to 5 mm, and for research purposes, may be taken at intervals of 1 mm or less.
인공삽입물의 형상정보를 얻기 위해서 필요한 뼈의 형태는 부위마다 다르지만(골반의 경우 약 20㎝ 정도), 3차원으로 모델링 하기 위해서는 일반적으로 대략 100장 이상의 이미지를 얻게 된다. 골형상은 복잡한 3차원적인 형상이기 때문에 척추에 수직적인 단면만을 제공하는 2차원 CT 이미지만으로는 환자의 골형상에 적합한 인공삽입물의 형상 설계 변수를 측정하기 어렵다. Although the shape of the bone required to obtain the shape information of the prosthesis varies from site to site (about 20 cm in the case of the pelvis), approximately 100 images are generally obtained to model in three dimensions. Since the bone shape is a complex three-dimensional shape, it is difficult to measure the shape design parameters of the prosthesis suitable for the bone shape of the patient only by the two-dimensional CT image that provides only the cross-section perpendicular to the spine.
이 때문에 2차원 CT 이미지를 이용하여 3차원으로 재구성하는 영상처리 기술이 활발히 연구되고 있다. For this reason, image processing techniques for reconstructing in three dimensions using two-dimensional CT images have been actively studied.
그러나 기존의 연구들은 CT 이미지에서 뼈의 윤곽선을 측정하고 이를 이용하여 3차원 재구성하는 방식이다. 이러한 방식은 골형상의 형태학적 정보만을 담고 있어 인공삽입물의 시술시 인공삽입물과 뼈와의 정합도는 높여줄 수 있지만 골밀도와 같은 생물학적인 정보를 표현해 주지 못하여 골고정 나사가 필요한 경우 충분한 고정력을 줄 수 없게 되어 재차 수술을 필요로 할 수도 있는 문제점이 있었다.Existing studies, however, have measured bone contours in CT images and used them to three-dimensional reconstruction. This method only contains the morphological information of the bone shape, which can increase the degree of coordination between the implant and the bone when the implant is being implanted, but it does not express biological information such as bone density, which gives sufficient fixation force when bone fixation screws are needed. There was a problem that could not be able to require surgery again.
또한, 고령의 환자나 골감소증, 골다공증 질환을 가지고 있는 환자의 경우 골밀도의 감소에 따라 뼈의 강도가 약해져서 인공삽입물의 고정력 확보에 어려움이 있다. 특히, 인공 고관절 재치환술의 경우 기존의 골반에 결합되어 있는 비구컵을 제거하고 새로운 비구컵으로 교체하게 되는데, 이때 비구컵과 골반 사이의 골 성장으로 붙어있던 골 조직이 비구컵의 교체와 함께 유실된다. 이로 인하여 골 고정부위의 제한과 약해진 뼈의 강도에 의해 고정력을 확보하기 힘들다. 따라서, 비구컵의 고정력을 높이기 위하여 골고정 나사를 이용하여 고정하게 되는데, 골고정 나사는 골반에 골밀도가 높은 부위에 고정되어야 한다. 하지만 환자마다 골밀도 정도가 다르므로, 맞춤형 비구컵 설계 시 고정 나사의 위치는 골밀도를 고려하여 설계되어야 하지만, 종래에는 어려움이 있었다.In addition, in elderly patients or patients with osteopenia, osteoporosis disease, the bone strength is weakened according to the decrease in bone density, which makes it difficult to secure the fixation force of the prosthesis. In particular, in the case of total hip arthroplasty, the acetabular cup coupled to the existing pelvis is removed and replaced with a new acetabular cup, wherein the bone tissue attached to the bone growth between the acetabular cup and the pelvis is lost with the replacement of the acetabular cup. Because of this, it is difficult to secure fixation force due to the limitation of the bone fixation and weakened bone strength. Therefore, in order to increase the fixing force of the acetabular cup is fixed using a bone fixation screw, the bone fixation screw should be fixed to the site of high bone density in the pelvis. However, since the degree of bone density varies from patient to patient, the location of the fixation screw should be designed in consideration of bone density when designing a customized acetabular cup, but there was a difficulty in the related art.
따라서 본 발명의 목적은 전술한 문제점을 해결하는 것을 목적으로 한다. Accordingly, an object of the present invention is to solve the above problems.
구체적으로, 본 발명은 골고정 나사에 충분한 고정력을 제공할 수 있도록 골밀도와 같은 생물학적인 정보를 이용하는 것을 목적으로 한다. Specifically, the present invention aims to use biological information such as bone density to provide sufficient fixation force to the bone fixation screw.
또한, 본 발명은 환자마다 다른 뼈 형상뿐만 아니라 골밀도 정도를 고려한 맞춤형 인공삽입물의 설계 변수 결정 및 골고정 나사의 위치 설정을 할 수 있도록 함에 있다.In addition, the present invention is to be able to determine the design parameters of the customized prosthesis considering the bone density as well as different bone shape for each patient and to position the bone fixation screw.
전술한 목적을 달성하기 위하여, 본 발명은 CT 이미지에서 측정할 수 있는 뼈의 윤곽선 정보를 이용하여, 3차원 영상을 재구성하고 골밀도 정보를 재구성한 영상의 색채화를 통하여, 맞춤형 인공삽입물의 설계 변수 및 골고정 나사의 위치를 결정한다. 다시 말해서, 본 발명은 CT 이미지를 이용하여 환자의 골형상 및 골밀도에 대한 3차원 영상을 구축하고, 환자의 3차원 골형상 및 골밀도 영상으로부터 인공삽입물의 설계 변수 및 골고정 나사 위치 정보를 결정한다.In order to achieve the above object, the present invention uses the contour information of the bone that can be measured in the CT image, reconstruct the three-dimensional image and colorize the image of the reconstructed bone density information, the design parameters of the custom implant And the position of the bone fixation screw. In other words, the present invention constructs a three-dimensional image of the bone shape and bone density of the patient using a CT image, and determines the design parameters and bone fixation screw position information of the implant from the three-dimensional bone shape and bone density image of the patient .
본 발명은 맞춤형 인공삽입물(Prosthesis)을 설계하는데 있어서 CT(Computed Tomography)로부터 획득한 뼈의 윤곽선 정보와 골밀도에 관한 정보를 이용하여 3차원 영상을 재구성(Reconstruction)하고, 구성된 영상에 골밀도를 입체적으로 가시화하여 시뮬레이션한다. The present invention reconstructs a three-dimensional image by using bone contour information and bone density information obtained from CT (Computed Tomography) in designing a customized prosthesis, and stereoscopically configures the bone density in the constructed image. Visualize and simulate.
또한, 본 발명에서는 맞춤형 인공삽입물의 설계 시, 뼈의 형태학적 정보뿐만 아니라 인공삽입물이 삽입될 위치의 골밀도 분포에 관한 정보를 포함하는 3차원 영상을 형성하여 인공삽입물의 위치, 크기 정보와 인공삽입물의 고정력을 높이기 위한 골 고정 나사의 위치정보를 포함하는 것을 특징으로 한다. In addition, in the present invention, when designing a custom implant, a three-dimensional image including not only the morphological information of the bone but also information on the bone density distribution of the position where the implant is to be inserted is formed to form the position, size information and the implant of the implant. Characterized in that it includes the position information of the bone fixing screw to increase the fixing force of.
또한, 본 발명에서는 윤곽선 추출 알고리즘을 이용하여 CT 이미지에서 뼈의 윤곽선 정보를 추출하고, 이에 대하여 3차원 보간을 통하여 3차원 영상을 형성하고, 다음으로 상기 3차원 영상에 골밀도 분포에 따른 색채맵을 설정한 후, 마지막으로 상기 영상에 골밀도 정보를 가시화시켜 인공 삽입물의 기본 설계변수 및 골밀도를 고려한 골 고정나사의 위치를 결정할 수 있도록 한다.Also, in the present invention, the contour information of the bone is extracted from the CT image by using the contour extraction algorithm, and a three-dimensional image is formed through the three-dimensional interpolation, and the color map according to the bone density distribution is then added to the three-dimensional image. Finally, the bone density information is visualized in the image to determine the position of the bone fixation screw considering the basic design parameters and the bone density of the artificial implant.
상기의 과제를 이루기 위한 본 발명에 따른 맞춤형 인공삽입물 설계 방법은 2D로 촬영된 뼈 이미지를 이용하여 3D(Three Dimensions) 영상을 생성하는 단계와; 상기 2D 뼈 이미지를 이용하여 골밀도 정보를 측정하는 단계와; 상기 골밀도 정보를 상기 3D 영상에 색채화하는 단계와; 상기 색채화된 3D 영상을 이용하여, 인공 삽입물을 설계하기 위한 인공삽입물의 크기, 위치, 골고정 나사 구멍의 위치 값을 산출하는 단계를 포함할 수 있다.According to an aspect of the present invention, there is provided a method of designing a customized prosthesis according to the present invention, comprising: generating a 3D (Three Dimensions) image using a bone image photographed in 2D; Measuring bone density information using the 2D bone image; Colorizing the bone density information on the 3D image; Using the colored 3D image, the method may include calculating a size, a position, and a position value of the bone fixation screw hole for designing the prosthesis.
상기 2D 뼈 이미지에서 골밀도를 측정하는 단계는 명암 차이로 표현되는 2D 영상의 명암값을 미리 설정된 경계값에 따라 측정하는 단계를 포함할 수 있다.Measuring the bone density in the 2D bone image may include measuring the contrast value of the 2D image represented by the contrast difference according to a preset boundary value.
상기 설계 방법은 상기 산출된 인공삽입물의 크기, 위치, 골고정 나사 위치 값을 이용하여 인공삽입물을 설계하는 단계를 더 포함할 수 있다.The design method may further include designing the prosthesis using the calculated size, position, and bone fixation screw position value of the prosthesis.
상기 색채화하는 단계는 골밀도 값에 따라 색채맵을 설정하는 단계와; 상기 설정된 색채맵 정보와, 포인트 좌표로 계산된 윤곽선 정보를 정합하는 단계와; 상기 생성된 3D 영상에 골밀도 색채맵을 색채화하는 단계를 포함할 수 있다.The coloring may include setting a color map according to a bone density value; Matching the set color map information with contour information calculated from point coordinates; The method may include colorizing the bone density color map on the generated 3D image.
상기 인공삽입물의 설계변수를 결정하는 단계는: 상기 색채화된 3D 영상에서, 인공삽입물 시술 부위의 중심, 크기, 위치 중 하나 이상을 추출할 수 있다.The determining of the design variable of the prosthesis may include: extracting at least one of a center, a size, and a position of the prosthesis procedure region from the colored 3D image.
상기 인공삽입물의 골고정 나사구멍의 위치 값을 산출하는 단계는: 상기 색채화된 3D 영상에서, 골밀도에 따라 정합된 색채맵을 이용하여, 골고정 나사구멍의 위치, 개수 중 적어도 하나 이상을 추출할 수 있다.The step of calculating a position value of the bone fixation screw hole of the prosthesis may include: extracting at least one of the position and the number of the bone fixation screw holes using the color map matched according to the bone density in the colored 3D image. can do.
본 발명에 따르면, CT 이미지에서 획득한 뼈의 윤곽선과 골밀도 정보를 이용하여 시각화된 골밀도 정보가 포함된 3차원 영상을 재구성하고 이를 이용하여 인공삽입물의 설계변수 및 골 고정나사의 위치정보를 제공함으로써, 환자의 골형상과 일치하는 인공삽입물의 설계가 가능하다. 이를 통하여 정합성을 높일 수 있고 골밀도가 높은 부위에 골 고정나사를 통한 고정이 이루어져 높은 고정력과 안정성을 통한 인공삽입물의 수명을 연장할 수 있으며, 시술의 성공률을 높일 수 있다.According to the present invention, by reconstructing the three-dimensional image including the bone density information visualized using the bone contour and bone density information obtained from the CT image and using this to provide the design variable of the prosthesis and the position information of the bone fixation screw For this reason, it is possible to design an implant that matches the bone shape of the patient. Through this, it is possible to increase the coherence and fixation through the bone fixation screw in the region of high bone density, prolong the life of the prosthesis through high fixation force and stability, and increase the success rate of the procedure.
특히 인공고관절 전치환술(Total hip replacement) 및 재치환술(Revision total hip replacement), 척추유합술(Spinal fusion), 본플레이트(Bone plate) 등의 시술시, 시술하고자 하는 위치에 대한 뼈의 형태학적 정보와 골밀도에 대한 생물학적 정보를 동시에 파악할 수 있음으로 인공삽입물의 정합성과 고정력을 향상시킬 수 있다.In particular, when the hip and total hip replacement, revision total hip replacement, spinal fusion, bone plate, etc. Simultaneous identification of biological information on bone density can improve the integrity and fixability of the prosthesis.
또한, 본 발명에서 제시하는 CT 이미지의 3차원 재구성을 통한 영상화와 골밀도의 시각화는 인공삽입물의 시술에 있어서 기존의 규격화된 인공삽입물의 사용으로 인하여 발생되는 골용해, 해리, 부정합 등의 문제점을 감소시킬 수 있다. 또한 골밀도를 고려한 골고정 나사의 설계가 가능하여 골지지대 부족, 고정부위의 제한으로 인공삽입물의 고정력 확보에 제한이 따르는 고령의 환자, 골감소증, 골다공증 환자의 경우에도 골밀도가 높은 부위를 찾아 골고정 나사로 고정할 수 있어서 인공삽입물의 고정력을 높여줄 수 있다.In addition, the imaging and visualization of bone density through the three-dimensional reconstruction of the CT image proposed in the present invention reduces problems such as osteolysis, dissociation, and misalignment caused by the use of the standardized prosthesis in the prosthesis procedure. You can. In addition, the bone fixation screw can be designed considering bone density, so the elderly patients, osteopenia, and osteoporosis patients who have limited bone support and limited fixation area can find a site with high bone density. Can be fixed to increase the fixing force of the prosthesis.
도 1은 본 발명의 일 실시예에 따라 인공삽입물의 설계 및 시뮬레이션 시스템의 전체 구성도를 나타낸 것이다.
도 2는 도 1에 도시된 영상 처리 장치 및 설계 장치를 구체적으로 나타낸 블록도이다.
도 3은 뼈 윤곽선이 추출된 이미지와 골밀도가 시각화된 이미지를 나타낸다.
도 4는 3차원 영상으로 재구성되는 것을 예시적으로 나타낸다.
도 5는 3차원으로 재구성된 영상에 골밀도 정보를 시각화하는 것을 나타낸다.
도 6은 인공 삽입물의 시술될 부위의 위치와 크기, 그리고 골고정 나사의 위치와 개수를 결정하는 것을 예시적으로 나타낸다.Figure 1 shows the overall configuration of the design and simulation system of the prosthesis according to an embodiment of the present invention.
FIG. 2 is a block diagram illustrating in detail the image processing apparatus and the design apparatus illustrated in FIG. 1.
3 shows an image in which bone contours are extracted and an image in which bone density is visualized.
4 exemplarily illustrates reconstruction into a 3D image.
5 illustrates visualization of bone density information in an image reconstructed in three dimensions.
6 exemplarily determines the location and size of the site to be treated and the location and number of bone fixation screws.
이하, 첨부된 도면을 참조하여 본 발명에 따른 바람직한 실시예를 상세히 설명하되, 도면 부호에 관계없이 동일하거나 유사한 구성 요소는 동일한 참조 번호를 부여하고 이에 대한 중복되는 설명은 생략하기로 한다. 이하의 설명에서 사용되는 구성요소에 대한 접미사 "부"는 명세서 작성의 용이함만이 고려되어 부여되거나 혼용되는 것으로서, 그 자체로 서로 구별되는 의미 또는 역할을 갖는 것은 아니다. 또한, 본 발명을 설명함에 있어서 관련된 공지 기술에 대한 구체적인 설명이 본 발명의 요지를 흐릴 수 있다고 판단되는 경우 그 상세한 설명을 생략한다. 또한, 첨부된 도면은 본 발명의 사상을 쉽게 이해할 수 있도록 하기 위한 것일 뿐, 첨부된 도면에 의해 본 발명의 사상이 제한되는 것으로 해석되어서는 아니 됨을 유의해야 한다.Hereinafter, exemplary embodiments of the present invention will be described in detail with reference to the accompanying drawings, wherein like reference numerals refer to like or similar elements throughout the several views, and redundant description thereof will be omitted. The suffix "part" for components used in the following description is given or mixed in consideration of ease of specification, and does not have meanings or roles that are distinguished from each other. In the following description, well-known functions or constructions are not described in detail since they would obscure the invention in unnecessary detail. It is to be noted that the accompanying drawings are only for the purpose of facilitating understanding of the present invention, and should not be construed as limiting the scope of the present invention with reference to the accompanying drawings.
이하, 도면을 참조하여 설명하면 다음과 같다.Hereinafter, with reference to the drawings as follows.
도 1은 본 발명의 일 실시예에 따라 인공삽입물의 설계 및 시뮬레이션 시스템의 전체 구성도를 나타낸 것이다.Figure 1 shows the overall configuration of the design and simulation system of the prosthesis according to an embodiment of the present invention.
도 1을 참조하여 알 수 있는 바와 같이, 본 발명의 일 실시예에 따른 시스템은 CT(Computed Tomography), MRI(Magnetic Resonance Imaging: 자기공명영상)와 같은 측정 장치(100)와, 3차원 영상 처리 장치(200)와, 설계 장치(300)를 포함한다. As can be seen with reference to Figure 1, the system according to an embodiment of the present invention is a
상기 CT, MRI와 같은 측정 장치(100)는 환자의 뼈를 단층 촬영하고, 촬영된 단층 이미지를 상기 3차원 영상 처리 장치(200)로 전달한다. 여기서 CT, MRI 영상은 뼈의 폭과 길이 형태를 정확히 보여주고 주변의 신경관이나 해부학적 구조의 위치를 파악할 수 있다는 장점이 있다. The
상기 3차원 영상 처리 장치(200)는 상기 단층 촬영 이미지를 기본으로 하여 각 섹션 별 골형상 단면의 윤곽선을 추출하고 이를 기반으로 3차원 이미지를 생성한다. 골형상은 기하학적으로 복잡하고 굴곡이 심하며, 개인적 편차가 존재하기 때문에 형상의 단순화에 따른 오차가 크게 발생한다. 따라서 2차원 CT영상을 기본으로 하여 각 섹션 별 골형상 단면의 윤곽선을 추출하고 이를 이용하여 3차원 보간을 통해 3차원 이미지를 생성한다.The three-dimensional
또한, 상기 3차원 영상 처리 장치(200)는 CT 이미지에서 골밀도를 측정해 시각화한다. In addition, the 3D
한편, 설계 장치(300)는 상기 완성된 3차원 뼈 영상에서 뼈의 형태학적 정보와 골밀도의 생물학적 정보를 이용하여 인공삽입물이 시술될 부위의 위치와 크기, 그리고 골고정 나사의 위치와 개수를 결정한다.Meanwhile, the
이하에서는, 도 2 내지 도 6을 참조하여 보다 구체적으로 설명하기로 한다.Hereinafter, with reference to FIGS. 2 to 6 will be described in more detail.
도 2는 도 1에 도시된 영상 처리 장치 및 설계 장치를 구체적으로 나타낸 블록도이다.FIG. 2 is a block diagram illustrating in detail the image processing apparatus and the design apparatus illustrated in FIG. 1.
도 2를 참조하여 알 수 있는 바와 같이, 상기 영상 처리 장치(200)는 CT 이미지 입력 단계(S210), 문턱값에 따른 윤곽선 검출 단계(S220), 골밀도 측정 및 상기 영상에 골밀도 분포에 따른 색채맵을 설정하는 단계(S230), 3차원 골 형상 재구성 단계(S240), 3차원 골 형상에 골밀도 시각화 단계(250)를 수행한다.As can be seen with reference to Figure 2, the
의료영상 정보는 DICOM format으로 표준화되어 그 정보를 기록하고 있다. CT 이미지 입력 단계(S210)는 DICOM reader를 이용하여 표준화된 DICOM 파일을 특별한 변환없이 불러올 수 있게 한다. Medical imaging information is standardized in the DICOM format and recorded. The CT image input step S210 allows a standardized DICOM file to be loaded without special conversion using a DICOM reader.
상기 윤곽선 검출 단계(S220)에 대해서 설명하면, 2차원 CT 영상은 256 단계의 그레이 스케일로 구성된 이미지이며, 윤곽선과 같은 의미 있는 데이터를 얻기 위해서는 영상처리 과정이 필요하다. 윤곽선을 추출하는 이미지 처리 과정은 512 x 512 픽셀의 2차원 CT 이미지에 행렬로 구성된 2차원 필터를 컨볼루션(Convolution) 시키고 기준값 이상의 출력값만을 윤곽선으로 추출하게 된다. 윤곽선을 검출하기 위해서는 함수의 변화분을 취할 수 있는 미분의 개념이 이용되며 Robert, Prewitt, Sobel 등의 알고리즘이 존재한다. 이러한 알고리즘들은 2차원 행렬로 구현되며 원 영상과 컨볼루션 시켜 필터의 특성에 따라 원하는 정보를 추출할 수 있다. 일반적으로 명도의 차이가 클수록 출력값은 커지게 되고, 반대로 명도의 차이가 두드러지지 않을수록 출력값은 작아짐으로 기준값 이상의 필터 출력값만을 골 윤곽선으로 검출할 수 있게 된다. Referring to the contour detection step (S220), the 2D CT image is an image composed of 256 gray scales, and an image processing process is required to obtain meaningful data such as an outline. In the image processing process of extracting the contour, a two-dimensional filter composed of a matrix is convolved into a two-dimensional CT image of 512 x 512 pixels, and only output values above the reference value are extracted as the contour. To detect the contour, the concept of differential that can take the variation of the function is used, and algorithms such as Robert, Prewitt, and Sobel exist. These algorithms are implemented in a two-dimensional matrix and can be convolved with the original image to extract the desired information according to the characteristics of the filter. In general, the greater the difference in brightness, the larger the output value. On the contrary, the less the difference in brightness, the smaller the output value, so that only the filter output value above the reference value can be detected as the bone outline.
다음으로, 골밀도 측정 및 상기 영상에 골밀도 분포에 따른 색채맵을 설정하는 단계(S230)에 대해서 설명하면 다음과 같다. Next, the step of measuring the bone density and setting the color map according to the bone density distribution in the image (S230) will be described.
골밀도의 정량적 측정을 위해 사용되는 방법은 방사선 흡수법, 이중에너지 방사선 측정법, 정량적 전산화 단층활영, 정량적 초음파와 정량적 자기공명영상 등이 있다. 이중 정량적 전산화 단층촬영 방법은 CT 촬영과 동일한 요건으로 피검자를 골밀도 측정용 플랫폼(예컨대 Phantom) 위에 누이고 촬영한 것으로 골밀도에 따라 영상에 명암 차이가 생기게 된다. 한편, 골밀도 측정을 위해 CT 이미지를 활용하는 방안에 대해 설명하면, 물체의 밀도차에 따라 조사된 x-ray가 투과된 정도가 달라져 밝고 어두운 부위로 나뉘어지며 이에 따라 골밀도가 높은 부위는 밝게 표시되고 골밀도가 낮은 부위는 어둡게 표시되므로, CT 영상으로부터 골밀도 분포에 관한 정보를 정량화하고 시각화할 수 있다.Methods used for the quantitative measurement of bone density include radiation absorption, dual energy radiometry, quantitative computed tomography, quantitative ultrasound and quantitative magnetic resonance imaging. Dual quantitative computed tomography (CT) has the same requirements as CT imaging, in which a subject is photographed while lying on a platform for measuring bone density (for example, a phantom). Meanwhile, if the CT image is used to measure the bone density, the degree of transmission of the irradiated x-ray varies depending on the density difference of the object, and thus the light is divided into bright and dark areas. Areas with low bone density are darkened so that information about bone density distribution can be quantified and visualized from CT images.
따라서, 본 발명에서 CT 이미지를 이용하는 골밀도를 측정 방법을 이용한다. 그리고, 골밀도의 분포를 색채로 표시함에 있어서, 최대 뼈 밀도 값부터 최소 뼈 밀도 값까지의 구간 별 색상을 지정한다. 여기서 골밀도가 높을수록 255에 가까우며 골밀도가 낮을수록 0에 가깝게 나타나게 되며, 0~255 사이의 명암 값의 일정한 경계치 값을 정하고 경계 값 마다 색을 정하여 2차원 윤곽선 검출된 영상에 색채화할 수 있다. 그러나, 골밀도가 매우 떨어지는 노인이나 골다공증 환자에 대해 선형적인 색채맵을 적용하면 낮은 골밀도의 색채로 표시되는 구간이 많아지게 되어 전체적으로 구분이 쉽지 않을 수 있는 단점이 있으므로, 각 환자의 골밀도 범위와 그 범위에서의 색상을 선택할 수 있도록 하는 것이 바림직하며, 이는 선형 보간되어 색상표로 나타내어질 수 있다.Therefore, in the present invention, a method for measuring bone density using a CT image is used. And, in displaying the distribution of bone density in color, the color for each section from the maximum bone density value to the minimum bone density value is specified. In this case, the higher the bone density, the closer to 255. The lower the bone density, the closer to 0. The constant threshold value of the contrast value between 0 and 255 is determined, and the color is determined for each boundary value to colorize the 2D contour detected image. However, when linear color maps are applied to elderly or osteoporosis patients with very low bone density, there are many disadvantages that can be difficult to distinguish as a whole because of the color section of low bone density. It is desirable to be able to select the color in the, which can be linearly interpolated and represented by a color table.
이와 같이 함으로써 3차원 골형상을 재구성하고(S240), 상기 구성된 영상에 골밀도 정보를 시각화한다(S250). 구체적으로, 상기 추출된 윤곽선 정보와 골밀도 정보는 3차원 모델링 가능한 포인트 좌표로 변환되고, 각 포인트 좌표들은 NURBS(Non-Uniform Rational B-Spline)방법을 이용하여 라인으로 보간되며 각 섹션 별 내/외곽 윤곽선으로 표현된다. 다음으로 각 섹션 별 윤곽선은 멀티 섹션 스윕(Multi section sweep) 방법으로 3차원 영상으로 재구성된다.In this way, the 3D bone shape is reconstructed (S240), and bone density information is visualized on the constructed image (S250). Specifically, the extracted contour information and bone density information are converted to point coordinates that can be modeled in three dimensions, and each point coordinates are interpolated into lines using a non-uniform Rational B-Spline (NURBS) method, and the inner and outer parts of each section are It is represented by an outline. Next, each section outline is reconstructed into a 3D image by a multi section sweep method.
한편, 상기 설계 장치(300)는 상기 재구성된 3차원 골 이미지를 통해 인공 삽입물의 크기와 삽입 위치를 결정하고(S310), 상기 시각화된 골밀도 이미지를 통해 골고정 나사 위치, 크기 개수를 결정할 수 있도록 한다(S320).On the other hand, the
그리고, 상기 설계 장치(300)는 상기 결정된 정보들을 바탕으로 환자의 골형상에 가장 정확하게 맞는 인공삽입물을 설계한다(S330).In addition, the
도 3은 뼈 윤곽선이 추출된 이미지와 골밀도가 시각화된 이미지를 나타낸다.3 shows an image in which bone contours are extracted and an image in which bone density is visualized.
도 3을 참조하여 알 수 있는 바와 같이, 256단계의 그레이스케일(Gray scale)로 구성된 CT 이미지에 일정한 경계치를 정하여 경계값 안에 들어오는 값들과 지정한 색을 매칭(Macthing)시켜 추출된 윤곽선 정보에 색채화를 통하여 나타내었다.As can be seen with reference to FIG. 3, a predetermined boundary value is determined on a CT image composed of 256 gray scales, color matching is performed on the contour information extracted by matching the values in the boundary value with a specified color. It is shown through.
도 4는 3차원 영상으로 재구성되는 것을 예시적으로 나타낸다.4 exemplarily illustrates reconstruction into a 3D image.
도 4를 참조하여 알 수 있는 바와 같이, 상기 추출된 윤곽선 정보는 멀티 섹션 스윕(Multi section sweep) 방법으로 3차원 영상으로 재구성한다. As can be seen with reference to FIG. 4, the extracted contour information is reconstructed into a 3D image by a multi section sweep method.
도 5는 3차원으로 재구성된 영상에 골밀도 정보를 시각화하는 것을 나타낸다.5 illustrates visualization of bone density information in an image reconstructed in three dimensions.
도 5를 참조하여 알 수 있는 바와 같이, 최종적으로 경계치의 명암값의 범위를 정하고, 3차원으로 재구성된 영상에 색채화를 통하여 골밀도 정보를 시각화한다.As can be seen with reference to FIG. 5, finally, the range of the contrast value of the boundary value is determined, and bone density information is visualized through colorization of the image reconstructed in three dimensions.
도 6은 인공 삽입물의 시술될 부위의 위치와 크기, 그리고 골고정 나사의 위치와 개수를 결정하는 것을 예시적으로 나타낸다.6 exemplarily determines the location and size of the site to be treated and the location and number of bone fixation screws.
인공삽입물에 사용되는 골고정 나사는 개수가 많을수록 높은 고정력을 얻게 되며 골밀도가 높은 부위에 위치하여야 한다. 이를 위하여 완성된 환자의 3차원 뼈 영상에서 뼈의 형태학적 정보와 골밀도의 생물학적 정보를 이용하여 인공삽입물이 시술될 부위의 위치와 크기, 그리고 골고정 나사의 위치와 개수를 결정한다. 결정된 인공삽입물의 설계 변수는 3차원 CAD 프로그램에서 불러들 수 있도록 변화되며, 이를 이용하여 골밀도를 고려한 맞춤형 인공삽입물을 설계한다. The larger the number of bone fixation screws used in the prosthesis, the higher the fixing force and the higher the bone density should be. For this purpose, the location and size of the implant site and the number and location of bone fixation screws are determined by using the morphological information of bone and biological information of bone density in the completed 3D bone image of the patient. The design parameters of the prosthesis are changed so that they can be loaded in the 3D CAD program. Using this, we design a custom prosthesis considering the bone density.
여기까지 설명된 본 발명에 따른 방법은 소프트웨어, 하드웨어, 또는 이들의 조합으로 구현될 수 있다. 예를 들어, 본 발명에 따른 방법은 소프트웨어 프로그램으로 구현되어, 상기 시뮬레이션 장치(200)의 저장 매체에 저장될 수 있고, 상기 시뮬레이션 장치(200)의 프로세서에 의해서 실행될 수 있다.The method according to the invention described thus far can be implemented in software, hardware, or a combination thereof. For example, the method according to the present invention may be implemented as a software program, stored in a storage medium of the
이상에서는 본 발명의 바람직한 실시예를 예시적으로 설명하였으나, 본 발명의 범위는 이와 같은 특정 실시예에만 한정되는 것은 아니므로, 본 발명은 본 발명의 사상 및 특허청구범위에 기재된 범주 내에서 다양한 형태로 수정, 변경, 또는 개선될 수 있다.While the present invention has been described with reference to exemplary embodiments, it is to be understood that the invention is not limited to the disclosed exemplary embodiments, but, on the contrary, May be modified, modified, or improved.
100: 측정 장치
200: 3차원 영상 처리 장치
300: 설계 장치100: measuring device
200: 3D image processing device
300: design device
Claims (6)
상기 2D 뼈 이미지를 이용하여 골밀도 정보를 측정하는 단계와;
상기 골밀도 정보를 상기 3D 영상에 색채화하는 단계와;
상기 색채화된 3D 영상을 이용하여, 인공 삽입물을 설계하기 위한 인공삽입물의 크기, 위치, 골고정 나사 구멍의 위치 값을 산출하는 단계를 포함하는 것을 특징으로 하는 맞춤형 인공삽입물 설계 방법.Generating a 3D (Three Dimensions) image using the bone image photographed in 2D;
Measuring bone density information using the 2D bone image;
Colorizing the bone density information on the 3D image;
Using the colored 3D image, calculating a size, a position, and a position value of a bone fixation screw hole for designing the prosthesis.
명암 차이로 표현되는 2D 영상의 명암값을 미리 설정된 경계값에 따라 측정하는 단계를 포함하는 것을 특징으로 하는 인공삽입물 설계 방법.The method of claim 1, wherein the step of measuring bone density in the 2D bone image
And measuring a contrast value of the 2D image represented by the contrast difference according to a preset boundary value.
골밀도 값에 따라 색채맵을 설정하는 단계와;
상기 설정된 색채맵 정보와, 포인트 좌표로 계산된 윤곽선 정보를 정합하는 단계와;
상기 생성된 3D 영상에 골밀도 색채맵을 색채화하는 단계를 포함하는 것을 특징으로 하는 인공삽입물 설계 방법.The method of claim 1, wherein the coloring step
Setting a color map according to the bone density value;
Matching the set color map information with contour information calculated from point coordinates;
And colorizing a bone density color map on the generated 3D image.
상기 인공 삽입물을 설계하기 위한 인공삽입물의 크기, 위치 값을 산출하는 단계는
상기 색채화된 3D 영상에서, 인공삽입물 시술 부위의 중심, 크기, 위치 중 하나 이상을 추출하는 것을 특징으로 하는 인공삽입물 설계 방법.The method of claim 1,
Computing the size and position of the prosthesis for designing the prosthesis is
The implant design method, characterized in that for extracting at least one of the center, the size, the position of the implant site in the colored 3D image.
상기 색채화된 3D 영상에서, 골밀도에 따라 정합된 색채맵을 이용하여, 골고정 나사구멍의 위치, 개수 중 적어도 하나 이상을 추출하는 것을 특징으로 하는 인공삽입물의 설계 방법.The method of claim 1, wherein calculating the position value of the bone fixation screw hole of the prosthesis
In the colored 3D image, using the color map matched according to the bone density, at least one or more of the position, the number of bone fixation screw holes is extracted.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020110012517A KR101179081B1 (en) | 2011-02-11 | 2011-02-11 | Apparatus for designing custom-made prosthesis based on bone density measurement |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| KR1020110012517A KR101179081B1 (en) | 2011-02-11 | 2011-02-11 | Apparatus for designing custom-made prosthesis based on bone density measurement |
Publications (2)
| Publication Number | Publication Date |
|---|---|
| KR20120092451A KR20120092451A (en) | 2012-08-21 |
| KR101179081B1 true KR101179081B1 (en) | 2012-09-03 |
Family
ID=46884478
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| KR1020110012517A KR101179081B1 (en) | 2011-02-11 | 2011-02-11 | Apparatus for designing custom-made prosthesis based on bone density measurement |
Country Status (1)
| Country | Link |
|---|---|
| KR (1) | KR101179081B1 (en) |
Families Citing this family (6)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| AU2014323515B2 (en) * | 2013-09-18 | 2018-08-23 | Stryker Corporation | Patient specific bone preparation for consistent effective fixation feature engagement |
| US20150119987A1 (en) * | 2013-10-28 | 2015-04-30 | Stryker Corporation | Implant design using heterogeneous bone properties and probabilistic tools to determine optimal geometries for fixation features |
| WO2017200785A1 (en) * | 2016-05-16 | 2017-11-23 | Think Surgical, Inc. | Implant design and computer assisted surgery |
| KR102313749B1 (en) * | 2019-12-23 | 2021-10-18 | 주식회사 메가젠임플란트 | Apparatus for Automatically Transforming Color of Computerized Tomography Images on Oral Cavity Based on Artificial Intelligence and Driving Method Thereof |
| KR102367095B1 (en) * | 2021-02-17 | 2022-02-23 | 김신재 | A method for preparing 3d bone model reflecting 2d image of target bone and a 3d bone model prepared thereby |
| KR102414601B1 (en) * | 2021-10-12 | 2022-07-04 | 에이아이다이콤 (주) | Method of Measuring Bone Mineral Density For Osteoporosis Diagnosis Based On Machine Learning And Bone Density Measurement Program Using The Same |
Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2719444B2 (en) | 1992-09-14 | 1998-02-25 | ルナー コーポレイション | Method and apparatus for automatically determining and analyzing bone morphology |
| JP2002523204A (en) | 1998-08-28 | 2002-07-30 | アーチ・デベロップメント・コーポレーション | Computer-based analysis method and apparatus for bone mass and structure |
-
2011
- 2011-02-11 KR KR1020110012517A patent/KR101179081B1/en active IP Right Grant
Patent Citations (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| JP2719444B2 (en) | 1992-09-14 | 1998-02-25 | ルナー コーポレイション | Method and apparatus for automatically determining and analyzing bone morphology |
| JP2002523204A (en) | 1998-08-28 | 2002-07-30 | アーチ・デベロップメント・コーポレーション | Computer-based analysis method and apparatus for bone mass and structure |
Also Published As
| Publication number | Publication date |
|---|---|
| KR20120092451A (en) | 2012-08-21 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US9681956B2 (en) | Acquiring and utilizing kinematic information for patient-adapted implants, tools and surgical procedures | |
| Helgason et al. | A modified method for assigning material properties to FE models of bones | |
| Taylor et al. | Accounting for patient variability in finite element analysis of the intact and implanted hip and knee: a review | |
| KR101179081B1 (en) | Apparatus for designing custom-made prosthesis based on bone density measurement | |
| Bryce et al. | Validation of three-dimensional models of in situ scapulae | |
| US20130188848A1 (en) | Automatic Implant Detection From Image Artifacts | |
| US20150250552A1 (en) | Advanced methods of modeling knee joint kinematics and designing surgical repair systems | |
| US20170258598A1 (en) | Advanced Methods of Modeling Knee Joint Kinematics and Designing Surgical Repair Systems | |
| Brodén et al. | Low-dose CT-based implant motion analysis is a precise tool for early migration measurements of hip cups: a clinical study of 24 patients | |
| ES2828728A1 (en) | METHOD TO OBTAIN AN IMAGE BIOMARKER THAT QUANTIFIES THE QUALITY OF THE TRABECULAR STRUCTURE OF THE BONES (Machine-translation by Google Translate, not legally binding) | |
| Durastanti et al. | Comparison of cartilage and bone morphological models of the ankle joint derived from different medical imaging technologies | |
| Koh et al. | Reconstruction of patient-specific femurs using X-ray and sparse CT images | |
| Galvez et al. | Error measurement between anatomical porcine spine, CT images, and 3D printing | |
| Svedmark et al. | A new CT method for assessing 3D movements in lumbar facet joints and vertebrae in patients before and after TDR | |
| Humbert et al. | 3D Evaluation of the acetabular coverage assessed by biplanar X-rays or single anteroposterior X-ray compared with CT-scan | |
| Boymans et al. | The femoral head center shifts in a mediocaudal direction during aging | |
| Oravec et al. | Digital tomosynthesis and high resolution computed tomography as clinical tools for vertebral endplate topography measurements: comparison with microcomputed tomography | |
| Parra et al. | Damage assessment of spinal bones due to prostate cancer | |
| Rooppakhun et al. | Advanced medical imaging and reverse engineering technologies in craniometric study | |
| Tho | Bone and joints modelling with individualized geometric and mechanical properties derived from medical images | |
| TWI693922B (en) | Preparation method of bone defect bone guide device | |
| Osman et al. | Contralateral breast volume measurement during chest CT for postmastectomy breast reconstruction | |
| Joseph et al. | Three dimensional reconstruction models for medical modalities: a comprehensive investigation and analysis | |
| CN110680353B (en) | Analysis device for femoral head integrity rate | |
| Crisco et al. | Three-dimensional joit kinematics using bone surface registration: A computer assisted approach with an application to the wrist joint in vivo |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| A201 | Request for examination | ||
| E902 | Notification of reason for refusal | ||
| E701 | Decision to grant or registration of patent right | ||
| GRNT | Written decision to grant | ||
| FPAY | Annual fee payment |
Payment date: 20150730 Year of fee payment: 4 |
|
| FPAY | Annual fee payment |
Payment date: 20160728 Year of fee payment: 5 |