Sparse Representation of Deformable 3D Organs with Spherical Harmonics and Structured Dictionary

This paper proposed a novel algorithm to sparsely represent a deformable surface (SRDS) with low dimensionality based on spherical harmonic decomposition (SHD) and orthogonal subspace pursuit (OSP). The key idea in SRDS method is to identify the subspaces from a training data set in the transformed spherical harmonic domain and then cluster each deformation into the best-fit subspace for fast and accurate representation. This algorithm is also generalized into applications of organs with both interior and exterior surfaces. To test the feasibility, we first use the computer models to demonstrate that the proposed approach matches the accuracy of complex mathematical modeling techniques and then both ex vivo and in vivo experiments are conducted using 3D magnetic resonance imaging (MRI) scans for verification in practical settings. All results demonstrated that the proposed algorithm features sparse representation of deformable surfaces with low dimensionality and high accuracy. Specifically, the precision evaluated as maximum error distance between the reconstructed surface and the MRI ground truth is better than 3 mm in real MRI experiments.

[1]  Fillia Makedon,et al.  A Spatio-Temporal Modeling Method for Shape Representation , 2006, Third International Symposium on 3D Data Processing, Visualization, and Transmission (3DPVT'06).

[2]  Dan Wang,et al.  In vivo tracking of 3D organs using spherical harmonics and subspace clustering , 2009, 2009 16th IEEE International Conference on Image Processing (ICIP).

[3]  Thomas Vetter,et al.  Modeling the Remaining Flexibility of Partially Fixed Statistical Shape Models , 2008 .

[4]  M. McPeek,et al.  Modeling Three-Dimensional Morphological Structures Using Spherical Harmonics , 2009, Evolution; international journal of organic evolution.

[5]  Sheng Chen,et al.  Orthogonal least squares methods and their application to non-linear system identification , 1989 .

[6]  Dmitry B. Goldgof,et al.  The Use of Three- and Four-Dimensional Surface Harmonics for Rigid and Nonrigid Shape Recovery and Representation , 1995, IEEE Trans. Pattern Anal. Mach. Intell..

[7]  Guido Gerig,et al.  Parametrization of Closed Surfaces for 3-D Shape Description , 1995, Comput. Vis. Image Underst..

[8]  Guillermo Sapiro,et al.  Sparse Representations for Three-Dimensional Range Data Restoration , 2011 .

[9]  Moo K. Chung,et al.  Tensor-Based Cortical Surface Morphometry via Weighted Spherical Harmonic Representation , 2008, IEEE Transactions on Medical Imaging.

[10]  Robin Sibson,et al.  What is projection pursuit , 1987 .

[11]  Moo K. Chung,et al.  Weighted Fourier Series Representation and Its Application to Quantifying the Amount of Gray Matter , 2007, IEEE Transactions on Medical Imaging.

[12]  Fillia Makedon,et al.  Spherical mapping for processing of 3D closed surfaces , 2006, Image Vis. Comput..

[13]  Shang-Liang Chen,et al.  Orthogonal least squares learning algorithm for radial basis function networks , 1991, IEEE Trans. Neural Networks.

[14]  Michel Couprie,et al.  An open, clinically-validated database of 3D+t cine-MR images of the left ventricle with associated manual and automated segmentation , 2007, The Insight Journal.

[15]  Joseph F. Murray,et al.  An improved FOCUSS-based learning algorithm for solving sparse linear inverse problems , 2001, Conference Record of Thirty-Fifth Asilomar Conference on Signals, Systems and Computers (Cat.No.01CH37256).

[16]  Kjersti Engan,et al.  Multi-frame compression: theory and design , 2000, Signal Process..

[17]  Stéphane Mallat,et al.  Matching pursuits with time-frequency dictionaries , 1993, IEEE Trans. Signal Process..

[18]  Dan Wang,et al.  Real time tracking of 3D organ surfaces using single MR image and limited optical viewing , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[19]  Horace Ho-Shing Ip,et al.  A multi-resolution statistical deformable model (MISTO) for soft-tissue organ reconstruction , 2009, Pattern Recognit..

[20]  Jean-Louis Dillenseger,et al.  Spherical Harmonics Based Intrasubject 3-D Kidney Modeling/Registration Technique Applied on Partial Information , 2006, IEEE Transactions on Biomedical Engineering.

[21]  Guoyan Zheng,et al.  A 2D/3D correspondence building method for reconstruction of a patient-specific 3D bone surface model using point distribution models and calibrated X-ray images , 2009, Medical Image Anal..

[22]  Fillia Makedon,et al.  Hemispherical Harmonic Surface Description and Applications to Medical Image Analysis , 2006, Third International Symposium on 3D Data Processing, Visualization, and Transmission (3DPVT'06).

[23]  Zhifeng Zhang,et al.  Adaptive time-frequency decompositions , 1994 .

[24]  Paul J. Besl,et al.  A Method for Registration of 3-D Shapes , 1992, IEEE Trans. Pattern Anal. Mach. Intell..

[25]  Dan Wang,et al.  Real time tracking of exterior and interior organ surfaces using sparse sampling of the exterior surfaces , 2010, 2010 IEEE International Conference on Acoustics, Speech and Signal Processing.

[26]  Bhaskar D. Rao,et al.  Sparse signal reconstruction from limited data using FOCUSS: a re-weighted minimum norm algorithm , 1997, IEEE Trans. Signal Process..

[27]  David J. Field,et al.  Sparse coding with an overcomplete basis set: A strategy employed by V1? , 1997, Vision Research.

[28]  Christopher C Thompson,et al.  Surgery Insight: natural orifice transluminal endoscopic surgery—an analysis of work to date , 2007, Nature Clinical Practice Gastroenterology &Hepatology.

[29]  Timothy F. Cootes,et al.  A minimum description length approach to statistical shape modeling , 2002, IEEE Transactions on Medical Imaging.

[30]  Guoyan Zheng,et al.  Statistical deformable bone models for robust 3D surface extrapolation from sparse data , 2007, Medical Image Anal..

[31]  A. Bruckstein,et al.  K-SVD : An Algorithm for Designing of Overcomplete Dictionaries for Sparse Representation , 2005 .

[32]  Long Quan,et al.  Accurate and Scalable Surface Representation and Reconstruction from Images , 2007, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[33]  Michael A. Saunders,et al.  Atomic Decomposition by Basis Pursuit , 1998, SIAM J. Sci. Comput..

[34]  M. Elad,et al.  $rm K$-SVD: An Algorithm for Designing Overcomplete Dictionaries for Sparse Representation , 2006, IEEE Transactions on Signal Processing.

[35]  Guoliang Xu,et al.  A general framework for surface modeling using geometric partial differential equations , 2008, Comput. Aided Geom. Des..

[36]  Mike E. Davies,et al.  Gradient Pursuits , 2008, IEEE Transactions on Signal Processing.

[37]  Ahmed H. Tewfik,et al.  Learning Sparse Representation Using Iterative Subspace Identification , 2010, IEEE Transactions on Signal Processing.