%?ANISODIFF?-?Anisotropic?diffusion.
%
%?Usage:
%??diff?=?anisodiff（im?niterations?kappa?lambda?option）
%
%?????????im??-?input?image
%?????????nititerations?
%?????????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
%
%?Reference:?
%?P.?Perona?and?J.?Malik.?
%?Scale-space?and?edge?detection?using?ansotropic?diffusion.
%?IEEE?Transactions?on?Pattern?Analysis?and?Machine?Intelligence?
%?12（7）:629-639?July?1990.
%
%?Author:?Peter?Kovesi???pk@cs.uwa.edu.au
%?Department?of?Computer?Science?&?Software?Engineering
%?The?University?of?Western?Australia
%
%?June?2000?????????

function?diff?=?anisodiff2（im?niterations?kappa?lambda?option）

im?=?double（im）;
[rowscols]?=?size（im）;
diff?=?im;
??
for?i?=?1:niterations
??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;

??%?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+exp（-（deltaN/kappa）.^2））;
????cS?=?1./（1+exp（-（deltaS/kappa）.^2））;
????cE?=?1./（1+exp（-（deltaE/kappa）.^2））;
????cW?=?1./（1+exp（-（deltaW/kappa）.^2））;
??end


??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í)間???名稱
-----------?---------??----------?-----??----

?????文件???????4155??2002-07-15?16:28??F0F1F2F3\ftrack.m

?????文件????????290??2007-05-24?20:27??F0F1F2F3\demo01.m

?????文件????????414??2007-05-19?11:10??F0F1F2F3\nanmax.m

?????文件????????741??2007-05-19?11:06??F0F1F2F3\nanmean.m

?????文件????????414??2007-05-19?11:09??F0F1F2F3\nanmin.m

?????文件?????812656??2002-03-11?14:27??F0F1F2F3\wioioi01.mat

?????文件???????1718??2002-06-18?01:47??F0F1F2F3\lpcsgram.m

?????文件???????1069??2002-03-11?15:20??F0F1F2F3\frmnts1.m

?????文件???????1928??2001-11-20?16:43??F0F1F2F3\anisodiff.m

?????目錄??????????0??2007-05-26?14:30??F0F1F2F3

-----------?---------??----------?-----??----

???????????????823603????????????????????11