I suspect this is a left over from aliasing squeeze into sq
The issue never appeared as this is part of some dead sub-function that never gets called: sub_create_NVR_file_HCP
$ cat NORDIC.m | grep sq\(
Q=sq(abs(KSP2(end-3,3:end-3,:,end))); Q1=mean(Q(:)); % take image noise from the bottom of the volume, where the signal is lowest - JIC
Q=sq(abs(KSP2(end-3,3,:,:))); Q1=mean(Q(:)); % take image noise from the bottom of the volume, where the signal is lowest - JIC
Q=sq(abs(KSP_NOISE(:,:,:,end))); Q11=mean(Q(:));
Q=sq((KSP2(end-3,2,:,:))); Q1a=mean(abs(Q(:))); % take image noise from the bottom of the volume, where the signal is lowest - JIC
Q=sq((KSP2(end-3,end-2,:,:))); Q1b=mean(abs(Q(:))); Q1=min(Q1a,Q1b);
Q=sq((KSP_NOISE(:,:,:,end))); Q11=mean(abs(Q(:)));
Q=sq((KSP2(end-3,2,:,:))); Q1a=mean(abs(Q(:))); % take image noise from the bottom of the volume, where the signal is lowest - JIC
Q=sq((KSP2(end-3,end-2,:,:))); Q1b=mean(abs(Q(:))); Q1=min(Q1a,Q1b);
Q=sq((KSP_NOISE(:,:,:,end))); Q11=mean(abs(Q(:)));
Q=sq((KSP2(end-3,2,:,:))); Q1a=mean(abs(Q(:))); % take image noise from the bottom of the volume, where the signal is lowest - JIC
Q=sq((KSP2(end-3,end-2,:,:))); Q1b=mean(abs(Q(:))); Q1=min(Q1a,Q1b);
Q=sq((KSP_NOISE(:,:,:,end))); Q11=mean(abs(Q(:)));
I suspect this is a left over from aliasing
squeeze
intosq
The issue never appeared as this is part of some dead sub-function that never gets called:
sub_create_NVR_file_HCP