CN103976754A - Automatic calculating method for three-dimensional acetabulum index - Google Patents
Automatic calculating method for three-dimensional acetabulum index Download PDFInfo
- Publication number
- CN103976754A CN103976754A CN201410206419.2A CN201410206419A CN103976754A CN 103976754 A CN103976754 A CN 103976754A CN 201410206419 A CN201410206419 A CN 201410206419A CN 103976754 A CN103976754 A CN 103976754A
- Authority
- CN
- China
- Prior art keywords
- dimensional
- index
- acetabular
- measured value
- bone
- 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.)
- Pending
Links
Abstract
The invention discloses an automatic calculating method for a three-dimensional acetabulum index. The method comprises the steps of expanding two-dimensional AI (artificial index) measurement to three-dimensional AI measurement, and allowing to obtain all AIs of the semispherical edge surrounding a whole acetabulum, so the lesion state of DDH (developmental dysplasia of the hip) can be more accurately represented. According to the automatic calculating method for the three-dimensional acetabulum AI, the calculation difference caused by the subjectivity of a clinical doctor in the traditional measuring process of the AI is avoided.
Description
Technical field
The invention belongs to hip dysplasia detection technique field, relate to a kind of automatic calculating method of three-dimensional acetabular index.
Background technology
Hip dysplasia (Developmental Dysplasia of the Hip, DDH) refers to because acetabular development defect causes acetabular bone bad to capital covering, is the commonly encountered diseases of child's cacomelia.Main manifestations is that the outer top of acetabular bone and the place ahead are damaged, and acetabular bone shoals, and hip joint moves outside center, also can be with hip joint subluxation in various degree, and dislocation or variations in skews.DDH is the common disease in current orthopaedics field, and different genetic background and not agnate sickness rate are different, and Neonatal Morbidity is 1.5～5 ‰ generally, at the average attack rate of U.S. neonate DDH roughly 1.5 ‰; At China's sickness rate, be 0.91 ‰～8.2 ‰, average out to 3.5 ‰.DDH patient mostly is infant and school age population, and such disease has a strong impact on patient's limb development and function and patient's quality of life.
Although there is multiple diagnosis and therapeutic modality for hip dysplasia, but because defect type and therapeutic modality complexity and children's's characteristics of lesion have individuality, only depend on traditional two dimension (2D) image and two-dimentional index, be difficult to accurately reflect locus and the anatomic information of patient's diseased region, therefore the acceptance rate not high (being about 60～70%) of performing the operation, femoral head is dislocated again, the incidence rate of the necrosis of femoral head epiphysis and the complication such as hip joint is stiff is higher.
Summary of the invention
The present invention realizes by following technical solution:
1, an automatic calculating method for three-dimensional acetabular index, is characterized in that, comprises following operation:
1) processor, according to the CT of hip joint or MRI data, extracts bone information based on PV-MAP method, and generates corresponding 3-D view;
2) in three dimensions, the central point of the Y shape cartilage of ilium both sides inwall is demarcated, then using the horizontal line Y-Y that connects two central points as rotating shaft, through rotating shaft, set up one with the crossing shear surface in acetabular bone edge; In shear surface, the angle θ by the straight line A-B of Y shape cartilage central point and acetabular bone marginal point formation and the horizontal line composition of Y-Y, carries out the acetabular index AI measured value that slope calculates this shear surface, obtains two-dimentional acetabular index subsequently;
3) in 3-D view space, when when rotating shaft rotates, to each the rotational shear face repeating step 2 forming) calculating, can obtain the AI measured value of 0 °～180 ° all rotational shear faces in rotating range, obtain three-dimensional acetabular index;
4) will calculate resulting AI measured value at different rotary shear surface and be mapped in two-dimensional direct angle coordinate system, obtain the AI distribution curve around acetabular bone hemisphere face edge, and by its output.
Describedly AI measured value is mapped to two-dimensional direct angle coordinate system is: X-axis is described the anglec of rotation from small to large successively; Y-axis is described each anglec of rotation and is calculated resulting three-dimensional acetabular index.
Compared with prior art, the present invention has following useful technique effect:
The present invention, by patient's hip joint image data is analyzed, proposes a kind of new three-dimensional acetabular index, and has set up corresponding three-dimensional quantitative measurement method, to reflect more accurately lesion locations and the degree of DDH, and the auxiliary detection that improves hip dysplasia.
The present invention expands to three-dimensional by two-dimentional AI index measurement, allows to obtain all AI indexes around whole acetabular bone hemisphere face edge, thereby can show more accurately its developmental condition feature.
The automatic account form of the 3D AI index that the present invention sets up, has avoided the calculated difference that in traditional 2D AI index measurement process, subjectivity brings, and the information that has more objectivity can be provided.
Accompanying drawing explanation
Fig. 1 is Y shape cartilage schematic diagram;
Fig. 2 is three-dimensional AI Index for Calculation flow chart;
Fig. 3 is the bone segmentation result based on PV-MAP algorithm.(a)-(c) be original CT image; (d)-(f) be bone segmentation result;
Fig. 4 is the quantitative measurement of three-dimensional AI index.(a) for Y-Y point is determined; (b) be the AI measured value of orthogonal plane; (c) be 3D AI exponential curve (upper curve represents left side exponential, and lower curve represents right side exponential);
Fig. 5 is the three-dimensional AI exponential schematic diagram (20 °-150 °) of normal person and DDH patient.(a) be normal group (N=13); (b) be abnormal group (N=15).
The specific embodiment
For the diagnosis of hip dysplasia, the quantitative measurement of acetabular index (Acetabulum Index, AI) is a very effective supplementary mode.For example: when patient's acetabular index is greater than 40 °, need to implement operative treatment.By quantitative analysis, can judge acetabular bone sophistication, also can provide foundation to the evaluation of aftertreatment effect simultaneously.At present under clinic diagnosis pattern, measure conventionally directly on patient's x-ray mating plate or CT faultage image and carry out, what obtain is only the acetabular index value of the two-dimentional coronalplane vertical with horizontal plane, as shown in Figure 1, is difficult to accurately damaged position and the degree of reflection.The three-dimensional configuration of considering acetabular bone changes, and adopts three-dimensional acetabular index can make doctor obtain all changes and the damaged situation around whole acetabular bone hemisphere face edge, thereby can show more accurately its developmental condition and anatomical features.Such as the 3D quantitative Analysis to AI, contribute to reflect the mutual relation between acetabular bone and femur, the relative distance that acetabular bone arrives nearly femoral head central point etc.Yet, lack clinically accurate definition and the measuring method to three-dimensional index at present, for its measurement has brought difficulty.The present invention is by measuring concept extension by two-dimentional AI, utilize projection variation and the coordinate translation technology of graphics to shine upon the organizational structure in three dimensions, based on CT or MRI faultage image, in conjunction with the mutual spatial relationship between acetabular bone and femur, a kind of new 3D AI index definition has been proposed, and corresponding automatic calculating method.
In three dimensions, as shown in FIG. 1 to 3, first the central point of the Y shape cartilage of ilium both sides inwall is demarcated.Then the horizontal line (Y-Y line) that connects two central points is defined as to rotating shaft, through rotating shaft, set up one with the crossing shear surface in acetabular bone edge.In shear surface, by bone central point and acetabular bone the marginal point straight line (A-B line) forming and the angle (θ) that Y-Y horizontal line forms, carry out slope calculating and obtain AI measured value subsequently.In 3-D view space, when when rotating shaft (Y-Y axle) rotates, each the rotational shear face forming is repeated to above-mentioned calculating, obtain all AI values in 0 °～180 ° rotating ranges.Finally these AI values are mapped in two-dimensional coordinate system, by the three-dimensional AI distribution curve obtaining around acetabular bone hemisphere face edge, as shown in Figure 4.
The automatic calculating side method of the 3D AI index that the present invention proposes, comprises following operation
1) read in CT or MRI data, based on PV-MAP algorithm, extract bone information;
Processor, according to the CT of hip joint or MRI data, extracts bone information based on PV-MAP method, and generates corresponding 3-D view;
2), in 3d space, establish Y-Y point; In O degree plane, according to the definition of AI index, calculate AI value;
In three dimensions, the central point of the Y shape cartilage of ilium both sides inwall is demarcated, then using the horizontal line Y-Y that connects two central points as rotating shaft, through rotating shaft, set up one with the crossing shear surface in acetabular bone edge; In shear surface, the angle θ by the straight line A-B of Y shape cartilage central point and acetabular bone marginal point formation and the horizontal line composition of Y-Y, carries out the acetabular index AI measured value that slope calculates this shear surface, obtains two-dimentional acetabular index subsequently;
3) take Y-Y as rotating shaft, each the rotational shear face forming is repeated to above-mentioned calculating, obtain 0 °～180 ° AI values;
In 3-D view space, when when rotating shaft rotates, to each the rotational shear face repeating step 2 forming) calculating, can obtain the AI measured value of 0 °～180 ° all rotational shear faces in rotating range, obtain three-dimensional acetabular index;
4) the AI value of generation is mapped in two-dimensional coordinate system, obtains three-dimensional AI distribution curve.
To calculate resulting AI measured value at different rotary shear surface and be mapped in two-dimensional direct angle coordinate system, obtain the AI distribution curve around acetabular bone hemisphere face edge, and by its output.
For the performance of assessment algorithm, in No.1 Hospital Attached to No.4 Army Medical College, select 4 volunteers, 10 patients to carry out the comparative research of Ipsilateral and strong side acetabular development situation.In 10 patients, abnormal 2 examples in left side, abnormal 3 examples in right side, abnormal 5 examples of bilateral; January to 42 years old age year old.In research process, use PQ6000 type multi-layer spiral CT to obtain sample data to patient's scanning, sweep limits is from the front lower ridge of ilium to lesser trochanter.Design parameter is: scanning bed thickness 3mm, rebuilds bed thickness 2mm, sweep parameter 120kV, 70-120mA (body size by patient determines).Experimental result shows:
1) doctor's manual measurement value and system automatic measurements are without marked difference, as shown in table 1;
2) as shown in Figure 5, contrast 2D AI value, 3D AI measured value provides the more information of horn of plenty, thereby can doctor improve the accuracy of diagnosis and treatment.
In the present invention, it is vital that skeleton in CT/MRI image is cut apart accurately, because its basis that to be quantitative Analysis realize, cuts apart quantitative analysis that error can carry out to the later stage and the simulation of surgical planning and brings gross error.For local body (the Partial Volume existing in image, PV) effect, the impact of especially weak edge effect, maximum posteriori criterion (the Maximum a Posteriori of employing based on PV mixed model, MAP) to the Image Segmentation Using obtaining, segmentation result as shown in Figure 3.
Table 1, positive side position doctor's manual measurement and system are measured AI value relatively automatically
Claims (2)
1. an automatic calculating method for three-dimensional acetabular index, is characterized in that, comprises following operation:
1) processor, according to the CT of hip joint or MRI data, extracts bone information based on PV-MAP method, and generates corresponding 3-D view;
2) in three dimensions, the central point of the Y shape cartilage of ilium both sides inwall is demarcated, then using the horizontal line Y-Y that connects two central points as rotating shaft, through rotating shaft, set up one with the crossing shear surface in acetabular bone edge; In shear surface, the angle θ by the straight line A-B of Y shape cartilage central point and acetabular bone marginal point formation and the horizontal line composition of Y-Y, carries out the acetabular index AI measured value that slope calculates this shear surface, obtains two-dimentional acetabular index subsequently;
3) in 3-D view space, when when rotating shaft rotates, to each the rotational shear face repeating step 2 forming) calculating, can obtain the AI measured value of 0 °～180 ° all rotational shear faces in rotating range, obtain three-dimensional acetabular index;
4) will calculate resulting AI measured value at different rotary shear surface and be mapped in two-dimensional direct angle coordinate system, obtain the AI distribution curve around acetabular bone hemisphere face edge, and by its output.
2. the automatic calculating method of three-dimensional acetabular index as claimed in claim 1, is characterized in that, describedly AI measured value is mapped to two-dimensional direct angle coordinate system is: X-axis is described the anglec of rotation from small to large successively; Y-axis is described each anglec of rotation and is calculated resulting three-dimensional acetabular index.
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410206419.2A CN103976754A (en) | 2014-05-15 | 2014-05-15 | Automatic calculating method for three-dimensional acetabulum index |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN201410206419.2A CN103976754A (en) | 2014-05-15 | 2014-05-15 | Automatic calculating method for three-dimensional acetabulum index |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN103976754A true CN103976754A (en) | 2014-08-13 |
Family
ID=51269102
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN201410206419.2A Pending CN103976754A (en) | 2014-05-15 | 2014-05-15 | Automatic calculating method for three-dimensional acetabulum index |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN103976754A (en) |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108348193A (en) * | 2015-08-05 | 2018-07-31 | F·贝特纳 | A technique for for measures perform the operation in cup top rake based on fluoroscopy |
| WO2023198118A1 (en) * | 2022-04-14 | 2023-10-19 | Beijing Zhuoyue Medical Technology Co., Ltd | Method and apparatus for determining acetabulum-to-femoral-head distance |
Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102204829A (en) * | 2011-05-09 | 2011-10-05 | 中国人民解放军第三军医大学第一附属医院 | Method for measuring central bearing diameter of hip joint |
-
2014
- 2014-05-15 CN CN201410206419.2A patent/CN103976754A/en active Pending
Patent Citations (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN102204829A (en) * | 2011-05-09 | 2011-10-05 | 中国人民解放军第三军医大学第一附属医院 | Method for measuring central bearing diameter of hip joint |
Non-Patent Citations (3)
| Title |
|---|
| 刘欣等: "骸关节发育不良的计算机辅助检测与治疗", 《中国生物医学工程进展—2007中国生物医学工程联合学术年会论文集（上册）》, 30 April 2007 (2007-04-30) * |
| 赵黎等: "骨性髋臼指数的三维测量与分析", 《中国矫形外科杂志》, vol. 18, no. 15, 31 August 2010 (2010-08-31) * |
| 颉强等: "三维CT 评估发育性髋关节脱位的骨性形态学特征", 《中国组织工程研究与临床康复》, vol. 11, no. 23, 10 June 2007 (2007-06-10) * |
Cited By (2)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN108348193A (en) * | 2015-08-05 | 2018-07-31 | F·贝特纳 | A technique for for measures perform the operation in cup top rake based on fluoroscopy |
| WO2023198118A1 (en) * | 2022-04-14 | 2023-10-19 | Beijing Zhuoyue Medical Technology Co., Ltd | Method and apparatus for determining acetabulum-to-femoral-head distance |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| JP6833912B2 (en) | Bone reconstruction and orthopedic implants | |
| Dandachli et al. | Analysis of cover of the femoral head in normal and dysplastic hips: new CT-based technique | |
| CN110522501B (en) | 3D printing personalized orthopedic implant construction and biomechanical optimization processing method | |
| CN108711187A (en) | The method that registration fusion CT and MRI signal establish human body lumbar vertebrae three-dimensional simulation model | |
| CN105832415B (en) | Navigation device of rotation and preparation method thereof and application method are precisely gone for femur | |
| Gelaude et al. | Quantitative computerized assessment of the degree of acetabular bone deficiency: total radial acetabular bone loss (TrABL) | |
| CN113069134A (en) | Three-dimensional measurement method for morphological parameters in femoral bone marrow cavity based on CT three-dimensional reconstruction image | |
| Chougule et al. | Development of patient specific implants for Minimum Invasive Spine Surgeries (MISS) from non-invasive imaging techniques by reverse engineering and additive manufacturing techniques | |
| CN113069135A (en) | Three-dimensional measurement method for morphological parameters of outer side wall of femur based on CT three-dimensional reconstruction image | |
| Wang et al. | Two-dimensional and three-dimensional cup coverage in total hip arthroplasty with developmental dysplasia of the hip | |
| CN103976754A (en) | Automatic calculating method for three-dimensional acetabulum index | |
| WO2023198119A1 (en) | System and method for determining acetabular morphology | |
| CN110811829B (en) | Construction method and system based on femoral rotation axis and varus analysis model | |
| WO2023198118A1 (en) | Method and apparatus for determining acetabulum-to-femoral-head distance | |
| Hsieh et al. | Three-dimensional hip morphology analysis using CT transverse sections to automate diagnoses and surgery managements | |
| Trajanović et al. | Morphometric analysis of the hip bone as the basis for reverse engineering | |
| Song et al. | Computer-aided modeling and morphological analysis of hip joint | |
| Albers et al. | Computer assisted diagnosis and treatment planning of femoroacetabular impingement (FAI) | |
| Shao et al. | Quantitative analysis of bony birth canal for periacetabular osteotomy patient by template fitting | |
| Chen et al. | A Modular Software System for 3D Measurement of Acetabulum in Arthritic Dysplasia Hips | |
| Donno | Statistical shape modelling of osteoarthritic knee joint | |
| Du et al. | 3D Automated Quantitative Calculations Based on CT Images of the Hip Joint | |
| Pradhan et al. | Role of Artificial Intelligence in 3-D Bone Image Reconstruction: A Review | |
| Pradhan et al. | Machine Learning in Radio Imaging | |
| Jayakar et al. | Cortical volumetry using 3D reconstruction of metacarpal bone from multi-view images |
Legal Events
| Date | Code | Title | Description |
|---|---|---|---|
| C06 | Publication | ||
| PB01 | Publication | ||
| C10 | Entry into substantive examination | ||
| SE01 | Entry into force of request for substantive examination | ||
| RJ01 | Rejection of invention patent application after publication |
Application publication date: 20140813 |
|
| RJ01 | Rejection of invention patent application after publication |