Print

Print


Hi Jesper

Perhaps it's a bit different. I acquired multishell DTI and the corresponding opposite blip at different TEs to tests this. The resulting field map with top up changed with every different TE and the unwarped resulting image had slightly different correction in shape and voxel values.

We are working in testing this in different levels, but with this experience, the DTI and opposite blip should have the same TE.

Do you think the differences in this case and your experience are not so important then?

Thanks

Eduardo 

El miércoles, 5 de octubre de 2016, Jesper Andersson <[log in to unmask]> escribió:
Dear Qiuyun (and Eduardo),

I’m trying to figure out how to preprocess a multi-shell diffusion dataset. For each shell, the diffusion weighted images (A>>P) and the b0 images (P>>A) were collected with the same echo time. Between shells, the echo time was not fixed (64, 76, 83, 94 ms respectively for each of the 4 shells).  But for all shells, the total EPI readout time was the same.
 
There’ll be different T2 weighting between different shells.  Is it possible to perform a joint eddy current correction with all shells?  If so, what’s the best way to carry out the TOPUP step, and is it possible to use the results as eddy input?

as long as there has been no re-shimming between the acquisition of the different shells, my experience is that this works well. I have done this many times on data acquired by the HCP during the piloting phase. The echo time doesn’t enter in to the models used by either topup or eddy.

Eduardo, it would be interesting to hear your experience and what it is that fails when trying to use all shells in the correction?

Qiuyun, there is no principal difference between this case and the case we describe on the EDDY help page. Run topup with your P->A b0 image and an A->P b0 image from one of the shells. I suggest using the first b0 from the shell that was acquired closest in time to the P->A b0. Make sure that the first volume in the 4D file you feed into topup is the A->P b0 image. Concatenate all the shells, doing this for data, bvals and bvecs alike. Make sure that the shell you took the b0 from for topup is the first shell in your concatenated data set. Then just run eddy in the usual way.

But please note that for your next step, the modelling of the diffusion signal, you need to address the problem of differing echo time for the different shells. Maybe Stam or Saad has some suggestions for how to do that?

Jesper


 
Any suggestions will be highly appreciated.  Thanks a lot in advance!
 
Best,
Qiuyun
 

The information in this e-mail is intended only for the person to whom it is
addressed. If you believe this e-mail was sent to you in error and the e-mail
contains patient information, please contact the Partners Compliance HelpLine at
http://www.partners.org/complianceline . If the e-mail was sent to you in error
but does not contain patient information, please contact the sender and properly
dispose of the e-mail.




--
Eduardo A. Garza-Villarreal, M.D., Ph.D.
  • CONACYT Research Fellow, Investigaciones Clínicas, Instituto Nacional de Psiquiatría, Mexico.
  • Visiting Senior Researcher, Center of Functionally Integrative Neuroscience, University of Aarhus, Denmark.