Streamlines from 4D-flow MRI have been used clinically for intracranial blood-flow tracking. However, deterministic and stochastic errors degrade streamline quality. The purpose of this study is to integrate displacement corrections, probabilistic streamlines, and novel fluid constraints to improve selective blood-flow tracking and emulate "virtual bolus injections." Both displacement artifacts (deterministic) and velocity noise (stochastic) inherently occur during phase-contrast MRI acquisitions. Here, two displacement correction methods, single-step and iterative, were tested in silico with simulated displacements and were compared with ground-truth velocity fields. Next, the effects of combining displacement corrections and constrained probabilistic streamlines were performed in 10 healthy volunteers using time-averaged 4D-flow data. Measures of streamline length and depth into vasculature were then compared with streamlines generated with no corrections and displacement correction alone using one-way repeated-measures analysis of variance and Friedman's tests. Finally, virtual injections with improved streamlines were generated for three intracranial pathology cases. Iterative displacement correction outperformed the single-step method in silico. In volunteers, the combination of displacement corrections and constrained probabilistic streamlines allowed for significant improvements in streamline length and increased the number of streamlines entering the circle of Willis relative to streamlines with no corrections and displacement correction alone. In the pathology cases, virtual injections with improved streamlines were qualitatively similar to dynamic arterial spin labeling images and allowed for forward/reverse selective flow tracking to characterize cerebrovascular malformations. Virtual injections with improved streamlines from 4D-flow MRI allow for flexible, robust, intracranial flow tracking.