Join the official 2019 Python Developers Survey

SSFP simulation

## SSFP

Simple steady-state free precession simulation. The goal is to provide a simple to use, pip-installable solution for simulating this wonderful pulse sequence.

In this package:

• bSSFP: bssfp()
• GS solution: gs_recon()

## Installation

Should be as easy as:

```pip install ssfp
```

## Usage

See ssfp.examples for typical usage. You can run examples like:

```python -m ssfp.examples.basic_bssfp
```

Balanced steady-state free precession can be called through bssfp(). This is an implementation of equations [1–2] in . These equations are based on the Ernst-Anderson derivation  where off-resonance is assumed to be subtracted as opposed to added (as in the Freeman-Hill derivation ). Hoff actually gets Mx and My flipped in the paper, so we fix that here. We also assume that the field map will be provided given the Freeman-Hill convention.

```from ssfp import bssfp

# Here's the simplest usage, see docstring for all the possible
# function arguments
sig = bssfp(T1, T2, TR, alpha)
```

We can also easily get the Geometric Solution to the elliptical signal model as described in  as follows:

```from ssfp import gs_recon
recon = gs_recon(phased_cycled_images, pc_axis=-1)

# Notice that we can specify the axis where the phase-cycles live
```

PLANET  is a method for simultaneous T1, T2 fitting for bSSFP phase-cycled data. Call like this:

```from ssfp import planet

# For a single pixel:
Meff, T1, T2 = planet(
phased_cycled_pixels, alpha, TR, T1_guess,
```

## Project details

This version 0.2.0 0.1.0 0.0.1