Inter- and Intra-Observer Variability of the Volume of Cervical Ossification of the Posterior Longitudinal Ligament Using Medical Image Processing Software

Article information

J Korean Neurosurg Soc. 2017;60(4):441-447
Publication date (electronic) : 2017 June 29
doi : https://doi.org/10.3340/jkns.2015.0708.014
1Department of Neurosurgery, Yonsei University College of Medicine, Seoul, Korea
2Spine and Joint Research Institute, Guro Cham Teun Teun Hospital, Seoul, Korea
3Department of Neurosurgery, Kangbuk Samsung Hospital, Sungkyunkwan University College of Medicine, Seoul, Korea
Address for reprints: Hyunchul Shin, M.D., Ph.D., Department of Neurosurgery, Kangbuk Samsung Hospital, Sungkyunkwan University College of Medicine, 29 Saemunan-ro, Jongno-gu, Seoul 03181, Korea, Tel: +82-2-2001-1791, Fax: +82-2-2001-2157, E-mail: nsdrshin@gmail.com
*These authors contributed equally to this work.
Received 2015 August 4; Revised 2016 January 27; Accepted 2016 February 26.

Abstract

Objective

Computed tomography (CT)-based method of three dimensional (3D) analysis (MIMICS®, Materialise, Leuven, Belgium) is reported as very useful software for evaluation of OPLL, but its reliability and reproducibility are obscure. This study was conducted to evaluate the accuracy of MIMICS® system, and inter- and intra-observer reliability in the measurement of OPLL.

Methods

Three neurosurgeons independently analyzed the randomly selected 10 OPLL cases with medical image processing software (MIMICS®) which create 3D model with Digital Imaging and Communication in Medicine (DICOM) data from CT images after brief explanation was given to examiners before the image construction steps. To assess the reliability of inter- and intra-examiner intraclass correlation coefficient (ICC), 3 examiners measured 4 parameters (volume, length, width, and length) in 10 cases 2 times with 1-week interval.

Results

The inter-examiner ICCs among 3 examiners were 0.996 (95% confidence interval [CI], 0.987–0.999) for volume measurement, 0.973 (95% CI, 0.907–0.978) for thickness, 0.969 (95% CI, 0.895–0.993) for width, and 0.995 (95% CI, 0.983–0.999) for length. The intra-examiner ICCs were 0.994 (range, 0.991–0.996) for volume, 0.996 (range, 0.944–0.998) for length, 0.930 (range, 0.873–0.947) for width, and 0.987 (range, 0.985–0.995) for length.

Conclusion

The medical image processing software (MIMICS®) provided detailed quantification OPLL volume with minimal error of inter- and intra-observer reliability in the measurement of OPLL.

INTRODUCTION

Ossification of the posterior longitudinal ligament (OPLL) of the cervical spine is one of the most important causes of cervical myelopathy, and it is more commonly observed in Asian population1,4,8,10,12,14,15,17,19,24). The size of OPLL is one of important factor relating to spinal canal stenosis and causing myelopathy16). So, evaluating the size and growth of OPLL is crucial to determine the timing of the operation and the risk factors for rapid progression of OPLL. Previous attempts have been made to measure the size of OPLL using the Tsuyama system, but, these methods provided limited information by involving two-dimensional (2D) imaging3,5,9,13,22,23). The technical improvements in computed tomography (CT)-based three-dimensional (3D) imaging analysis have made accurate 3D measurement of OPLL16,20). Recently, a study reported that novel CT-based method of 3D analysis (MIMICS®, Materialise, Leuven, Belgium) provided the detailed OPLL classification and quantification in the ossified volume11). This method appeared to be very useful for evaluation of OPLL, but accidental errors could be occurred because the identification of ossification was not completely automatic. So, its reliability and reproducibility are obscure. Hence, the present study was designed to evaluate the accuracy of MIMICS® system to evaluate the OPLL, and inter- and intra-observer reliability of this method in the volume of OPLL were measured.

MATERIALS AND METHODS

Study design and subjects enrolment

A total of 10 cases (5 male and 5 female) that underwent surgery for cervical OPLL from 1999 to 2013 was randomly selected from the OPLL database in this analysis. This study was approved by the institutional review board of university hospital and the informed consent was obtained from all patients. Three neurosurgeons participated in this study as examiners. They had no prior experience or training of MIMICS® software which create 3D model with digital imaging and communications in medicine (DICOM) data from CT images. Only the brief explanation was given to theses examiners just before the image construction steps.

Image acquisition and data transfer

CT images of the cervical spine were obtained using a CT scanner (Brivo; GE Healthcare, Waukesha, WI, USA) under uniform conditions (slice thickness, 1.5 mm; voltage, 120 kV; current, 75 mA). The images were saved in the DICOM format and transmitted into a workstation (Z820; HP, Palo Alto, CA, USA) via PACS. Then medical image processing software (MIMICS®) imported the data for 3D reconstruction and measurement. The time was calculated from the step of importing CT data into the software to the end of measurement.

3D reconstruction and measurement

The process consisted of 4 steps (Fig. 1). First, the spine was segmented by threshold which is the function of including only those pixels of the image with a value higher than or equal to the threshold value (segmentation). Second, the segmented object was further separated into the vertebra and the OPLL by making another mask on the OPLL and removing floating pixels (region growing). Whereas initial segmentation is easily achieved by threshold function on a set threshold value, further separation needs diligent work of removing the vertebra attached to the OPLL. The separation was difficult due to obscure boundary between the vertebra and the OPLL. Third, a 3D object of OPLL was reconstructed (3D reconstruction) (Fig. 2). Fourth, the parameters were measured on the processed images (measurement) (Fig. 3). The width was defined as the maximum distance along the axis of the coronal plane of OPLL. The thickness was defined as the maximum distance along the axis of the sagittal plane of OPLL. The length was defined as the longest dimension of OPLL. The measurement of width and thickness was made on the axial CT images. The length was measured on the 3D reconstructed objects. The volume was calculated automatically by the software. A total of 240 measurements (4 parameters for each of the 10 patients by 3 examiners with 2 times) were performed.

Fig. 1

Steps for 3D reconstruction and measurement. A: CT scanning with 3 mm thickness. B: The obtained images are transmitted digitally via PACS. C: 3D software imports CT data in the DICOM format. D: The segmented object (visualized by a colored-mask) is obtained by including only those pixels of the image with a value higher than or equal to the threshold value (thresholding). E: The segmented object is further split by making another mask on OPLL and removing floating pixels (region growing). F: 3D representation is created and the volume is calculated automatically by the software. CT: computed tomography, DICOM: Digital Imaging and Communications in Medicine, OPLL: ossification of the posterior longitudinal ligament.

Fig. 2

Reconstructed 3D images of OPLL after the calculation of the 3D. A: Anterior view. B: Lateral view. C: Posterior view. D: A volumetric mesh with the material assignment for finite element analysis. OPLL: ossification of the posterior longitudinal ligament.

Fig. 3

Measurements. A: The width was defined as the maximum distance along the axis of the coronal plane of OPLL and the thickness was defined as the maximum distance along the axis of the sagittal plane of OPLL. B: The length was defined as the longest dimension of OPLL. C: The volume was calculated automatically by the software in the 3D-reconstructed image. OPLL: ossification of the posterior longitudinal ligament.

Validation of the measurement

To assess the reliability of the measurement, 3 examiners measured 4 parameters in 10 cases 2 times with 1 week interval. Inter-examiner and intra-examiner intraclass correlation coefficient (ICC) and 95% CI were assessed. Inter-examiner reliability was assessed by the ICC of data obtained for the measurement of each parameter among all three examiners. Intra-examiner reliability was also assessed between the first and second measurements by all three examiners.

Statistical analysis

A 95% confidence level was used for all tests. All data were analyzed by SPSS ver. 18 (SPSS Inc., Chicago, IL, USA). Values are reported as the mean±standard deviation. Each of statistical method is described at the end of each result. A p-value of <0.05 was considered statistically significant.

RESULTS

Measurement process

The sample consisted of 10 subjects (5 male and 5 female). The mean age was 49.5±10.9 years (range, 36–63 years). The mean time for measuring each case was 23.7±4.3 minutes (range, 15–30 minutes). The measurement speed became quicker from 24.7±4.1 minutes initially to 22.7±4.3 minutes afterwards (p=0.003). The rate-limiting step was the separation step to split the vertebra and OPLL. The mean volume was 1915.2±1027.3 mm3 (range, 780.2–4123.1 mm3). The mean width was 11.0±2.9 mm (range, 5.4–15.9 mm). The mean thickness was 3.6±1.3 mm (range, 2.0–6.3 mm). The mean length was 27.6±10.4 mm (range, 11.2–49.5 mm). There was no significant differences in volume, width, thickness, and length among the examiners (p=0.969; p=0.539; p=0.712; p=0.878).

Validation of the measurement

The ICC among three examiners for the measurement of volume was 0.996 (95% CI, 0.987–0.999). The ICC for thickness was 0.973, and the CI was 0.907–0.978. The ICC for width was 0.969, and the CI was 0.895–0.993. The ICC for length was 0.995, and the CI was 0.983–0.999. The ICC between the two measurement by the same examiner was 0.994 (range, 0.991–0.996) for volume, 0.996 (range, 0.944–0.998) for length, 0.930 (range, 0.873–0.947) for width, and 0.987 (range, 0.985–0.995) for length, indicating high intra-examiner reliability.

DISCUSSION

Measurements of ossification length and thickness using 2D X-ray images are possible to a certain extent3,5,13,21), but, 2D evaluation with plain lateral radiographs has some limitations. First, its accuracy is reduced due to errors associated with the imaging procedure, because the ossification width and volume are not measurable. Ossification not identifiable by X-ray examination can be evaluated on CT images, and therefore, the classification discrimination may differ from that obtained by X-ray imaging. Second, the growth can occur in any direction, only craniocaudal and ventrodorsal progression are depicted on plain radiographs. Growth in the oblique direction is projected only orthographically onto a radiograph. Indeed, multidimensional evaluation of ossification can be achieved with CT images, but it is difficult to evaluate continuity on 2D X-ray images. Other disadvantages of X-ray analysis include inability to identify small-area ossifications and to accurately locate the ossification.

Compared with radiography, CT scanning has several advantages. CT scanning is the most sensitive diagnostic method for detecting small ossifications or calcifications of the ligament, which are likely to be missed on radiographs18). However, con ventional CT scanning also has some limitation of slice thickness. Measurement of OPLL on conventional radiographs is unreliable because the slice is likely to be thick, owing to the limited number of films used. It is also difficult to slice OPLL on the same level and angle to compare past and present ossification. However, these limitations can be overcome using helical scanning with multi-detector CT and digital viewers, but, its reliability and reproducibility are obscure2).

Three-dimensional evaluation with CT scanning is better suited to measuring OPLL growth. OPLL could be masked by the shoulder girdle shadows in the lower cervical spine (Fig. 4), and might be less distinct because the ossification is less densely calcified in young patients7). CT based method of 3D analysis can precisely evaluate the volume of OPLL at one point and provide information about the volume of OPLL with minimal analytical error. Evaluation of continuity and classification of ossification are comparatively simple in 3D images2). With CT images, changes in ossification form and thickness on the caudal side can be evaluated in detail in patients with ossifications growing toward the cranial side. Slight graininess of the images associated with slice thickness is seen at ossification boundaries. Hence, the accurate superimposition by voxel based registration could provide comparison of OPLL before and after surgery. Indeed, it is also possible to derive an accurate numerical estimate of both overall absolute value of the volume of ossification and its rate of progression over time, and to compare these volumes with those derived from subsequent scans by calculating exact volumes.

Fig. 4

Representative images obtained in a 47-year-old woman with cervical OPLL. A: Lateral X-ray showing an ossified lesion at the level of C3–4, however, the cervico-thoracic junction is not visualized adequately due to shoulder shadow. B: Sagittal-reconstruction CT showing mixed type OPLL. C: Axial CT image at the level of C4 demonstrating central OPLL. OPLL: ossification of the posterior longitudinal ligament, CT: computed tomography.

Izumi et al.11) previously reported CT-based 3D analysis method to measure volume changes in OPLL. The subjects were 20 OPLL patients who were being followed using the MIMICS® software (slice thickness, 1.25 mm; voltage, 120 kV; current, 178 mA or slice thickness, 0.5–1.0 mm; voltage, 120 kV; current, 75 mA) to calculate the volume. They concluded that this 3D method used allowed detailed OPLL classification and quantification of change in the ossified volume and appears to be very useful for quantitative evaluation of OPLL with only minimal measurement error. The ICC calculated in this study was obviously high, suggesting that the evaluation of ossification volume was accurate and valid, but, they also commented that accidental errors might have occurred because this method used for identification of ossification was not completely automatic.

In the study to evaluate the inter- and intra-observer reliability of lateral radiograph, axial CT, 2D and 3D reconstructed CT images based OPLL classification, the inter- and intra-observer ICCs were only 0.51 and 0.67 for the lateral radiograph, even in combination with the axial CT images, 0.70 and 0.85 for 2D CT images, and 0.76 and 0.86 for 3D CT images, respectively2). These values showed a good to excellent range for the 2D and 3D reconstructed CT images while those of the lateral radiograph indicated a fair range. These ICCs could improve the correlation by using software system (Visualization Toolkit (http://www.vtk.org/) as 0.968 (95% CI, 0.880–0.992) for inter-observer reproducibility and 0.987 (95% CI, 0.953–0.997) for intra-observer reproducibility of volume measurement6). Recent CT-based method of 3D analysis using MIMICS® software also provided the detailed OPLL quantification in the ossified volume with excellent inter- (>0.969) and intra-observer reproducibility (>0.930) although the examiners were not professional users. Similarly result was also reported in the study using MIMICS® with 0.856 of inter-observer ICC and 0.999 of intra-observer ICC. Therefore, this measurement method seems to be universally applicable method using similar viewers with equivalent function.

CONCLUSION

The medical image processing software (MIMICS®) provided detailed quantification of OPLL volume with minimal error of inter- and intra-observer reliability in the measurement of OPLL.

Inter-examiner reliability of each measurement

Intra-examiner reliability of each measurement

References

1. Bakay L, Cares HL, Smith RJ. Ossification in the region of the posterior longitudinal ligament as a cause of cervical myelopathy. J Neurol Neurosurg Psychiatry 33:263–268. 1970;
2. Chang H, Kong CG, Won HY, Kim JH, Park JB. Inter- and intra-observer variability of a cervical OPLL classification using reconstructed CT images. Clin Orthop Surg 2:8–12. 2010;
3. Chiba K, Kato Y, Tsuzuki N, Nagata K, Toyama Y, Iwasaki M, et al. Computer-assisted measurement of the size of ossification in patients with ossification of the posterior longitudinal ligament in the cervical spine. J Orthop Sci 10:451–456. 2005;
4. Chiba K, Ogawa Y, Ishii K, Takaishi H, Nakamura M, Maruiwa H, et al. Long-term results of expansive open-door laminoplasty for cervical myelopathy--average 14-year follow-up study. Spine (Phila Pa 1976) 31:2998–3005. 2006;
5. Chiba K, Yamamoto I, Hirabayashi H, Iwasaki M, Goto H, Yonenobu K, et al. Multicenter study investigating the postoperative progression of ossification of the posterior longitudinal ligament in the cervical spine: a new computer-assisted measurement. J Neurosurg Spine 3:17–23. 2005;
6. Choi JH, Shin JJ, Kim TH, Shin HS, Hwang YS, Park SK. Does intra-medullary signal intensity on MRI affect the surgical outcomes of patients with ossification of posterior longitudinal ligament? J Korean Neurosurg Soc 56:121–129. 2014;
7. Fujimori T, Iwasaki M, Nagamoto Y, Ishii T, Sakaura H, Kashii M, et al. Three-dimensional measurement of growth of ossification of the posterior longitudinal ligament. J Neurosurg Spine 16:289–295. 2012;
8. Han IH, Suh SH, Kuh SU, Chin DK, Kim KS. Types and prevalence of coexisting spine lesions on whole spine sagittal MR images in surgical degenerative spinal diseases. Yonsei Med J 51:414–420. 2010;
9. Hori T, Kawaguchi Y, Kimura T. How does the ossification area of the posterior longitudinal ligament progress after cervical laminoplasty? Spine (Phila Pa 1976) 31:2807–2812. 2006;
10. Iwasaki M, Okuda S, Miyauchi A, Sakaura H, Mukai Y, Yonenobu K, et al. Surgical strategy for cervical myelopathy due to ossification of the posterior longitudinal ligament: Part 2: Advantages of anterior decompression and fusion over laminoplasty. Spine (Phila Pa 1976) 32:654–660. 2007;
11. Izumi T, Hirano T, Watanabe K, Sano A, Ito T, Endo N. Three-dimensional evaluation of volume change in ossification of the posterior longitudinal ligament of the cervical spine using computed tomography. Eur Spine J 22:2569–2574. 2013;
12. Jayakumar PN, Kolluri VR, Vasudev MK, Srikanth SG. Ossification of the posterior longitudinal ligament of the cervical spine in Asian Indians--a multiracial comparison. Clin Neurol Neurosurg 98:142–148. 1996;
13. Kawaguchi Y, Kanamori M, Ishihara H, Nakamura H, Sugimori K, Tsuji H, et al. Progression of ossification of the posterior longitudinal ligament following en bloc cervical laminoplasty. J Bone Joint Surg Am 83-A:1798–1802. 2001;
14. Lee T, Chacha PB, Khoo J. Ossification of posterior longitudinal ligament of the cervical spine in non-Japanese Asians. Surg Neurol 35:40–44. 1991;
15. Matsunaga S, Kukita M, Hayashi K, Shinkura R, Koriyama C, Sakou T, et al. Pathogenesis of myelopathy in patients with ossification of the posterior longitudinal ligament. J Neurosurg 96(2 Suppl):168–172. 2002;
16. Matsunaga S, Nakamura K, Seichi A, Yokoyama T, Toh S, Ichimura S, et al. Radiographic predictors for the development of myelopathy in patients with ossification of the posterior longitudinal ligament: a multicenter cohort study. Spine (Phila Pa 1976) 33:2648–2650. 2008;
17. Morio Y, Nagashima H, Teshima R, Nawata K. Radiological pathogenesis of cervical myelopathy in 60 consecutive patients with cervical ossification of the posterior longitudinal ligament. Spinal Cord 37:853–857. 1999;
18. Nagata K, Sato K. Diagnostic imaging of cervical ossification of the posterior longitudinal ligament. In : Yonenobu K, Nakamura K, Toyama Y, eds. OPLL: Ossification of the Posterior Longitudinal Ligament ed 2th ed. Tokyo: Springer; 2006. p. 127–143.
19. Ogata N, Kawaguchi H. Ossitication of the posterior longitudinal ligament of spine (OPLL). Clin Calcium 14:42–48. 2004;
20. Seichi A. Updates on ossification of posterior longitudinal ligament. Image diagnosis of ossification of posterior longitudinal ligament and associated diseases. Clin Calcium 19:1426–1434. 2009;
21. Takatsu T, Ishida Y, Suzuki K, Inoue H. Radiological study of cervical ossification of the posterior longitudinal ligament. J Spinal Disord 12:271–273. 1999;
22. Tsuyama N. Ossification of the posterior longitudinal ligament of the spine. Clin Orthop Relat Res (184):71–84. 1984;
23. Yonenobu K, Tsuzuki N, Nagata K, Toyama Y, Kato Y, Iwasaki M. Computer-assisted measurement of ossified lesion in ossification of the posterior longitudinal ligament of the cervical spine. Bone 16:283–286. 2002;
24. Zhong ZM, Chen JT. Phenotypic characterization of ligamentum flavum cells from patients with ossification of ligamentum flavum. Yonsei Med J 50:375–379. 2009;

Article information Continued

Fig. 1

Steps for 3D reconstruction and measurement. A: CT scanning with 3 mm thickness. B: The obtained images are transmitted digitally via PACS. C: 3D software imports CT data in the DICOM format. D: The segmented object (visualized by a colored-mask) is obtained by including only those pixels of the image with a value higher than or equal to the threshold value (thresholding). E: The segmented object is further split by making another mask on OPLL and removing floating pixels (region growing). F: 3D representation is created and the volume is calculated automatically by the software. CT: computed tomography, DICOM: Digital Imaging and Communications in Medicine, OPLL: ossification of the posterior longitudinal ligament.

Fig. 2

Reconstructed 3D images of OPLL after the calculation of the 3D. A: Anterior view. B: Lateral view. C: Posterior view. D: A volumetric mesh with the material assignment for finite element analysis. OPLL: ossification of the posterior longitudinal ligament.

Fig. 3

Measurements. A: The width was defined as the maximum distance along the axis of the coronal plane of OPLL and the thickness was defined as the maximum distance along the axis of the sagittal plane of OPLL. B: The length was defined as the longest dimension of OPLL. C: The volume was calculated automatically by the software in the 3D-reconstructed image. OPLL: ossification of the posterior longitudinal ligament.

Fig. 4

Representative images obtained in a 47-year-old woman with cervical OPLL. A: Lateral X-ray showing an ossified lesion at the level of C3–4, however, the cervico-thoracic junction is not visualized adequately due to shoulder shadow. B: Sagittal-reconstruction CT showing mixed type OPLL. C: Axial CT image at the level of C4 demonstrating central OPLL. OPLL: ossification of the posterior longitudinal ligament, CT: computed tomography.

Table 1

Inter-examiner reliability of each measurement

Data set ICC 95% CI
Volume 0.996 0.987–0.999
Thickness 0.973 0.907–0.994
Width 0.969 0.895–0.993
Length 0.995 0.983–0.999

ICC: intraclass correlation coefficient, CI: confidence interval

Table 2

Intra-examiner reliability of each measurement

Examiner Volume Thickness Width Length
SHC 0.996 0.998 0.947 0.995
CPG 0.995 0.944 0.873 0.985
PSH 0.991 0.994 0.913 0.978
Range 0.991–0.996 0.944–0.998 0.873–0.947 0.985–0.995