Skip to content

Commit

Permalink
removed 1000
Browse files Browse the repository at this point in the history
  • Loading branch information
scott committed May 2, 2017
1 parent db648c6 commit bbf57f9
Show file tree
Hide file tree
Showing 2 changed files with 11 additions and 11 deletions.
10 changes: 5 additions & 5 deletions svca4_extractRefGui.m
Original file line number Diff line number Diff line change
Expand Up @@ -183,7 +183,7 @@ function extract_Callback(hObject, eventdata, handles)
end

% add times for saving
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes GRAYtac(s,:)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes GRAYtac(s,:)'];

% - Write SVCA reference TAC to txt file
if handles.save_txt.Value == 1 && handles.remCereb.Value == 0
Expand All @@ -209,7 +209,7 @@ function extract_Callback(hObject, eventdata, handles)
end

figure;
plot(svca4.PET_standardEndTimes,GRAYtac(s,:)*1000)
plot(svca4.PET_standardEndTimes,GRAYtac(s,:))
figure;
imagesc(squeeze(PET(:,45,:,20)))
figure;
Expand Down Expand Up @@ -238,7 +238,7 @@ function extract_Callback(hObject, eventdata, handles)
% fclose(fid);

if handles.save_mean.Value == 1 && handles.remCereb.Value == 0
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(GRAYtac)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(GRAYtac)'];

if ifeedback == 0
fname = sprintf('%s/TACs/mean_svcaRef_G%.2fW%.2fB%.2fT%.2f.txt', svca4.outputPath, svca4.quantiles);
Expand All @@ -249,9 +249,9 @@ function extract_Callback(hObject, eventdata, handles)
fprintf(fid,'%.1f\t%.1f\t%.4f\n', myGRAY_TAC');
fclose(fid);
figure;
plot(svca4.PET_standardEndTimes,mean(GRAYtac)*1000)
plot(svca4.PET_standardEndTimes,mean(GRAYtac))
elseif handles.save_mean.Value == 1 && handles.remCereb.Value == 1
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(GRAYtac)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(GRAYtac)'];

if ifeedback == 0
fname = sprintf('%s/TACs/mean_noCB_svcaRef_G%.2fW%.2fB%.2fT%.2f.txt', svca4.outputPath, svca4.quantiles);
Expand Down
12 changes: 6 additions & 6 deletions svca4_voiTACGui.m
Original file line number Diff line number Diff line change
Expand Up @@ -167,7 +167,7 @@ function create_TAC_Callback(hObject, eventdata, handles)
vals(s,t) = mean(PET_t(indsVOI));
end
if handles.save_txt.Value
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes vals(s,:)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes vals(s,:)'];
fname = sprintf('%s/TACs/%s_%s', svca4.outputPath, subj{s},handles.txt_name.String);
fid = fopen(fname, 'w');
fprintf(fid,'start[seconds]\tend[seconds]\tTAC[1/1]\n');
Expand Down Expand Up @@ -199,7 +199,7 @@ function create_TAC_Callback(hObject, eventdata, handles)
end

if handles.save_txt.Value
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes vals(s,:)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes vals(s,:)'];
fname = sprintf('%s/TACs/Norm_%s_%s', svca4.outputPath, subj{s},handles.txt_name.String);
fid = fopen(fname, 'w');
fprintf(fid,'start[seconds]\tend[seconds]\tTAC[1/1]\n');
Expand All @@ -215,7 +215,7 @@ function create_TAC_Callback(hObject, eventdata, handles)
plot(svca4.PET_standardEndTimes,mean(vals),'LineWidth',2)
title('Mean VOI TAC')
if handles.meanTAC.Value == 1
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(vals)'*1000];
myGRAY_TAC = [svca4.PET_standardStartTimes svca4.PET_standardEndTimes mean(vals)'];
if handles.norm.Value
fname = sprintf('%s/TACs/%s_%s', svca4.outputPath, 'Norm_mean',handles.txt_name.String);
else fname = sprintf('%s/TACs/%s_%s', svca4.outputPath, 'mean',handles.txt_name.String);
Expand All @@ -227,17 +227,17 @@ function create_TAC_Callback(hObject, eventdata, handles)
handles.meanTAC.Value = 0;
end

else plot(svca4.PET_standardEndTimes,vals*1000,'LineWidth',2)
else plot(svca4.PET_standardEndTimes,vals,'LineWidth',2)
title('VOI TAC')

end
xlabel('Time after injection (sec)')

switch handles.norm.Value
case 0
ylabel('kbq/cc')
ylabel('1/1')
case 1
ylabel('normalized kbq/cc')
ylabel('normalized 1/1')
end


Expand Down

0 comments on commit bbf57f9

Please sign in to comment.