Qualitative Simulation of Photon Transport in Free Space Based on Monte Carlo Method and Its Parallel Implementation

Algorithm 1

Uniform framework for MC-based photon transport model.

Require: Set the initial parameters of the phantom, light source and optical imaging system.

1: Take out the preset photon package number  𝑛  and initialize detection matrix  𝑃(𝐫_𝑑)=0.

2: for  k = 1 to 𝑛  do

3:  Sample the light source and obtain the information of the kth photon package.

4:   Trace the transport process of the kth photon package in tissues.

5:   if  photon_dead = 0  then

6:     Record the parameters  𝐫,𝐬and𝑃(𝐫).

7:     Compute the visibility factor  𝜀(𝐫,𝐫_𝑑)  with (6).

8:     if  𝜀(𝐫,𝐫_𝑑)=1 then

9:          Compute the imaging spot or the imaging facula center  𝐫_𝑑.

10:        Compute the power  𝑃^′(𝐫_𝑑)  with (9).

11:        Update the power at the point  𝐫_𝑑  :  𝑃(𝐫_𝑑)=𝑃(𝐫_𝑑)+𝑃^′(𝐫_𝑑).

12:     else

13:        Break the current loop.

14:    end if

15:     else

16:       Break the current loop.

17:     end if

18:  end for