曲波变换核心代码及其实现

曲波变换核心代码及其实现，以及经过本人实践开发

clc;clear all;
img = double（imread（‘Lena.jpg‘））;
n=512;
sigma = 20;
is_real = 1;
noise=randn（n）;
noisy_img = img + sigma*noise;

X = randn（n）;
C = fdct_wrapping（X02）;
xxx=C;

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;
  end
end  

C = fdct_wrapping（noisy_img12）; 
yyy=C;
Ct = C;
for s = 2:length（C）
  thresh = 3*sigma + sigma*（s == length（C））;
  for w = 1:length（C{s}）
      T=thresh*E{s}{w};
      yhx=abs（C{s}{w}）;
    Ct{s}{w} = C{s}{w}.* （yhx> T）;
  end   %T=k*sigma*Ek取3或4thresh=k*sigma
end
restored_img = real（ifdct_wrapping（Ct1））; 
subplot（131）; imagesc（img）; colormap gray; axis（‘image‘）;
subplot（132）; imagesc（noisy_img）; colormap gray; axis（‘image‘）;
subplot（133）; imagesc（restored_img）; colormap gray; axis（‘image‘）;

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

     文件      34011  2014-06-19 15:45  曲波wrap\a2.jpg

     文件      85851  2014-06-19 15:54  曲波wrap\aa.jpg

     文件     627952  2014-09-14 13:39  曲波wrap\atlas_athlon[4].dll.zip

     文件      44640  2014-03-15 15:15  曲波wrap\barbara512.jpg

     文件        922  2015-07-12 21:03  曲波wrap\Cuv_denoise.m

     文件        932  2014-08-30 15:20  曲波wrap\Cuv_denoise2.m

     文件       2076  2015-07-12 21:13  曲波wrap\denoise.m

     文件      15750  2014-07-09 16:20  曲波wrap\fdct_wrapping.m

     文件        966  2014-07-05 20:03  曲波wrap\fdct_wrapping_demo_basic.m

     文件       4021  2014-06-26 13:02  曲波wrap\fdct_wrapping_demo_denoise.m

     文件       8415  2014-07-03 12:14  曲波wrap\fdct_wrapping_demo_denoise_enhanced.m

     文件        895  2014-07-05 20:48  曲波wrap\fdct_wrapping_demo_disp.m

     文件        932  2014-07-05 21:16  曲波wrap\fdct_wrapping_demo_recon.m

     文件       2063  2014-03-15 15:15  曲波wrap\fdct_wrapping_demo_wave.m

     文件       1919  2014-03-15 15:15  曲波wrap\fdct_wrapping_dispcoef.m

     文件       8577  2014-05-12 13:21  曲波wrap\fdct_wrapping_param.m

     文件       8647  2014-03-15 15:15  曲波wrap\fdct_wrapping_param.m~

     文件        832  2014-07-14 15:46  曲波wrap\fdct_wrapping_pos2idx.m

     文件        901  2014-07-09 15:37  曲波wrap\fdct_wrapping_window.m

     文件      16096  2014-04-06 18:03  曲波wrap\ifdct_wrapping.m

     文件      32087  2014-03-15 15:15  曲波wrap\Lena.jpg

     文件      97974  2014-03-15 15:15  曲波wrap\LenaCombined.jpg

     文件     190989  2014-06-08 14:14  曲波wrap\QQ.jpg

     文件       1856  2014-07-02 15:49  曲波wrap\seismodel.m

     目录          0  2015-07-12 21:17  曲波wrap

     文件       2076  2015-07-13 09:24  曲波wrap\denoise.asv

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

              1191380                    26