function?[kappa]?=?curvature（phi?dx?dy）
%
%?function?[kappa]?=?curvature（phi?dx?dy）
%
%?Calculates?the?curvature?of?the?function?phi（xy）
%
%?Author:?Baris?Sumengen??sumengen@ece.ucsb.edu
%?http://vision.ece.ucsb.edu/~sumengen/
%



dx2?=?dx*dx;
dy2?=?dy*dy;

phi_xp1?=?zeros（size（phi））;
phi_xm1?=?zeros（size（phi））;
phi_yp1?=?zeros（size（phi））;
phi_ym1?=?zeros（size（phi））;

phi_xp1（:2:end）?=?phi（:1:end-1）;
phi_xp1（:1）?=?2*phi_xp1（:2）-phi_xp1（:3）;
phi_xm1（:1:end-1）?=?phi（:2:end）;
phi_xm1（:end）?=?2*phi_xm1（:end-1）-phi_xm1（:end-2）;

phi_yp1（2:end:）?=?phi（1:end-1:）;
phi_yp1（1:）?=?2*phi_yp1（2:）-phi_yp1（3:）;
phi_ym1（1:end-1:）?=?phi（2:end:）;
phi_ym1（end:）?=?2*phi_ym1（end-1:）-phi_ym1（end-2:）;

phi_x?=?（phi_xm1?-?phi_xp1）/（2*dx）;
phi_y?=?（phi_ym1?-?phi_yp1）/（2*dy）;
phi_xx?=?（phi_xm1?-?2*phi?+?phi_xp1）/dx2;
phi_yy?=?（phi_ym1?-?2*phi?+?phi_yp1）/dy2;

dummy1?=?circshift（phi?[11]）;
dummy1（12:end）?=?2*dummy1（22:end）-dummy1（32:end）;
dummy1（:1）?=?2*dummy1（:2）-dummy1（:3）;

dummy2?=?circshift（phi?[-1-1]）;
dummy2（end1:end-1）?=?2*dummy2（end-11:end-1）-dummy2（end-21:end-1）;
dummy2（:end）?=?2*dummy2（:end-1）-dummy2（:end-2）;

dummy3?=?circshift（phi?[1-1]）;
dummy3（11:end-1）?=?2*dummy3（21:end-1）-dummy3（31:end-1）;
dummy3（:end）?=?2*dummy3（:end-1）-dummy3（:end-2）;

dummy4?=?circshift（phi?[-11]）;
dummy4（end2:end）?=?2*dummy4（end-12:end）-dummy4（end-22:end）;
dummy4（:1）?=?2*dummy4（:2）-dummy4（:3）;

phi_xy?=?（dummy1?+?dummy2?-?dummy3?-?dummy4）/（4*dx*dy）;

abs_grad_phi_sq?=?（phi_x.*phi_x?+?phi_y.*phi_y）;
abs_grad_phi_cube?=?abs_grad_phi_sq.^1.5;

kappa?=?（phi_xx.*phi_y.*phi_y?-?2.*phi_y.*phi_x.*phi_xy?...
????+?phi_yy.*phi_x.*phi_x）./（abs_grad_phi_cube?+?（abs_grad_phi_cube?==?0））;

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

????.......??????1730??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_ENO1_SD.m

????.......??????1302??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_ENO1.m

????.......??????1463??2005-05-05?00:00??LevelSetMethods\evolve_normal_ENO3.m

????.......??????1461??2005-05-05?00:00??LevelSetMethods\evolve_normal_ENO2.m

????.......??????1463??2005-05-05?00:00??LevelSetMethods\evolve_normal_ENO1.m

????.......??????1749??2005-05-05?00:00??LevelSetMethods\evolve_kappa.m

????.......?????11010??2005-05-12?07:10??LevelSetMethods\evolve2D.m

????.......??????1600??2005-05-05?00:00??LevelSetMethods\der_WENO_plus.m

????.......??????1602??2005-05-05?00:00??LevelSetMethods\der_WENO_minus.m

????.......??????1521??2005-05-05?00:00??LevelSetMethods\der_ENO3_plus.m

????.......??????1525??2005-05-05?00:00??LevelSetMethods\der_ENO3_minus.m

????.......??????1246??2005-05-05?00:00??LevelSetMethods\der_ENO2_plus.m

????.......??????1250??2005-05-05?00:00??LevelSetMethods\der_ENO2_minus.m

????.......???????618??2005-05-06?00:00??LevelSetMethods\der_ENO1_plus.m

????.......???????624??2005-05-05?00:00??LevelSetMethods\der_ENO1_minus.m

????.......??????1786??2005-05-05?00:00??LevelSetMethods\curvature.m

????.......???????408??2005-05-05?00:00??LevelSetMethods\get_dt_normal_vector_kappa.m

????.......???????363??2005-05-05?00:00??LevelSetMethods\get_dt_normal_vector.m

????.......???????399??2005-05-05?00:00??LevelSetMethods\get_dt_normal_kappa.m

????.......???????356??2005-05-05?00:00??LevelSetMethods\get_dt_normal.m

????.......???????360??2005-05-05?00:00??LevelSetMethods\get_dt_kappa.m

????.......???????673??2005-05-05?00:00??LevelSetMethods\evolve_vector_WENO.m

????.......???????672??2005-05-05?00:00??LevelSetMethods\evolve_vector_ENO3.m

????.......???????672??2005-05-05?00:00??LevelSetMethods\evolve_vector_ENO2.m

????.......???????674??2005-05-05?00:00??LevelSetMethods\evolve_vector_ENO1.m

????.......??????1465??2005-05-05?00:00??LevelSetMethods\evolve_normal_WENO.m

????.......??????1640??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_WENO_SD.m

????.......??????1302??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_WENO.m

????.......??????1642??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_ENO3_SD.m

????.......??????1304??2005-05-05?00:00??LevelSetMethods\evolve_normal_vector_ENO3.m

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