基于Curvelet变换的去噪

用Curvelet变换进行图像去噪，包括去噪算法和软阈值的提取算法

% fdct_usfft_demo_denoise.m -- Image denoising via curvelet thresholding
 
% img = double（imread（‘Lena.jpg‘））;
%oimg=double（rgb2gray（imread（‘SARfield.bmp‘）））;
oimg=double（rgb2gray（imread（‘race.png‘）））;
% figure;
% imshow（uint8（img））;
img=log（oimg+1e-5）;
n = size（img1）;
sigma =0;        
is_real = 1;
 
noisy_img = img+sigma*randn（n）;
 
disp（‘Compute all thresholds‘）;
X = randn（n）;
tic; C = fdct_usfft（Xis_real）; toc; 
 
% Compute norm of curvelets （Monte Carlo） 
E = cell（size（C））;
 for s=1:length（C）
   E{s} = cell（size（C{s}））;
   for w=1:length（C{s}）
     A = C{s}{w};
     E{s}{w} = median（abs（A（:） - median（A（:））））/.6745; % Estimate noise level with robust estimator
   end
 end
 
% Take curvelet transform
disp（‘ ‘）;
disp（‘Take curvelet transform: fdct_usfft‘）;

 属性            大小     日期    时间   名称
----------- ---------  ---------- -----  ----

     文件       1535  2008-07-28 15:52  curvelet去噪程序\denoise.m

     文件       1474  2008-07-28 15:52  curvelet去噪程序\softthresholded.m

     目录          0  2008-07-28 15:54  curvelet去噪程序

     文件        218  2007-06-05 03:14  www.pudn.com.txt

----------- ---------  ---------- -----  ----

                 3227                    4