Bloch-Siegert B1 mapping module

[mathieuboudreau]

The idea for this module started off of a closed thread (#135 (comment))

I'll repeat the details and resources available that I mentioned there here.

I released the code for my B1 mapping paper recently, and the Bloch Siegert function is here: https://github.com/mathieuboudreau/B1Paper_Analysis/blob/master/src/b1/fitDataBSB1.m

• It'll need some refactoring, and probably some explanations/thinking about what some of the code does. My apologies for having "magic numbers" in there, most of them are are constants that I remember why they are there (e.g. 4096 was because of the datatype of the phase data I used being signed and having a range of +- 4096 to represent the phase angle).

The biggest hurdle I've always seen for this being a qMRLab module is that the AmpInt and Kbs parameters need to be known ahead of time, or calculated from the RF pulse curve. I've dug in a bit deeper in my old hard drives to look for some of the code I used to generate the RF pulse curves into a *.pta format for Siemens and also calculate the AmpInt and Kbs values from those *.pta files. I just found them, and uploaded them on GitHub: https://github.com/mathieuboudreau/bloch-sieg-pulse-tools.

• Again, these are not as well commented as I would like since they're from my early coding/PhD days, but at least they're here now!

I've also taken the liberty of extracting the relevant section of my PhD thesis where I discuss Bloch-Siegert mapping, which is available here: https://drive.google.com/file/d/17PKUa9eag66USoDjCn2aKDsni971zdcl/view?usp=sharing

• My entire thesis should also be available for free on the McGill website here: https://escholarship.mcgill.ca/concern/theses/73666714q?locale=en

I also uploaded one of my Bloch-Siegert datasets on OSF: https://osf.io/c9mqj/.