1. A method performed by a nuclear magnetic resonance imaging (MRI) system, comprising: applying, at a scanner comprising a magnet, gradient coils and a radio frequency system, a gradient and a radio frequency pulse sequence to an imaging target to partition a three-dimensional (3D) Cartesian k-space into half projections in a ky-kz plane, the half projections divided into interleaves, each interleave comprising one or more of the half projections and providing an asymmetric radial sampling of the ky-kz plane; acquiring, at the scanner, magnetic resonance signals to sample phase and slice encodings along each of the half projections; and processing, at a data processing system in communication with the scanner, the acquired magnetic resonance signals to reconstruct at least one dynamic frame of the imaging target to provide 3D volumetric imaging of contrast dynamics.