tedana.model.monoexponential¶
Functions to estimate S0 and T2* from multi-echo data.
Functions
fit_decay (data, tes, mask, masksum) |
Fit voxel-wise monoexponential decay models to data |
fit_decay_ts (data, tes, mask, masksum) |
Fit voxel- and timepoint-wise monoexponential decay models to data |
-
fit_decay
(data, tes, mask, masksum)[source]¶ Fit voxel-wise monoexponential decay models to data
Parameters: - data ((S x E [x T]) array_like) – Multi-echo data array, where S is samples, E is echos, and T is time
- tes ((E,)
list
) – Echo times - mask ((S,) array_like) – Boolean array indicating samples that are consistently (i.e., across time AND echoes) non-zero
- masksum ((S,) array_like) – Valued array indicating number of echos that have sufficient signal in given sample
Returns: - t2s_limited ((S,)
numpy.ndarray
) – Limited T2* map. The limited map only keeps the T2* values for data where there are at least two echos with good signal. - s0_limited ((S,)
numpy.ndarray
) – Limited S0 map. The limited map only keeps the S0 values for data where there are at least two echos with good signal. - t2ss ((S x E-1)
numpy.ndarray
) – Voxel-wise T2* estimates using ascending numbers of echoes, starting with 2. - s0vs ((S x E-1)
numpy.ndarray
) – Voxel-wise S0 estimates using ascending numbers of echoes, starting with 2. - t2s_full ((S,)
numpy.ndarray
) – Full T2* map. For voxels affected by dropout, with good signal from only one echo, the full map uses the T2* estimate from the first two echoes. - s0_full ((S,)
numpy.ndarray
) – Full S0 map. For voxels affected by dropout, with good signal from only one echo, the full map uses the S0 estimate from the first two echoes.
Notes
Fit monoexponential decay function to all values for a given voxel across TRs, per TE, to estimate voxel-wise and :
Replace infinite values in map with 500 and NaN values in map with 0.
Generate limited and maps by doing something.
-
fit_decay_ts
(data, tes, mask, masksum)[source]¶ Fit voxel- and timepoint-wise monoexponential decay models to data
Parameters: - data ((S x E x T) array_like) – Multi-echo data array, where S is samples, E is echos, and T is time
- tes ((E,)
list
) – Echo times - mask ((S,) array_like) – Boolean array indicating samples that are consistently (i.e., across time AND echoes) non-zero
- masksum ((S,) array_like) – Valued array indicating number of echos that have sufficient signal in given sample
Returns: - t2s_limited_ts ((S x T)
numpy.ndarray
) – Limited T2* map. The limited map only keeps the T2* values for data where there are at least two echos with good signal. - s0_limited_ts ((S x T)
numpy.ndarray
) – Limited S0 map. The limited map only keeps the S0 values for data where there are at least two echos with good signal. - t2s_full_ts ((S x T)
numpy.ndarray
) – Full T2* timeseries. For voxels affected by dropout, with good signal from only one echo, the full timeseries uses the single echo’s value at that voxel/volume. - s0_full_ts ((S x T)
numpy.ndarray
) – Full S0 timeseries. For voxels affected by dropout, with good signal from only one echo, the full timeseries uses the single echo’s value at that voxel/volume.