THESE/codes/snake/basic_code/ExternalForceImage3D.m

   1 function Eextern = ExternalForceImage3D(I,Wline, Wedge,Sigma)\r
   2 % Eextern = ExternalForceImage3D(I,Wline, Wedge,Sigma)\r
   3 % \r
   4 % inputs, \r
   5 %  I : The image\r
   6 %  Sigma : Sigma used to calculated image derivatives \r
   7 %  Wline : Attraction to lines, if negative to black lines otherwise white\r
   8 %          lines\r
   9 %  Wterm : Attraction to terminations of lines (end points) and corners\r
  10 %\r
  11 % outputs,\r
  12 %  Eextern : The energy function described by the image\r
  13 %\r
  14 % Function is written by D.Kroon University of Twente (July 2010)\r
  15 \r
  16 Ix=ImageDerivatives3D(I,Sigma,'x');\r
  17 Iy=ImageDerivatives3D(I,Sigma,'y');\r
  18 Iz=ImageDerivatives3D(I,Sigma,'z');\r
  19 \r
  20 Eline =  imgaussian(I,Sigma);\r
  21 Eedge = sqrt(Ix.^2 + Iy.^2 + Iz.^2); \r
  22 \r
  23 Eextern= (Wline*Eline - Wedge*Eedge); \r
  24 \r