User Tools

Site Tools


biac:analysis:topup_correction

Differences

This shows you the differences between two versions of the page.

Link to this comparison view

Both sides previous revision Previous revision
Next revision
Previous revision
Next revision Both sides next revision
biac:analysis:topup_correction [2021/08/24 18:22]
cmp12
biac:analysis:topup_correction [2021/08/24 18:26]
cmp12
Line 8: Line 8:
  
 You need to calculate readout time in seconds of the PEPOLAR images ( 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 of the PEPOLAR images ( the physical time it takes to get the acquisition matrix of a single slice ) and get the polarity direction ( phase encode direction ).
 +
  
 the readout time in seconds for the parameter file will be: the readout time in seconds for the parameter file will be:
-<code> 
 **For singleband images :** **For singleband images :**
 +<code>
 readout = (echospacing * (acquisitionmatrix[0] * (percentsampling/100))) / 1e6  readout = (echospacing * (acquisitionmatrix[0] * (percentsampling/100))) / 1e6 
  
 echospacing in the BXH header is in microseconds echospacing in the BXH header is in microseconds
 +</code>
  
  
 **For mutli-band images the readout calculation is more complex:**  **For mutli-band images the readout calculation is more complex:** 
 +<code>
 readout = ( ( ceil ((1/Round_factor) * AcquisitionMatrixPE / Asset_R_factor ) * Round_factor) - 1 ) * EchoSpacing * 0.000001 readout = ( ( ceil ((1/Round_factor) * AcquisitionMatrixPE / Asset_R_factor ) * Round_factor) - 1 ) * EchoSpacing * 0.000001
  
Line 26: Line 29:
 Round_factor = 4 if partial fourier ( PFF is in scanoptions ), 2 if full fourier.  most all scans will be 4 Round_factor = 4 if partial fourier ( PFF is in scanoptions ), 2 if full fourier.  most all scans will be 4
 Asset_R_factor = the reciprocal of the first value of dcm tag (0043,1083) Asset_R_factor = the reciprocal of the first value of dcm tag (0043,1083)
- 
 </code> </code>
  
Line 111: Line 113:
 ==== Create your acq_params.txt file === ==== Create your acq_params.txt file ===
  
-the readout time in seconds for the parameter file will be: +calculate your readout time from above
-<code>readout = (echospacing * (acquisitionmatrix[0] * (percentsampling/100))) / 1e6  +
- +
-echospacing in BXH header is in microseconds +
-</code>+
  
 the regular B0s will get the "1" and the reverse phase B0s will get the "-1" the regular B0s will get the "1" and the reverse phase B0s will get the "-1"
Line 206: Line 204:
 If you are running dwidenoise, do it BEFORE dwipreproc.  If you are running dwidenoise, do it BEFORE dwipreproc. 
  
-==== Calculate your readout time ==== +==== Calculate your readout time as mentioned above ====
-<code>readout = (echospacing * (acquisitionmatrix[0] * (percentsampling/100))) / 1e6  +
- +
-echospacing in BXH header is in microseconds +
-</code>+
  
 ==== Prepare your input datasets ==== ==== Prepare your input datasets ====
biac/analysis/topup_correction.txt · Last modified: 2023/02/23 18:43 (external edit)