This shows you the differences between two versions of the page.
Both sides previous revision Previous revision Next revision | Previous revision Next revision Both sides next revision | ||
biac:analysis:topup_correction [2019/02/22 14:31] cmp12 [Run dwipreproc -rpe_all] |
biac:analysis:topup_correction [2019/09/04 17:45] cmp12 |
||
---|---|---|---|
Line 7: | Line 7: | ||
There will be two series of single timepoint EPI images. | There will be two series of single timepoint EPI images. | ||
- | You need to calculate readout time in seconds ( the physical time it takes to get the acquisition matrix of a single slice ) and get the polarity direction ( phase encode direction ). | + | You need to calculate readout time in seconds |
the readout time in seconds for the parameter file will be: | the readout time in seconds for the parameter file will be: | ||
Line 15: | Line 15: | ||
</ | </ | ||
+ | unfortunately at this point the polarity of the images will have to be determined from visual inspection. we aren't provided enough information in the metadata to give an entry into the BXH file ( yet ). | ||
+ | |||
+ | Here is a rough guide to help with inspection. | ||
+ | {{: | ||
- | the polarity for the entry will have to be determined from the **seriesdescription**, | ||
reverse will be " | reverse will be " | ||
Line 62: | Line 65: | ||
applytopup --imain=../ | applytopup --imain=../ | ||
applytopup --imain=../ | applytopup --imain=../ | ||
+ | </ | ||
+ | |||
+ | If you have the scenario where the functional data was acquired with acceleration, | ||
+ | |||
+ | < | ||
+ | echospacing in the BXH header is in microseconds | ||
</ | </ | ||
Line 192: | Line 201: | ||
there are 2 scenarios that typically apply here | there are 2 scenarios that typically apply here | ||
- | 1) a short acquisition with RPE B0s | + | 1) DWI and a short acquisition with RPE B0s |
- | 2) an entire | + | 2) 2 entire |
In both scenarios it is important to create your data with normal phase encoding direction first, followed by reversed. | In both scenarios it is important to create your data with normal phase encoding direction first, followed by reversed. | ||
Line 264: | Line 273: | ||
< | < | ||
#concat the 2 series ( regular, reversed ) | #concat the 2 series ( regular, reversed ) | ||
- | bxh_concat | + | bxh_concat |
+ | 7_LAS.bxh bia6_00260_008_LAS.bxh both | ||
#extract the gradients | #extract the gradients | ||
- | extractdiffdirs --fsl both bvecs bvals | + | extractdiffdirs --fsl both.bxh bvecs bvals |
#convert to mif | #convert to mif | ||
Line 289: | Line 299: | ||
dwipreproc dwi.mif dwi_corr.mif -rpe_all -pe_dir AP -readout_time 0.10656 -debug | dwipreproc dwi.mif dwi_corr.mif -rpe_all -pe_dir AP -readout_time 0.10656 -debug | ||
- | -rpe_all signals that you've replicated A | + | -rpe_all signals that you've replicated |
+ | A | ||
LL the directions with a rpe acquisition | LL the directions with a rpe acquisition | ||
-pe_dir is the phase encode direction of your regular acquisition | -pe_dir is the phase encode direction of your regular acquisition |