From 801cd9e5d190bea3d5a2050e7de078a540437819 Mon Sep 17 00:00:00 2001 From: scott Date: Fri, 5 May 2017 12:09:34 +0200 Subject: [PATCH] vals std nan --- svca4_calculate.m | 4 +++- svca4_calculate_classes.m | 6 ++++-- svca4_iterateGui.m | 8 ++++++-- 3 files changed, 13 insertions(+), 5 deletions(-) diff --git a/svca4_calculate.m b/svca4_calculate.m index 8b1781d..4371f36 100644 --- 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 index 5c80eb9..2c663d9 100644 --- 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 index bf32045..48bebb8 100644 --- 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;