clear?all
clc
pfilt?=‘maxflat‘;
dfilt?=?‘dmaxflat7‘;
%nlevs?=?[0?0?4?4?5];????%?Number?of?levels?for?DFB?at?each?pyramidal?level
nlevs?=?[123];
th?=?3;?????????????????????%?lead?to?3*sigma?threshold?denoising
rho?=?3;????????????????????%?noise?level

%?Test?image:?the?usual?suspect...
im?=?imread（‘camera.tif‘）;
im?=?double（im）;
%im?=?double（im）/?256;
[mn]=size（im）;
f=?wgn（mn10）;
var=std（f（:））;
nim=im+n;
nim=nim/256;
%?Generate?noisy?image.?
sig?=?std（im（:））;
sigma?=?sig?/?rho;
sig?=?std（im（:））;
sigma?=?sig?/?rho;
%nim?=?im?+?sigma?*?randn（size（im））;
nn=nim;
%nim=[12;34];

%%%%%?Contourlet?Denoising?%%%%%
%?Contourlet?transform
y?=?nsctdec（nimnlevsdfilt?pfilt）;
%y1=y{11};%?64x64??double?array?（低頻系數）
%y2=y{12};%64x128??64x128?
%y3=y{13};%128x128x4?
%y4=y{14};%128x256x4???256x128x4?
??????????%y5=y{15};%[32x128?double]x8??[128x32?double]x8
??????????%y6=y{16};%[32x256?double]x16??[256x32?double]x16

%相關去噪

for?x=2:3
EC=zeros（1length（nlevs）+1）;
EP=zeros（1length（nlevs）+1）;
for?k=1:2.^（x-1）
????nimm=y{1x}{1k};
????[mn]=size（nimm）;
nnimm（::）=zeros（m+2n+2）;
nnimm（11）=0;nnimm（1end）=0;nnimm（end1）=0;nnimm（endend）=0;
????????nnimm（12:end-1）=nimm（1:）;%復制邊界象素點
????????nnimm（2:end-11）=nimm（:1）;
????????nnimm（m+22:end-1）=nimm（m:）;
????????nnimm（2:end-1n+2）=nimm（:n）;
for?i=1:m
????for?j=1:n
??????nnimm（i+1j+1）=nimm（ij）;
????end
end
for?i=1:m
????for?j=1:n
???????Ec1=0;Ec2=0;
????????for?a=1:3?%窗口大小為3
????????????for?b=1:3
????????????????%c=nnim（i+a-1j+b-1）;
????????????????c1=nnimm（i+a-1j+b-1）;
????????????????Ec1=c1+Ec1;
???????????????%?c2=（nnim（i+a-1j+b-1）-Ec1/9）^2;
????????????????%Ec2=c2+Ec2;
???????????????%?noise（ij）=sum（（nnim（i+a-1j+b-1）-sum（nnim（i+a-1j+b-1））/（a*b））^2）;
????????????????
????????????end
????????end
????????c2=（nimm（ij）-Ec1/9）^2;
????????Ec2=c2+Ec2;
???????noise（ij）=Ec2;?
????end
end
????for?i=1:m
????????for?j=1:n??????
??%相關量?
?while?1
??P_coeffs{1x}{1k}=?y{1x}{1k}.*?y{1x+1}{1k};
??EC_coeffs{1x}{1k}=y{1x}{1k}.^2;
??%歸一化?
???EC（1k）=EC（1k）+EC_coeffs{1x}{1k}（ij）;
???EP_coeffs{1x}{1k}=?P_coeffs{1x}{1k}.^2;
???EP（1k）=EP（1k）+EP_coeffs{1x}{1k}（ij）;
?????Pav_coeffs{1x}{1k}=?sqrt（abs（P_coeffs{1x}{1k}）.*sqrt（EP（1k）/EC（1k）））;
?????
??????if（abs（?Pav_coeffs{1x}{1k}（ij））>abs（y{1x}{1k}（ij）））
??????????%ynew{1x}{1k}（ij）=Pav_coeffs{1x}{1k}（ij）;
??????????ynew{1x}{1k}（ij）=sgn（y{1x}{1k}（ij））.*max（abs（y{1x}{1k}（ij）-noise（ij）/var）0）;
??????????y{1x}{1k}（ij）=0;
??????????P_coeffs{1x}{1k}（ij）=0;
??????
??????end
???if（abs（EC_coeffs{1x}{1k}（ij））-min（noise（:））>1e-20）%?do?while
?????break
????end
??end
????????end
???????disp（i）;
???????disp（j）;
????end
end
???????disp（x）;
???????disp（k）;
end
ynew{11}=y{11};ynew{14}=y{14};
m?=?nsctrec（ynew?dfilt?pfilt）;

%XIAOBO
y?=?pdfbdec（nim?‘9-7‘‘pkva‘?zeros（length（nlevs）?1））;
[c?s]?=?pdfb2vec（

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----

?????文件????????336??2008-11-09?23:23??Nonsubsampled?Contourlet?Transform\code\noscdenoise.m

?????文件???????2969??2004-10-17?16:42??Nonsubsampled?Contourlet?Transform\code\nsct\atrousc.c

?????文件???????7168??2004-10-17?16:45??Nonsubsampled?Contourlet?Transform\code\nsct\atrousc.dll

?????文件???????8664??2005-01-24?21:51??Nonsubsampled?Contourlet?Transform\code\nsct\atrousc.mexmac

?????文件???????1234??2005-11-03?09:10??Nonsubsampled?Contourlet?Transform\code\nsct\atrousdec.m

?????文件??????19406??2004-09-09?15:01??Nonsubsampled?Contourlet?Transform\code\nsct\atrousfilters.m

?????文件???????1012??2004-10-17?16:45??Nonsubsampled?Contourlet?Transform\code\nsct\atrousrec.m

?????文件???????3214??2005-01-21?15:59??Nonsubsampled?Contourlet?Transform\code\nsct\decdemo.m

?????文件???????1495??2004-12-22?22:48??Nonsubsampled?Contourlet?Transform\code\nsct\dfbdecdemo.m

?????文件??????13234??2004-10-18?12:28??Nonsubsampled?Contourlet?Transform\code\nsct\dfilters.m

?????文件???????2607??2004-12-22?21:26??Nonsubsampled?Contourlet?Transform\code\nsct\dmaxflat.m

?????文件???????1103??2003-04-10?11:30??Nonsubsampled?Contourlet?Transform\code\nsct\efilter2.m

?????文件???????1861??2003-11-06?00:25??Nonsubsampled?Contourlet?Transform\code\nsct\extend2.m

?????文件????????856??2003-04-10?11:30??Nonsubsampled?Contourlet?Transform\code\nsct\ld2quin.m

?????文件????????662??2003-11-05?10:17??Nonsubsampled?Contourlet?Transform\code\nsct\ldfilter.m

?????文件????????854??2003-11-05?10:19??Nonsubsampled?Contourlet?Transform\code\nsct\mctrans.m

?????文件????????775??2003-04-10?11:30??Nonsubsampled?Contourlet?Transform\code\nsct\modulate2.m

?????文件???????4662??2005-11-03?10:47??Nonsubsampled?Contourlet?Transform\code\nsct\nsctdec.m

?????文件???????3257??2005-11-03?10:33??Nonsubsampled?Contourlet?Transform\code\nsct\nsctrec.m

?????文件???????5046??2005-01-20?23:13??Nonsubsampled?Contourlet?Transform\code\nsct\nsdfbdec.m

?????文件???????4931??2005-01-20?23:14??Nonsubsampled?Contourlet?Transform\code\nsct\nsdfbrec.m

?????文件????????889??2005-12-19?11:12??Nonsubsampled?Contourlet?Transform\code\nsct\nsfbdec.m

?????文件????????821??2004-10-21?07:11??Nonsubsampled?Contourlet?Transform\code\nsct\nsfbrec.m

?????文件???????2424??2005-01-20?23:51??Nonsubsampled?Contourlet?Transform\code\nsct\nssfbdec.m

?????文件???????2785??2005-01-20?23:50??Nonsubsampled?Contourlet?Transform\code\nsct\nssfbrec.m

?????文件???????1552??2004-08-07?18:43??Nonsubsampled?Contourlet?Transform\code\nsct\parafilters.m

?????文件???????1161??2003-11-07?19:58??Nonsubsampled?Contourlet?Transform\code\nsct\qupz.m

?????文件???????1557??2005-01-20?23:45??Nonsubsampled?Contourlet?Transform\code\nsct\README.txt

?????文件???????1981??2003-04-10?11:30??Nonsubsampled?Contourlet?Transform\code\nsct\resampz.m

?????文件???????1001??2005-01-20?23:27??Nonsubsampled?Contourlet?Transform\code\nsct\shownsct.m

............此處省略64個文件信息