A phase difference enhanced imaging method comprising: a first step of creating a phase difference image PD(x), x denoting a position within the image, using a first complex image including a magnitude image and a phase image obtained from a magnetic resonance signal and a second complex image obtained by low-pass filtering the first complex image; a second step of creating an enhanced image through selecting a phase corresponding to an objective tissue in the phase difference image PD(x) taking account of change in distribution of phase difference due to the filtering, and selecting an enhancement function w(θ) for enhancing the selected phase, and calculating the enhancement image as w(PD(x)); and a third step of creating a phase difference enhanced image by masking the magnitude image with the enhanced image w(PD(x)), characterized in that the enhancement function w(θ) is expressed by the following expressions: w(θ)=1···(-σ≤θ≤σ) w(θ)=exp(-a*(Abs(θ)-σ) b)···(θ is out of the above rage) a, b, σ: parameter having a value of a real number Abs(θ): absolute value of θ, wherein a and b are determined based on the filter size of the low pass filter, and to maximize contrast between the objective tissue in the phase difference enhanced image and a background thereof or a contrast-to-noise ratio which is based on said contrast, and σ is chosen to reduce noise in the phase difference image PD(x). a first step of creating a phase difference image PD(x), x denoting a position within the image, using a first complex image including a magnitude image and a phase image obtained from a magnetic resonance signal and a second complex image obtained by low-pass filtering the first complex image; a second step of creating an enhanced image through selecting a phase corresponding to an objective tissue in the phase difference image PD(x) taking account of change in distribution of phase difference due to the filtering, and selecting an enhancement function w(θ) for enhancing the selected phase, and calculating the enhancement image as w(PD(x)); and a third step of creating a phase difference enhanced image by masking the magnitude image with the enhanced image w(PD(x)), characterized in that the enhancement function w(θ) is expressed by the following expressions: w(θ)=1···(-σ≤θ≤σ) w(θ)=exp(-a*(Abs(θ)-σ) b)···(θ is out of the above rage) a, b, σ: parameter having a value of a real number Abs(θ): absolute value of θ, wherein a and b are determined based on the filter size of the low pass filter, and to maximize contrast between the objective tissue in the phase difference enhanced image and a background thereof or a contrast-to-noise ratio which is based on said contrast, and σ is chosen to reduce noise in the phase difference image PD(x). w(θ)=1···(-σ≤θ≤σ) w(θ)=exp(-a*(Abs(θ)-σ) b)···(θ is out of the above rage) a, b, σ: parameter having a value of a real number Abs(θ): absolute value of θ, wherein a and b are determined based on the filter size of the low pass filter, and to maximize contrast between the objective tissue in the phase difference enhanced image and a background thereof or a contrast-to-noise ratio which is based on said contrast, and σ is chosen to reduce noise in the phase difference image PD(x).