2009 WinawerDoughertyWandell: Difference between revisions

From Psych 221 Image Systems Engineering
Jump to navigation Jump to search
imported>Winawer
No edit summary
imported>Winawer
No edit summary
Line 104: Line 104:


:<math>x_n = \frac{1}{N} \sum_{k=0}^{N-1} X_k e^{\frac{2\pi i}{N} k n} \quad \quad n = 0,\dots,N-1.</math>
:<math>x_n = \frac{1}{N} \sum_{k=0}^{N-1} X_k e^{\frac{2\pi i}{N} k n} \quad \quad n = 0,\dots,N-1.</math>
A simple description of these equations is that the complex numbers <math>X_k</math> represent the amplitude and phase of the different sinusoidal components of the input "signal" <math>x_n</math>.  The DFT computes the <math>X_k</math> from the <math>x_n</math>, while the IDFT shows how to compute the <math>x_n</math> as a sum of sinusoidal components <math>(1/N) X_k e^{\frac{2\pi i}{N}k n}</math> with [[frequency]] <math>k/N</math> cycles per sample.  By writing the equations in this form, we are making extensive use of [[Euler's formula]] to express sinusoids in terms of complex exponentials, which are much easier to manipulate. In the same way, by writing <math>X_k</math> in [[Complex_number#Polar_form|polar form]], we immediately obtain the sinusoid amplitude <math>A_k</math> and phase <math>\phi_k</math> from the [[complex argument|complex modulus and argument]] of <math>X_k</math>, respectively:
:<math>A_k = |X_k| = \sqrt{\operatorname{Re}(X_k)^2 + \operatorname{Im}(X_k)^2},</math>
:<math>\varphi_k = \arg(X_k) = \operatorname{atan2}\big( \operatorname{Im}(X_k), \operatorname{Re}(X_k) \big).</math>
Note that the normalization factor multiplying the DFT and IDFT (here 1 and 1/''N'') and the signs of the exponents are merely conventions, and differ in some treatments. The only requirements of these conventions are that the DFT and IDFT have opposite-sign exponents and that the product of their normalization factors be 1/''N''.  A normalization of <math>1/\sqrt{N}</math> for both the DFT and IDFT makes the transforms [[unitary matrix|unitary]], which has some theoretical advantages, but it is often more practical in numerical computation to perform the scaling all at once as above (and a unit scaling can be convenient in other ways).
(The convention of a negative sign in the exponent is often convenient because it means that <math>X_k</math> is the amplitude of a "positive frequency" <math>2\pi k/N</math>. Equivalently, the DFT is often thought of as a [[matched filter]]: when looking for a frequency of +1, one correlates the incoming signal with a frequency of −1.)
In the following discussion the terms "sequence" and "vector" will be considered interchangeable.

Revision as of 17:06, 8 December 2009

Back to Psych 204 Projects 2009

Project Title - Retinotopic maps in MNI space

Your project overview goes here. For example: Much of the visual cortex is organized into visual field maps: nearby neurons have receptive fields at nearby locations in the image. These maps are usually identified in individual subjects. The precise location of each map may be different in different brains. For this project, we asked how the quality of the maps would compare using (a) standard retinoptic methods on individual brains or (b) group-averaged brains projected into MNI space.
Note that this is a project template. Other styles are possible. For example, you could use a multiple page format.

Background

Retinotopic maps

You can use subsections if you like. Below is an example of a retinotopic map. Or, to be precise, below will be an example of a retinotopic map once the image is uploaded. To add an image, simply put text like this inside double brackets 'MyFile.jpg | My figure caption'. When you save this text and click on the link, the wiki will ask you for the figure.
Figure 1

Below is another example of a reinotopic map in a different subject.
Figure 2

Once you upload the images, they look like this. Note that you can control many features of the images, like whether to show a thumbnail, and the display resolution.

Figure 3


MNI space

MNI is an abbreviation for Montreal Neurological Institute.

Methods

Measuring retinotopic maps

Retinotopic maps were obtained in 5 subjects using Population Receptive Field mapping methods Dumoulin and Wandell (2008). These data were collected for another research project in the Wandell lab. We re-analyzed the data for this project, as described below.

Subjects

Subjects were 5 healthy volunteers.

MR acquisition

Data were obtained on a GE scanner. Et cetera.

MR Analysis

The MR data was analyzed using mrVista software tools.

Pre-processing

All data were slice-time corrected, motion corrected, and repeated scans were averaged together to create a single average scan for each subject. Et cetera.

PRF model fits

PRF models were fit with a 2-gaussian model.

MNI space

After a pRF model was solved for each subject, the model was trasnformed into MNI template space. This was done by first aligning the high resolution t1-weighted anatomical scan from each subject to an MNI template. Since the pRF model was coregistered to the t1-anatomical scan, the same alignment matrix could then be applied to the pRF model.
Once each pRF model was aligned to MNI space, 4 model parameters - x, y, sigma, and r^2 - were averaged across each of the 6 subjects in each voxel.

Et cetera.


Results - What you found

Retinotopic models in native space

Some text. Some analysis. Some figures.

Retinotopic models in individual subjects transformed into MNI space

Some text. Some analysis. Some figures.

Retinotopic models in group-averaged data on the MNI template brain

Some text. Some analysis. Some figures.


Retinotopic models in group-averaged data projected back into native space

Some text. Some analysis. Some figures.

Conclusions

Here is where you say what your results mean.

References - Resources and related work

References

Software

Appendix I - Code and Data

Code

zip file with code

Data

zip file with my data

Appendix II - Work partition (if a group project)

Brian and Bob gave the lectures. Jon mucked around on the wiki.


Test Equations

This is a test of equation use on our wiki. The text below is copied and pasted from wikipedia's article on the DFT

Definition

The sequence of N complex numbers x0, ..., xN−1 is transformed into the sequence of N complex numbers X0, ..., XN−1 by the DFT according to the formula:

where i is the imaginary unit and is a primitive N'th root of unity. (This expression can also be written in terms of a DFT matrix; when scaled appropriately it becomes a unitary matrix and the Xk can thus be viewed as coefficients of x in an orthonormal basis.)

The transform is sometimes denoted by the symbol , as in or or .

The inverse discrete Fourier transform (IDFT) is given by