1 / 33

# Statistical Parametric Mapping - PowerPoint PPT Presentation

Statistical Parametric Mapping. Lecture 11 - Chapter 13 Head motion and correction. Textbook : Functional MRI an introduction to methods , Peter Jezzard, Paul Matthews, and Stephen Smith. Many thanks to those that share their MRI slides online. http://www.fmrib.ox.ac.uk/fsl/.

I am the owner, or an agent authorized to act on behalf of the owner, of the copyrighted work described.

## PowerPoint Slideshow about ' Statistical Parametric Mapping' - ormand

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.While downloading, if for some reason you are not able to download a presentation, the publisher may have deleted the file from their server.

- - - - - - - - - - - - - - - - - - - - - - - - - - E N D - - - - - - - - - - - - - - - - - - - - - - - - - -
Presentation Transcript

### Statistical Parametric Mapping

Lecture 11 - Chapter 13

Textbook: Functional MRI an introduction to methods, Peter Jezzard, Paul Matthews, and Stephen Smith

Many thanks to those that share their MRI slides online

• FSL software for Motion Correction using FMRIB’s Linear Image Registration Tool - FLIRT.

2

3

4

m14, m24, and m34 are x, y, and z translations

3 each rotations, scales, shears

Let R represent (x, y, z) column vectors and M the 4x4 transform matrix

R’ = M R when you want to calculate primed location from unprimed one

R = M-1 R’ when you want to calculate unprimed location from primed one

4x4 translation matrix T (3 parameters)

M= [S][Rx][Ry][Rz][T]

4x4 scale matrix S (3 parameters)

rotation matrix Rx for rotation about x-axis

(1 parameter)

rotation matrix Ry for rotation about y-axis

(1 parameter)

rotation matrix Rz for rotation about z-axis

(1 parameter)

Order of matrices important. Above ordering does translations first to match origins, then rotations about new origin, and finally scaling of the aligned image. For MC of fMRI no scale, i.e. use on 3 translation and 3 rotation parameters.

2

3

4

Nearest neighbor leads to similar image histograms in transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

Mango provides these three interpolation options when saving an image with transform applied.

1 transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

2

3

4

Same Modality transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

Motion Correction

1 transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

2

3

4

1 transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

2

3

Talariach's Coordinate System transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

Z = +1 mm

• AC-PC line

• AC as origin

• Bounding Box

• 136 x 172 x 118 mm

• Right-handed system

Origin

(AC)

Averages of Anatomy & Function transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

Average MRI - 16 Subjects

Average PET - 15 Subjects

( FDG )

CC transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

TN

SC

CB

Fitted AC-PC Line

Manual SN with High-Resolution MRI

Mango has a plugin for fitting high resolution MRI to the Talairach standard.

CT transformed images but least accurate in terms of interpolation values. Sinc interpolation is considered to be the most accurate. Trilinear is intermediate in quality but relative fast.

PET

MRI

SPECT

Multi-modality registration to a standard brain space (i.e. a particular brain atlas).

Using FLIRT to fit a T1W MR image to the MNI305 3-D average brain template (template brain feature outline indicated by red lines).

Note the large rotation about the y-axis indicated in the left image.

Atlas Based Registration of fMRI brain template (template brain feature outline indicated by red lines).

EPI

structural

MNI152

M1

M2

M= [M2][M1]

MNI152 structural brain template is average of 152 3-D T1W images after affine registration.

What is a 7 DOF transform?

All three scales are identical to correct for possible differences in spatial calibrations.

Global (affine) vs. Regional (warping) to target brain brain template (template brain feature outline indicated by red lines).

4x4 affine transform registration to target brain

High degree of freedom warp to register to target brain

(~2x105 deformation vectors)

Target Brain

This could be stage 2 in fMRI alignment to atlas.

1 brain template (template brain feature outline indicated by red lines).

2

3

4

1 brain template (template brain feature outline indicated by red lines).

2

R’ = MR

For 6 DOFs M includes rotations about x, y, and z axes and translations along x, y, and z.

Total motion is calculated from distance = {(x’-x)2+ (y’-y)2 +(z’-z)2}1/2 calculated for each voxel within the brain.

Figure 13.1. brain template (template brain feature outline indicated by red lines). Motion correction in a finger tapping experiment. Fig 1a (upper figure) shows the motion component ( z translation) exhibiting the strongest stimulus-correlated motion. The experimental “on” periods are shown at the bottom of the figure. Fig 1b (lower figure) shows a set of superior axial images with no correction (top row), after realignment only (middle row) and following full correction( bottom row).

Full correction models signal loss due to movement.

Figure 13.2 brain template (template brain feature outline indicated by red lines).. Motion correction in a visual stimulation experiment. Fig 2a (upper figure) shows the motion component ( z translation) exhibiting the strongest stimulus-correlated motion. The experimental “on” periods are shown at the bottom. Fig 2b (lower figure) shows a set of superior axial images with no correction (top row), after realignment only (middle row) and following full correction( bottom row)

Figure 13.3 brain template (template brain feature outline indicated by red lines).. Two slices from a group activation map ( n=6) in schizophrenic patients before ( upper row) and after (lower row) group correction for subject motion.

visual brain template (template brain feature outline indicated by red lines).

No MC

auditory

visual

MC

auditory

http://www.fmrib.ox.ac.uk/fsl/ brain template (template brain feature outline indicated by red lines).