Described here are systems and methods for producing an image that depicts blood flow stasis using magnetic resonance imaging (MRI), Doppler echocardiography, or other medical instruments for measuring flow velocities in a human body. A time series of three-dimensional (3D) image volumes is provided, where this time series of 3D image volumes contains flow velocity information at voxel locations in a 3D volume in a subject. One or more regions-of-interest are then segmented from the 3D image volumes. For each voxel in the regions-of-interest, velocity magnitudes are calculated. Using the velocity magnitudes, a flow stasis volume is produced by computing a relative stasis value for each voxel location in the corresponding region-of-interest. This flow stasis volume can be provided as a 3D flow stasis image, or a flow stasis map can be produced by projecting the flow stasis volume onto a two-dimensional (2D) plane.
