Skip to content
New issue

Have a question about this project? Sign up for a free GitHub account to open an issue and contact its maintainers and the community.

By clicking “Sign up for GitHub”, you agree to our terms of service and privacy statement. We’ll occasionally send you account related emails.

Already on GitHub? Sign in to your account

Convert the flip angle maps to B1+ efficiency maps [nT/V] #18

Closed
evaalonsoortiz opened this issue Mar 9, 2024 · 4 comments
Closed

Convert the flip angle maps to B1+ efficiency maps [nT/V] #18

evaalonsoortiz opened this issue Mar 9, 2024 · 4 comments

Comments

@evaalonsoortiz
Copy link
Contributor

The notebook currently uses the same code that is used in the RF shimming project: https://github.com/shimming-toolbox/rf-shimming-7t/blob/main/content/index.ipynb to convert flip angle maps [degrees] to B1+ efficiency maps [nT/V].

The code code contains the following line:
requested_fa = 90 # saturation flip angle -- hard-coded in sequence

Could someone @jcohenadad @dpapp86 confirm that the 90degree flip angle is not hard-coded and that it should be whatever flip angle that was used for the acquisition? I am unsure if this 90 degree flip angle pertains to a pulse that is run before the actual image acquisition itself. @dpapp86 I seem to recall you mentioning something like that.

The turboFLASH scans acquired at VE sites all have a 90 degree flip angle. However, those acquired at the VB17 site have different flip angles (which can be found in the .json). Should our code be retrieving the flip angle recorded in the .json file at all times?

Moreover, DREAM flip angle maps are not acquired with a 90 degree flip angle at VE or VB17 sites. Does this flip-angle-to-B1+-efficiency code still apply for DREAM maps? And if so, should the flip angle be adjusted accordingly?

@dpapp86
Copy link
Contributor

dpapp86 commented Mar 10, 2024

Hello!

Wrt the saturation pulse:

That flip angle, the saturation flip angle (the pulse that is used to probe B1+ efficiency, and not the excitation pulse), is hard-coded into the sequence on VE (to the best of my knowledge). On VB17 (Mt Sinai, so far), it is set in the "Sat Flip Angle" field under Sequence/Special.

For our acquisitions, this has always been set to 80 deg in Mt Sinai (which has caused some workflow issues with the SOP, as Alan documented in the scan notes and you picked up on here: #13 (comment)

For VE scanner (all but Mt Sinai), this 90 degree does not show up anywhere in the default json files. The "FlipAngle" value refers to the excitation FA, 10 degrees.

For VB, however, this field, "FlipAngle", describes the saturation flip angle, being either 80 or 60 or so, as per the name of the sequence.

Needless to say, this is going to make processing annoying. I suggest you treat all VE sites the same, with an assumption of Saturation FA=90, and parse the Mt Sinai site separately, using the FlipAngle field of the JSON file.

WRT DREAM:

How to deal with the DREAM sequences and how to convert them into nT/V maps is described in my original processing code, here: https://colab.research.google.com/drive/1HZYovXC4X-skmH-PcvumpXXXcOWESb0L#scrollTo=621d4320

Luckily, the FlipAngle in the JSON files of the DREAM sequences is the saturation flip angle, so it is easy to exctract it, and my code does it.

@evaalonsoortiz
Copy link
Contributor Author

Thanks @dpapp86 this is very helpful. If I'm following your code correctly, the conversion for FA-to-nT/V seems to be the same equation for both TFL and DREAM?

My last comment is that, the "flip angle" field in the .json files I've generated seems to be 90 degrees for all TFL scans, and not 10 degrees ... Should I trust this to be the saturation flip angle?

Screenshot 2024-03-10 at 9 58 50 PM

@dpapp86
Copy link
Contributor

dpapp86 commented Mar 11, 2024 via email

@evaalonsoortiz
Copy link
Contributor Author

This issue has been resolved. The code now looks for the saturation flip angle in the json files.

Sign up for free to join this conversation on GitHub. Already have an account? Sign in to comment
Labels
None yet
Projects
None yet
Development

No branches or pull requests

2 participants