A method for controlling a magnetic resonance imaging diagnostic apparatus for flow-sensitive black blood imaging of a region of interest of a subject using a gradient echo pulse sequence, the method comprising: generating a slice gradient magnetic field, a phase-encode gradient magnetic field, and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis, and a read-out axis, respectively; generating a gradient echo pulse sequence, which includes a flow-weighting dephase gradient-magnetic-field pulse in at least one axis selected from the slice axis, the phase-encode axis, and the read-out axis; and acquiring three-dimensional nuclear magnetic resonance black blood image data of the region of interest using the gradient echo pulse sequence, wherein the dephase gradient-magnetic field pulse corresponds to a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in the region of interest of the subject, wherein the method further comprises setting an echo time of the gradient echo pulse sequence in such a manner that a difference in susceptibility between flows in the arteries and veins and the surrounding tissue is weighted in the contrast of the black blood image data such that both the effect of flow and the effect of susceptibility are reflected in the black blood image data; creating an amplitude image and a phase image from the black blood image data; synthesizing the amplitude image and the phase image to create a phase mask image that represents a part greatly phase-shifted and a part slightly phase-shifted; performing a three-dimensional process by multiplying the phase mask image and the amplitude image, thereby creating a three-dimensional black blood image; and performing a minimum intensity projection of the three-dimensional black blood image thereby visualizing both the veins and the arteries clearly and accurately. generating a slice gradient magnetic field, a phase-encode gradient magnetic field, and a read-out gradient magnetic field that extend in a slice axis, a phase-encode axis, and a read-out axis, respectively; generating a gradient echo pulse sequence, which includes a flow-weighting dephase gradient-magnetic-field pulse in at least one axis selected from the slice axis, the phase-encode axis, and the read-out axis; and acquiring three-dimensional nuclear magnetic resonance black blood image data of the region of interest using the gradient echo pulse sequence, wherein the dephase gradient-magnetic field pulse corresponds to a dephase amount for weighting a signal-level decrease resulting from flows in the arteries and veins present in the region of interest of the subject, wherein the method further comprises setting an echo time of the gradient echo pulse sequence in such a manner that a difference in susceptibility between flows in the arteries and veins and the surrounding tissue is weighted in the contrast of the black blood image data such that both the effect of flow and the effect of susceptibility are reflected in the black blood image data; creating an amplitude image and a phase image from the black blood image data; synthesizing the amplitude image and the phase image to create a phase mask image that represents a part greatly phase-shifted and a part slightly phase-shifted; performing a three-dimensional process by multiplying the phase mask image and the amplitude image, thereby creating a three-dimensional black blood image; and performing a minimum intensity projection of the three-dimensional black blood image thereby visualizing both the veins and the arteries clearly and accurately.