The present invention relates to a method for estimating perfusion indices (PI) in a mammal (200), e.g. a human. Data (DAT, DAT′) representative of a contrast agent concentration as a function of time of an injected contrast agent is obtained from a medical imaging system (100). Perfusion indices (PI) are found by applying a parametric model (PM) for capillary transit time distributions as a function of time, and a minorize-maximization (MM) type procedure, such as an expectation-maximization (EM) type procedure, with regularization. The minorize-maximization type procedure has an exact analytical expression for the variance of an observation error of the contrast agent (Cε) in a non-linear observation model for the contrast agent concentration used in the maximization step. Clinical tests performed for 7 patients show improved MTT mapping as compared to singular value decomposition (SVD), and reduced sensitivity to delay.La présente invention concerne un procédé d'estimation d'indices de perfusion (PI) chez un mammifère (200), par exemple un humain. Des données (DAT, DAT') représentant une concentration en agent de contraste en fonction du temps d'un agent de contraste injecté sont obtenues dans un système d'imagerie médicale (100). Des indices de perfusion (PI) sont trouvés par application d'un modèle paramétrique (PM) de distributions du temps de transit capillaire en fonction du temps, et une procédure de type minorisation-maximisation (MM), comme une procédure de type attente-maximisation (EM), avec régularisation. La procédure de type minorisation-maximisation a une expression analytique exacte pour la variance d'une erreur d'observation de l'agent de contraste (Cε) dans un modèle d'observation non linéaire pour la concentration en agent de contraste utilisée dans l'étape de maximisation. Des tests cliniques réalisés pour 7 patients ont présenté un mappage MTT amélioré comparé à une décomposition en valeurs singulières (SVD), et une sensibilité réduite au reta
