352c352
<     if ( m_LabelImage->GetPixel( neighIndex ) != AlivePoint )
---
>     if ( m_LabelImage->GetPixel( neighIndex ) == FarPoint )
362c362
<     if ( m_LabelImage->GetPixel( neighIndex ) != AlivePoint )
---
>     if ( m_LabelImage->GetPixel( neighIndex ) == FarPoint )

Index: itkFastMarchingImageFilter.h
===================================================================
RCS file: /cvsroot/Insight/Insight/Code/Algorithms/itkFastMarchingImageFilter.h,v
retrieving revision 1.40
diff -u -r1.40 itkFastMarchingImageFilter.h
--- itkFastMarchingImageFilter.h	23 Apr 2009 03:53:35 -0000	1.40
+++ itkFastMarchingImageFilter.h	13 May 2009 17:29:12 -0000
@@ -162,7 +162,7 @@
    * away points; TrialPoints represent points within a narrowband of the
    * propagating front; and AlivePoints represent points which have already
    * been processed. */
-  enum LabelType { FarPoint, AlivePoint, TrialPoint };
+  enum LabelType { FarPoint, AlivePoint, TrialPoint, InitialTrialPoint };
 
   /** LabelImage typedef support. */
   typedef Image<unsigned char, itkGetStaticConstMacro(SetDimension)> LabelImageType;
Index: itkFastMarchingImageFilter.txx
===================================================================
RCS file: /cvsroot/Insight/Insight/Code/Algorithms/itkFastMarchingImageFilter.txx,v
retrieving revision 1.52
diff -u -r1.52 itkFastMarchingImageFilter.txx
--- itkFastMarchingImageFilter.txx	21 Dec 2008 19:13:11 -0000	1.52
+++ itkFastMarchingImageFilter.txx	13 May 2009 17:28:42 -0000
@@ -239,7 +239,7 @@
         }
 
       // make this a trial point
-      m_LabelImage->SetPixel( node.GetIndex(), TrialPoint );
+      m_LabelImage->SetPixel( node.GetIndex(), InitialTrialPoint );
 
       outputPixel = node.GetValue();
       output->SetPixel( node.GetIndex(), outputPixel );
@@ -289,7 +289,8 @@
       } 
 
     // is this node already alive ?
-    if ( m_LabelImage->GetPixel( node.GetIndex() ) != TrialPoint ) 
+    double label = m_LabelImage->GetPixel( node.GetIndex() );
+    if ( label == AlivePoint )
       {
       continue;
       }
@@ -349,7 +350,8 @@
       {
       neighIndex[j] = index[j] - 1;
       }
-    if ( m_LabelImage->GetPixel( neighIndex ) != AlivePoint )
+    double label = m_LabelImage->GetPixel( neighIndex );
+    if ( label != AlivePoint && label != InitialTrialPoint )
       {
       this->UpdateValue( neighIndex, speedImage, output );
       }
@@ -359,7 +361,8 @@
       {
       neighIndex[j] = index[j] + 1;
       }
-    if ( m_LabelImage->GetPixel( neighIndex ) != AlivePoint )
+    label = m_LabelImage->GetPixel( neighIndex );
+    if ( label != AlivePoint && label != InitialTrialPoint )
       {
       this->UpdateValue( neighIndex, speedImage, output );
       }