Abstract: During the interpolation process of Partial Volume Distribution, the random perturbation relating to the nine grid points around the perturbing point is taken into account, and the perturbation is evatuated to calculate the weights in PV interpolation. Experimental results indicate this algorithm can avoid the local extremum on grid points and non-grid points, and effectively smooth the object function through mutual information.