An expectation-maximisation approach for simultaneous pixel classification and tracer kinetic modelling in dynamic contrast enhanced-magnetic resonance imaging

• Published in: Medical & biological engineering & computing Vol. 49; no. 4; pp. 485 - 495
• Main Authors: Sansone, Mario, Fusco, Roberta, Petrillo, Antonella, Petrillo, Mario, Bracale, Marcello
• Format: Journal Article
• Language: English
• Published: Berlin/Heidelberg Springer-Verlag 01.04.2011; Springer Nature B.V
• Subjects: Accuracy; Algorithms; Bioengineering; Biomedical and Life Sciences; Biomedical Engineering and Bioengineering; Biomedicine; Computer Applications; Contrast Media - pharmacokinetics; Human Physiology; Humans; Image Interpretation, Computer-Assisted - methods; Image processing systems; Imaging; Magnetic resonance imaging; Magnetic Resonance Imaging - methods; Missing data; Models, Biological; Neoplasms - diagnosis; NMR; Nuclear magnetic resonance; Original Article; Permeability; Plasma; Radiology; Simulation; Studies; Tumors; Expectation-maximisation; DCE-MRI; Tracer kinetic modelling; Classification
• ISSN: 0140-0118; 1741-0444; 1741-0444
• DOI: 10.1007/s11517-010-0695-x

Traditionally, tracer kinetic modelling and pixel classification of DCE-MRI studies are accomplished separately, although they could greatly benefit from each other. In this article, we propose an expectation-maximisation scheme for simultaneous pixel classification and compartmental modelling of DCE-MRI studies. The key point in the proposed scheme is the estimation of the kinetic parameters ( K trans and K ep ) of the two-compartmental model. Typically, they are estimated via nonlinear least-squares fitting. In our scheme, by exploiting the iterative nature of the EM algorithm, we use instead a Taylor expansion of the modelling equation. We developed the theoretical framework for the particular case of two classes and evaluated the performances of the algorithm by means of simulations. Results indicate that the accuracy of the proposed method supersedes the traditional pixel-by-pixel scheme and approaches the theoretical lower bound imposed by the Cramer–Rao theorem. Preliminary results on real data were also reported.