%?本程序用于測試超聲圖像的去噪效果
%?方法包括中值濾波、自適應維納濾波、各向異性擴散、小波變換以及contourlet變換
%?白培瑞?20110423?調試通過
%?有幾個問題需要驗證：1）信噪比的計算公式是否認可？
%2）各種去噪方法的參數選擇對信噪比的影響？
clear?all
close?all
%**************************************************************************
%?讀入原始圖像
%im=phantom（512）+randn（512）;???????%原程序所處理的圖像
im=imread（‘ph055test-2.JPG‘）;?????%phantom055的圖像，儀器是ALOKA的a10
im=rgb2gray（uint8（im））;?
im=double（im）;??
figure?
subplot（321）?
imshow（im[]）?
title（‘原始圖像‘）?
%**************************************************************************
%?中值濾波效果
im_medfilt=medfilt2（im[5?5]）;
subplot（322）?
imshow（im_medfilt[]）?
title（sprintf（‘中值濾波去噪?（SNR?=?%.2f?dB）‘SNR（imim_medfilt））?‘FontSize‘?10）;
%**************************************************************************
%?維納濾波效果
im_wiener=wiener2（im[5?5]）;
subplot（323）?
imshow（im_wiener[]）?
title（sprintf（‘維納濾波去噪?（SNR?=?%.2f?dB）‘SNR（imim_wiener））?‘FontSize‘?10）;
%*************************************************************************
%?各向異性擴散去噪效果
num_iter=15;
delta_t=1/7;
kappa=30;
option=2;
%?ad=anisodiff2D（snum_iterdelta_tkappaoption）;
%Convert?input?image?to?double.
im=double（im）;???%輸入圖像和輸出圖像必須是雙精度的
%?PDE?（partial?differential?equation）?initial?condition.
diff_im?=im;
%?Center?pixel?distances.
dx=1;
dy=1;
dd=sqrt（2）;
%?2D?convolution?masks?-?finite?differences.
hN=[0?1?0;?0?-1?0;?0?0?0];
hS=[0?0?0;?0?-1?0;?0?1?0];
hE=[0?0?0;?0?-1?1;?0?0?0];
hW=[0?0?0;?1?-1?0;?0?0?0];
hNE=[0?0?1;?0?-1?0;?0?0?0];
hSE=[0?0?0;?0?-1?0;?0?0?1];
hSW=[0?0?0;?0?-1?0;?1?0?0];
hNW=[1?0?0;?0?-1?0;?0?0?0];
%Anisotropic?diffusion.?各向異性擴散
for?t=1:num_iter
?%?Finite?differences.?[imfilter（..‘conv‘）?can?be?replaced?by?conv2（..‘same‘）]
????????nablaN=imfilter（diff_imhN‘conv‘）;
????????nablaS=imfilter（diff_imhS‘conv‘）;???
????????nablaW=imfilter（diff_imhW‘conv‘）;
????????nablaE=imfilter（diff_imhE‘conv‘）;???
????????nablaNE=imfilter（diff_imhNE‘conv‘）;
????????nablaSE=imfilter（diff_imhSE‘conv‘）;???
????????nablaSW=imfilter（diff_imhSW‘conv‘）;
????????nablaNW=imfilter（diff_imhNW‘conv‘）;?
??%?Diffusion?function.
????????if?option==1
????????????cN=exp（-（nablaN/kappa）.^2）;
????????????cS=exp（-（nablaS/kappa）.^2）;
????????????cW=exp（-（nablaW/kappa）.^2）;
????????????cE=exp（-（nablaE/kappa）.^2）;
????????????cNE=exp（-（nablaNE/kappa）.^2）;
????????????cSE=exp（-（nablaSE/kappa）.^2）;
????????????cSW=exp（-（nablaSW/kappa）.^2）;
????????????cNW=exp（-（nablaNW/kappa）.^2）;
????????elseif?option==2
????????????cN=1./（1+（nablaN/kappa）.^2）;
????????????cS=1./（1+（nablaS/kappa）.^2）;
????????????cW=1./（1+（nablaW/kappa）.^2）;
????????????cE=1./（1+（nablaE/kappa）.^2）;
????????????cNE=1./（1+（nablaNE/kappa）.^2）;
????????????cSE=1./（1+（nablaSE/kappa）.^2）;
????????????cSW=1./（1+（nablaSW/kappa）.^2）;
????????????cNW=1./（1+（nablaNW/kappa）.^2）;
????????end
?%?Discrete?PDE?solution.
????????diff_im=diff_im?+?...
??????????????????delta_t*（...
??????????????????（1/（dy^2））*cN.*nabl

?屬性????????????大小?????日期????時間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2011-09-05?18:09??non_rein_20110706\
?????文件???????21670??2011-03-30?10:42??non_rein_20110706\00000010-1.jpg
?????文件???????41344??2011-03-28?09:29??non_rein_20110706\00000010.jpg
?????文件????????4761??2011-05-11?10:42??non_rein_20110706\bpr_denosing.asv
?????文件????????4753??2011-05-11?10:43??non_rein_20110706\bpr_denosing.m
?????文件????????3359??2010-04-20?16:34??non_rein_20110706\Cell_area.m
?????目錄???????????0??2011-09-20?10:39??non_rein_20110706\Chan-Vese\
?????文件???????40154??2011-03-28?09:29??non_rein_20110706\Chan-Vese\00000001.jpg
?????文件???????44478??2011-03-28?09:30??non_rein_20110706\Chan-Vese\00000006.jpg
?????文件???????41344??2011-03-28?09:29??non_rein_20110706\Chan-Vese\00000010.jpg
?????文件???????16663??2009-03-26?19:20??non_rein_20110706\Chan-Vese\4colors.JPG
?????文件???????24291??2009-03-26?19:20??non_rein_20110706\Chan-Vese\anti-mass.jpg
?????文件????????2879??2009-03-26?19:20??non_rein_20110706\Chan-Vese\brain.jpg
?????文件?????????848??2009-03-26?19:20??non_rein_20110706\Chan-Vese\checkstop.m
?????文件???????14373??2009-03-26?19:41??non_rein_20110706\Chan-Vese\chenvese.m
?????文件???????39889??2011-03-27?17:23??non_rein_20110706\Chan-Vese\cyst_sim.gif
?????文件????????1973??2011-03-11?10:31??non_rein_20110706\Chan-Vese\demo?chenvese.m
?????文件????????2021??2011-03-27?12:12??non_rein_20110706\Chan-Vese\demochenvese.m
?????文件??????136210??2009-03-26?19:20??non_rein_20110706\Chan-Vese\flowers.jpg
?????文件?????????389??2009-03-26?19:20??non_rein_20110706\Chan-Vese\Heaviside.m
?????文件??????809765??2009-03-26?19:20??non_rein_20110706\Chan-Vese\HELP.pdf
?????文件?????????767??2009-03-26?19:20??non_rein_20110706\Chan-Vese\kappa.m
?????文件????????1327??2009-10-05?23:09??non_rein_20110706\Chan-Vese\license.txt
?????文件???????23792??2011-03-27?09:32??non_rein_20110706\Chan-Vese\long?axis.jpg
?????文件????????8687??2011-03-27?15:21??non_rein_20110706\Chan-Vese\long_axis.jpg
?????文件????????2101??2009-03-26?19:20??non_rein_20110706\Chan-Vese\maskcircle2.m
?????文件????????4052??2009-03-26?19:01??non_rein_20110706\Chan-Vese\READ?ME.txt
?????文件?????????855??2009-03-26?19:20??non_rein_20110706\Chan-Vese\reinitialization.m
?????文件?????????706??2009-03-26?19:20??non_rein_20110706\Chan-Vese\showphi.m
?????文件????????5281??2011-03-28?14:17??non_rein_20110706\Chan-Vese\ultrasoundimageprocess.asv
?????目錄???????????0??2011-09-05?18:09??non_rein_20110706\Chan-Vese\參考文獻\
............此處省略38個文件信息