Hi,
Looks like you remove the baseline signal from the data matrix, but not from the aif matrix.
%subtraction of baseline signal
sbase=mean(data(1:t1,:,:,:),1);
tmp=data((t1+1):t2,:,:,:)-repmat(mean(data(1:t1,:,:,:),1),[t2-t1 1 1 1]);
data=tmp;
%% Added code here
masktemp = mask(2:end,:,:,:);
clear mask
mask = masktemp;
Can you check that this is correct?
All the best,
James
Hi,
Looks like you remove the baseline signal from the data matrix, but not from the aif matrix.
%subtraction of baseline signal
sbase=mean(data(1:t1,:,:,:),1);
tmp=data((t1+1):t2,:,:,:)-repmat(mean(data(1:t1,:,:,:),1),[t2-t1 1 1 1]);
data=tmp;
%% Added code here
masktemp = mask(2:end,:,:,:);
clear mask
mask = masktemp;
Can you check that this is correct?
All the best,
James