- This topic has 7 replies, 5 voices, and was last updated 3 years, 2 months ago by
mblesac.
-
AuthorPosts
-
-
August 28, 2019 at 1:43 am #1877
jblomm6
ParticipantDear Brainsuite experts,
I am currently using your toolbox for distortion correction (described in Bhushan et al. 2015) on a multishell diffusion dataset.
However, I am not interested in calculating the DTI model, as I will be using CSD tractography. Is it possible to tell the bdp command to only apply distortion correction? It would save me a lot of computational time and memory.
Thanks in advance,
Jeroen
-
August 29, 2019 at 10:35 am #1878
sychoi
ParticipantThank you for using BrainSuite.
The short answer is no.
If a diffusion modeling flag is not indicated then BDP will calculate tensors by default. -
December 16, 2019 at 8:49 am #2017
miguel_rivas
ParticipantDear BrainSuite members,
We are dealing with b0 distortion correction in diffusion and bold EPI images. Unfortunately we neither have the FIELDMAP acquisition nor the EPI with opposite phase encoding direction (anterior-posterior, posterior-anterior). In this context, we would like to use the toolbox for distortion correction (described in Bhushan et al. 2015) in order to correct the B0 in EPI images. However we would like to know two questions about the toolbox:
1- Is there any incompatibility if I use the EPI corrected images (output of the BDP) in other softwares?
2- Is possible to use the toolbox in bold EPI images? We are conducting a fMRI study and I would like to know if it is possible to correct this artifact in Bold EPI images with this toolbox.
Thanks in advance.
Best regards,
-
December 16, 2019 at 2:12 pm #2018
sychoi
ParticipantHey Miguel,
1. Generally you should be fine to port it into other softwares. However, we can’t claim complete knowledge or responsibility over what is being done by other tools outside of BrainSuite so you will have to use your best judgement to figure out the right steps to take for your study.
2. It may be possible to use our registration-based distortion correction on BOLD images as the method is based on the physics of EPI distortion. However, applying this method on BOLD images has not been properly tested or examined by BrainSuite so we can’t make any guarantees on the performance. It’s absolutely worth giving it a try. The full description of the methods can be found here: Bhushan, Chitresh, et al. “Co-registration and distortion correction of diffusion and anatomical images based on inverse contrast normalization.” Neuroimage 115 (2015): 269-280.
You can find the source code for BDP here: https://forums.brainsuite.org/download/
or you can try to hack the command line tool by inputting dummy bvec and bval files, just making sure you have the right number of values to match the number of volumes in your BOLD images. Also at least one of the bvalues should equal to 0 since it the b=0 image is used as the reference image for the registration template.If you give it a try, please let me know how it goes by replying on this thread. Very interested!
Good luck!
-
December 17, 2019 at 12:09 am #2019
miguel_rivas
ParticipantDear sychoi,
Thank your for your reply, I tried to run the bdp.sh script using as input the T1 (skullstripped and bias field corrected) and a 4D image with the 246 bold volumes following (I think) your recommendations. As you said, for the bvals and bvecs files I created:
bvals file: matrix with 3 rows and 246 columns
bvecs file: matrix with 1 row and 246 columnsWhen I ran the script everything worked fine but stopped with the following error:
rmnlab1:bdp neucoga$ bash bdp.sh /Applications/BrainSuite19a/bdp/sujetos/T1_image.bfc.nii.gz –tensor –nii /Applications/BrainSuite19a/bdp/sujetos/fmri.nii.gz -g /Applications/BrainSuite19a/bdp/sujetos/bvecs -b /Applications/BrainSuite19a/bdp/sujetos/bvals
BDP Version: 19a (build #0074), released 2019-02-14
================================================================================
Setting up dataset and inputs
================================================================================
Reading input flags…Checking input files…
BDP could not find any .mask.nii.gz file. BDP will use input bfc file itself as
brain mask for registration and statistics. All voxels with intensity>0 in bfc
will be treated as voxels in brain. You can specify a custom brain mask by using
flag –t1-mask <maskfile_name>. The custom mask must overlay correctly with
input BFC image in BrainSuite.Processing data with fileprefix:
/Applications/BrainSuite19a/bdp/sujetos/T1_image================================================================================
Co-registration and Distortion Correction
================================================================================
Reading the input parameters for co-registration…
Total physical memory found: 64.00GBChecking orientation information…Done
Extracting 0-diffusion (b=0) image from input DWIs…DWI mask is not defined in input flags. BDP will try to estimate a (pseudo) mask
from 0-diffusion (b=0) image. Automatic mask estimation may not be accurate in
some sitations and can affect overall quality of co-registration. In case
co-registration is not accurate, you can define a DWI mask by using flag
–dwi-mask <mask_filename>. The mask can be generated and hand edited in
BrainSuite interface. This mask would be used only for registration purposes
(and not for statistics computation).
Saved (pseudo) masks: /Applications/BrainSuite19a/bdp/sujetos/T1_image.dwi.RAS.mask.nii.gzStarting Registration based distortion Correction…
Loading data…Starting rigid registration of input images…
Reading input data…
Setting/generating masks…
Matching centroids (approx. align)…Warning: Matrix is singular to working precision.
> In register_files_affine (line 304)
In EPI_correct_files_registration_INVERSION (line 109)
In coregister_diffusion_mprage_pipeline (line 219)
In BrainSuite_Diffusion_pipeline (line 26)
Warning: Matrix is singular to working precision.
> In register_files_affine (line 305)
In EPI_correct_files_registration_INVERSION (line 109)
In coregister_diffusion_mprage_pipeline (line 219)
In BrainSuite_Diffusion_pipeline (line 26)Matching resolution of the images…Error using eig
Input to EIG must not contain NaN or Inf.Error in get_original_grid_data (line 44)
Error in register_files_affine (line 313)
Error in EPI_correct_files_registration_INVERSION (line 109)
Error in coregister_diffusion_mprage_pipeline (line 219)
Error in BrainSuite_Diffusion_pipeline (line 26)
MATLAB:eig:matrixWithNaNInf
What could be the error? I have installed the Matlab 2015b runtime…
Thanks for your help.
Best regards,
-
December 17, 2019 at 2:45 am #2021
miguel_rivas
ParticipantDear sychoi,
I have successfully run the distortion correction on the BOLD images. The error was that the bvals and bvecs matrices (bvals: 3×246 and bvecs: 1×246) had to have only a zero in the first column and different values to 0 in the rest. However, at this moment I have two important questions:
1- Will it have any effect on the image if there are values other than zero in the rest of the columns? I used the value 1000. I think that this will not be a problem because b=0 will be used as the reference image. That’s right?
2. I would like to do the fMRI analysis in SPM12. In which preprocessing step should I include the EPI distortion correction with brainsuite? Maybe after realignment for motion correction? Pleas any insight about this would be very appreciated. This would be my pipeline
– Slice timing
– Realign: Estimate and Reslice
– Coregister: Estimate
– Segment
– Dartel: create templates
– Dartel: Normalise to MNI structural scans.
– Dartel: Normalise to MNI functional scans. -
December 19, 2019 at 11:02 am #2032
Breogan
ParticipantDear Brain Suit experts
I’m a beginner in DWI data preprocessing. I am interested in performing Registration-based distortion correction. I managed to perform this step correctly but would like to know if the Bvecs for the resulting image are altered, if so, how could I know the new values associated with the nifti output.
Thanks in advance.
best regards
-
July 10, 2021 at 6:56 am #2878
mblesac
ParticipantDear experts,
It is my first time using BrainSuite. I would like to use the Registration-based distortion correction in my neonatal data (T1w + dMRI). I believe I would have to input everything manually (especially the masks, as the standard methods don’t work with this data). Could you provide me some advice about how to launc the script and the correct flags to use? Thanks in advance.
Best regards,
Manuel
-
-
AuthorPosts
- You must be logged in to reply to this topic.