add parameter to phase-based-masking

korbinian90 · Jul 6, 2022 · 3545d04 · 3545d04 · korbinian90 · Jul 11, 2022
1 parent 80454b1
commit 3545d04
Show file tree

Hide file tree

Showing 2 changed files with 3 additions and 3 deletions.
diff --git a/Project.toml b/Project.toml
@@ -1,7 +1,7 @@
 name = "MriResearchTools"
 uuid = "557dad86-b9bd-4533-8525-1a39684d020f"
 authors = ["Korbinian Eckstein"]
-version = "1.0.0"
+version = "1.0.1"
 
 [deps]
 DataStructures = "864edb3b-99cc-5e75-8d2d-829cb0a9cfe8"

diff --git a/src/masking.jl b/src/masking.jl
@@ -69,11 +69,11 @@ Original MATLAB algorithm:
     PB=imclose(PB,se);
     mask{2}=round(imopen(PB,se));
 """
-function phase_based_mask(phase; filter=true)
+function phase_based_mask(phase; filter=true, threshold=1.0)
     strel = sphere(6, ndims(phase))
     laplacian = imfilter(sign.(phase), Kernel.Laplacian(1:ndims(phase), ndims(phase)))
     test = imfilter(abs.(laplacian), strel)
-    PB = test .< (500 * 6)
+    PB = test .< (500 * 6 * threshold)
     if filter
         PB = LocalFilters.closing(PB, strel)
         PB = LocalFilters.opening(PB, strel)