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