Statistical (or stochastic) modelling techniques model individual photon trajectories and have the advantage that the Poisson error is incorporated into the model in a natural and elegant way. The most commonly used statistical technique in diffuse optics, and that which is often regarded as the gold standard to which other techniques are compared, is the Monte Carlo method. The geometry of the model is defined in terms of µa, µs, and the refractive index, and the trajectories of photons, or packets of photons, are followed until they either escape from the object under study or are absorbed. By continuing until the required counting statistics are obtained, data with arbitrarily low statistical errors can be simulated. In optical imaging, Monte Carlo techniques are commonly used to calculate light propagation in non-diffusive domains where the diffusion approximation does not hold (Boas et al. 2002, Okada and Delpy 2003, Hayashi et al. 2003), or to validate results obtained using other, faster, methods (Schweiger et al. 1995, Chernomordik et al. 2002b, Dehghani et al. 2003a).
http://medicalimaging.wikia.com/wiki/Forward_problem_in_diff...