In an embodiment, a method is disclosed which includes: carrying out interpolation by using a slice away from a slab boundary, and substituting a slice having slab boundary artifacts existing in the slab boundary, to obtain an interpolated image; carrying out Fourier transform on the interpolated image to generate first K-space data; carrying out Fourier transform on the original image to generate second K-space data; merging the first K-space data with the second K-space data, wherein the weight of the first K-space data is greater than that of the second K-space data in the middle of the K-space, and the weight of the second K-space data is greater than that of the first K-space data at the edge of the K-space; and carrying out inverse Fourier transform on the merged K-space data.