Fix disconnected control point input voxels appearing in perimeter file

layerfMRI · Mar 14, 2022 · 7365b03 · 7365b03
1 parent c505112
commit 7365b03
Showing 1 changed file with 5 additions and 1 deletion.
diff --git a/src/LN2_MULTILATERATE.cpp b/src/LN2_MULTILATERATE.cpp
@@ -681,7 +681,11 @@ int main(int argc, char*  argv[]) {
         // Translate 0 crossing
         for (uint32_t ii = 0; ii != nr_voi; ++ii) {
             i = *(voi_id + ii);
-            *(flood_dist_data + i) -= thr_radius;
+            if (*(flood_dist_data + i) != 0) {
+                *(flood_dist_data + i) -= thr_radius;
+            } else {  // These voxels are disconnected to centroid cluster
+                *(flood_dist_data + i) = 666;
+            }
         }
 
         // Mark voxels inside perimeter