untouch

scott-love · Apr 4, 2017 · 81e0560 · 81e0560
1 parent 40140be
commit 81e0560
Show file tree

Hide file tree

Showing 2 changed files with 9 additions and 9 deletions.
diff --git a/svca4_PlotClassGui.m b/svca4_PlotClassGui.m
@@ -33,10 +33,10 @@ function svca4_PlotGui_OpeningFcn(hObject, eventdata, handles, varargin)
 % This sets up the initial plot - only do when we are invisible
 % so window can get raised using svca4_plotclassgui.
 if strcmp(get(hObject,'Visible'),'off')
-    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classes_it00(svca4.GMWM_sel,1,:))),'-b','LineWidth',2); hold on
-    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classes_it00(svca4.GMWM_sel,2,:))),'-g','LineWidth',2)
-    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classes_it00(svca4.BLOOD_sel,3,:))),'-r','LineWidth',2)
-    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classes_it00(svca4.TSPO_sel,4,:))),'-k','LineWidth',2)
+    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classesnn_it00(svca4.GMWM_sel,1,:))),'-b','LineWidth',2); hold on
+    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classesnn_it00(svca4.GMWM_sel,2,:))),'-g','LineWidth',2)
+    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classesnn_it00(svca4.BLOOD_sel,3,:))),'-r','LineWidth',2)
+    plot(svca4.PET_standardEndTimes,nanmean(squeeze(svca4.classesnn_it00(svca4.TSPO_sel,4,:))),'-k','LineWidth',2)
     legend('Grey','White','Blood','TSPO')
     xlabel('Time (sec)')
     set(gca,'FontSize',14)

diff --git a/svca4_calculate_classes.m b/svca4_calculate_classes.m
@@ -35,12 +35,12 @@
     svca4.PET_durations = svca4.PET_standardDurations;
 
     %%% load brain mask %%%
-    MASK_struct = load_nii(fullfile(svca4.MASK_dir, svca4.MASK_list{fi}));
+    MASK_struct = load_untouch_nii(fullfile(svca4.MASK_dir, svca4.MASK_list{fi}));
     MASK = single(MASK_struct.img);
     %clear MASK_struct
 
     %%% load PET image %%%
-    PET_struct = load_nii(fullfile(svca4.PET_dir, svca4.PET_list{fi}));
+    PET_struct = load_untouch_nii(fullfile(svca4.PET_dir, svca4.PET_list{fi}));
     PET = single(PET_struct.img);
     svca4.Res = PET_struct.hdr.dime.pixdim([2 4 3]); %
     xDim = size(PET,1);
@@ -66,7 +66,7 @@
     if isBLOOD
         if isfield(svca4,'BANANA_list')
             BANANA_fname = fullfile(svca4.BANANA_dir, svca4.BANANA_list{fi});
-            BANANA_struct = load_nii(BANANA_fname);
+            BANANA_struct = load_untouch_nii(BANANA_fname);
             BANANA = single(BANANA_struct.img); clear BANANA_struct
             BM4D = repmat(BANANA, [1 1 1 numel(svca4.BLOOD_frames)]);
             firstFrames = PET_norm(:,:,:,svca4.BLOOD_frames).*single(BM4D);
@@ -92,7 +92,7 @@
     isGMWM = any(svca4.GMWM_sel==fi);
     if isGMWM
         SEG_fname = fullfile(svca4.SEG_dir, svca4.SEG_list{fi});
-        SEG_struct = load_nii(SEG_fname);
+        SEG_struct = load_untouch_nii(SEG_fname);
         GM = single(SEG_struct.img).*MASK; clear SEG_struct;
         WM = GM; WM(WM~=2) = 0; WM(WM==2)=1;
         GM(GM~=1) = 0;
@@ -124,7 +124,7 @@
     isINF = any(svca4.TSPO_sel==fi);
     if isINF
         INF_fname = fullfile(svca4.INF_dir, svca4.INF_list{fi});
-        INF_struct = load_nii(INF_fname);
+        INF_struct = load_untouch_nii(INF_fname);
         INF = single(INF_struct.img); clear INF_struct;
         % it's not here we want to dilate!
 %         if svca4.TSPODilateParameter