Combining edge-gradient information to improve adaptive discontinuity-preserving smoothing of multispectral images