English
- function thr = Donoho(x)
- %用Donoho通用阈值公式计算阈值 x为要进行处理的图像
- % thr = delta * sqrt( 2 * log(n))
- % n为信号的长度或尺寸
- % delta = MAD / 0.6745 -经验公式，其中MAD为小波分解后高子带系数的中值
- n = prod( size(x) ); %图像尺寸
- %计算delta
- [C, S] = wavedec2(x, 1, 'db1'); %小波分解
- d = C( prod( S(1,:) ) + 2 * prod( S(2,:) ) + 1 : end); %HH子带系数
- delta = median( abs(d) ) / 0.6745;
- %计算阈值
- thr = delta * sqrt(2*log(n));
