Self calibration iso-pathlength point in cylindrical tissue geometry : solution of steady-state photon diffusion based on the extrapolated zero-boundary

Near-infrared optical techniques permit tissue diagnosis by surface measurement. However, the geometrical shape of this interface profiles the intensity of the surface measurement, which is found to have an iso-pathlength (IPL) point allowing for absorption identification independent of tissue scattering. The IPL point was projected in Monte Carlo (MC) simulation, validated experimentally in cylindrical tissues, but remains underappreciated through analytical approaches. In this work, we present an analytical solution of an IPL point for steady-state diffusion based on the extrapolated zero-boundary condition. The same IPL points were found when comparing this solution to 3-D MC simulations for a tissue radius range of 5-8mm. © 2018 Optical Society of America under the terms of the OSA Open Access Publishing Agreement


Introduction
Light-tissue interactions are commonly studied for near infra-red (NIR) spectroscopy for imaging and diagnosis.While imaging is bound to the surface due to the high tissue scattering [1][2][3][4], many diagnosis methods, such as the photoplethysmography (PPG) and pulse oximetery, focus on sensing instead of imaging, requiring rapid forward estimation of the surface measurement in response to preferably a homogeneous optical entity of the medium [5,6].
There are two inversion approaches for the investigation of light-tissue interactions: one numerical and the other analytical.The most robust numerical method is the Monte Carlo (MC) simulation [7], which is a study of statistical photon migration based on the optical properties of the different tissue regions.However, the accuracy of such statistical approaches is limited when speed is concerned.Analytical methods based on the diffusion approximation to radiative transfer allow fast forward computation [8].These methods apply mathematical tools to identify Green functions that fulfill the boundary conditions of the problem in a specific geometry.The common boundary condition of zero fluence rate is set at the tissue-air interface, and two inputs are defined; an isotropic source inside the tissue and an additional image source outside the tissue.The idea is that since photons have an initial direction of entering the tissue normal to the boundary, the diffuse photon propagation in a tissue with a reduced scattering coefficient of μ s ', having units of cm −1 , can be attributed to a boundaryaffected isotropic source after they scatter and completely lose their memory at r = R 0 -R a (Fig. 1), where: a R =1/ '.In this method, a solution based on Bessel and modified Bessel functions was found for cylindrical geometries [9].However, a more accurate extrapolated zero-boundary condition yields a solution based on modified Bessel functions of the first and second kind [10][11][12][13].In this paper we will use the solution of steady-state photon diffusion based on the extrapolated zero-boundary and compare it to MC simulations in cylindrical geometry.

Solution of steady-state diffusion in cylindrical concave geometry
In the extrapolated zero-boundary condition method, the Fresnel reflection due to the refractive index change on the tissue-detector interface is accounted for by placing the boundary of zero fluence rate away from the physical tissue-air interface.This is accomplished by defining an extrapolated boundary located a distance R b = 2A•D outside the turbid medium, meaning at r = R 0 + R b (Fig. 1), where D is the diffusion coefficient to be defined below and A accounts for the effect of the refractive index mismatch as A = (1 + R eff )/(1-R eff ), with R eff being the effective reflection coefficient [14].The image source, that sets the fluence rate to zero on the extrapolated boundary, is now located at r = R 0 + R a + 2R b .
In steady-state, photon diffusion in a homogeneous medium is described by [5,6]: where Ψ is the photon fluence rate in units of W•cm −2 , S is the illumination source intensity in units of W•cm −3 , µ a is the absorption coefficient of the tissue, having units of cm −1 , and the diffusion coefficient D, having units of cm [15], is defined by µ a and µ s ' to be: Fourier analysis of Eq. ( 2) derives the Green functions that fulfill the boundary condition.For the cylindrical geometry of the source and detector on the same azimuthal plane as specified in Fig. 1, the following solution to Eq. (2) was found [13]: where This azimuthal distribution of the photon fluence rate correlates to the previously defined full scattering profile (FSP) [16].

MC simulation of cylindrical geometry
MC simulations revealed that the FSP has a fixed intensity point identified as the isopathlength (IPL) point, which does not depend on the reduced scattering coefficient.The location of this point depends solely on the geometry of the tissue-air interface.These findings were validated by tissue mimicking phantoms [17], as well as in vivo human finger measurements [18].The MC simulation is based on the assumption that photons enter the tissue at a certain known direction and interact with the tissue via absorption or scattering processes in accordance to its optical properties.Once the photon enters the tissue interface, its probability to scatter after a distance dr is given by: ( ) where s  is a random number bounded within [-1, 1].Hence, the subsequent scattering angle is uniformly distributed when the reduced scattering property is concerned.This process is repeated until the photon reencounters the tissue-air interface.The photons' locations in this tissue-air interface assemble the FSP.

IPL point in MC simulation and steady-state diffusion theory
The FSPs of a 6mm radius cylinder was simulated by such a 3-D MC simulation with an absorption coefficient of µ a = 0.025cm −1 and different reduced scattering coefficients.The FSPs from the MC simulations are shown in Fig. 2(a) (gray asterisks, black squares, and red triangles correspond respectively to a reduced scattering coefficient of 10, 16 and 26cm −1 ).The FSPs reveal an IPL point at 171°.The FSPs based on the steady-state diffusion theory are developed by assuming a spatially impulse light source injected to the tissue at an initial direction normal to the tissue surface (Eq.( 4)), R eff = 0.43, 150 orders of the Bessel and modified Bessel functions (m = 0-150), and 50 k eff values.The results are presented in Fig. 2(b) (black dotted, blue solid, green dotdashed and red dashed lines represent respectively a reduced scattering coefficient of 16, 18, 20 and 26 cm −1 ).In the diffusion theory, lower reduced scattering coefficients are less relevant and thus not as accurate since they yield a large transport mean free path (MFP'), which is too far from the light source.Hence, the lowest reduced scattering coefficient was chosen to be 16cm −1 , corresponding to MFP' = 0.625mm, in order to maintain an order of magnitude from the radius.The FSPs present an IPL point at an identical position of 171°.
Next, we compare both methods (Fig. 2 (c)).The FSPs of the MC simulation for different reduced scattering coefficients (black squares and red triangles represent respectively a reduced scattering coefficient of 16 and 26 cm -1 ) are used for examining the FSPs of the diffusion theory (black dotted and red dashed lines represent respectively a reduced scattering coefficient of 16 and 26 cm −1 ).These two methods, although producing very different FSPs, converge to the same IPL point.The existence of a single point matching IPL over the diverse FSPs from the numerically straightforward method of diffusion theory when compared to the computationally intensive MC method is encouraging.
It is well known that diffusion theory-based solutions are inaccurate near the light source in addition to low scattering condition, when in comparison to the more precise MC methods.The inaccuracy of the diffusion-based solution presented in Eq. ( 4) at a near-the-entry distance is also apparent from the IPL point's location in various tissue radiuses (Fig. 3).A comparison between MC simulation and diffusion theory (red circles and blue squares in Fig. 3) reveals this divergence in the low radius range since the IPL point is too close to the source.A recent work has demonstrated a simple approach to enhancing the diffusion-based estimation of the photon fluence rate at distances near the point-of-entry for surface measurement [19].That approach, when extended to the analysis of the FSP in cylindrical geometry, is expected to enhance the match between analytical prediction and MC measurement of the IPL.

Discussion and conclusions
In real-time sensing applications, the reading from sensors implemented with IPL principles may be conducted when the sensor concave geometries change in adapting to tissue-air interface.In such situations, any change of the sensor geometry, in terms of the radius of the concave shape, must be accounted for rapidly for accurate inversion of the tissue optical properties.Analytical knowledge of the IPL point will be important for real-time IPL-based identification of tissue absorption properties for cylindrical tissue geometries.
For endoscopic applications involving a convex tissue geometry that can be simplified by a cylinder of air surrounded by tissue of homogeneous properties, the solution of photon diffusion as a counterpart of Eq. ( 4) is the following: A simulation of the FSPs for the convex geometry (Fig. 4) presents results similar to that in a concave geometry of the same radius (Fig. 2).This finding was expected, since we have projected in the past that the IPL point phenomena in different geometries can be predicted by the distance between the light source and the detector [20].The simulation results will encourage further experimental studies of IPL behaviors in such a convex geometry, both in MC simulation and using phantom materials, for gauging the solution of steady steady-state photon diffusion.Furthermore, the geometry where the source and detector are not on the same azimuthal plane is also possible.The evaluation of the IPL point along the longitudinal direction should be investigated as well.
The comparison of the IPL points given by MC and by diffusion theory for the longitudinal direction, whether that will give the same IPL points comparing to the presented azimuthal directions or not, could give additional support to this work.
To conclude, in this study we derived the FSP from the steady-state diffusion theory based on the extrapolated zero-boundary condition.A comparison of different reduced scattering coefficients revealed the IPL point, where the intensity remains constant regardless of the scattering property of the tissue.The position of the IPL according to this theoretical method is in agreement with the MC simulations performed under the same geometrical parameters and optical properties [21,22].

Fig. 1 .
Fig. 1.Details of a cylindrical geometry corresponding to a concave tissue medium.The red arrow indicates the incident light.The positions of the image source, the equivalent isotropic source, and the extrapolated zero-boundary (dashed line) are also illustrated.
new direction (θ new ) with respect to the immediate previous (old) direction (θ old ) is amended by the anisotropy g:

Fig. 2 .
Fig. 2. FSPs and the IPL point resulted from 3-D MC simulation and diffusion theory for steady-state photon illumination into a concave cylindrical tissue domain.Tissue radius is R 0 = 6mm.The FSP has an IPL point which is constant for different reduced scattering coefficients: (a) In the MC simulation gray asterisks, black squares and red triangles correspond respectively to μ s ' = 10, 16 and 26cm −1 .(b) In the diffusion theory black dotted, blue solid, green dot-dashed and red dashed lines represent respectively μ s ' = 16, 18, 20 and 26 cm −1 .(c) Comparison of 3-D MC simulation and diffusion theory reveals a common IPL point at 171°.

Fig. 4 .
Fig.4.The FSPs of steady-state photon diffusion for a convex geometry with a radius of 6mm, where a cylinder of air is surrounded by tissue (Eq.(7).Different reduced scattering coefficients (μ s ' = 16 and 26 cm −1 corresponding to black dotted and red dashed lines) have a common IPL point at 171°.