|Deletions are marked like this.||Additions are marked like this.|
|Line 7:||Line 7:|
|Note: there are also other problems with this subject, but for the purpose of this tutorial, we will focus only on the skull strip error.|
Fixing a bad skull strip
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.
As you can see, the skull strip for this subject is very poor for both timepoints MR1 and MR2, therefore affecting the results of the base and the longitudinals downstream. Note: there are also other problems with this subject, but for the purpose of this tutorial, we will focus only on the skull strip error.
The best way to fix a problem like this is to adjust the watershed parameters. This can only be done on the cross-sectionals, and you'll have to use the edited cross-sectionals to recreate the base and the longitudinals afterwards.
First, adjust the watershed parameters for both MR1 and MR2 using the following commands:
recon-all -skullstrip -wsthresh 20 -clean-bm -no-wsgcaatlas -subjid OAS2_0004_MR1 recon-all -skullstrip -wsthresh 20 -clean-bm -no-wsgcaatlas -subjid OAS2_0004_MR2
Take a look at your output volume and verify the result of the new skull stripping is correct.
tkmedit OAS2_0004_MR1 brainmask.mgz -aux T1.mgz -surfs tkmedit OAS2_0004_MR2 brainmask.mgz -aux T1.mgz -surfs
After you are satisfied with the results, remember to run the rest of the cross sectional recon-all steps for all timepoints that you made an adjustment to:
recon-all -autorecon2 -autorecon3 -subjid OAS2_0004_MR1 recon-all -autorecon2 -autorecon3 -subjid OAS2_0004_MR2
To save time, you can check in with the recon-all-status.log in the scripts directory of these subjects to see if the normalization step is complete. Once that step is complete, you can start recreating your base.
recon-all -base OAS2_0004 -tp OAS2_0004_MR1 -tp OAS2_0004_MR2 -all
Then, the final step is to recreate the longitudinals once base and cross sectional runs are complete.
recon-all -long OAS2_0004_MR1 OAS2_0004 -all recon-all -long OAS2_0004_MR2 OAS2_0004 -all
**NOTE: When recreating the base and longitudinals, it is recommended to remove or rename the existing base and longitudinal directories and start those processes from scratch.
You can examine the corrected version of the base and longitudinals by using the following commands in separate terminals.
tkmedit OAS2_0004 brainmask.mgz -aux T1.mgz -surfs tkmedit OAS2_0004_MR1.long.OAS2_0004_fixed brainmask.mgz -aux T1.mgz -surfs tkmedit OAS2_0004_MR2.long.OAS2_0004_fixed brainmask.mgz -aux T1.mgz -surfs
For more details on skull strip errors and different ways to fix them, take a look here.