function Output = denoiseByTv(Ima, Thr) ImaTv = imfilter(Ima, 'Laplacian');ImaTh = (ImaTv< Thr)*ImaTv;output = Ima - ImaTh;end