Bayesian spatiotemporal segmentation of combined PET-CT data using a bivariate poisson mixture model

This paper presents an unsupervised algorithm for the joint segmentation of 4-D PET-CT images. The proposed method is based on a bivariate-Poisson mixture model to represent the bimodal data. A Bayesian framework is developed to label the voxels as well as jointly estimate the parameters of the mixture model. A generalized four-dimensional Potts-Markov Random Field (MRF) has been incorporated into the method to represent the spatio-temporal coherence of the mixture components. The method is successfully applied to 4-D registered PET-CT data of a patient with lung cancer. Results show that the proposed model fits accurately the data and allows the segmentation of different tissues and the identification of tumors in temporal series.

[1]  Thomas Beyer,et al.  Non-small cell lung cancer: dual-modality PET/CT in preoperative staging. , 2003, Radiology.

[2]  Jayaram K. Udupa,et al.  Co-segmentation of Functional and Anatomical Images , 2012, MICCAI.

[3]  Hadj Batatia,et al.  4D CT image reconstruction based on interpolated optical flow fields , 2010, 2010 IEEE International Conference on Image Processing.

[4]  M. Girolami,et al.  Riemann manifold Langevin and Hamiltonian Monte Carlo methods , 2011, Journal of the Royal Statistical Society: Series B (Statistical Methodology).

[5]  Anand Rangarajan,et al.  A Bayesian Joint Mixture Framework for the Integration of Anatomical Information in Functional Image Reconstruction , 2000, Journal of Mathematical Imaging and Vision.

[6]  Byung-Tae Kim,et al.  Non-small cell lung cancer: prospective comparison of integrated FDG PET/CT and CT alone for preoperative staging. , 2005, Radiology.

[7]  J. Buatti,et al.  Globally Optimal Tumor Segmentation in PET-CT Images: A Graph-Based Co-segmentation Method , 2011, IPMI.

[8]  Isabelle Bloch,et al.  Joint variational segmentation of CT-PET data for tumoral lesions , 2010, 2010 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[9]  O Caselles,et al.  The retrospective binning method improves the consistency of phase binning in respiratory-gated PET/CT , 2012, Physics in medicine and biology.

[10]  Sylvia Richardson,et al.  Markov chain concepts related to sampling algorithms , 1995 .

[11]  Nelly Pustelnik,et al.  A wavelet-based quadratic extension method for image deconvolution in the presence of poisson noise , 2009, 2009 IEEE International Conference on Acoustics, Speech and Signal Processing.

[12]  Hugh Gribben,et al.  MAP-MRF segmentation of lung tumours in PET/CT images , 2009, 2009 IEEE International Symposium on Biomedical Imaging: From Nano to Macro.

[13]  N. Dobigeon,et al.  Bayesian segmentation of chest tumors in pet scans using a Poisson-Gamma mixture model , 2011, 2011 IEEE Statistical Signal Processing Workshop (SSP).

[14]  David Dagan Feng,et al.  Lung tumor delineation in PET-CT images using a downhill region growing and a Gaussian mixture model , 2011, 2011 18th IEEE International Conference on Image Processing.

[15]  Junjie Bai,et al.  Globally optimal lung tumor co-segmentation of 4D CT and PET images , 2013, Medical Imaging.

[16]  Jeffrey A. Fessler,et al.  Statistical image reconstruction for polyenergetic X-ray computed tomography , 2002, IEEE Transactions on Medical Imaging.

[17]  Christian P. Robert,et al.  Monte Carlo Statistical Methods , 2005, Springer Texts in Statistics.