gdimap (version 0.0-3)

simul.fandtasia: Simulation of Curved Fibre Bundles for von Mises-Fisher Fibre Orientation Mapping

Description

The synthesized field of diffusion profiles generated by simul.fandtasia are used to reconstruct ODF profiles using the GQI method. ODF profiles and fibre directions are estimated by relying on von Mises-Fisher (vMF) distributions for directional mapping.

Usage

simul.fandtasia(gridsz=32, b=4000, depth=3, sigma=0.01,
 threshold=0.5, showglyph=FALSE, savedir=tempdir(), snapshot=FALSE,
 pngfig="fandtasia")

Arguments

gridsz

dimension of squared grid to use in simulation (default 32)

b

strength of the magnetic diffusion gradient (default b-value=4000).

depth

sampling densities on the hemisphere used in simulation (default N=321; depth=3).

sigma

Rician noise level used in simulation; (default 0.01).

threshold

thresholding generalized fractional anisotropy (GFA) value at each voxel (default: 0.5).

showglyph

logical variable controlling visualization of voxel glyphs (default: FALSE).

savedir

directory for saving/loading processed results (default: tempdir().

snapshot

logical variable controlling line-maps snapshots (default: FALSE).

pngfig

base-name of the PNG file when snapshot is TRUE (default "fandtasia").

Value

• simul.fandtasia returns a field of 32x32 diffusion profiles in NIfTI format.

concept

• Simulation
• Glyph field mapping

Details

Noisy profiles may be simulated by adding Rician noise to the simulated diffusion profile, with a user defined standard deviation level specified as $\sigma$ (SNR=1/$\sigma$). The procedure is adapted from Barmpoutis' code to generate synthetic tensor diffusion-weighted MRI fields. The procedure is very time intensive for grids of size 32x32.

References

Ferreira da Silva, A. R. Facing the Challenge of Estimating Human Brain White Matter Pathways. In Proc. of the 4th International Joint Conference on Computational Intelligence (Oct. 2012), K. Madani, J. Kacprzyk, and J. Filipe, Eds., SciTePress, pp. 709-714.

Hornik, K., and Gruen, B. movMF: Mixtures of von Mises-Fisher Distributions, 2012. R package version 0.1-0.

Barmpoutis, A. Tutorial on Diffusion Tensor MRI using Matlab. Electronic Edition, University of Florida, 2010, http://www.mathworks.com/matlabcentral/fileexchange/file_infos/26997-fandtasia-toolbox.

See Also

simul.fandtasiaSignal, simulglyph.vmf, simul.simplefield

Examples

## simulation with a new generated field of profiles,
    ## of size 16x16 (for speed), with added noise 
    simul.fandtasia(gridsz=16, sigma=0.01)
    ## same as before, but showing crossing-fibre glyphs 
    simul.fandtasia(gridsz=16, sigma=0.01, showglyph=TRUE)
    ## using a 32x32 data field as in the original reference
    ## Warning: time-consuming example
    simul.fandtasia()