From ebe5661476f1a0688d05e645a271c864b860d6b4 Mon Sep 17 00:00:00 2001
From: "eduardo.davila@creatis.insa-lyon.fr"
 <eduardo.davila@creatis.insa-lyon.fr>
Date: Fri, 27 Feb 2026 14:35:40 +0100
Subject: [PATCH] Clean code

---
 ...bbcreaVtkBifurcationsInSkeletonization.cxx | 56 ++++++++++---------
 1 file changed, 31 insertions(+), 25 deletions(-)

diff --git a/bbtk_creaVtk_PKG/src/bbcreaVtkBifurcationsInSkeletonization.cxx b/bbtk_creaVtk_PKG/src/bbcreaVtkBifurcationsInSkeletonization.cxx
index 8b05524..477f9e3 100644
--- a/bbtk_creaVtk_PKG/src/bbcreaVtkBifurcationsInSkeletonization.cxx
+++ b/bbtk_creaVtk_PKG/src/bbcreaVtkBifurcationsInSkeletonization.cxx
@@ -48,37 +48,43 @@ void BifurcationsInSkeletonization::Process()
             tmpLstX.erase( tmpLstX.begin() );
             tmpLstY.erase( tmpLstY.begin() );
             tmpLstZ.erase( tmpLstZ.begin() );
-            w   = w + 1;
-            tmpImage->SetScalarComponentFromDouble( px,py,pz,0,w);
-            ok=false;
-            for ( i=px-1 ; i<=px+1 ; i++)
+            if (tmpImage->GetScalarComponentAsDouble(px,py,pz,0)<=1)
             {
-                for ( j=py-1 ; j<=py+1;j++)
+                w   = w + 1;
+                tmpImage->SetScalarComponentFromDouble( px,py,pz,0, w );
+                ok=false;
+                for ( i=px-1 ; i<=px+1 ; i++)
                 {
-                    for ( k=pz-1 ; k<=pz+1 ; k++ )
+                    for ( j=py-1 ; j<=py+1;j++)
                     {
-                        if (!((px==i) && (py==j) && (pz==k)))
+                        for ( k=pz-1 ; k<=pz+1 ; k++ )
                         {
-                            if ((i>=0) &&(j>=0) &&(k>=0) &&(i<dim[0]) &&(j<dim[1]) &&(k<dim[2]) )
+                            if (!((px==i) && (py==j) && (pz==k)))
                             {
-                                if (1==tmpImage->GetScalarComponentAsDouble(i,j,k,0))
+                                if ((i>=0) &&(j>=0) &&(k>=0) &&(i<dim[0]) &&(j<dim[1]) &&(k<dim[2]) )
                                 {
-                                    ok=true;
-                                    tmpLstX.push_back(i);
-                                    tmpLstY.push_back(j);
-                                    tmpLstZ.push_back(k);
-                                } // if graylevel == 1
-                            } // if ijk>0  &&  ijk<dim
-                        } // if px,py,pz != i,j,k
-                    } // for k
-                } // for j
-            } // for i
-            if (ok==false)
-            {
-                leavesLstX.push_back(px);
-                leavesLstY.push_back(py);
-                leavesLstZ.push_back(pz);
-            } //
+                                    if (1==tmpImage->GetScalarComponentAsDouble(i,j,k,0))
+                                    {
+                                        ok=true;
+//                                        tmpLstX.push_back(i);
+//                                        tmpLstY.push_back(j);
+//                                        tmpLstZ.push_back(k);
+                                                                                tmpLstX.insert(tmpLstX.begin() , i);
+                                                                                tmpLstY.insert(tmpLstY.begin() , j);
+                                                                                tmpLstZ.insert(tmpLstZ.begin() , k);
+                                    } // if graylevel == 1
+                                } // if ijk>0  &&  ijk<dim
+                            } // if px,py,pz != i,j,k
+                        } // for k
+                    } // for j
+                } // for i
+                if (ok==false)
+                {
+                    leavesLstX.push_back(px);
+                    leavesLstY.push_back(py);
+                    leavesLstZ.push_back(pz);
+                } //
+            }// if value > 1
         } // while tmpLstX.size()!=0
         
         size = leavesLstX.size();
-- 
2.53.0