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) model can be applied, whereby a means for estimating factor curves in the myocardium that are physiologically meaningful is provided. Additional optional aspects include performing a discretization to transform continuous-time K-FADS model into a discretetime K-FADS model and application of an iterative Voxel-Resolution myocardial blood flow (V-MBF) algorithm. A V-MBF algorithm based on a model that accounts for the fact that the shape of TACs due to ischemic and normal tissue are different can be included.