more experiments...

[FrontAlgorithms.git] / lib / fpa / Image / ExtractEndPointsAndBifurcationsFromMinimumSpanningTree.hxx

diff --git a/lib/fpa/Image/ExtractEndPointsAndBifurcationsFromMinimumSpanningTree.hxx b/lib/fpa/Image/ExtractEndPointsAndBifurcationsFromMinimumSpanningTree.hxx
index d2bc8287b1b1c0beab591afc7d6a24ddc1d75651..16d03e83ce3aa2e2bb5551f98b5feb0fc186cbc7 100644 (file)
--- a/lib/fpa/Image/ExtractEndPointsAndBifurcationsFromMinimumSpanningTree.hxx
+++ b/lib/fpa/Image/ExtractEndPointsAndBifurcationsFromMinimumSpanningTree.hxx
@@ -116,6 +116,8 @@ _MarkSphere(
   const unsigned long& l
   )
 {
+  double rr = r * double( 1.5 );
+
   // Get marking region
   auto rreg = this->m_MarkImage->GetRequestedRegion( );
   auto spac = this->m_MarkImage->GetSpacing( );
@@ -126,7 +128,7 @@ _MarkSphere(
   bool all = false;
   for( unsigned int d = 0; d < _TImage::ImageDimension; ++d )
   {
-    unsigned long s = std::ceil( r / double( spac[ d ] ) );
+    unsigned long s = std::ceil( rr / double( spac[ d ] ) );
     if( s < 3 )
     {
       s = 3;
@@ -158,9 +160,9 @@ _MarkSphere(
   {
     typename _TImage::PointType pnt;
     this->m_MarkImage->TransformIndexToPhysicalPoint( spIt.GetIndex( ), pnt );
-    // TODO: if( double( pnt.EuclideanDistanceTo( c ) ) <= r || all )
-    // TODO: if( spIt.Get( ) == 0 )
-      spIt.Set( l );
+    if( double( pnt.EuclideanDistanceTo( c ) ) <= rr || all )
+      if( spIt.Get( ) == 0 )
+        spIt.Set( l );
 
   } // rof
 }