Compatible with Pulseq 1.4.0
🚨🚨🚨 NOTE: This is the dev
branch which hosts the bleeding edge version. For the most recent, stable release,
switch to the master
branch by clicking here. 🚨🚨🚨
- 👥 Contributors
- 📚 Citations
- 🔨 Installation
- ⚡ Lightning-start - PyPulseq in your browser!
- 🏃♂ Quickstart - example scripts
- 🤿 Deep dive - custom pulse sequences
- 👥 Contributing and Community guidelines
- 📖 References
- 📃 API documentation
Pulse sequence design is a significant component of MRI research. However, multi-vendor studies require researchers to be acquainted with each hardware platform's programming environment.
PyPulseq enables vendor-neutral pulse sequence design in Python [1,2]. The pulse sequences can be
exported as a .seq
file to be run on Siemens/GE/Bruker hardware by leveraging their respective
Pulseq interpreters. This tool is targeted at MRI pulse sequence designers, researchers, students and other interested
users. It is a translation of the Pulseq framework originally written in Matlab [3].
👉 Currently, PyPulseq is compatible with Pulseq 1.4.0. 👈
It is strongly recommended to first read the Pulseq specification before proceeding. The specification document defines the concepts required for pulse sequence design using PyPulseq.
If you use PyPulseq in your work, cite the following publications:
Ravi, Keerthi, Sairam Geethanath, and John Vaughan. "PyPulseq: A Python Package for MRI Pulse Sequence Design." Journal
of Open Source Software 4.42 (2019): 1725.
Ravi, Keerthi Sravan, et al. "Pulseq-Graphical Programming Interface: Open source visual environment for prototyping
pulse sequences and integrated magnetic resonance imaging algorithm development." Magnetic resonance imaging 52 (2018):
9-15.
Design pulse sequences using pypulseq
in your browser! Check out the ⚡ Lightning-start section to
learn how!
- @bilal-tasdelen
- @sairamgeethanath
- @schuenke
- @tonggehua
Please email/submit PR/open an issue if any contributors are missing.
2. 📚 Citations (reverse chronological)
- Hennig, J., Barghoorn, A., Zhang, S. and Zaitsev, M., 2022. Single shot spiral TSE with annulated segmentation. Magnetic Resonance in Medicine.
- Niso, G., Botvinik-Nezer, R., Appelhoff, S., De La Vega, A., Esteban, O., Etzel, J.A., Finc, K., Ganz, M., Gau, R., Halchenko, Y.O. and Herholz, P., 2022. Open and reproducible neuroimaging: from study inception to publication.
- Tong, G., Gaspar, A.S., Qian, E., Ravi, K.S., Vaughan, J.T., Nunes, R.G. and Geethanath, S., 2022. Open-source magnetic resonance imaging acquisition: Data and documentation for two validated pulse sequences. Data in Brief, 42, p.108105.
- Tong, G., Gaspar, A.S., Qian, E., Ravi, K.S., Vaughan Jr, J.T., Nunes, R.G. and Geethanath, S., 2022. A framework for validating open-source pulse sequences. Magnetic resonance imaging, 87, pp.7-18.
- Karakuzu, A., Appelhoff, S., Auer, T., Boudreau, M., Feingold, F., Khan, A.R., Lazari, A., Markiewicz, C., Mulder, M.J., Phillips, C. and Salo, T., 2021. qMRI-BIDS: an extension to the brain imaging data structure for quantitative magnetic resonance imaging data. medRxiv.
- Karakuzu, A., Biswas, L., Cohen‐Adad, J. and Stikov, N., 2021. Vendor‐neutral sequences and fully transparent workflows improve inter‐vendor reproducibility of quantitative MRI. Magnetic Resonance in Medicine.
- Geethanath, S., Single echo reconstruction for rapid and silent MRI. (ISMRM) (2021).
- Qian, E. and Geethanath, S., Open source Magnetic rEsonance fingerprinting pAckage (OMEGA). (ISMRM) (2021).
- Ravi, K.S., O'Reilly, T., Vaughan Jr, J.T., Webb, A. and Geethanath, S., Seq2prospa: translating PyPulseq for low-field imaging. (ISMRM) (2021).
- Ravi, K.S., Vaughan Jr, J.T. and Geethanath, S., PyPulseq in a web browser: a zero footprint tool for collaborative and vendor-neutral pulse sequence development. (ISMRM) (2021).
- Ravi, K.S. and Geethanath, S., 2020. Autonomous magnetic resonance imaging. Magnetic Resonance Imaging, 73, pp.177-185.
- Nunes, Rita G., et al. "Implementation of a Diffusion-Weighted Echo Planar Imaging sequence using the Open Source Hardware-Independent PyPulseq Tool." ISMRM & SMRT Virtual Conference & Exhibition, International Society for Magnetic Resonance in Medicine (ISMRM) (2020).
- Loktyushin, Alexander, et al. "MRzero--Fully automated invention of MRI sequences using supervised learning." arXiv preprint arXiv:2002.04265 (2020).
- Jimeno, Marina Manso, et al. "Cross-vendor implementation of a Stack-of-spirals PRESTO BOLD fMRI sequence using TOPPE and Pulseq." ISMRM & SMRT Virtual Conference & Exhibition, International Society for Magnetic Resonance in Medicine (ISMRM) (2020).
- Clarke, William T., et al. "Multi-site harmonization of 7 tesla MRI neuroimaging protocols." NeuroImage 206 (2020): 116335.
- Geethanath, Sairam, and John Thomas Vaughan Jr. "Accessible magnetic resonance imaging: a review." Journal of Magnetic Resonance Imaging 49.7 (2019): e65-e77.
- Tong, Gehua, et al. "Virtual Scanner: MRI on a Browser." Journal of Open Source Software 4.43 (2019): 1637.
- Archipovas, Saulius, et al. "A prototype of a fully integrated environment for a collaborative work in MR sequence development for a reproducible research." ISMRM 27th Annual Meeting & Exhibition, International Society for Magnetic Resonance in Medicine (ISMRM) (2019).
- Pizetta, Daniel Cosmo. PyMR: a framework for programming magnetic resonance systems. Diss. Universidade de São Paulo (2018).
>=Python 3.6, virtual environment recommended:
pip install pypulseq
- Create a new notebook on Google Colab
- Install PyPulseq
- Get going!
Or, explore an example notebook:
- Copy URL of an example notebook from [here][section-notebook-examples]
- On Google Colab, insert the copied link to get started
Every example script creates a pulse sequence, plots the pulse timing diagram and writes a .seq
file to disk.
Getting started with pulse sequence design using PyPulseq
is simple:
- Install PyPulseq
- First, define system limits in
Opts
and then create aSequence
object with it:import pypulseq as pp system = pp.Opts(max_grad=32, grad_unit='mT/m', max_slew=130, slew_unit='mT/m/s') seq = pp.Sequence(system=system)
- Then, design gradient, RF or ADC pulse sequence events:
Nx, Ny = 256, 256 # matrix size fov = 220e-3 # field of view delta_k = fov / Nx # RF sinc pulse with a 90 degree flip angle rf90, _, _ = pp.make_sinc_pulse(flip_angle=90, duration=2e-3, system=system, slice_thickness=5e-3, apodization=0.5, time_bw_product=4) # Frequency encode, trapezoidal event gx = pp.make_trapezoid(channel='x', flat_area=Nx * delta_k, flat_time=6.4e-3, system=system) # ADC readout adc = pp.make_adc(num_samples=Nx, duration=gx.flat_time, delay=gx.rise_time, system=system)
- Add these pulse sequence events to the
Sequence
object from step 2. One or more events can be executed simultaneously, simply pass them all to theadd_block()
method. For example, thegx
andadc
pulse sequence events need to be executed simultaneously:seq.add_block(rf90) seq.add_block(gx, adc)
- Visualize plots:
seq.plot()
- Generate a
.seq
file to be executed on a real MR scanner:seq.write('demo.seq')
Get in touch regarding running the .seq
files on your Siemens/GE/Bruker scanner.
PyPulseq
adheres to a code of conduct adapted from the Contributor Covenant code of conduct.
Contributing guidelines can be found here.
- Ravi, Keerthi, Sairam Geethanath, and John Vaughan. "PyPulseq: A Python Package for MRI Pulse Sequence Design." Journal of Open Source Software 4.42 (2019): 1725.
- Ravi, Keerthi Sravan, et al. "Pulseq-Graphical Programming Interface: Open source visual environment for prototyping pulse sequences and integrated magnetic resonance imaging algorithm development." Magnetic resonance imaging 52 (2018): 9-15.
- Layton, Kelvin J., et al. "Pulseq: a rapid and hardware‐independent pulse sequence prototyping framework." Magnetic resonance in medicine 77.4 (2017): 1544-1552.