kpnl:reu_dti_demo
Differences
This shows you the differences between two versions of the page.
Both sides previous revisionPrevious revisionNext revision | Previous revision | ||
kpnl:reu_dti_demo [2016/06/27 15:27] – [Laboratory Report] admin | kpnl:reu_dti_demo [2016/06/27 17:06] (current) – [Visualizing Tracts with TrackVis] admin | ||
---|---|---|---|
Line 3: | Line 3: | ||
====== Diffusion Weighted Imaging ====== | ====== Diffusion Weighted Imaging ====== | ||
- | A white matter fiber connects the output of one neuron to one or more target neurons. There are well known bundles or **tracts** of fibers that run through the brain and spinal cord and connect distant brain regions. Here are a few examples of major tracts: | + | A white matter fiber connects the output of one neuron to one or more target neurons. There are well known bundles or **tracts** of fibers that run through the brain and spinal cord and connect distant brain regions. Here are a few examples of major tracts |
- | * Corticospinal tract. | + | * Corticospinal tract. |
- | * Corpus | + | * See DTI image [[http://www.ajnr.org/content/25/ |
- | * Superior longitudinal fasciculus.It is discussed | + | * Cingulum. A tract that runs along the top of the corpus callsoum connecting frontal, temporal, and parietal regions |
- | * Inferior longitudinal fasciculus. | + | * See DTI image [[http:// |
+ | * Corpus | ||
+ | * See DTI image [[http:// | ||
+ | * Corona Radiata. Widely spreading (fan-like) projections that connect the basal ganglia | ||
+ | * See DTI image [[http://www.ajnr.org/content/25/ | ||
+ | * Superior longitudinal fasciculus. | ||
+ | * See DTI image [[http://www.ajnr.org/ | ||
+ | * Inferior longitudinal fasciculus. Connects the occipital and temporal lobes. | ||
+ | * See DTI image [[http://www.ajnr.org/content/25/ | ||
+ | * Forceps major (aka posterior forceps). Connects the occipital lobes across the splenium of the corpus callosum. | ||
+ | * See DTI image [[http:// | ||
+ | * Forceps minor (aka anterior forceps). Connects the frontal lobes across the genu of the corpus callosum. | ||
+ | * See DTI image [[http:// | ||
- | A review of tracts using DTI can be found [[http:// | + | ===== Tractography ===== |
- | + | ||
- | In this lab we will investigate relatively new methods in MRI to visualize the integrity and direction of white matter tracts. | + | |
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | + | ||
- | ===== Data used in this lab ===== | + | |
- | * | + | |
- | + | ||
- | + | ||
- | ===== Part 1: Introduction to Diffusion and White Matter Tracts ===== | + | |
- | + | ||
- | + | ||
- | + | ||
- | ==== White Matter Tracts ==== | + | |
- | + | ||
- | An important part of the lab will be knowing your way around the white-matter and being able to identity some white-matter tracts, particularly when trying to identity locations of crossing fibers. Let's learn about some tracts now. | + | |
- | + | ||
- | === FSLView Setup === | + | |
- | <WRAP center round tip 80%> | + | |
- | The info below requires knowledge of fslview and using the 'Atlas Tools' | + | |
- | </ | + | |
- | + | ||
- | **1.** Start up '' | + | |
- | + | ||
- | **2.** Click on **Tools** => **Toolbars** => **Atlas tools** | + | |
- | + | ||
- | **3.** Toggle your view layout to have bigger slices by clicking the following button below. | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | You want your view to look this below. | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | **4.** Click on the '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | **5.** Then select the '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | Move this window to the side so you can see all three slices in FSLView while being able to change the selected tract in this atlas. | + | |
- | + | ||
- | + | ||
- | === Learning about White Matter Tracts === | + | |
- | + | ||
- | <WRAP center round important 100%> | + | |
- | <WRAP centeralign> | + | |
- | <wrap em>LAB REPORT Part 1 - # | + | |
- | </ | + | |
- | We will now identify various white matter structures. | + | |
- | * For each of the four white matter tracts indicated below, please take a screenshot and add it along with the tract label to your lab report. | + | |
- | </ | + | |
- | + | ||
- | **6.** Identify each of the following structures one by one. Note we are only going through the left hemisphere tracts here to conserve time. For the following four, make sure to take a screenshot. | + | |
- | + | ||
- | * **Corticospinal tract L**: nerves within the corticospinal tract are involved in movement of muscles of the body. | + | |
- | * **Cingulum (cingulate gyrus) L**: white matter fibers that project from the cingulate gyrus to the entorhinal cortex. | + | |
- | * **Forceps minor**: connects the lateral and medial surfaces of the frontal lobes and crosses the midline via the genu of the corpus callosum. | + | |
- | * **Inferior fronto-occipital fasciculus L**: passes backward from the frontal lobe to the occipital and temporal lobes. Also note the proximity of this tract to the forceps minor and to the corticospinal tract (this might be important later when discussing crossing fibers)! | + | |
- | + | ||
- | Feel free to explore some of the other fiber tracts on your own. Do you notice any other tracts that might either cross or come very close to each other? | + | |
- | + | ||
- | === Crossing Fibers === | + | |
- | + | ||
- | **7.** To get an idea of this possibility (of adjacent or crossing fibers), load several of the white matter tracts you just went through at the same time. You can do so by running the following command below in the terminal. (remember, you can just copy and paste it) | + | |
- | + | ||
- | <code bash> | + | |
- | cd ~/ | + | |
- | + | ||
- | fslview ${FSLDIR}/ | + | |
- | cingulum_L.nii.gz -l ' | + | |
- | corticospinal_tract_L.nii.gz -l ' | + | |
- | forceps_minor.nii.gz -l ' | + | |
- | inferior_fronto-occipital_fasciculus_L.nii.gz -l ' | + | |
- | </ | + | |
- | + | ||
- | We have just opened four tracts primarily in the left hemisphere: | + | |
- | + | ||
- | * **Corticospinal tract** in blue indicating it mainly goes between superior and inferior | + | |
- | * **Cingulum** in green indicating it mainly goes between anterior and posterior | + | |
- | * **Forceps minor** in red indicating it mainly goes between lateral and medial | + | |
- | * **Inferior fronto-occipital fasciculus** in green like the cingulum | + | |
- | + | ||
- | Feel free to adjust the transparency of each of the tracts via the slider at the bottom or to hide/show each of the tracts by double-clicking the eye icon next to the name of the tract. Note that the values for each tract represent the probability that tract was found for a subject in that dataset. Here, we are displaying tracts where a value was found for any of the subjects. | + | |
- | + | ||
- | <WRAP center round important 100%> | + | |
- | <WRAP centeralign> | + | |
- | <wrap em>LAB REPORT Part 1 - # | + | |
- | </ | + | |
- | * Take one screenshot of an area where some of these tracts intersect or come close to each other. | + | |
- | * What issues might you face in trying to identity the different tracts using diffusion weighted imaging? | + | |
- | </ | + | |
- | + | ||
- | <WRAP centeralign> | + | |
- | <wrap hi> | + | |
- | **//Before proceeding you should close all open '' | + | |
- | </ | + | |
- | </ | + | |
- | ===== Part 2: Computing FA and DT images using FSL/FDT ===== | + | |
- | + | ||
- | In this section we will compute FA and DT maps using the [[http:// | + | |
- | + | ||
- | <WRAP center round alert 100%> | + | |
- | Note, the '' | + | |
- | + | ||
- | <wrap em> Use the '' | + | |
- | </ | + | |
- | + | ||
- | Open up a new terminal window and enter in the following: | + | |
- | + | ||
- | <code bash> | + | |
- | mkdir -p ~/ | + | |
- | cp ~/ | + | |
- | cp ~/ | + | |
- | cp ~/ | + | |
- | cd ~/ | + | |
- | </ | + | |
- | + | ||
- | In the same terminal window you used above: | + | |
- | + | ||
- | * Type **'' | + | |
- | * Because you launched fsl within a dti subject folder, this will simplify navigation when selecting the input files. | + | |
- | + | ||
- | As this exercise will focus upon white matter tracts, note that '' | + | |
- | ==== WM Atlas ==== | + | |
- | + | ||
- | In this section, we will continue to use the **JHU White-Matter Tractography Atlas** in '' | + | |
- | + | ||
- | - Start '' | + | |
- | - Open a standard brain (**File** => **Open Standard**), | + | |
- | - Reveal the atlas tools: **Tools** => **Toolbars** => **Atlas Tools** | + | |
- | - Once the atlas toolbar appears click the '' | + | |
- | - Select the **JHU White-Matter Tractography Atlas** and **JHU ICBM DTI-81 White Matter Labels** | + | |
- | - // | + | |
- | - Click '' | + | |
- | + | ||
- | Now in the atlas toolbar, you should see your two new white matter atlases. If you move your cursor around on the brain, information in this toolbar will update with the probability that the selected voxel is a given white matter label or tract. | + | |
- | + | ||
- | To also visualize the probability distributions for a specific tract (as in the previous [[: | + | |
- | + | ||
- | - Select the **Structures** button in the atlas toolbar | + | |
- | - Select the **JHU White-Matter Tractography Atlas** in the dropdown | + | |
- | - Select the check boxes for **Locate selected structure** and **Preview selected structure' | + | |
- | - Now you can select a particular tract in the list and it will display on your standard brain. Feel free to move this structure list to the side so you can view the whole brain. | + | |
- | + | ||
- | <WRAP center round tip 70%> | + | |
- | Leave this FSLView window open as you'll want to refer back to it later in the lab. | + | |
- | </ | + | |
- | + | ||
- | + | ||
- | ==== Compensate for eddy currents ==== | + | |
- | + | ||
- | Applying strong gradients such as those used in diffusion imaging creates compensatory magnetic fields that oppose the gradients. These are called 'eddy currents' | + | |
- | + | ||
- | The problem in MRI is that spatial encoding is exquisitely dependent upon the strength of the magnetic field gradient. The spatial location of an hydrogen nucleus (proton) is encoded by applying a magnetic field of a precise strength to a location in space. If the magnetic field strength is augmented or diminished by eddy currents, there will be a change in the perceived location of the proton. Thus, the image will be sheared or otherwise distorted. | + | |
- | + | ||
- | There are algorithms that attempt to compensate for eddy current distortions, | + | |
- | + | ||
- | <WRAP center round info 80%> | + | |
- | To save time during class, **the eddy current corrections have been calculated beforehand**. Do not re-run the eddy current corrections. The eddy current corrected data is called **33470_LAS_eddy.nii.gz** (remember, .nii.gz means that it is in compressed (gzipped) NIFTI-1 format) | + | |
- | </ | + | |
- | + | ||
- | ==== Strip your skulls! ==== | + | |
- | + | ||
- | <WRAP center round tip 80%> | + | |
- | Need a refresher on using BET? **[[psyc410_s15: | + | |
- | </ | + | |
- | + | ||
- | You will need a brain mask for this exercise. A brain mask is a binary image that contains ones wherever there is brain and zeros everywhere else (e.g., skull, scalp, outside of brain). Such a mask can then be used to remove from computation any region that is not brain. Brain masks can be created with FSL's Brain Extraction Tool (BET). Be sure to **specify your eddy corrected images as input**. | + | |
- | + | ||
- | The '' | + | |
- | + | ||
- | <WRAP center round tip 80%> | + | |
- | You might find using the option for lowering the '' | + | |
- | </ | + | |
- | + | ||
- | Also, because we need a brain mask, you need to check the option (under '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | Before proceeding, make sure you have a good brain mask and a good skull stripped brain. Quality is important. Make a note of their names, as this will need to be provided in the next step. | + | |
- | ==== Compute DTI and FA ==== | + | |
- | + | ||
- | Choose the '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | Now //(note the numbered list below corresponds to the numbers in the image below)//: | + | |
- | + | ||
- | - Check the '' | + | |
- | - Be sure to specify the eddy-current corrected diffusion data '' | + | |
- | - Specify your brain mask from the skull stripping step | + | |
- | - FSL will automatically set the output file | + | |
- | - The gradient directions are provided in the '' | + | |
- | - The b vals are the '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | Press '' | + | |
- | + | ||
- | ==== Examine the DTI and FA results ==== | + | |
- | + | ||
- | You will need '' | + | |
- | + | ||
- | **1.** Examine your FA (fractional anisotropy) image. | + | |
- | + | ||
- | * Click **File** => **Open** | + | |
- | * Open the '' | + | |
- | + | ||
- | <WRAP center round help 80%> | + | |
- | Can you see what tissue types have high FA, and what tissue types have low FA? If the answers are not obvious - make sure to check with me. | + | |
- | </ | + | |
- | + | ||
- | **2.** Examine the DTI images for the 3 orthogonal axes that describe the tensor at each voxel. | + | |
- | + | ||
- | * Click **File** => **Add** | + | |
- | * Add the '' | + | |
- | + | ||
- | <WRAP center round info 80%> | + | |
- | Note that the V1 image ('' | + | |
- | </ | + | |
- | + | ||
- | You can appreciate this by plotting the color coded principal diffusion axes. | + | |
- | + | ||
- | - First, let's hide '' | + | |
- | - Highlight '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | This brings up a box of options. Choose DTI display options and then select '' | + | |
- | + | ||
- | {{ : | + | |
- | + | ||
- | <WRAP center round info 80%> | + | |
- | The lines point in the direction of the principal direction of diffusion. | + | |
- | * RED = Left/ | + | |
- | * BLUE = Superior/ | + | |
- | * GREEN = Anterior/ | + | |
- | </ | + | |
- | + | ||
- | <WRAP center round tip 60%> | + | |
- | To get larger images of a particular orthogonal slice of the brain (axial, sagittal, or coronal) you can switch to single view. | + | |
- | * Select '' | + | |
- | * Close the window with the three views | + | |
- | * To cycle through the three cardinal axes, press the '' | + | |
- | {{ : | + | |
- | </ | + | |
- | + | ||
- | + | ||
- | <WRAP center round important 100%> | + | |
- | <WRAP centeralign> | + | |
- | <wrap em>LAB REPORT Part 2 - # | + | |
- | </ | + | |
- | * Use your new found knowledge of major white matter tracts to identify the structures below in your V1 image. Include labeled screenshots in your report. | + | |
- | * corpus callosum | + | |
- | * corticospinal tract | + | |
- | + | ||
- | //Hint: you can refer back to the JHU white matter atlas in your previous instance of '' | + | |
- | </ | + | |
- | ===== Part 3: Tractography ===== | + | |
==== Visualizing Tracts with TrackVis ==== | ==== Visualizing Tracts with TrackVis ==== | ||
- | <WRAP center round info 80%> | + | <WRAP center round info 90%> |
- | The sample data for this exercise are located in '' | + | The sample data for this exercise are located in '' |
</ | </ | ||
- | In this section we will use [[http:// | + | We will use [[http:// |
**1.** Open '' | **1.** Open '' | ||
- | <WRAP center round tip 80%> | + | <WRAP center round tip 90%> |
If you don't see the icon in your dock, you can open it from your '' | If you don't see the icon in your dock, you can open it from your '' | ||
</ | </ | ||
Line 299: | Line 47: | ||
'' | '' | ||
- | <WRAP center round alert 80%> | + | <WRAP center round alert 90%> |
Before proceeding, watch the [[http:// | Before proceeding, watch the [[http:// | ||
</ | </ | ||
Line 308: | Line 56: | ||
* In the lower-right '' | * In the lower-right '' | ||
- | <WRAP center round todo 60%> | + | {{ : |
- | Add some screenshots of the track and peroperty windows and options. | + | |
- | + | ||
- | Show an example of how changing the length affects which fiber tracts are shown. | + | |
- | </ | + | |
Play with the controls to get a sense of how you can navigate through space, rotate the brain, adjust which fiber tracts are displayed, etc. | Play with the controls to get a sense of how you can navigate through space, rotate the brain, adjust which fiber tracts are displayed, etc. | ||
- | <WRAP center round tip 80%> | + | <WRAP center round tip 90%> |
It'll take a little bit of time playing around before you get the hang of it. Your goal isn't to become an expert user, but you should get a feel for the basic image/track manipulation options. And if things go totally off the rails, you can always close and reopen the program to start fresh. | It'll take a little bit of time playing around before you get the hang of it. Your goal isn't to become an expert user, but you should get a feel for the basic image/track manipulation options. And if things go totally off the rails, you can always close and reopen the program to start fresh. | ||
* Only view the long tracks and then only view the short tracks | * Only view the long tracks and then only view the short tracks | ||
Line 333: | Line 76: | ||
Let's create a ball-shaped region of interest (ROI) to select some tracts. This will help us identify only the tracks that run through areas we're interested in (our ROI) | Let's create a ball-shaped region of interest (ROI) to select some tracts. This will help us identify only the tracks that run through areas we're interested in (our ROI) | ||
- | <WRAP center round tip 80%> | + | <WRAP center round tip 90%> |
To clear your workspace, you might consider turning off, or hiding, the current slice-based view. To do so, right-click on '' | To clear your workspace, you might consider turning off, or hiding, the current slice-based view. To do so, right-click on '' | ||
</ | </ | ||
Line 345: | Line 88: | ||
* Double-click on the '' | * Double-click on the '' | ||
- | You can drag the ball around within your three orthogonal slice windows along the bottom of the screen. | + | {{ : |
- | <WRAP center round todo 60%> | ||
- | Include a sceenshot of what I mean by " | ||
- | </ | ||
+ | You can drag the ball around within your three orthogonal slice windows along the bottom of the screen. | ||
**3.** If your three orthogonal slices are not moving as you move your sphere around, click on the '' | **3.** If your three orthogonal slices are not moving as you move your sphere around, click on the '' | ||
Line 356: | Line 97: | ||
**4** You can add a second ROI sphere by simply repeating steps 1-3. | **4** You can add a second ROI sphere by simply repeating steps 1-3. | ||
- | <WRAP center round tip 80%> | + | <WRAP center round tip 90%> |
You can similarly create a TrackGroup based upon a hand-traced ROI like in the [[http:// | You can similarly create a TrackGroup based upon a hand-traced ROI like in the [[http:// | ||
</ | </ | ||
- | **Find Some Tracts** | + | **Find Some Tracts!!** |
- | + | ||
- | Use these methods, or others you discover, to isolate the following tracts. | + | |
- | + | ||
- | <WRAP center round todo 60%> | + | |
- | Include links and pictures to help them identify these images | + | |
- | </ | + | |
- | * minor forceps | + | Use these methods, or others you discover, to isolate some of the tracts described at the [[http:// |
- | * major forceps | + | |
- | * corticospinal tract | + | |
- | * cingulum | + | |
kpnl/reu_dti_demo.1467059254.txt.gz · Last modified: 2016/06/27 15:27 by admin