Diseased Region Detection of Longitudinal Knee MRI Data

Statistical analysis of longitudinal cartilage changes in osteoarthritis (OA) is of great importance and still a challenge in knee MRI data analysis. A major challenge is to establish a reliable correspondence across subjects within the same latent subpopulations. We develop a novel Gaussian hidden Markov model (GHMM) to establish spatial correspondence of cartilage thinning across both time and subjects within the same latent subpopulations and make statistical inference on the detection of diseased regions in each OA patient. A hidden Markov random field (HMRF) is proposed to extract such latent subpopulation structure. The EM algorithm and pseudo-likelihood method are both considered in making statistical inference. The proposed model can effectively detect diseased regions and present a localized analysis of longitudinal cartilage thickness within each latent subpopulation. Simulation studies and diseased region detection on 2D thickness maps extracted from full 3D longitudinal knee MRI Data for Pfizer Longitudinal Dataset are performed, which show that our proposed model outperforms standard voxel-based analysis.

[1]  Paul F. Whelan,et al.  Segmentation of the Left Ventricle of the Heart in 3-D+t MRI Data Using an Optimized Nonrigid Temporal Model , 2008, IEEE Transactions on Medical Imaging.

[2]  Stephen M. Smith,et al.  Segmentation of brain MR images through a hidden Markov random field model and the expectation-maximization algorithm , 2001, IEEE Transactions on Medical Imaging.

[3]  T. Spector,et al.  Osteoarthritis: New Insights. Part 1: The Disease and Its Risk Factors , 2000, Annals of Internal Medicine.

[4]  M. Cugmas,et al.  On comparing partitions , 2015 .

[5]  William M. Rand,et al.  Objective Criteria for the Evaluation of Clustering Methods , 1971 .

[6]  J. Besag.  Statistical Analysis of Non-Lattice Data , 1975 .

[7]  Stuart Geman,et al.  Markov Random Field Image Models and Their Applications to Computer Vision , 2010 .

[8]  F. Eckstein,et al.  Magnetic resonance imaging (MRI) of articular cartilage in knee osteoarthritis (OA): morphological assessment. , 2006, Osteoarthritis and cartilage.

[9]  Jh. Kellgren.  Radiological assessment of osteoarthritis , 1957 .

[10]  D. Rubin,et al.  Maximum likelihood from incomplete data via the EM - algorithm plus discussions on the paper , 1977 .

[11]  Felix Eckstein,et al.  A Technique for Regional Analysis of Femorotibial Cartilage Thickness Based on Quantitative Magnetic Resonance Imaging , 2008, IEEE Transactions on Medical Imaging.

[12]  Felix Eckstein,et al.  Does the use of ordered values of subregional change in cartilage thickness improve the detection of disease progression in longitudinal studies of osteoarthritis? , 2009, Arthritis and rheumatism.

[13]  W. Qian,et al.  Estimation of parameters in hidden Markov models , 1991, Philosophical Transactions of the Royal Society of London. Series A: Physical and Engineering Sciences.

[14]  A. Woolf,et al.  Burden of major musculoskeletal conditions. , 2003, Bulletin of the World Health Organization.

[15]  J. Kellgren,et al.  Radiological Assessment of Osteo-Arthrosis , 1957, Annals of the rheumatic diseases.

[16]  F. Cicuttini,et al.  Comparison of conventional standing knee radiographs and magnetic resonance imaging in assessing progression of tibiofemoral joint osteoarthritis. , 2005, Osteoarthritis and cartilage.

[17]  Jean-Pierre Raynauld,et al.  Quantitative magnetic resonance imaging of articular cartilage in knee osteoarthritis , 2003, Current opinion in rheumatology.

[18]  J. Besag.  On the Statistical Analysis of Dirty Pictures , 1986 .

[19]  Trevor Hastie,et al.  Risk factors for progressive cartilage loss in the knee: a longitudinal magnetic resonance imaging study in forty-three patients. , 2002, Arthritis and rheumatism.