Evaluation of osseointegration using image analysis and visualization of 2D and 3D image data

Computerized image analysis, the discipline of using computers to automatically extract information from digital images, is a powerful tool for automating time consuming analysis tasks. In this thesis, image analysis and visualization methods are developed to facilitate the evaluation of osseointegration, i.e., the biological integration of a load-carrying implant in living bone. Adequate osseointegration is essential in patients who are in need of implant treatment. New implant types, with variations in bulk material and surface structural parameters, are continuously being developed. The main goal is to improve and speed up the osseointegration and thereby enhance patient well-being. The level of osseointegration can be evaluated by quantifying the bone tissue in proximity to the implant in e.g., light microscopy images of thin cross sections of bone implant samples extracted from humans or animals. This operator dependent quantitative analysis is cumbersome, time consuming and subjective. Furthermore, the thin sections represent only a small region of the whole sample. In this thesis work, computerized image analysis methods are developed to automate the quantification step. An image segmentation method is proposed for classifying the pixels of the images as bone tissue, non-bone tissue or implant. Subsequently, bone area and bone implant contact length in regions of interest are quantified. To achieve an accurate classification, the segmentation is based on both intensity and spatial information of the pixels. The automated method speeds up and facilitates the evaluation of osseointegration in the research laboratories. Another aim of this thesis is extending the 2D analysis to 3D and presenting methods for visualization of the 3D image volumes. To get a complete picture, information from the whole sample should be considered, rather than thin sections only. As a first step, 3D imaging of the implant samples is evaluated. 3D analysis methods, which follow the helix shaped implant thread and collects quantified features along the path, are presented. Additionally, methods for finding the position of the 2D section in the corresponding 3D image volume, i.e., image registration, are presented, enabling a direct comparison of the data from the two modalities. These novel and unique 3D quantification and visualization methods support the biomaterial researchers with improved tools for gaining a wider insight into the osseointegration process, with the ultimate goal of improved quality of life for the patients.

[1]  Kentaro Uesugi,et al.  Bone structure and mineralization demonstrated using synchrotron radiation computed tomography (SR-CT) in animal models: preliminary findings , 2003, Journal of Bone and Mineral Metabolism.

[2]  Avinash C. Kak,et al.  Principles of computerized tomographic imaging , 2001, Classics in applied mathematics.

[3]  N. Otsu A threshold selection method from gray level histograms , 1979 .

[4]  Jayaram K. Udupa,et al.  A framework for evaluating image segmentation algorithms , 2006, Comput. Medical Imaging Graph..

[5]  M Gavrilovic,et al.  Quantification of colocalization and cross‐talk based on spectral angles , 2009, Journal of microscopy.

[6]  Joakim Lindblad Surface Area Estimation of Digitized Planes Using Weighted Local Configurations , 2003, DGCI.

[7]  M. J. D. Powell,et al.  An efficient method for finding the minimum of a function of several variables without calculating derivatives , 1964, Comput. J..

[8]  Roger Lundqvist Atlas-Based Fusion of Medical Brain Images : Methods and Applications , 2001 .

[9]  Joakim Lindblad,et al.  Quantification of Bone Remodeling in SRmuCT Images of Implants , 2009, SCIA.

[10]  David P. Dobkin,et al.  The quickhull algorithm for convex hulls , 1996, TOMS.

[11]  P. Rüegsegger,et al.  A microtomographic system for the nondestructive evaluation of bone architecture , 2006, Calcified Tissue International.

[12]  David L Cochran,et al.  The evidence for immediate loading of implants. , 2006, The journal of evidence-based dental practice.

[13]  خالد بالطو,et al.  Quantification of periapical bone destruction in mice by microcomputed tomography. , 2000 .

[14]  A Robb Richard Biomedical Imaging: Past, Present and Predictions , 2006 .

[15]  Jayaram K. Udupa,et al.  Iterative relative fuzzy connectedness for multiple objects with multiple seeds , 2007, Comput. Vis. Image Underst..

[16]  Felix Beckmann,et al.  3D analysis of bone formation around titanium implants using micro-computed tomography (μCT) , 2006, SPIE Optics + Photonics.

[17]  Gunilla Borgefors,et al.  Hierarchical Chamfer Matching: A Parametric Edge Matching Algorithm , 1988, IEEE Trans. Pattern Anal. Mach. Intell..

[18]  D. Williams,et al.  The Williams Dictionary of Biomaterials: L , 1999 .

[19]  William Becker,et al.  A descriptive study on retrieved non-threaded and threaded implant designs. , 2005, Clinical oral implants research.

[20]  Françoise Peyrin,et al.  Synchrotron radiation /spl mu/CT of biological tissue , 2002, Proceedings IEEE International Symposium on Biomedical Imaging.

[21]  Lucia Ballerini,et al.  Comparison of histomorphometrical data obtained with two different image analysis methods , 2007, Journal of materials science. Materials in medicine.

[22]  Leo Joskowicz,et al.  Effective Intensity-Based 2D/3D Rigid Registration between Fluoroscopic X-Ray and CT , 2003, MICCAI.

[23]  Max A. Viergever,et al.  Mutual-information-based registration of medical images: a survey , 2003, IEEE Transactions on Medical Imaging.

[24]  Hiroshi Nakada,et al.  Micro-CT Analysis of Rabbit Cancellous Bone Around Implants , 2007 .

[25]  Börje Svensson,et al.  Single-dose irradiation followed by implant insertion in rat bone. An investigative study to find a critical level for osseointegration , 2010 .

[26]  Joakim Lindblad,et al.  High-Precision Boundary Length Estimation by Utilizing Gray-Level Information , 2009, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[27]  T Albrektsson,et al.  Quantitative and qualitative investigations of surface enlarged titanium and titanium alloy implants. , 1998, Clinical oral implants research.

[28]  J. Koplowitz,et al.  Design Of Perimeter Estimators For Digitized Planar Shapes , 1988, Other Conferences.

[29]  Calle Lejdfors High-Level GPU Programming: Domain-Specific Optimization and Inference , 2008 .

[30]  Robert C. Bolles,et al.  Parametric Correspondence and Chamfer Matching: Two New Techniques for Image Matching , 1977, IJCAI.

[31]  L. Grodzins,et al.  Optimum energies for x-ray transmission tomography of small samples. Applications of synchrotron radiation to computerized tomography I , 1983 .

[32]  G. S. Leventhal,et al.  Titanium, a metal for surgery. , 1951, The Journal of bone and joint surgery. American volume.

[33]  John A Jansen,et al.  The quantitative assessment of peri-implant bone responses using histomorphometry and micro-computed tomography. , 2009, Biomaterials.

[34]  D. Hill,et al.  Medical image registration , 2001, Physics in medicine and biology.

[35]  Bruno De Man,et al.  An outlook on x-ray CT research and development. , 2008, Medical physics.

[36]  Colin Studholme,et al.  An overlap invariant entropy measure of 3D medical image alignment , 1999, Pattern Recognit..

[37]  G. Harry van Lenthe,et al.  CT-based visualization and quantification of bone microstructure in vivo , 2008 .

[38]  Michael White,et al.  The pope and the heretic : a true story of Giordano Bruno, the man who dared to defy the Roman Inquisition , 2002 .

[39]  D Chappard,et al.  Synchrotron X-ray microtomography (on a micron scale) provides three-dimensional imaging representation of bone ingrowth in calcium phosphate biomaterials. , 2003, Biomaterials.

[40]  Milan Sonka,et al.  Image Processing, Analysis and Machine Vision , 1993, Springer US.

[41]  P Morberg,et al.  Cutting directions of bone with biomaterials in situ does influence the outcome of histomorphometrical quantifications. , 1995, Biomaterials.

[42]  Julia F. Barrett,et al.  Artifacts in CT: recognition and avoidance. , 2004, Radiographics : a review publication of the Radiological Society of North America, Inc.

[43]  Torsten Rohlfing,et al.  Fast intensity-based 2D-3D image registration of clinical data using light , 2003, Proceedings Ninth IEEE International Conference on Computer Vision.

[44]  Leonardo Bocchi,et al.  Image Segmentation by a Genetic Fuzzy c-Means Algorithm Using Color and Spatial Information , 2004, EvoWorkshops.

[45]  J. Chu,et al.  CT and PET lung image registration and fusion in radiotherapy treatment planning using the chamfer-matching method. , 1999, International journal of radiation oncology, biology, physics.

[46]  Carina B. Johansson,et al.  On tissue reactions to metal implants , 1991 .

[47]  P. Thurner,et al.  Comparison of microfocus- and synchrotron X-ray tomography for the analysis of osteointegration around Ti6Al4V implants. , 2004, European cells & materials.

[48]  S. Bradbury,et al.  Contrast Techniques in Light Microscopy , 1996, Microscopy Today.

[49]  J. Hickman,et al.  An investigation into the correlation between the electrical potentials of metals and their behaviour in biological fluids. , 1953, The Journal of bone and joint surgery. British volume.

[50]  Rodney A. Kennedy,et al.  A Survey of Medical Image Registration on Multicore and the GPU , 2010, IEEE Signal Processing Magazine.

[51]  Ann Wennerberg,et al.  The impact of oral implants - past and future, 1966-2042. , 2005, Journal.

[52]  J C Southam,et al.  Structural changes around screws used in the treatment of fractured human mandibles. , 1971, The British journal of oral surgery.

[53]  W. Eric L. Grimson,et al.  2D-3D rigid registration of X-ray fluoroscopy and CT images using mutual information and sparsely sampled histogram estimators , 2001, Proceedings of the 2001 IEEE Computer Society Conference on Computer Vision and Pattern Recognition. CVPR 2001.

[54]  F Peyrin,et al.  Bulk and interface investigations of scaffolds and tissue-engineered bones by X-ray microtomography and X-ray microdiffraction. , 2007, Biomaterials.

[55]  Daniel van Steenberghe,et al.  Histologic evaluation of Brånemark clinic oral implants retrieved from grafted sites. , 2006, Clinical implant dentistry and related research.

[56]  Tim A. J. Hopper,et al.  Quantitative MRI and Micro-CT of bone architecture : applications and limitations in orthopaedics , 2005 .

[57]  Goldberg,et al.  Genetic algorithms , 1993, Robust Control Systems with Genetic Algorithms.

[58]  Punam K. Saha,et al.  Perimeter and Area Estimations of Digitized Objects with Fuzzy Borders , 2003, DGCI.

[59]  H. Hansson,et al.  Osseointegrated titanium implants. Requirements for ensuring a long-lasting, direct bone-to-implant anchorage in man. , 1981, Acta orthopaedica Scandinavica.

[60]  Michel Adamek Creating a Graphical User InterfaceTemplate for Izolde : The complete design process, focusing on usability and design , 2010 .

[61]  TOR Hildebrand,et al.  Quantification of Bone Microarchitecture with the Structure Model Index. , 1997, Computer methods in biomechanics and biomedical engineering.

[62]  C. Johansson,et al.  Resonance frequency measurements in vivo and related surface properties of magnesium-incorporated, micropatterned and magnesium-incorporated TiUnite, Osseotite, SLA and TiOblast implants. , 2009, Clinical oral implants research.

[63]  John A. Jansen,et al.  ErratumErratum to “Osteoconductive modifications of Ti-implants in a goat defect model: Characterization of bone growth with SR μCT and histology”: [Biomaterials 26 (2005) 3009–3019] , 2006 .

[64]  R. Fisher THE USE OF MULTIPLE MEASUREMENTS IN TAXONOMIC PROBLEMS , 1936 .

[65]  Robert G. Paden,et al.  Metal artifact reduction image reconstruction algorithm for CT of implanted metal orthopedic devices: a work in progress , 2009, Skeletal Radiology.

[66]  Supun Samarasekera,et al.  Fuzzy Connectedness and Object Definition: Theory, Algorithms, and Applications in Image Segmentation , 1996, CVGIP Graph. Model. Image Process..

[67]  A. Eriksson,et al.  Heat-induced bone tissue injury : an in vivo investigation of heat tolerance of bone tissue and temperature rise in the drilling of cortical bone , 1984 .

[68]  Victoria Franke Stenport,et al.  Evaluations of bone tissue integration to pure and alloyed titanium implants. , 2008, Clinical implant dentistry and related research.

[69]  T. Martin,et al.  Recent developments in X-ray imaging with micrometer spatial resolution. , 2006, Journal of synchrotron radiation.

[70]  Fumihiko Ino,et al.  A GPGPU Approach for Accelerating 2-D/3-D Rigid Registration of Medical Images , 2006, ISPA.

[71]  Rafael C. González,et al.  Digital image processing, 3rd Edition , 2008 .

[72]  Jean-Louis Coatrieux,et al.  A moment-based three-dimensional edge operator , 1993 .

[73]  P Morberg,et al.  Importance of ground section thickness for reliable histomorphometrical results. , 1995, Biomaterials.

[74]  Frank Deinzer,et al.  Extended Global Optimization Strategy for Rigid 2D/3D Image Registration , 2007, CAIP.