Hole-Filling by Rank Sparsity Tensor Decomposition for Medical Imaging

Surface integrity of 3D medical data is crucial for surgery simulation or virtual diagnoses. However, undesirable holes often exist due to external damage on bodies or accessibility limitation on scanners. To bridge the gap, hole-filling for medical imaging is a popular research topic in recent years. Considering that medical image, e.g. CT or MRI, has the natural form of tensor, we recognize the problem of medical hole-filling as the extension of PCP problem from matrix case to tensor case. Since the new problem in tensor case is much more difficult than the matrix case, we design an efficient algorithm for the extension by relaxation technique. The most significant feature of our algorithm is that unlike traditional methods which follow a strictly local approach, our method fixes the hole by the global structure in the specific medical data. Another important difference to previous algorithm is that our algorithm is able to automatically separate the completed data from the hole in an implicit manner. Our experiments demonstrate that the proposed method can lead to satisfactory result.

[1]  Yi Ma,et al.  Robust principal component analysis? , 2009, JACM.

[2]  J Vander Sloten,et al.  Design and manufacturing of cranioplasty implants by 3-axis cnc milling. , 2002, Technology and health care : official journal of the European Society for Engineering and Medicine.

[3]  C. Gold,et al.  Delete and insert operations in Voronoi/Delaunay methods and applications , 2003 .

[4]  Xiaodong Li,et al.  Stable Principal Component Pursuit , 2010, 2010 IEEE International Symposium on Information Theory.

[5]  John Wright,et al.  Robust Principal Component Analysis: Exact Recovery of Corrupted Low-Rank Matrices via Convex Optimization , 2009, NIPS.

[6]  Kecheng Liu,et al.  Shape recovery algorithms using level sets in 2-D/3-D medical imagery: a state-of-the-art review , 2002, IEEE Transactions on Information Technology in Biomedicine.

[7]  James F. O'Brien,et al.  Spectral surface reconstruction from noisy point clouds , 2004, SGP '04.

[8]  Emmanuel J. Candès,et al.  A Singular Value Thresholding Algorithm for Matrix Completion , 2008, SIAM J. Optim..

[9]  Brent J Liu,et al.  Utilizing data grid architecture for the backup and recovery of clinical image data. , 2005, Computerized medical imaging and graphics : the official journal of the Computerized Medical Imaging Society.

[10]  P. Tseng Convergence of a Block Coordinate Descent Method for Nondifferentiable Minimization , 2001 .

[11]  Joos Vandewalle,et al.  A Multilinear Singular Value Decomposition , 2000, SIAM J. Matrix Anal. Appl..

[12]  Xiangjian He,et al.  Tensor error correction for corrupted values in visual data , 2010, 2010 IEEE International Conference on Image Processing.

[13]  Mark A. Ganter,et al.  Real-time finite element modeling for surgery simulation: an application to virtual suturing , 2004, IEEE Transactions on Visualization and Computer Graphics.

[14]  Demetri Terzopoulos,et al.  Multilinear Analysis of Image Ensembles: TensorFaces , 2002, ECCV.

[15]  Jieping Ye,et al.  Tensor Completion for Estimating Missing Values in Visual Data , 2009, IEEE Transactions on Pattern Analysis and Machine Intelligence.

[16]  Martin Klein,et al.  Individual CAD/CAM Fabricated Glass-Bioceramic Implants in Reconstructive Surgery of the Bony Orbital Floor , 2006, Plastic and reconstructive surgery.

[17]  Yin Zhang,et al.  Fixed-Point Continuation for l1-Minimization: Methodology and Convergence , 2008, SIAM J. Optim..

[18]  Robert D. Nowak,et al.  Platelets: a multiscale approach for recovering edges and surfaces in photon-limited medical imaging , 2003, IEEE Transactions on Medical Imaging.