vals std nan

scott-love · May 5, 2017 · 801cd9e · 801cd9e
1 parent 07c0801
commit 801cd9e
Show file tree

Hide file tree

Showing 3 changed files with 13 additions and 5 deletions.
diff --git a/svca4_calculate.m b/svca4_calculate.m
@@ -34,7 +34,9 @@
     for t=1:svca4.nFrames
         PET_t = TARGET(:,:,:,t);
         vals  = PET_t(indMASK) - mean(PET_t(indMASK));
-        vals = vals/std(vals(:));
+        if std(vals(:)) ~= 0
+            vals = vals/std(vals(:));
+        end
         PET_t_norm = PET_norm(:,:,:,t);
         PET_t_norm(indMASK) = vals;
         PET_norm(:,:,:,t) = PET_t_norm;

diff --git a/svca4_calculate_classes.m b/svca4_calculate_classes.m
@@ -1,7 +1,7 @@
 function [TAC_TABLE,TAC_TABLEnn] = svca4_calculate_classes(svca4)
 
 clear TAC_TABLE TAC_TABLEnn
-global svca4
+%global svca4
 
 svca4
 %%% %%%
@@ -54,7 +54,9 @@
     for t=1:svca4.nFrames
         PET_t = PET(:,:,:,t);
         vals  = PET_t(indMASK) - mean(PET_t(indMASK));
-        vals = vals/std(vals(:));
+        if std(vals(:)) ~= 0
+            vals = vals/std(vals(:));
+        end
         PET_t_norm = PET_norm(:,:,:,t);
         PET_t_norm(indMASK) = vals;
         PET_norm(:,:,:,t) = PET_t_norm;

diff --git a/svca4_iterateGui.m b/svca4_iterateGui.m
@@ -128,7 +128,9 @@ function iterate_Callback(hObject, eventdata, handles)
         for t=1:svca4.nFrames
             PET_t = PET(:,:,:,t);
             vals  = PET_t(indMASK) - mean(PET_t(indMASK));
-            vals = vals/std(vals(:));
+            if std(vals(:)) ~= 0
+                vals = vals/std(vals(:));
+            end
             PET_t_norm = PET_norm(:,:,:,t);
             PET_t_norm(indMASK) = vals;
             PET_norm(:,:,:,t) = PET_t_norm;
@@ -271,7 +273,9 @@ function iterate_Callback(hObject, eventdata, handles)
         for t=1:svca4.nFrames
             PET_t = TARGET(:,:,:,t);
             vals  = PET_t(indMASK) - mean(PET_t(indMASK));
-            vals = vals/std(vals(:));
+            if std(vals(:)) ~= 0
+                vals = vals/std(vals(:));
+            end
             PET_t_norm = PET_norm(:,:,:,t);
             PET_t_norm(indMASK) = vals;
             PET_norm(:,:,:,t) = PET_t_norm;