# Forward model analysis in diffuse optical tomography using RTE and FEM.

INTRODUCTIONDiffuse optical tomography (DOT) [1, 9] is a method of imaging technique using near-infrared (NIR) region for bio-medical applications that has potential to provide the three-dimensional images of the optical parameters such as absorption coefficient and the scattering coefficient [12, 16]. The light transport of photons through a biological tissue is established during diffusion equation. This type of photons is scattered or absorbed in turbid media. Since the high level of scattering in a most biological tissue medium image reconstruction in DOT is suffering from ill-posed nonlinear problem [7, 19]. It has so many applications in a medical imaging. The biomedical applications of DOT include early stage detection of breast cancer, measurement of brain function and others in bio-measurement fields.

The forward model frequently utilized the light transport models it contains stochastic models such as Monte-Carlo [8] simulation method or deterministic models such as RTE or diffusion equation (DE). In the previous techniques, photon transport in a biological tissue can be equivalently modeled in numbers with Monte-Carlo simulation [5, 13] method or analytically by RTE. It is referred to as the Boltzmann transport equation (BTE), which is the term of the balance of energy within a boundary volume element of the scattering medium. The analytical results to the RTE [4, 18] for a very easy method but for more practical media; with very complex multiple scattering effects and the numerical methods are required. It has some disadvantages such as difficult to solve without introducing approximation and it is computationally very expensive. So the DOT uses the diffusion equation (DE), which is describing the light transport problem in a biological tissue. The solution to the diffusion equations for photon transports are more computationally efficient and can be analyzed deeply.

The numerical method contains finite difference method (FDM), finite element method (FEM) [1], finite volume method (FVM) and boundary element methods (BEM). The forward model is to predict the optical flux density at the detectors given in a geometric model with optical parameters as source-detector location and functionality. The forward model problem was developed by using the finite element method (FEM) for solving the DE [15] to the RTE. FEM is a numerical method [11] for calculating the approximation solutions to the boundary value problem and to generate the optical properties and optical flux for a given distributions of an absorption and diffusion coefficient. The robin boundary conditions are incorporated at the boundary values. It has the true representation of complex geometry, simple representation of the total solution, generates the spatially resolved images and develops the low-resolution [17,24] of the functional image with high-resolution complementary structures. In the inverse model, the solution of the forward model problem is used for calculating the jacobian matrix and the simultaneous equation is solved using a conjugate gradient function.

Mohammad Ali Ansari et al [20] has provided a volume fraction of blood tissue by light absorption using boundary element method (BEM) with a suitable procession and accuracy. Tissue depth is the main parameter in diffuse optical tomography which BEM has a good ability to reconstruct. However, in wide tissue layers its accuracy decrease due to the small value of data received from reflectance. Young Sik Jun and Woon Sik Baek [21] has developed the light transport in a highly scattered tissue medium. The biological tissue is modeled by using the diffusion equation that is Boltzmann transport equation (BTE). However, it has a large amount of variables for reflection and to obtain the result of BTE is very difficult. Gao and Zhao [22] analyzed the numerical solver for radiative transport equation. The mathematical calculation of RTE is challenging due to its high dimension. Bo Bi and Bo Han [18] analyzed the forward model design of DOT using radiative transfer equation (RTE), which provides accurate and efficient results within biological tissue medium and models the inverse problem using regularization method. In this method, we can get the most accurate results from inner measurements. In an execution time, we cannot obtain the inner measurement data frequently. Hence, the RTE cannot reconstruct inclusions with complicated internal structure accurately with a small amount of boundary angular-averaged measurements. Hence, we have evaluated the finite element method numerically. It is the most flexible one in terms of dealing with complex geometry and difficult boundary conditions. Table I shows the quantities and dimensions of the diffusion equation.

The review of the paper is as follows in section II, we had discussed about the forward modeling. The forward model was described by the diffusion equation (DE) governing boundary value problems; the solution was obtained by radiative transfer equation (RTE) and finite element method (FEM). In section III, results and discussion based on the forward model were analyzed for normal and cancer affected people. In section IV, we discuss in the conclusion that the diffusion equation modeled by FEM performance was better compared to RTE.

II. Forward Model:

The forward model of DOT is to find out the outgoing current on the detectors when the incident impulse and the absorption coefficients and scattering coefficients are known. It describes photon propagation in tissue [14]. In an experimental method of getting potential measurements is as follows. First, a set of laser source [10] devices and detectors on the boundary of the object guess distribute the optical properties absorption and scattering coefficient. The incident impulse launched from the laser source and the resulting measurement from all the detectors was recorded.

A. Radiative transfer equation:

The radiative transfer equation (RTE) [18, 1] describing the photon propagation in tissue medium. The RTE is solved by analytically in a biological tissue; the RTE is commonly accepted model for light propagation in a biological tissue. The time domain in an RTE is given by,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (1)

Where, the radiance is denoted by L(r, t, [??]) , q(r, t, [??]) is the source of inside and the scattering phase function is denoted by f ([??],[??]). The energy transfer per unit time is denoted by L and the photons in a solid angle is [d.sup.2] [??] during an elemental area da is specified by its outgoing normal unit n at position [??] is set by,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (2)

The exitance [GAMMA] with a unit area perpendicular to [??] is achieved by the above integrating equation. Over all angle is obtained by,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (3)

In RTE, both spatial location and the angular direction of the discretization are implemented and the equations are derived by analytically or stochastically.

B. Diffusion equation:

The forward model frequently handles the light transport model such as a diffusion equation (DE). In DOT, the diffusion equation in a frequency domain [4,23] can be obtained by,

(- [nabla].k(r)[nabla] + ([[mu].sub.a](r) + jw/c))[phi](r,w)

= [q.sub.0](r,w) (4)

Where w = 2[pi]f , k(r) is a diffusion coefficient, [phi](r) is a photon flux density at the position r, c is the speed of light in a medium and [[mu].sub.a] (r) is the absorption coefficient. The diffusion coefficient k (r) is given by,

k(r) = [1/3([[mu].sub.a](r)+[[mu]'.sub.s](r)] (5)

[[mu].sub.a] and [[mu]'.sub.s] illustrate the absorption and the reduced scattering coefficient ([[mu].sub.a] [much less than] [[mu]'.sub.s]) [13] respectively. The reduced scattering coefficient is described by [[mu]'.sub.s] = [[mu].sub.s] (1 - g) with [[mu].sub.s] represents the scattering coefficient and the g is an anisotropy factor, with the typical range of 0.7 to 0.9. The photon input is from a source of constant intensity located at r = [r.sub.0].

The output optical flux on the boundary is given by,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (6)

Where, the normal unit vector to the surface [partial derivative][OMEGA], is denoted by [??]. Type III robin boundary condition [12] is used for diffusion equation, as it models refractive index mismatch at the boundary of the soft tissue. The robin boundary condition is solved using FEM and [phi](r) is obtained.

C. Solution of DE using Finite element method:

The light diffusion equation in a frequency domain is numerically calculated by using the finite element method (FEM) [3, 6]. In the FEM, we have applied the final technique in the form well-known as the galerkin formulation method. The Galerkin formulation method [2, 4] proceeds solution of the optical flux [phi],

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (7)

Where, [s.sub.0] (r, t) is the source term. If it is also a solution of,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (8)

Where, [psi](r)is a function from a suitable test space satisfying the same boundary conditions as [phi]. Integration by the parts,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (9)

For this problem to be solvable, [psi] and all first derivatives must be integral over [OMEGA]. The solution [phi] to the diffusion equation (4) is approximated by a continuous function,

[[phi].sup.h] (r, t) = [[summation].sup.D.sub.i=1] [[phi].sub.i](t)[[psi].sub.i](r)[v.sup.h] (10)

such that,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (11)

This can be expressed in matrix notation as,

[K(k)+C([[mu].sub.a]c)][phi]+M [partial derivative]/[partial derivative]t = Q + [beta] (12)

Where,

[MATHEMATICAL EXPRESSION NOT REPRODUCIBLE IN ASCII] (13)

We assume the boundary condition [[psi].sub.i] = 0 on[partial derivative][OMEGA]. The D x D matrices K, C, and B are a light while, they will have non-zero entries only when [N.sub.j] and [N.sub.i] are vertices of the same element.

RESULTS AND DISCUSSION

The radiative transfer equation (RTE) and the finite element method (FEM) solution to the forward model design, we consider a two-dimensional circular image of human tissue has 20mm radius with six sources and detectors placed at the boundary of the head. In a simulation method, we execute the programs on a 3.0 GHz window 7 PC with 8 GB RAM in MATLAB 2013 a. The optical properties of the forward model were set as 0.01 [cm.sup.-1] absorption and 1.0 [cm.sup.-1] scattering coefficients for normal people.

[FIGURE 1 OMITTED]

[FIGURE 2 OMITTED]

[FIGURE 3 OMITTED]

[FIGURE 4 OMITTED]

Figure 1 displays the angular discretization of the 2-D image and in this simulation which was obtained by placing the six detectors at the different angles on the tissue boundary. These phasor angles are at 0[degrees] or 360[degrees], 45[degrees], 135[degrees], 180[degrees], 225[degrees] and 315[degrees] respectively. Figure 2 display the 2-D circular mesh image for RTE. This type of circular mesh image is formed by the triangular element in our calculation with 1097 nodes and 9600 elements. In this method, the number of nodes in a forward problem for RTE is always greater than the number of nodes in the reconstructed mesh. It is used to improve our spatial resolution. Figure 3 shows the real measurement values of absorption and scattering coefficients for normal people ranges from 0.01 [cm.sup.-1] to 1.0 [cm.sup.-1]. Figure 4 shows the real measurement values of absorption and scattering coefficient for cancer affected people, it ranges from 2.4 [cm.sup.-1] to 5.5 [cm.sup.-1] respectively.

[FIGURE 5 OMITTED]

[FIGURE 6 OMITTED]

[FIGURE 7 OMITTED]

Figure 5 represents the two-dimensional circular mesh diagram formed by the triangular element for the FEM solution of the forward model. Figure 6 shows the angular discretization of the two-dimensional image for FEM using diffusion equation. Figure 7 represents the magnitude of the optical properties such as absorption and scattering coefficient for FEM.

[FIGURE 8 OMITTED]

Table II and Table III presents the absorption and scattering coefficient values along with optical flux reading for the normal people and cancer affected people as observed from the reconstructed image. We can distinguish the normal people and the affected people since normal people have absorption coefficient less than 2 [cm.sup.-1] and scattering coefficient values less than 5 [cm.sup.-1]. In the case of cancer affected people the absorption coefficient, exceeds 2 [cm.sup.-1] as well as scattering coefficient exceeds 5 [cm.sup.-1]. Figure 8 represents the optical parameters such as absorption coefficient in the range of 0.01 [cm.sup.-1] to 1.4 [cm.sup.-1] and the scattering coefficient in the range of 1 [cm.sup.-1] to 3.8 [cm.sup.-1] for normal people.

[FIGURE 9 OMITTED]

[FIGURE 10 OMITTED]

Figure 9 represents the optical parameters such as absorption coefficient in the range of 2.4 [cm.sup.-1] to 5.4 [cm.sup.-1] and the scattering coefficient in the range of 5.7 [cm.sup.-1] to 12.18 [cm.sup.-1] in tumor patient. Figure 10 shows the optical flux value for normal and cancer affected people. The optical flux (f) is energy per unit time that is radiated from a source over optical wavelengths in nm. The above graph shows the relationship between the optical flux and wavelengths in the range of 750 nm to 840 nm. The normal people will have lower flux value compared to cancer affected people. We can find the resolution of the reconstructed images by determining the signal to noise ratio (SNR) and CPU time for RTE and FEM. The SNR ratio is given as,

SNR = 10[log.sub.10][[parallel]signal[parallel].sub.2]/[[parallel]noise[parallel].sub.2] (14)

In, table IV parameters values are compared that is SNR and CPU time. The reconstructed image of an existing system provides low SNR ratio and higher CPU time about 44.231. The RTE gives 37.0645 CPU time and also provides higher SNR ratio compared than the existing system. Then, we compared to RTE, the FEM provides the higher SNR ratio and very less CPU time. Therefore, FEM is a most important method for improving system performance.

Conclusion:

By using algorithms in diffuse optical tomography (DOT) with NIR light ranging from 780 nm to 850 nm, we had performed the numerical simulation using the radiative transfer equation (RTE) and diffusion equation solved using finite element method (FEM) to find out the result of the soft tissue optical parameters that is absorption and scattering coefficient explicitly for normal and cancer affected people. The optical flux was also evaluated based on the forward model numerically. We had evaluated the performance of a system by varying the position of source and detectors, thereby calculating the geometrical symmetry of the tissue in a turbid medium. FEM is more efficient and flexible in implementing different boundary conditions. The parameter values such as SNR and CPU time are compared for both algorithms. FEM provides better performance and high rate of SNR value than RTE. By analyzing these values, Cancer affected tissue was predicted from the forward model, which in turn avoids the reconstruction of the image. It is proved that the FEM gives efficient results compared to RTE.

REFERENCES

[1.] Arridge, S.R., 1999. "Optical tomography in medical imaging," Inverse Problems, 15(2): R41-R93.

[2.] Arridge, S.R., M. Schweiger, M. Hiraoka and D.T. Delpy, 1993. "Finite Element Approach for Modeling Photon Transport in Tissue," Medical Physics, 20: 299-309.

[3.] Schweiger, M., S.R. Arridge and I. Nissil, 2005. "Gauss-Newton method for image reconstruction in diffuse optical tomography, ''Physics in Medicine and Biology, 50(10): 2365-2386.

[4.] Biswas, S.K., K. Rajan, R.M. Vasu and D. Roy, 2011."Accelerated Gradient Based Diffuse Optical Tomographic Image Reconstruction," Medical Physics, 38: 539.

[5.] Boas, D.A., J.P. Culver, J.J. Stott and A.K. Dunn, 2002. "Three-Dimensional Monte Carlo Code for Photon Migration through Complex Heterogeneous Media Including the Adult Human Head," Optics Express, 10(3): 159-170.

[6.] Schweiger, M., S.R. Arridge and D.T. Delpy, 1993. "Application of the Finite-Element Method for the Forward and Inverse Models in Optical Tomography", Journal of Mathematical Imaging and vision, 3: 263283.

[7.] Arridge, S.R., M. Schweiger, M. Hiraoka and D.T. Delpy, 1995. "The Finite Element model for propagation of light in scattering media: Boundary and source conditions," Med. Phys., 22: 1779-1792.

[8.] Sassaroli, 2011. "Fast Perturbation Monte Carlo Method for Photon Migration in Heterogeneous Turbid Media," Optics Letters, 36(11): Opt. http://dx.doi.org/10.1364/OL.36.002095

[9.] Gibson and H. Dehghani, 2009. "Diffuse optical imaging," Phil. Trans. R. Soc., A 367, 3055-3072.

[10.] Abdoulaev, G.S. and A.H. Hielscher, 2003. "Three-Dimensional Optical Tomography with the Equation of Radiative Transfer," Journal of Electronic Imaging, 12(4): 594-601.

[11.] Yalavarthy, P.K., K. Karlekar, H.S. Patel, R M. Vasu, M. Pramanik, P.C. Mathias, B. Jain and P.K. Gupta, 2005. "Experimental Investigation of Perturbation Monte-Ca Based Derivative Estimation for Imaging Low scattering Tissue," Optics Express, 13(3): 985-988.

[12.] Arridge, S.R., J.C. Schotland, 2009. "Optical Tomography: forward and inverse problems," Inverse Problems, 25: 123010-123069.

[13.] Groenhuis, R.A.J., H.A. Ferwerda and J.J. Ten Bosch, 1983. "Scattering and absorption of turbid materials determined from reflection measurements (parts 1 and 2)," Appl. Opt., 22: 2456-2467.

[14.] Gao, H., S. Osher and H.K. Zhao, 2012. "Quantitative photo acoustic tomography," in Mathematical Modeling in Biomedical Imaging II, Lecture Notes in Mathematics, pp: 131-158, Springer, Berlin, Germany.

[15.] Samir Kumar Biswas, Rajan Kanhirodan, Ram Mohan Vasu, 2013. "Models and Algorithms for Diffuse Optical Tomographic System," Int. J. Communications, Network and System Sciences, 6: 489-496.

[16.] Gibson, A.P., J.C. Hebden and S.R. Arridge, "Recent advances in diffuse optical imaging," Physics in Medicine and Biology, 50: 4, article R1.

[17.] Sukanyadevi, R., K. Umamaheswari, S. Sathiyamoorthy, 2013. Resolution improvement in Diffuse Optical Tomography, International Journal of Computer Applications ICIIIOES (9): 37-41.

[18.] Bo Bi, Bo Han, Weimin Han, Jinping Tang and Li Li, 2015."Image Reconstruction for Diffuse Optical Tomography Based on Radiative Transfer Equation" computational and mathematical methods in Medicine, Article ID 286161.

[19.] Umamaheswari, K., S. Sathiyamoorthy, 2015. Stein's Unbiased Risk Estimate Regularization (SURE) for Diffuse Optical Tomography (DOT) System Enhances Image Reconstruction with High Contrast to Noise Ratio (CNR), International Journal of Applied Engineering Research, 10(24): 21186-21191.

[20.] Ansari, M.A., M. Erfanzadeh, Z. Hosseini, E. Mohajerani, 2014. "Diffuse Optical Tomography: Image Reconstruction and Verification," J Lasers Med Sci., 5(1): 13-8.

[21.] Young sik Jun and woon sik baek, 2010. "Experimental Reconstruction Images of Tissue Phantom by Diffuse Optical Tomography," Journal of Physics: Conference Series, 224: 012146.

[22.] Gao, H. and H.K. Zhao, 2013. "Analysis of a numerical solver for radiative transport equation, "Mathematics of Computation, 82(281): 153-172.

[23.] Ebraheem Sultan, Laleh Najafizadeh, Amir H. Gandjbakhche, Kambiz Pourrezaei and Afshin Daryousha, 2013. "Accurate optical parameter extraction procedure for broadband near-infrared spectroscopy of brain matter," Journal of Biomedical Optics, 18(1): 017008.

[24.] Sultan, E., K. Manseta, A. Khwaja, L. Najafizadeh, A. Gandjbakhche, K. Pourrezaei, A.S. Daryoush, 2011. "Modeling and Tissue Parameter Extraction Challenges for Free Space Broadband fNIR Brain Imaging Systems," Proc. of SPIE 7902,790223

(1) K. Uma Maheswari, (2) Dr. S. Sathiyamoorthy, (3) J. Hilda Amala

(1,2,3) Department of ECE, J.J. College of Engineering & Technology Trichy. India

Received 25 April 2016; Accepted 28 May 2016; Available 2 June 2016

Address For Correspondence:

K. Uma Maheswari, Department of ECE, J.J. College of Engineering & Technology Trichy. India

Table 1: Description of quantities S. No Quantity Description Dimension 1. [bar.f] Average cosine of None (direction independent) phase function 2. k (r) Diffusion ([m.sup.2] coefficient [s.sup.-1]) 3. [??] normal unit vector (mm) 4. [s.sub.0] (r, t) isotropic component ([m.sup.-3]) of the source 5. P radial distance (mm) 6. [[mu].sub.a] (r) absorption ([m.sup.-1]) coefficient at position r 7. [[mu].sub.s](r) scattering ([m.sup.-1]) coefficient at position r 8. [[mu].sub.s]'(r) reduced scattering ([m.sup.-3]) coefficient 9. [phi](r, t) photon density at ([m.sup.-3]) point r, at time t 10. c speed of light in a (m) medium Table II: Absorption and scattering coefficient values for normal people S. No Absorption Scattering Optical flux coefficient coefficient [phi] a.u. ([[mu].sub.a]]) ([[mu].sub.s]]) [cm.sup.-1] [cm.sup.-1] 1 0.01 1 1.18E-015 2 0.5 2.15 2.22E-016 3 0.8 2.7 1.97E-015 4 1.4 3.8 1.32E-015 Table III: Absorption And Scattering Coefficient Values For Cancer Affected People S. No Absorption Scattering Optical flux coefficient coefficient [phi] a.u. ([[mu].sub.a]) ([[mu].sub.s]) [cm.sup.-1] [cm.sup.-1] 1 2.4 5.7 3.25E-015 2 3.2 7 2.96E-015 3 4.7 10.3 2.74E-015 4 5.4 12.18 1.88E-015 Table IV: Performance Analysis Of Forward Model Design Parameters Existing Methods system Radiative Finite transfer element equation method (RTE) (FEM) SNR 1.3321 1.9148 2.5039 CPU time 44.231 37.0645 35.7215

Printer friendly Cite/link Email Feedback | |

Author: | Maheswari, K. Uma; Sathiyamoorthy, S.; Amala, J. Hilda |
---|---|

Publication: | Advances in Natural and Applied Sciences |

Date: | Jun 1, 2016 |

Words: | 3509 |

Previous Article: | An integrated approach for deaf speech enhancement using modified SOM clustering clustering for deaf speech enhancement. |

Next Article: | An effective clustering method with search structure for large multidimensional dynamic indexes. |

Topics: |