@@ -674,18 +674,16 @@ def _bin_geodesic_shells(
674674 v = mesh .vertices .view (np .ndarray )
675675 e_m = float (mesh .edges_unique_length .mean ()) # mean mesh-edge length
676676
677- c_soma = soma .centre
678- soma_verts = set () if soma .verts is None else set (map (int , soma .verts ))
679- soma_vids = np .fromiter (soma_verts , dtype = np .int64 )
677+ c_soma = soma .centre
678+ soma_verts = set () if soma .verts is None else set (map (int , soma .verts ))
679+ soma_vids = np .fromiter (soma_verts , dtype = np .int64 )
680680
681681 # ------------------------------------------------------------------
682682 # build a vertex list for every connected surface patch
683683 # ------------------------------------------------------------------
684684 comp_vertices = [
685685 np .asarray (c , dtype = np .int64 ) for c in gsurf .components ()
686686 ]
687-
688- soma_vids = np .fromiter (soma_verts , dtype = np .int64 )
689687 all_shells : List [List [np .ndarray ]] = []
690688
691689 for cid , verts in enumerate (comp_vertices ):
@@ -748,7 +746,6 @@ def _bin_geodesic_shells(
748746 comp_idx = np .fromiter (
749747 (inner [i ] for i in comp ), dtype = np .int64
750748 )
751- # comps.append(comp_idx)
752749 if split_elongated_shells and len (comp ) < 1500 : # hard-coded for now, if too large, might be a soma
753750 for part in _split_comp_if_elongated (comp_idx , v ,
754751 aspect_thr = split_aspect_thr ,
@@ -1660,10 +1657,9 @@ def skeletonize(
16601657
16611658 The algorithm proceeds in eight conceptual stages:
16621659
1663- 0. optional pre-skeletonization soma detection>
16641660 1. geodesic shell binning of every connected surface patch
16651661 2. cluster each shell ⇒ interior node with local radius
1666- 3. optional post-skeletonization soma detection>
1662+ 3. optional post-skeletonization soma detection
16671663 4. project mesh edges ⇒ graph edges between nodes
16681664 5. optional collapsing of soma-like/fat nodes near the centroid
16691665 6. optional bridging of disconnected components
0 commit comments