Abstract: A multiscale total variation algorithm is proposed. By taking advantage of the multiscale structure after wavelet transform, the proposed algorithm applies different regularization parameters to low frequency and high frequency blocks. Numerical experiments show that the proposed algorithm eliminates the noise better while preserving the meaningful edges comparing to using Wavelet method or total variation method alone.