lszxs4
2012-03-25, 22:57
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));
我想问一下这个语句是什么意思 谢谢了。
d = C( prod( S(1,:) ) + 2 * prod( S(2,:) ) + 1: end); %HH子带系数
%用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));
我想问一下这个语句是什么意思 谢谢了。
d = C( prod( S(1,:) ) + 2 * prod( S(2,:) ) + 1: end); %HH子带系数