An object of the invention is to obtain a magnetic resonance spectroscopic image to which the MAC summation is applied with high accuracy and in short time, even though a phase characteristic distribution of the MAC has a spatial non-uniformity, in the MRSI measurement using a magnetic resonance imaging apparatus provided with a MAC. A snapshot marker is inserted into a cache, data continues to be flushed to disk until the snapshot marker is flushed from the cache, and a snapshot is performed. The method further includes a cell classification step of classifying each of the cells into an internal cell positioned inside or outside the target object or a boundary cell including the boundary data, and a cut point determination step of determining cut points of edges of the boundary cell based on the boundary data. In the method, the first substrate is formed.