Abstract
This study aimed at developing a fully automated bone segmentation method for the human knee (femur and tibia) from magnetic resonance (MR) images. MR imaging was acquired on a whole body 1.5T scanner with a gradient echo fat suppressed sequence using an extremity coil. The method was based on the Ray Casting technique which relies on the decomposition of the MR images into multiple surface layers to localize the boundaries of the bones and several partial segmentation objects being automatically merged to obtain the final complete segmentation of the bones. Validation analyses were performed on 161 MR images from knee osteoarthritis patients, comparing the developed fully automated to a validated semi-automated segmentation method, using the average surface distance (ASD), volume correlation coefficient, and Dice similarity coefficient (DSC). For both femur and tibia, respectively, data showed excellent bone surface ASD (0.50 ± 0.12 mm; 0.37 ± 0.09 mm), average oriented distance between bone surfaces within the cartilage domain (0.02 ± 0.07 mm; −0.05 ± 0.10 mm), and bone volume DSC (0.94 ± 0.05; 0.92 ± 0.07). This newly developed fully automated bone segmentation method will enable large scale studies to be conducted within shorter time durations, as well as increase stability in the reading of pathological bone.






Similar content being viewed by others
References
Berthiaume MJ, Raynauld JP, Martel-Pelletier J, Labonté F, Beaudoin G, Bloch DA, Choquette D, Haraoui B, Altman RD, Hochberg M, Meyer JM, Cline GA, Pelletier JP (2005) Meniscal tear and extrusion are strongly associated with the progression of knee osteoarthritis as assessed by quantitative magnetic resonance imaging. Ann Rheum Dis 64:556–563
Dalvi R, Abugharbieh R, Wilson DC, Wilson D (2007) Multi-contrast MR for enhanced bone imaging and segmentation. Conf Proc IEEE Eng Med Biol Soc 2007:5620–5623
Ding C, Cicuttini F, Jones G (2007) Tibial subchondral bone size and knee cartilage defects: relevance to knee osteoarthritis. Osteoarthr Cartil 15:479–486
Dogdas B, Shattuck D, Leahy RM (2002) Segmentation of the skull in 3D human MR images using mathematical morphology. Proceedings of the SPIE. Med Imaging 4684:1553–1562
Felson DT, Niu J, Guermazi A, Roemer F, Aliabadi P, Clancy M, Torner J, Lewis CE, Nevitt MC (2007) Correlation of the development of knee pain with enlarging bone marrow lesions on magnetic resonance imaging. Arthritis Rheum 56:2986–2992
Folkesson J, Dam E, Olsen OF, Pettersen P, Christiansen C (2005) Automatic segmentation of the articular cartilage in knee MRI using a hierarchical multi-class classification scheme. Med Image Comput Comput Assist Interv 8:327–334
Folkesson J, Olsen OF, Pettersen P, Dam E, Christiansen C (2005) Combining binary classifiers for automatic cartilage segmentation in knee MRI. In: Liu Y, Jiang T, Zhang C (eds) Lecture notes in computer science. Springer-Verlag, Berlin, pp 230–239
Fripp J, Bourgeat P, Crozier S, Ourselin S (2007) Segmentation of the bones in MRI of the knee using phase, magnitude and shape information. Acad Radiol 10:1201–1208
Fripp J, Crozier S, Warfield SK, Ourselin S (2007) Automatic segmentation of the bone and extraction of the bone-cartilage interface from magnetic resonance images of the knee. Phys Med Biol 52:1617–1631
Fripp J, Warfield SK, Crozier S and Ourselin S. (2007) Automatic segmentation of articular cartilage in magnetic resonance images of the knee. Lecture notes in computer science, medical image computing and computer-assisted intervention 10:186–194
Fripp J, Crozier S, Warfield SK, Ourselin S (2010) Automatic segmentation and quantitative analysis of the articular cartilages from magnetic resonance images of the knee. IEEE Trans Med Imaging 29:55–64
Gambini J, Mejail M, Jacobo J and Delrieux C. (2004) SAR image segmentation through B-spline deformable contours and fractal dimension. Proceedings of the international society for photogrammetry and remote sensing (ISPRS) July 15–23, Istanbul, Turkey
Gerig G, Jomier M and Chakos M. (2001) VALMET: A new validation tool for assessing and improving 3D object segmentation. Proceedings of the international society and conference series on medical image computing and computer-assisted intervention 2208:516–523
Grau V, Mewes AUJ, Alcaniz M, Kikinis R, Warfield SK (2004) Improved watershed transform for medical image segmentation using prior information. IEEE Trans Med Imaging 23:447–457
Hamarneh G, Li X (2009) Watershed segmentation using prior shape and appearance knowledge. Image Vis Comput 27:59–68
Heinze P, Meister D, Kober R, Raczkowsky J, Worn H (2002) Atlas-based segmentation of pathological knee joints. Stud Health Technol Inform 85:198–203
Kapur T, Beardsley PA, Gibson SF, Grimson WEL and Wells WM. (1998) Model based segmentation of clinical knee MRI. Proceedings of the IEEE international workshop on model-based 3D image analysis, Bombay, India:97–106
Kauffmann C, Gravel P, Godbout B, Gravel A, Beaudoin G, Raynauld J-P, Martel-Pelletier J, Pelletier J-P, DeGuise JA (2003) Computer-aided method for quantification of cartilage thickness and volume changes using MRI: validation study using a synthetic model. IEEE Trans Biomed Eng 50:978–988
Kubassova O, Boyle RD and Pyatnizkiy M. (2005) Bone segmentation in metacarpophalangeal MR data. Proceedings of the third international conference on advances in pattern recognition (ICAPR):726–735
Levoy M (1990) Efficient ray tracing of volume data. ACM Trans Gr 9:245–261
Loeuille D, Chary-Valckenaere I, Champigneulle J, Rat AC, Toussaint F, Pinzano-Watrin A, Goebel JC, Mainard D, Blum A, Pourel J, Netter P, Gillet P (2005) Macroscopic and microscopic features of synovial membrane inflammation in the osteoarthritic knee: correlating magnetic resonance imaging findings with disease severity. Arthritis Rheum 52:3492–3501
Lorenz C, von Berg J (2005) Fast automated object detection by recursive casting of search rays. Int Congr Ser 1281:230–235
Lorigo LM, Faugeras O, Grimson WEL, Keriven R and Kikinis R. (1998) Segmentation of bone in clinical knee MRI using texture-based geodesic active contours. Lecture notes in computer science, proceedings of the first international conference on medical image computing and computer-assisted intervention 1496:1195–1204
Mayunga SD, Coleman DJ, Zhang Y (2007) A semi-automated approach for extracting buildings from QuickBird imagery applied to informal settlement mapping. Int J Remote Sens 28:2343–2357
Otsu N (1979) A threshold selection method from gray-level histograms. IEEE Trans Syst Man Cybern 9:62–66
Pelletier JP, Raynauld JP, Abram F, Haraoui B, Choquette D, Martel-Pelletier J (2008) A new non-invasive method to assess synovitis severity in relation to symptoms and cartilage volume loss in knee osteoarthritis patients using MRI. Osteoarthr Cartil 16(Suppl 3):S8–S13
Raynauld JP, Kauffmann C, Beaudoin G, Berthiaume MJ, de Guise JA, Bloch DA, Camacho F, Godbout B, Altman RD, Hochberg M, Meyer JM, Cline G, Pelletier JP, Martel-Pelletier J (2003) Reliability of a quantification imaging system using magnetic resonance images to measure cartilage thickness and volume in human normal and osteoarthritic knees. Osteoarthr Cartil 11:351–360
Raynauld JP, Martel-Pelletier J, Berthiaume MJ, Labonté F, Beaudoin G, de Guise JA, Bloch DA, Choquette D, Haraoui B, Altman RD, Hochberg M, Meyer JM, Cline G, Pelletier JP (2004) Quantitative magnetic resonance imaging evaluation of knee osteoarthritis progression over two years and correlation with clinical symptoms and radiologic changes. Arthritis Rheum 50:476–487
Raynauld JP, Martel-Pelletier J, Bias P, Laufer S, Haraoui B, Choquette D, Beaulieu AD, Abram F, Dorais M, Vignon E, Pelletier JP (2009) Protective effects of licofelone, a 5-lipoxygenase and cyclo-oxygenase inhibitor, versus naproxen on cartilage loss in knee osteoarthritis: a first multicentre clinical trial using quantitative MRI. Ann Rheum Dis 68:938–947
Roemer FW, Frobell R, Hunter DJ, Crema MD, Fischer W, Bohndorf K, Guermazi A (2009) MRI-detected subchondral bone marrow signal alterations of the knee joint: terminology, imaging appearance, relevance and radiological differential diagnosis. Osteoarthr Cartil 17:1115–1131
Wolf M, Welerich P, Niemann H (1997) Automatic segmentation and 3D-registration of a femoral bone in MR images of the knee. Pattern Recognit Image Anal 7:152–165
Yuille AL, Poggio TA (1986) Scaling theorems for zero crossings. IEEE Trans Pattern Anal Mach Intell 8:15–25
Acknowledgments
The authors would like to thank Yvan Ross for his involvement in the execution of the computation needed for the validation protocol, Françoys Labonté, PhD and François Guéritaud, PhD for their critical review and comments, and Virginia Wallis for her assistance with the manuscript preparation.
Author information
Authors and Affiliations
Corresponding author
Appendix: Bone Localization Procedure
Appendix: Bone Localization Procedure
The Ray Casting technique requires some prior knowledge including an approximation of the bone location to set the observers in each image. This appendix describes the autonomous procedure.
The analysis of the intensity histogram H shows a mixture of two probability densities, one centered in the low values, and the other in the high values. Because of its stability and robustness, the Otsu’s algorithm [25] was chosen to analyze the histogram, giving a decomposition of the histogram as a sum of two Gaussian distributions N (μ d , σ d ) and \( N(\mu_{{_{b} }} ,\sigma_{{_{b} }} ) \). The histogram’s dark tissue peak is μ d and bright tissue peak is \( \mu_{{_{b} }} \). Values inferior to the dark peak are called very dark and values superior to the bright peak, very bright. Otsu’s algorithm also provides a decision threshold s of separation between dark and bright, and mixture parameters α d and α b , allowing the histogram to be written:
An algorithm was designed to capture IFemur and ITibia on the MR image where the bones, femur and tibia (dark), are surrounded by cartilage and muscle (very bright), further surrounded by fat (dark), then skin (bright), and finally the background (very dark).
One has to start from Ω and evaluate the following set in order of decreasing size, where Leg is the leg, and MCart is the cartilage and muscle set:
Let us define H convex (X) the convex hull of a finite set X. Considering the decision threshold s, the anatomical convexity of the leg in each axial slice, and the brightness of the tissues, each axial subset of the leg can be approximated by
In order to have an appropriate cartilage and muscle set, surrounding femur and tibia, the selection relies on bright intensity and on convexity properties. Thus, MCart can be written as:
Figure 7a, b shows two sagittal slices of the intermediate results, the Leg z=t set (Fig. 7a) and the MCart z=t set (Fig. 7b).
a, b A sagittal slice of the image restricted to the Leg z=t set (a) and restricted to the smaller MCart z=t set (b). The black parts of the slices denote the outside of each set while parts with image content denote the inside of each set. c A representative sagittal slice and d the sets Leg z=t , MCart z=t , IFemurandITibia, respectively, in dark, medium, light and bright gray
Finally, the bone part inside this set can be identified by the very dark tissue inside MCart such that:
Femur and tibia sets are easily defined by separating the two largest components of FemurTibia as shown in Fig. 7c and d, where the largest is the set Femur and the second largest is the set Tibia. To facilitate the decomposition, an opening morphology operator ∘ is used [12] with structuring disk D of small diameter, e.g. five pixels. The set FemurTibia admits the following decomposition:
Rights and permissions
About this article
Cite this article
Dodin, P., Martel-Pelletier, J., Pelletier, JP. et al. A fully automated human knee 3D MRI bone segmentation using the ray casting technique. Med Biol Eng Comput 49, 1413–1424 (2011). https://doi.org/10.1007/s11517-011-0838-8
Received:
Accepted:
Published:
Issue Date:
DOI: https://doi.org/10.1007/s11517-011-0838-8