A myocardial blood flow analysis scan includes incorporating a pharmacological kinetic model with the standard factor analysis model where each time activity curve is assumed to be a linear combination of factor curves. Pharmacological kinetics based factor analysis of dynamic structures (K-FADS-II) model can be applied, whereby estimating factor curves in the myocardium can be physiologically meaningful is provided. Additional optional aspects include performing a discretization to transform continuous-time K-FADS-II model into a discrete-time K-FADS-II model and application of an iterative Improved Voxel-Resolution myocardial blood flow (IV-MBF) algorithm. Where the model is applied without assumption that a right ventricle tissue curve and a left ventricle tissue curve obey a particular mathematical relationship, a least squares objective function can be applied to obtain estimates for parameters of the pharmacological kinetic model by applying a majorize-minimize optimization technique to iteratively estimate the curves.