rmHotellingT2.m: added option to average across slices

julienbesle · Oct 27, 2016 · 4205d10 · 4205d10
1 parent c05d35c
commit 4205d10
Showing 1 changed file with 8 additions and 2 deletions.
diff --git a/mrLoadRet/Plugin/GLM_v2/combineTransformOverlayFunctions/rmHotellingT2.m b/mrLoadRet/Plugin/GLM_v2/combineTransformOverlayFunctions/rmHotellingT2.m
@@ -9,7 +9,8 @@
 
 function [P,N,T2] = rmHotellingT2(varargin)
 
-justDoCentralSlice=true;
+justDoCentralSlice=false;
+average5CentralSlices = false;
 
 amplitudePhase=false;
 
@@ -51,6 +52,11 @@
   oldKkk=kkk;
   kkk=1;
   X = X(:,:,ceil(oldKkk/2),:,:);
+elseif average5CentralSlices
+  mrWarnDlg('(rmHotellingT2) Averaging 5 central slices and restricting analysis to this');
+  oldKkk=kkk;
+  kkk=1;
+  X = nanmean(X(:,:,ceil(oldKkk/2)+(-2:2),:,:),3);
 end
 
 X=reshape(X,[k*kk*kkk n p]);
@@ -96,7 +102,7 @@
 P=reshape(P,[k kk kkk]);
 N=reshape(N,[k kk kkk]);
 
-if justDoCentralSlice
+if justDoCentralSlice || average5CentralSlices
     T2=repmat(T2,[1 1 oldKkk]);
     T2(:,:,[1:ceil(oldKkk/2)-1 ceil(oldKkk/2)+1:11])=NaN;
     P=repmat(P,[1 1 oldKkk]);