A New Level-Set-Based Protocol for Accurate Bone Segmentation From CT Imaging

A new medical image segmentation pipeline for accurate bone segmentation from computed tomography (CT) imaging is proposed in this paper. It is a two-step methodology, with a pre-segmentation step and a segmentation refinement step, as follows. First, the user performs a rough segmenting of the desired region of interest. Second, a fully automatic refinement step is applied to the pre-segmented data. The automatic segmentation refinement is composed of several sub-steps, namely, image deconvolution, image cropping, and interpolation. The user-defined pre-segmentation is then refined over the deconvolved, cropped, and up-sampled version of the image. The performance of the proposed algorithm is exemplified with the segmentation of CT images of a composite femur bone, reconstructed with different reconstruction protocols. Segmentation outcomes are validated against a gold standard model, obtained using the coordinate measuring machine Nikon Metris LK V20 with a digital line scanner LC60-D and a resolution of $28~\mu \text{m}$ . High sub-pixel accuracy models are obtained for all tested data sets, with a maximum average deviation of 0.178 mm from the gold standard. The algorithm is able to produce high quality segmentation of the composite femur regardless of the surface meshing strategy used.

[1]  Tony F. Chan,et al.  A Multiphase Level Set Framework for Image Segmentation Using the Mumford and Shah Model , 2002, International Journal of Computer Vision.

[2]  M W Vannier,et al.  Validation study of skull three-dimensional computerized tomography measurements. , 1990, American journal of physical anthropology.

[3]  T. Albrektsson,et al.  Implant fixation improved by close fit. Cylindrical implant-bone interface studied in rabbits. , 1988, Acta orthopaedica Scandinavica.

[4]  E L Nickoloff,et al.  A simplified approach for modulation transfer function determinations in computed tomography. , 1985, Medical physics.

[5]  O. Harrysson,et al.  Custom-designed orthopedic implants evaluated using finite element analysis of patient-specific computed tomography data: femoral-component case study , 2007, BMC musculoskeletal disorders.

[6]  Joanne B. Adams,et al.  Patient-specific approach in total knee arthroplasty. , 2008, Orthopedics.

[7]  David E. Breen,et al.  Contour-based surface reconstruction using MPU implicit models , 2007, Graph. Model..

[8]  William J Hozack,et al.  Custom-fit total knee arthroplasty (OtisKnee) results in malalignment. , 2008, The Journal of arthroplasty.

[9]  W. Kalender,et al.  Accuracy limits for the determination of cortical width and density: the influence of object size and CT imaging parameters. , 1999, Physics in medicine and biology.

[10]  J Lund-Larsen,et al.  Determination of Hounsfield value for CT-based design of custom femoral stems. , 1999, The Journal of bone and joint surgery. British volume.

[11]  Andrew H. Gee,et al.  High resolution cortical bone thickness measurement from clinical CT data , 2010, Medical Image Anal..

[12]  Bernhard Preim,et al.  Visualization in Medicine: Theory, Algorithms, and Applications , 2007 .

[13]  Hideki Yoshikawa,et al.  Limits on the accuracy of 3-D thickness measurement in magnetic resonance images-effects of voxel anisotropy , 2003, IEEE Transactions on Medical Imaging.

[14]  Asmaa Maloul,et al.  The Impact of Voxel Size-Based Inaccuracies on the Mechanical Behavior of Thin Bone Structures , 2011, Annals of Biomedical Engineering.

[15]  D. Okada,et al.  Digital Image Processing for Medical Applications , 2009 .

[16]  Lee W Goldman,et al.  Principles of CT: Radiation Dose and Image Quality* , 2007, Journal of Nuclear Medicine Technology.

[17]  Elena De Momi,et al.  Does computer-assisted surgery benefit leg length restoration in total hip replacement? Navigation versus conventional freehand , 2009, International Orthopaedics.

[18]  David A. Steinman,et al.  An image-based modeling framework for patient-specific computational hemodynamics , 2008, Medical & Biological Engineering & Computing.

[19]  J. Geleijns,et al.  Evaluation of 4 Multisection CT Systems in Postoperative Imaging of a Cochlear Implant: A Human Cadaver and Phantom Study , 2008, American Journal of Neuroradiology.

[20]  Toby P. Breckon,et al.  Fundamentals of Digital Image Processing: A Practical Approach with Examples in Matlab , 2011 .

[21]  Toru Matsumoto,et al.  Imaging of small spherical structures in CT: simulation study using measured point spread function , 2007, Medical & Biological Engineering & Computing.

[22]  T. Hangartner,et al.  Evaluation of cortical bone by computed tomography , 1996, Journal of bone and mineral research : the official journal of the American Society for Bone and Mineral Research.

[23]  R L Wixson,et al.  The rationale, design characteristics, and preliminary results of a primary custom total hip prosthesis. , 1989, Clinical orthopaedics and related research.

[24]  J. Parvizi,et al.  Biology of implant osseointegration. , 2009, Journal of musculoskeletal & neuronal interactions.

[25]  Yan Kang,et al.  A new accurate and precise 3-D segmentation method for skeletal structures in volumetric CT data , 2003, IEEE Transactions on Medical Imaging.

[26]  Richard Szeliski,et al.  PSF estimation using sharp edge prediction , 2008, 2008 IEEE Conference on Computer Vision and Pattern Recognition.

[27]  P S Walker,et al.  Strains and micromotions of press-fit femoral stem prostheses. , 1987, Journal of biomechanics.

[28]  Beat Schmutz,et al.  Effects of CT image segmentation methods on the accuracy of long bone 3D reconstructions. , 2011, Medical engineering & physics.

[29]  M. Unser,et al.  Interpolation revisited [medical images application] , 2000, IEEE Transactions on Medical Imaging.

[30]  David L. Donoho,et al.  De-noising by soft-thresholding , 1995, IEEE Trans. Inf. Theory.

[31]  Ming Ye,et al.  Precision of cortical bone reconstruction based on 3D CT scans , 2009, Comput. Medical Imaging Graph..

[32]  Eric Keppel,et al.  Approximating Complex Surfaces by Triangulation of Contour Lines , 1975, IBM J. Res. Dev..

[33]  K. Egiazarian,et al.  Blind image deconvolution , 2007 .

[34]  J. Weiss,et al.  Subject-specific finite element model of the pelvis: development, validation and sensitivity studies. , 2005, Journal of biomechanical engineering.

[35]  Francisco Martins,et al.  TOWARDS A MEDICAL IMAGE ENHANCEMENT PROCEDURE BASED ON WAVELET TRANSFORM ANALYSIS , 2012 .

[36]  Tony F. Chan,et al.  Active contours without edges , 2001, IEEE Trans. Image Process..

[37]  I A Cunningham,et al.  Signal and noise in modulation transfer function determinations using the slit, wire, and edge techniques. , 1992, Medical physics.

[38]  M Besser,et al.  The utility of histological examination of tissue removed during elective joint replacement. A preliminary assessment. , 1999, The Journal of bone and joint surgery. American volume.

[39]  Kazunori Yoshiura,et al.  The accuracy of 3-dimensional magnetic resonance 3D vibe images of the mandible: an in vitro comparison of magnetic resonance imaging and computed tomography. , 2007, Oral surgery, oral medicine, oral pathology, oral radiology, and endodontics.

[40]  Steven W. Smith,et al.  The Scientist and Engineer's Guide to Digital Signal Processing , 1997 .

[41]  James V. Miller,et al.  Bias in the Localization of Curved Edges , 2004, ECCV.

[42]  J. C. Dainty,et al.  Iterative blind deconvolution method and its applications , 1988 .

[43]  Russell H. Taylor,et al.  Development of a Surgical Robot for Cementless Total Hip Arthroplasty , 1992, Clinical orthopaedics and related research.

[44]  P S Walker,et al.  Design and fabrication of cementless hip stems. , 1988, Clinical orthopaedics and related research.

[45]  P S Walker,et al.  Design of intramedullary femoral stems using computer graphics. , 1985, Engineering in medicine.

[46]  F. Brossa,et al.  The relevance of hydroxyapatite and spongious titanium coatings in fixation of cementless stems , 1998, Archives of Orthopaedic and Trauma Surgery.

[47]  William H. Richardson,et al.  Bayesian-Based Iterative Method of Image Restoration , 1972 .

[48]  Aldo von Wangenheim,et al.  Influence of Tomographic Slice Thickness and Field of View Variation on the Reproduction of Thin Bone Structures for Rapid Prototyping Purposes —An in Vitro Study , 2013 .

[49]  Stuart B Goodman,et al.  A conical-collared intramedullary stem can improve stress transfer and limit micromotion. , 2004, Clinical biomechanics.

[50]  B Notarberardino,et al.  An efficient approach to converting three-dimensional image data into highly accurate computational models , 2008, Philosophical Transactions of the Royal Society A: Mathematical, Physical and Engineering Sciences.

[51]  Purang Abolmaesumi,et al.  An estimation/correction algorithm for detecting bone edges in CT images , 2005, IEEE Transactions on Medical Imaging.

[52]  Jitendra Malik,et al.  Scale-Space and Edge Detection Using Anisotropic Diffusion , 1990, IEEE Trans. Pattern Anal. Mach. Intell..

[53]  J. L. Alves,et al.  3D reconstruction of a spinal motion segment from 2D medical images: Objective quantification of the geometric accuracy of the FE mesh generation procedure , 2013, 2013 IEEE 3rd Portuguese Meeting in Bioengineering (ENBENG).

[54]  I. Johnstone,et al.  Adapting to Unknown Smoothness via Wavelet Shrinkage , 1995 .

[55]  Y. Goussard,et al.  Experimental determination of CT point spread function anisotropy and shift-variance , 1997, Proceedings of the 19th Annual International Conference of the IEEE Engineering in Medicine and Biology Society. 'Magnificent Milestones and Emerging Opportunities in Medical Engineering' (Cat. No.97CH36136).

[56]  Jean-Luc Starck,et al.  Deconvolution and Blind Deconvolution in Astronomy , 2007 .

[57]  Roger L. Easton,et al.  Fourier Methods in Imaging: Easton/Fourier Methods in Imaging , 2010 .

[58]  H. Yoshikawa,et al.  Accuracy analysis of three‐dimensional bone surface models of the forearm constructed from multidetector computed tomography data , 2009, The international journal of medical robotics + computer assisted surgery : MRCAS.

[59]  Anna Vilanova,et al.  Correction for the dislocation of curved surfaces caused by the PSF in 2D and 3D CT images , 2005, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[60]  Rik Huiskes,et al.  The accuracy of criteria for automatic 3-D graphics reconstruction of bone from computer tomography , 1987 .

[61]  J Vander Sloten,et al.  Accuracy assessment of CT-based outer surface femur meshes , 2008, Computer aided surgery : official journal of the International Society for Computer Aided Surgery.

[62]  E. Samei,et al.  A method for measuring the presampled MTF of digital radiographic systems using an edge test device. , 1998, Medical physics.

[63]  R. Jacobs,et al.  A comparative evaluation of Cone Beam Computed Tomography (CBCT) and Multi-Slice CT (MSCT). Part II: On 3D model accuracy. , 2010, European journal of radiology.

[64]  W. Bargar,et al.  Shape the implant to the patient. A rationale for the use of custom-fit cementless total hip implants. , 1989, Clinical orthopaedics and related research.

[65]  Roger L. Easton,et al.  Fourier Methods in Imaging , 2011, J. Electronic Imaging.

[66]  Issei Mori,et al.  Deriving the modulation transfer function of CT from extremely noisy edge profiles , 2009, Radiological physics and technology.

[67]  Normand Robert,et al.  Generalized method for computation of true thickness and x-ray intensity information in highly blurred sub-millimeter bone features in clinical CT images , 2012, Physics in medicine and biology.

[68]  A Fenster,et al.  A method for modulation transfer function determination from edge profiles with correction for finite-element differentiation. , 1987, Medical physics.

[69]  Lucas J. van Vliet,et al.  On the Location Error of Curved Edges in Low-Pass Filtered 2-D and 3-D Images , 1994, IEEE Trans. Pattern Anal. Mach. Intell..