%?ANISODIFF?-?Anisotropic?diffusion.
%
%?
%??diff?=?anisodiff（im?niter?kappa?lambda?option）
%
%?
%?????????im?????-?input?image
%?????????niter??-?number?of?iterations.
%?????????kappa??-?conduction?coefficient?20-100??
%?????????lambda?-?max?value?of?.25?for?stability
%?????????option?-?1?Perona?Malik?diffusion?equation?No?1
%??????????????????2?Perona?Malik?diffusion?equation?No?2
%
%?Return
%?????????diff???-?diffused?image.
%
%?kappa?controls?conduction?as?a?function?of?gradient.??If?kappa?is?low
%?then?mall?intensity?gradients?are?able?to?block?conduction?and?hence?diffusion
%?across?step?edges.??A?large?value?reduces?the?influence?of?intensity
%?gradients?on?conduction.
%
%?lambda?controls?speed?of?diffusion?（you?usually?want?it?at?a?maximum?of
%?0.25）
%
%?Diffusion?equation?1?preserve?high?contrast?edges?over?low?contrast?ones.
%?Diffusion?equation?2?favours?wide?regions?over?smaller?ones.

%?Reference:?
%?P.?Perona?and?J.?Malik.?
%?Scale-space?and?edge?detection?using?anisotropic?diffusion.
%?IEEE?Transactions?on?Pattern?Analysis?and?Machine?Intelligence?
%?12（7）:629-639?July?1990.


function?diff?=?anisodiff（im?niter?kappa?lambda?option）

if?ndims（im）==3
??error（‘Anisodiff?only?operates?on?2D?grey-scale?images‘）;
end

im?=?double（im）;
[rowscols]?=?size（im）;
diff?=?im;

%{
var?=?2;
x?=?（-4:4）;
g?=?exp（-x.*x/（2*var））;?g??=?g/sum（g）;

blurred?=?conv2（img‘same‘）;
im_b?=?conv2（blurredg‘‘same‘）;?

%}

for?i?=?1:niter
?%?fprintf（‘\rIteration?%d‘i）;

??%?Construct?diffl?which?is?the?same?as?diff?but
??%?has?an?extra?padding?of?zeros?around?it.
??diffl?=?zeros（rows+2?cols+2）;
??diffl（2:rows+1?2:cols+1）?=?diff;

??%?North?South?East?and?West?differences
??deltaN?=?diffl（1:rows2:cols+1）???-?diff;
??deltaS?=?diffl（3:rows+22:cols+1）?-?diff;
??deltaE?=?diffl（2:rows+13:cols+2）?-?diff;
??deltaW?=?diffl（2:rows+11:cols）???-?diff;
??%deltaN?=?diff;deltaW;
??

??%?Conduction

??if?option?==?1
????cN?=?exp（-（deltaN/kappa）.^2）;
????cS?=?exp（-（deltaS/kappa）.^2）;
????cE?=?exp（-（deltaE/kappa）.^2）;
????cW?=?exp（-（deltaW/kappa）.^2）;
??elseif?option?==?2
????cN?=?1./（1?+?（deltaN/kappa）.^2）;
????cS?=?1./（1?+?（deltaS/kappa）.^2）;
????cE?=?1./（1?+?（deltaE/kappa）.^2）;
????cW?=?1./（1?+?（deltaW/kappa）.^2）;
??end

??%?APPLYING?FOUR-POINT-TEMPLETE?FOR?numerical?solution?of?DIFFUSION?P.D.E.
??
??diff?=?diff?+?lambda*（cN.*deltaN?+?cS.*deltaS?+?cE.*deltaE?+?cW.*deltaW）;

%??Uncomment?the?following?to?see?a?progression?of?images
%??subplot（ceil（sqrt（niterations））ceil（sqrt（niterations））?i）
%??imagesc（diff）?colormap（gray）?axis?image

end
%fprintf（‘\n‘）;

?屬性????????????大小?????日期????時(shí)間???名稱
-----------?---------??----------?-----??----
?????目錄???????????0??2012-11-13?15:25??Anisotropic?Diffusion\
?????目錄???????????0??2012-11-13?15:24??Anisotropic?Diffusion\DPAD\
?????文件?????3891320??2012-09-18?18:42??Anisotropic?Diffusion\DPAD\1?On?the?Estimation?of?the?coefficient?of?variation?for?anisotropic?diffusion?speckle?filtering.pdf
?????文件????????9580??2012-08-27?15:10??Anisotropic?Diffusion\DPAD\DPAD_toolbox.zip
?????目錄???????????0??2012-11-13?15:24??Anisotropic?Diffusion\OSRAD\
?????文件?????2616381??2012-09-10?08:34??Anisotropic?Diffusion\OSRAD\1?Oriented?speckle?reducing?anisotropic?diffusion.pdf
?????文件?????????761??2012-11-13?15:23??Anisotropic?Diffusion\OSRAD\說(shuō)明.txt
?????目錄???????????0??2012-11-13?15:24??Anisotropic?Diffusion\PM_AD\
?????文件?????1332481??2012-10-21?22:29??Anisotropic?Diffusion\PM_AD\Scale-space?and?edge?detection?using?anisotropic?diffusion.pdf
?????文件????????2694??2006-02-13?15:32??Anisotropic?Diffusion\PM_AD\anisodiff.m
?????目錄???????????0??2012-11-13?15:11??Anisotropic?Diffusion\SRAD\
?????文件????????5286??2006-04-10?15:42??Anisotropic?Diffusion\SRAD\SRAD.m
?????文件??????574071??2012-09-05?11:52??Anisotropic?Diffusion\SRAD\SRAD.zip
?????文件????????9307??2006-04-10?15:42??Anisotropic?Diffusion\SRAD\SRAD_c.c
?????文件????????8704??2006-04-10?15:42??Anisotropic?Diffusion\SRAD\SRAD_c.dll
?????文件??????768150??2012-10-13?10:38??Anisotropic?Diffusion\SRAD\Speckle?reducing?anisotropic?diffusion.pdf
?????文件????????1290??2006-04-10?15:42??Anisotropic?Diffusion\SRAD\ex_SRAD.m
?????文件???????18011??2006-03-08?17:33??Anisotropic?Diffusion\SRAD\licence.txt
?????文件???????67898??2006-03-08?14:18??Anisotropic?Diffusion\SRAD\ultrasound?phantom.bmp