In the medical field, the hardening of tissues is one of important informations used in diagnosis or understanding progress of disease, a quantitative measuring method of hardening is important for objective diagnosis. It has been proposed MRE(Magnetic Resonance Elastography) method that measures an index of hardening, viscoelastic properties in a noninvasive. Because the S/N ratio of MRE images go down when measuring viscoelastic properties from local wavelength and local damping factor of a propagating wave in MRE method, methods using multiple phase MRE images have been examined to decrease the effect of noise. We propose a method measuring viscoelastic properties after Fitting a function for multiple phase MRE images in this research. This proposed method has a advantage to set up arbitrarily the variation rate of a space direction of viscoelastic properties or the spatial resolution of measuring values according to changing of the noise included in images, though it applies viscoelastic wave for multiple phase MRE images. We confirmed the effectiveness of a proposed method by experiment using simulation images and experiment using silicone-gel phantom.