Dear Darren
> Note the following lines of code from spm_sections.m. Intervening lines are
> marked by .....
>
> % extract data from SPM{t} [sagittal {Ts} coronal {Tc} transverse {Tt}]
> %----------------------------------------------------------------------------
> Q = find(abs(XYZ(1,:) - L(1)) < V(4)/2);
> .....
> .....
>
> Q = find(abs(XYZ(2,:) - L(2)) < V(5));
> .....
> .....
>
> Q = find(abs(XYZ(3,:) - L(3)) < V(6));
> .....
> .....
>
> My understanding of each of these lines is that it finds data within each
> particular plane based on the location of the ">" markers. My question
> concerns why it is that for the sagittal plane data is found within V(4)/2
> (i.e. half voxel), but for the coronal and transverse planes data is found
> within V(5) and V(6) respectively (i.e. a full voxel). Should not all the
> planes be V(x)/2?
I think you are absolutely right. I cannot remember whether there was
any rationale for this anomaly or not. It will be corrected in SPM98
Many thanks - Karl
%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%%
|