Differences between revisions 4 and 5
Deletions are marked like this. Additions are marked like this.
Line 10: Line 10:
freeview -v wm1_edits_before/mri/wm.mgz:colormap=heat:opacity=0.4 wm1_edits_before/mri/brainmask.mgz \
-f wm1_edits_before/surf/lh.white:edgecolor=blue wm1_edits_before/surf/lh.pial:edgecolor=red \
wm1_edits_before/surf/rh.white:edgecolor=blue wm1_edits_before/surf/rh.pial:edgecolor=red \
freeview -v wm1_edits_before/mri/brainmask.mgz \
wm1_edits_before/mri/wm.mgz:colormap=heat:opacity=0.4 \
-f wm1_edits_before/surf/lh.white:edgecolor=blue \
wm1_edits_before/surf/lh.pial:edgecolor=red \
wm1_edits_before/surf/rh.white:edgecolor=blue \
wm1_edits_before/surf/rh.pial:edgecolor=red \
Line 35: Line 38:
To begin editing voxels, make '''sure''' the wm.mgz is highlighted in the volume menu, and on top on the brainmask.mgz. click on the edit voxels option ( {{attachment:voxel_edit.gif}} ) and enter 255 into the brush value box. your menu should look like this:<<BR>>
{{attachment:menu1.jpeg}}<<BR>>
'''Remember:''' the wm.mgz must be on top and and the brush value 255 in order to edit voxels in the wm.mgz, if it is not visible on top you will be editing the brainmask.mgz. You should only be editing the brainmask for pial edits and not wm edits.
In the Voxel edit menu that pops up, click on the first icon "freehand" and select a Brush size of 1 or 2. Make sure the Reference is the wm volume. This menu should look like this:<<BR>>
{{attachment:menu2.jpeg}}<<BR>>
To begin editing voxels, make '''sure''' the wm.mgz is '''highlighted''' in the volume menu, and on top on the brainmask.mgz. click on the edit voxels option ( {{attachment:voxel_edit.gif}} ) and enter a brush size of either 1 or 2 and click on the Recon editing checkbox (setting the brush value to 255 and the Eraser value to 1). In order to be ready for editing, both the Volumes tab and Voxel Edit window should look exactly like this:<<BR>>
{{attachment:wmeditsview.jpeg}}<<BR>>
'''Remember:''' the wm.mgz must be on top and and the 'Recon editing' box checked (brush value set to 255) in order to edit voxels in the wm.mgz, if it is not visible on top you will be editing the brainmask.mgz. You should only be editing the brainmask for pial edits and not wm edits.
<<BR>>
Line 45: Line 47:
freeview -v wm1_edits_after/mri/wm.mgz:colormap=heat:opacity=0.4 wm1_edits_after/mri/brainmask.mgz \
-f wm1_edits_after/surf/lh.white:edgecolor=blue wm1_edits_after/surf/lh.pial:edgecolor=red \
wm1_edits_after/surf/rh.white:edgecolor=blue wm1_edits_after/surf/rh.pial:edgecolor=red
freeview -v wm1_edits_after/mri/brainmask.mgz \
wm1_edits_after/mri/wm.mgz:colormap=heat:opacity=0.4 \
-f wm1_edits_after/surf/lh.white:edgecolor=blue \
wm1_edits_after/surf/lh.pial:edgecolor=red \
wm1_edits_after/surf/rh.white:edgecolor=blue \
wm1_edits_after/surf/rh.pial:edgecolor=red \

top | previous

Making Edits to the White Matter

  • To follow this exercise exactly be sure you've downloaded the tutorial data set before you begin. If you choose not to download the data set you can follow these instructions on your own data, but you will have to substitute your own specific paths and subject names.

Sometimes the white matter is not segmented correctly: sometimes voxels that should be white matter are excluded, and other times voxels that should not be white matter are included in error. Either of these occurences can be fixed with simple manual edits. Below you will find examples of a few of the common problems and how to fix them. For some people it is easiest to see these problems when looking at the inflated surface in tksurfer. Others can pick them out while viewing the volumes in tkmedit. Either method is fine for identifying these problems. However they can only be fixed using tkmedit with the wm.mgz volume.

Geometric inaccuracy due to brain lesion

wm1_edits_before is an example of white matter being excluded from the wm.mgz volume due to the presence of a brain lesion. This page will walk you through fixing this specific example as well as other common problems with the wm.mgz volume. First, make sure you have wm1_edits_before loaded in tkmedit:


freeview -v wm1_edits_before/mri/brainmask.mgz \
wm1_edits_before/mri/wm.mgz:colormap=heat:opacity=0.4 \
-f wm1_edits_before/surf/lh.white:edgecolor=blue \
wm1_edits_before/surf/lh.pial:edgecolor=red \
wm1_edits_before/surf/rh.white:edgecolor=blue \
wm1_edits_before/surf/rh.pial:edgecolor=red \


Again, this will bring up the brainmask.mgz volume, the wm.mgz volume, and the surfaces for both hemispheres. Click in the main viewing area and use Alt-C or the arrows arrows.jpeg to flip between the wm.mgz and brainmask.mgz on top of each other. with the wm.mgz loaded on top of the brainmask.mgz volume, it should look like this:
view.gif.jpeg
If you look at coronal slice 160 you will notice an area on the right hemisphere where the white surface (yellow line) does not follow the surface of the brain, but in fact cuts into it. This geometric inaccuracy is caused by a lesion where white matter has been marked as non-white matter.
wm_defect_vol.gif.png
You can see this as a dimple or hole on the inflated surface in tksurfer too (see next image). If not already open, in a second terminal window, open the rh surface in tksurfer and you can see the dimple or hole on the medial side:


tksurfer wm1_edits_before rh inflated



lesiontksurfer.jpg
With the wm.mgz volume loaded into tkmedit you can see that this area has been left out of the wm volume completely.
wm_defect.gif.png

To fix this problem you will need to fill in the missing voxels in the wm.mgz volume. First, it's a good idea to use the page up and down keys to scroll through the individual slices in the volume until you get a good idea of where the problem starts and ends. You will want to start filling in voxels when the inaccuracy appears, and keep filling them in slice by slice until the problem is no longer visible. The coronal view and a brush radius of one or two are good settings for painting in voxels. Zoom In and Out either by using the scroll button on your mouse or holding down shift and the up or down keys. To begin editing voxels, make sure the wm.mgz is highlighted in the volume menu, and on top on the brainmask.mgz. click on the edit voxels option ( voxel_edit.gif ) and enter a brush size of either 1 or 2 and click on the Recon editing checkbox (setting the brush value to 255 and the Eraser value to 1). In order to be ready for editing, both the Volumes tab and Voxel Edit window should look exactly like this:
wmeditsview.jpeg
Remember: the wm.mgz must be on top and and the 'Recon editing' box checked (brush value set to 255) in order to edit voxels in the wm.mgz, if it is not visible on top you will be editing the brainmask.mgz. You should only be editing the brainmask for pial edits and not wm edits.
You are now ready to edit. click back on the volume, and using your left mouse button fill in the wm hole with voxels. If you fill in too many, you can press shift while using your left mouse button to erase voxels. You can keep flipping back between the brainmask and the wm volume using your Alt-C combination if you want to check that what you are filling in is indeed white matter. At any time you can use the CTRL-Z combination to undo as many actions as you like. When you're done editing the volume should look like this. wm_defect_filled.jpeg
While filling the hole in slice 160, you will notice that the hole is not completely enclosed. To see where the boundary should be, switch back and forth between the brainmask.mgz and wm.mgz volume with Alt-C. Use the brainmask.mgz volume as a guide, and make the boundary of the voxels that you paint in correspond as closely as possible to the boundaries suggested in the brain volume. Go to the next slice, and fill in this slice in the same way. You will need to continue filling in slices in this fashion until the region is completely filled in. If you would like to compare your edits to ones that have been done for you you can open a second freeview window with the wm1_edits_after wm.mgz volume, and for each slice check your results with the corresponding slice in the wm1_edits_after volume. At any time, you can save the changes you've made by clicking on the wm.mgz, making sure it is highlighted in the menu to the left side of the viewing area and click the save volume button save_volume.gif.jpeg or save in the file menu. You can verify your results by viewing the changes made to the wm1_edits_after wm.mgz volume To view this corrected subject (which was corrected by following the detailed instructions provided to you):

freeview -v wm1_edits_after/mri/brainmask.mgz \
wm1_edits_after/mri/wm.mgz:colormap=heat:opacity=0.4 \ 
-f wm1_edits_after/surf/lh.white:edgecolor=blue \
wm1_edits_after/surf/lh.pial:edgecolor=red \
wm1_edits_after/surf/rh.white:edgecolor=blue \
wm1_edits_after/surf/rh.pial:edgecolor=red \

Tips:

  • Sometimes regions may appear enclosed in some slices (i.e. appear as holes), open in subsequent slices (i.e. no longer appear as holes), then enclosed again as you scroll forward. The rule of thumb when editing these is to keep filling until you reach the slice where they finally open up and are no longer enclosed.
  • When editing geometric inaccuracies in the future, you may find it easier to edit the first slice and last slice first, thereby 'capping' the slices that need filling. Then you just need to fill all the slices between the endcaps.
  • Sometimes there is variation between a subject's left and right hemispheres, so that in a particular slice one hemisphere's region will finally 'open up', but the other hemisphere's region is still enclosed. In such cases, continue to fill the enclosed region only, even though the other is open. This will address the topological problem.
  • Problems like the one shown here, due to brain lesions, should always be fixed with edits to the wm.mgz volume. Do not use control points to try and automatically adjust the intensity in these areas.

After you have saved all of your edits, you could recreate the final surfaces with the command:

recon-all -autorecon2-wm -subjid wm1_edits_before

Do not run this command if you are conducting the tutorial! This step will take a long time and there is no need to run it for the tutorial purposes.

FsTutorial/WhiteMatterEdits_freeview (last edited 2018-01-12 12:48:45 by EmmaBoyd)