1 / 39

SPM Pre-Processing

SPM Pre-Processing. Oli Gearing + Jack Kelly Methods for Dummies 07-04-2004. Talk Outline. Slice Timing Realignment Coregistration Normalisation Smoothing. Part 1 - Jack. Part 2 - Oli. Talk Outline. Slice Timing Realignment Coregistration Normalisation Smoothing. Henri Lartigue.

Download Presentation

SPM Pre-Processing

An Image/Link below is provided (as is) to download presentation Download Policy: Content on the Website is provided to you AS IS for your information and personal use and may not be sold / licensed / shared on other websites without getting consent from its author. Content is provided to you AS IS for your information and personal use only. Download presentation by click this link. While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server. During download, if you can't get a presentation, the file might be deleted by the publisher.

E N D

Presentation Transcript


  1. SPM Pre-Processing Oli Gearing + Jack Kelly Methods for Dummies 07-04-2004

  2. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing Part 1 - Jack Part 2 - Oli

  3. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing

  4. Henri Lartigue

  5. slit-scan photography exposed last exposed first

  6. slit-scan photography exposed last exposed first

  7. slit-scan photography exposed last exposed first

  8. Slice Timing SPM assumes each scan is ‘instantaneous’ RAW CORRECTED TR TA 0 3 6 time in seconds

  9. Slice timing Only needed if: • Temporal dynamics of evoked responses are important and if • TR is sufficiently small to permit interpolation ( <3 seconds ) • BioPhysical latency is on the order of seconds • Usually unnecessary if latency differences are modelled in SPM analysis “proper” using temporal derivatives.

  10. Slice timing Output: afilename.hdr afilename.img afilename.mat

  11. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing

  12. Subject motion: front line defence

  13. Realignment of subject motion • Why bother? • Subsequent analysis assumes that voxel = bit of brain • (e.g. Subtraction and averaging) • haemodynamic response is small compared to signal from movement • increase sensitivity of T-Test (movement contributes to variance)

  14. Realignment of subject motion • When to do it? • Must be before Normalization • Can be either before or after slice time correction (disadvantages to both options) • For interleaved acquisitions it’s recommended to slice time correct first • For sequential acquisitions it’s recommended to realign first

  15. Realignment of subject motion filename.img filename.mat reslice coregister rfilename.img

  16. Realignment • Realignment involves: • “coregister” - Estimate 6 parameters (3x translations, 3x rotations) of an affine rigid body transform. • Aim: minimize COST FUNCTION computed between each successive scan and a reference scan

  17. Realignment Realignment involves: 2. “reslice” Apply transformation by re-sampling the data

  18. fMRI adjustment In extreme cases, up to 90% of the variance in fMRI time-series can be accounted for by effects of movement AFTER realignment: • subject movement between slice acquisition • interpolation artefacts • nonlinear distortion due to magnetic field inhomogeneities (EPI distortion) • spin-excitation history (especially if TR approaches T1)

  19. fMRI adjustment • Adjustment can be carried out as either: • Part of the pre-processing step or • Embodied in model estimation during the ‘real’ analysis.

  20. EPI undistortion • EPI images are distorted relative to the structural scans • Bigger magnet = more distortion

  21. EPI undistortion Different tissues have different magnetic susceptibilities Magnetic field warps at tissue boundaries But the field gradient encodes position!

  22. EPI undistortion • It is possible to directly measure the magnetic field across the head, and then use this information to undistort the EPI images after reconstruction.

  23. EPI undistortion new feature in SPM2:

  24. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing

  25. Coregistration • Align different modalities (eg PET & MRI) • Align functional (EPI) with structural (T1) • Optimize parameters describing rigid body transformation to match functional with structural • SPM99’s 3-step coregistration procedure is replaced by an “information theoretic objective function” in SPM2

  26. Part 2 • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing Part 1 - Jack Part 2 - Oli

  27. Normalisation and smoothing • The story so far… • fMRI time data set • Movement between scans has been corrected for (realignment) • Functional data has been overlaid onto the high resolution anatomical data (co-registration) • What next…?

  28. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing

  29. Normalisation What do we want from fMRI? Analysis within subject data Analysis between subjects But how do we compare 2 different brains? Squash the subjects data into a common 3D brain space.

  30. The Talairach brain template1 1Talairach and Tournoux, 1988

  31. How is the data warped? • Either anatomical scan or functional data is used to estimate warping parameters, using one of the following models: • 12 parameter affine transformation • Low frequency basis spatial functions • Vector field specifying the mapping for each voxel • For Dummies: find the most probable warping parameters given the data

  32. Template Normalised Image

  33. Problems with normalisation • Structural alignment does not mean functional alignment • Differences in gyral anatomy and physiology lead to non-perfect fit • Strict warping to template will create non-existent features • Brain pathology may confuse the normalising procedure

  34. Talk Outline • Slice Timing • Realignment • Coregistration • Normalisation • Smoothing

  35. Smoothing How? Intensity value of a voxel is replaced by a weighted average of the neighbouring voxels Why smooth? Render the errors more normal in their distribution (I.e. Gaussian) For inter-subject analyses Increase signal-noise ratio

  36. Summary • Realignment - (adjust for movement between slices) • Co-registration - (link functional scans to anatomical scan) • Normalisation - (warp functional data into template space) • Smoothing - (to increase signal to noise ratio)

  37. Any questions • Ask Dan / Will / Lucy

  38. Resources SPM99 manual www.fil.ion.ucl.ac.uk/spm/course/manual/spatial.htm What’s new in SPM2 www.fil.ion.ucl.ac.uk/spm/spm2.html#New Cambridge’s CBU imaging www.mrc-cbu.cam.ac.uk/Imaging/Common/fmridefaults.shtml General info (on everything) www.WikiPedia.org

More Related