A pulse sequence for time-of-flight (TOF) magnetic resonance angiography (MRA) includes a fatsat segment, a magnetization transfer segment, and a spatial saturation segment that are applied by an MR apparatus to acquire MR data for image reconstruction with improved image quality. The pulse sequence is constructed such that at the beginning of each iteration of the inner loop of a 3D acquisition, a fatsat pulse is applied. After the fatsat pulse, MR data is acquired in a series of imaging segments with well-suppressed fat signal. Effective fat suppression is achieved by sampling central k-space data first, before signal from fat relaxes back to a pre-saturation level. Each imaging segment is immediately preceded by one of a MT pulse or a spatial saturation pulse and immediately followed by the other one of the MT pulse or the spatial saturation pulse.