Solve issue about calling the wrong mask function

Author: P-VS

mems_fmri_run.m

@@ -111,7 +111,7 @@
 
     case 4 %T2* weigghted (wi(t)=TEi * exp(-TEi/T2*(t)),sum(TEi * exp(-TEi/T2*(t))))
 
-        mask = my_spmbatch_mask(tefuncdat(:,:,:,:,1));
+        mask = MEHB_mask(tefuncdat(:,:,:,:,1));
         mask_ind = find(mask>0);
 
         %based on https://github.com/jsheunis/fMRwhy/tree/master
@@ -168,7 +168,7 @@
 
     case 5 %T2* mapping
 
-        mask = my_spmbatch_mask(tefuncdat(:,:,:,:,1));
+        mask = MEHB_mask(tefuncdat(:,:,:,:,1));
         mask_ind = find(mask>0);
 
         %based on https://github.com/jsheunis/fMRwhy/tree/master
@@ -225,4 +225,4 @@
     Vout = MEHB_write_vol_4d(Vout,funcdat(:,:,:,ti));
 
     out(ti) = {Vout.fname};
-end
+end