A novel image denoising method is proposed based on multiscale wavelet thresholding (WT) and bilateral filtering (BF). First, the image is decomposed into multiscale subbands by wavelet transform. Then, from the top scale to the bottom scale, we apply BF to the approximation subbands and WT to the detail subbands. The filtered subbands are reconstructed back to approximation subbands of the lower scale. Finally, subbands are reconstructed in all the scales, and in this way the denoised image is formed. Different from conventional methods such as WT and BF, it can smooth the low-frequency noise efficiently. Experiment results on the image Lena and Rice show that the peak signal-to-noise ratio (PSNR) is improved by at least 3 dB and 0.7 dB compared with using the WT and BF, respectively. In addition, the computational time of the proposed method is almost comparable with that of WT but much less than that of BF.