fixed bug in kmeansSegmentation

ANTsX · Nov 23, 2018 · 8113542 · 8113542
1 parent ab2e7cf
commit 8113542
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/DESCRIPTION b/DESCRIPTION
@@ -1,7 +1,7 @@
 Package: ANTsRCore
 Type: Package
 Title: Core Software Infrastructure for 'ANTsR'
-Version: 0.6.2.2
+Version: 0.6.3
 Date: 2018-11-02
 Authors@R: c(
     person(given = "Brian B.", family = "Avants", role = "cre", 

diff --git a/R/kmeansSegmentation.R b/R/kmeansSegmentation.R
@@ -46,16 +46,16 @@ kmeansSegmentation <- function(img, k, kmask = NA, mrf = 0.1, verbose=FALSE,
                                ...) {
   img = check_ants(img)
   dim <- img@dimension
+  kmimg = iMath(img, "Normalize")
   if (is.na(kmask)) {
-    kmimg = iMath(img, "Normalize")
     kmask <- getMask(kmimg, 0.01, 1, cleanup = 2)
   } else {
     kmask = check_ants(kmask)
   }
   kmask = iMath(kmask, "FillHoles") %>% thresholdImage(1,2)
   nhood <- paste(rep(1, dim), collapse = "x")
   mrf <- paste("[", mrf, ",", nhood, "]", sep = "")
-  kmimg <- atropos( a = kmimg, m = mrf, c = "[5,0]",
+  kmimg <- atropos( a = img, m = mrf, c = "[5,0]",
     i = paste("kmeans[",k, "]", sep = ""),
     verbose = verbose,
     x = kmask,