Reliable fusion of knee bone laser scans to establish ground truth for cartilage thickness measurement

We are interested in establishing ground truth data for validating morphology measurements of human knee cartilage from MR imaging. One promising approach is to compare the high-accuracy 3D laser scans of dissected cadaver knees before and after the dissolution of their cartilage. This requires an accurate and reliable method to fuse the individual laser scans from multiple views of the cadaver knees. Unfortunately existing methods using Iterative Closest Point (ICP) algorithm from off-the-shell packages often yield unreliable fusion results. We identify two major sources of variation: (i) the noise in depth measurements of the laser scans is significantly high and (ii) the use of point-to-point correspondence in ICP is not suitable due to sampling variation in the laser scans. We resolve the first problem by performing adaptive Gaussian smoothing on each individual laser scans prior to the fusion. For the second problem, we construct a surface mesh from the point cloud of each scan and adopt a point-to-mesh ICP scheme for pairwise alignment. The complete surface mesh is constructed by fusing all the scans in the order maximizing mutual overlaps. In experiments on 6 repeated scanning trials of a cadaver knee, our approach reduced the alignment error of point-to-point ICP by 30% and reduced coefficient of variation (CV) of cartilage thickness measurements from 5% down to 1.4%, significantly improving the method's repeatability.

[1]  V C Mow,et al.  Quantitation of articular surface topography and cartilage thickness in knee joints using stereophotogrammetry. , 1991, Journal of biomechanics.

[2]  Thomas H. Berquist,et al.  Imaging of articular pathology: MRI, CT, arthrography , 1997 .

[3]  André Guéziec 'Meshsweeper': Dynamic Point-to-Polygonal-Mesh Distance and Applications , 2001, IEEE Trans. Vis. Comput. Graph..

[4]  Felix Eckstein,et al.  Measuring cartilage morphology with quantitative magnetic resonance imaging. , 2004, Seminars in musculoskeletal radiology.

[5]  Yangsheng Wang,et al.  Mesh Smoothing via Adaptive Bilateral Filtering , 2005, International Conference on Computational Science.

[6]  Benjamin B. Kimia,et al.  Surface Reconstruction from Point Clouds by Transforming the Medial Scaffold , 2007, Sixth International Conference on 3-D Digital Imaging and Modeling (3DIM 2007).

[7]  Kuo-Young Cheng,et al.  A sharpness dependent filter for mesh smoothing , 2005, Comput. Aided Geom. Des..

[8]  Benjamin B. Kimia,et al.  Accurate Measurement of Cartilage Morphology Using a 3D Laser Scanner , 2006, CVAMIA.

[9]  Marc Levoy,et al.  A volumetric method for building complex models from range images , 1996, SIGGRAPH.

[10]  Marc Levoy,et al.  Efficient variants of the ICP algorithm , 2001, Proceedings Third International Conference on 3-D Digital Imaging and Modeling.

[11]  T P Andriacchi,et al.  Considerations in measuring cartilage thickness using MRI: factors influencing reproducibility and accuracy. , 2005, Osteoarthritis and cartilage.

[12]  Kari Pulli,et al.  Multiview registration for large data sets , 1999, Second International Conference on 3-D Digital Imaging and Modeling (Cat. No.PR00062).

[13]  P A Dieppe,et al.  No loss of cartilage volume over three years in patients with knee osteoarthritis as assessed by magnetic resonance imaging. , 2002, Osteoarthritis and cartilage.

[14]  Marc Levoy,et al.  Zippered polygon meshes from range images , 1994, SIGGRAPH.

[15]  Daniel Cohen-Or,et al.  Bilateral mesh denoising , 2003 .

[16]  Muhammad Hussain Fast Decimation of Polygonal Models , 2008, ISVC.

[17]  J. Reginster,et al.  The prevalence and burden of arthritis. , 2002, Rheumatology.