gdimap (version 0.1-9)

simul.fandtasiaSignal: Simulation of Crossing-Fibre Diffusion Profiles

Description

simul.fandtasiaSignal generates a synthesized field of diffusion profiles following Barmpoutis's algorithm.

Usage

simul.fandtasiaSignal(g, gridsz=32, b=4000, sigma=NULL, savedir=tempdir())

Arguments

g

matrix of 3D points on the S2 shell used in simulation.

gridsz

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

b

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

sigma

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

savedir

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

Value

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

Details

The diffusion field simulates a field of dimension gridsz x gridsz of diffusion profiles for testing crossing fibre orientation mapping. The procedure simul.fandtasiaSignal is an adaptation of Barmpoutis' MATLAB code to generate synthetic tensor diffusion-weighted MRI fields.

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.0-2.

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.fandtasia, simul.simplefield

Examples

## Not run: 
#     ## shell grid
#     s2 <- s2tessel.zorder(depth=3) 
#     simul.fandtasiaSignal(g=s2$pc, gridsz=16, sigma=0.01)
# ## End(Not run)