Images from a canine study acquired with sequence B, showing the presence of artifacts and steps of the proposed method. Left: a systolic image (22nd frame in 55 total frames) of the study, showing optimal BOLD sensitivity and significant flow artifacts. Middle: The IR image as defined in the methods indicating that artifact information is retained (only pixels in air are shown as found by the level-set segmentation): the arrows point to artifacts. Right: The distribution (H) of the IR values in air: bracket indicates IR values from ghost artifact regions (see arrows in middle image) that act as outliers. Due to the presence of these outliers the excess kurtosis for this stack was y=QK=21.