1. Introduction {#s0005}
===============

Dural hematoma often occurs after traumatic brain injury. Inability to make accurate and timely diagnosis for a reasonable treatment plan can result in irreversible brain damage and endanger the life of the patient. Therefore, non-invasive detection of traumatic dural hematoma is always a research focus in the biomedical engineering field ([@b0050]).

Due to the high absorption of 650--900 nm near-infrared light by the hemoglobin molecules within the tissue, the optical properties of the brain tissue tested can be obtained by analyzing the emitting light, so that rapid and non-invasive detection of dural hematoma can be achieved ([@b0125], [@b0035]). Near-infrared spectroscopy technology has been widely used in clinical applications such as functional neuroimaging ([@b0085]), brain tumor imaging ([@b0075], [@b0130]), cerebral blood flow measurement ([@b0070]), and brain hematoma detection ([@b0020]). Britton Chance at the University of Pennsylvania in the US has proposed the determination of the presence of cerebral hematoma by comparing the difference between the optical contrast between of two positions in the brain ([@b0090]). This research further extends this method. Based on the correlation between the differences in near-infrared light density (ΔOD), the degree of hematoma can be analyzed. However, the thicknesses of scalp and skull vary among different individuals due to different growing-up environments, age, race, gender, and other factors. As the thickness of scalp and skull varies, the location of the hematoma changes ([@b0040]). Therefore, rational selection of the position and number of detectors is essential for reducing the impact of inter-subject variation and improving the accuracy of the model.

Some researchers have used the tMCimg method to simulate the relationship between the effective depth and the distance between detector and the light source. They concluded that the effective detection depth is twice the distance between the light source and the detector ([@b0120]). Strangman et al. analyzed the effect of the changes in the thickness of scalp and skull on the selection of distances between the light source and the detector. Methods such as the Uninformative Variable Elimination (UVE), Synergy Interval Regression and Genetic Algorithm have wide applications in variable selection. The variable importance in the projection (VIP) method examines the importance of independent variables in modeling. This method has been widely applied in fields like epidemiological analysis ([@b0015]), remote sensing ([@b0135], [@b0105]), biochemical analysis ([@b0025]), and blood component testing ([@b0045]). In this study, it is noted that signals detected by the detectors have multiple linear dependencies. Thus, by using the VIP analysis technique, the detectors with strong predictive and interpretative abilities for hematoma can be selected. Through the VIP analysis, the number of detectors was reduced from 30 to 4 in this study and the model's prediction capability was improved ([@b0095]). This study introduced a novel idea for the selection of distance between the near-infrared light source and the detectors, upon which clinical application of hematoma prediction can be applied based on differential near-infrared optical density.

2. Materials and methods {#s0010}
========================

2.1. VIP analysis based on the partial least-squares regression {#s0015}
---------------------------------------------------------------

VIP is an auxiliary analysis technology based on the partial least squares method. It can be used for the determination and selection of important independent variables. When the correlation between variables is strong, it can describe the explanatory power of independent variables to dependent variables through the comprehensive principal components of the relevant independent variables. It also can be used to screen the independent variables based on their VIP abilities and to analyze the explanatory power of each independent variable on the dependent variables after preliminary PLS analysis. For example, assuming there is a dependent variable y and the independent variables are x~1~, x~2~, ..., xk. Eq. [(1)](#e0005){ref-type="disp-formula"} gives the VIP value of independent variable number j:$$\text{VIP}_{\text{j}} = \sqrt{\frac{\text{k}}{\sum\limits_{h = 1}^{m}r^{2}(y\text{,}c_{h})}\sum_{h = 1}^{m}r^{2}(y\text{,}c_{h})w_{\mathit{hj}}^{2}}$$where k is the number of independent variables; c~h~ is the principal component extracted from the relevant independent variables; r (y, c~h~) is the correlation coefficient of the dependent variable and the principal component, representing the explanatory power of the principal component to y; and w~hj~ is the weight of the independent variable on principal component. Since the explanatory power of x~j~ to y is conveyed through the principal component c~h~, if the explanatory power of c~h~ to y is very strong and the effect of x~j~ to c~h~ is very significant, it can be considered that the explanatory power of x~j~ to y is strong.

The VIP value of each independent variable represents the explanatory power of the independent variable on the dependent variable, also known as the importance level of the independent variable for the modeling. If a reasonable VIP threshold is set for selecting high VIP-value independent variables and excluding low-valued independent variables, the accuracy of the model would not be affected much but the complexity of the model would be greatly reduced. This will a have great impact on the system speed, ease of use and cost ([@b0125]).

2.2. Differential absorption dural hematoma detection {#s0020}
-----------------------------------------------------

[Fig. 1](#f0005){ref-type="fig"} shows the diagram of brain hematoma detection using multi-channel differential absorption spectrophotometry. In the diagram, S and S\' are the incident light source. D~1~-D~n~ and D~1~′-D~n~\' are the light detectors placed at equal interval space respectively. The detection was conducted on symmetrical positions of the human head ([@b0100]). First, the detection was conducted on the right side of the head (R site) once, and the light intensities detected by n detectors were I~1~-I~n~ respectively. Then, detection was performed on the left side of the head (L site) and the light intensities detected by the n detectors were I~1~′-I~n~\' respectively. I~0~ was the incident light intensity from both the left and the right sides. By using Eqs. [(2)](#e0010){ref-type="disp-formula"}, [(3)](#e0015){ref-type="disp-formula"}, the light absorbance at each location OD~1~-OD~n~ and OD~1~′-OD~n~\'can be obtained respectively.$$\mathit{OD}_{i} = \lg\left( \frac{I_{0}}{I_{i}} \right)\text{,}\quad(i = 1\text{,}2\text{,}\ldots\text{,n})$$$$\mathit{OD}_{i}^{\prime} = \lg\left( \frac{I_{0}}{I_{i}^{\prime}} \right)\text{,}\quad(i = 1\text{,}2\text{,}\ldots\text{,n})$$Fig. 1The structure of brain model.

By dividing the light absorbance from symmetrical locations, the differential absorbance can be calculated.$$\Delta\mathit{OD}_{i} = \frac{\mathit{OD}_{i}}{\mathit{OD}_{i}^{\prime}} = \frac{\lg\frac{I_{0}}{I_{i}}}{\lg\frac{I_{0}}{I_{i}^{\prime}}} = \lg\frac{I_{i}^{\prime}}{I_{i}}\text{,}\quad(i = 1\text{,}2\ldots\text{,}5)$$

According to the anatomy of the human brain, the brain model consists of 5 layers, namely, the layers of scalp, skull, cerebrospinal fluid (CSF), grey matter, and white matter ([Fig. 1](#f0005){ref-type="fig"}). The brain structure and optical parameters at 840 nm wavelength are shown in [Table 1](#t0005){ref-type="table"} ([@b0055], [@b0115]). The definitions of the brain parameters are as follows: n is the refractive index; *μ*~a~ (cm^−1^) is the absorption coefficient; *μ*~s~ (cm^−1^) is the scattering coefficient; g is the anisotropic factor; d (cm) is the tissue thickness.Table 1Brain model optical parameters.Brain Layer*nμ~a~* (cm^−1^)*μ~s~′* (cm^−1^)*gd* (cm)Scalp1.450.2118.10.9d1Skull1.450.1915.20.9d2CSF1.450.052.30.90.2Grey matter1.450.4220.90.90.4White matter1.450.1786.50.93.4

Monte Carlo simulation was used to establish the brain model in this study. Monte Carlo simulation can describe the transmission path of photons within any tissue structure. It is considered as the simulation most close to reality and is known as the "golden standard" in describing photon transmission trajectory in biological tissue ([@b0065]). The simulated photon number in this study was set as 10^8^. According to the study by Strangman et al., the thickness of the human scalp is 6.9 ± 3.6 mm and the thickness of the skull is 6.0 ± 1.9 mm. In this study, the thickness of the scalp and skull were set as 1.3 cm, of which the thickness of the scalp was 7.0 mm and the thickness of skull was 6.0 mm ([@b0110]).

Clinical cerebral hematoma can be categorized as subdural hematoma, epidural hematoma, intracerebral hematoma and subarachnoid hemorrhage. Most hematomas caused by traumatic brain injury are positioned extradurally or intradurally under the skull, which is located at the third level, i.e. the cerebrospinal fluid level in the brain model ([@b0030], [@b0005]). When traumatic brain hematoma occurs, the absorption coefficient of the dural position will increase significantly. Clinical studies have shown that the hematoma absorption coefficient increases over 10 times, while the normal dural position absorption coefficient is 0.05 cm^−1^. In order to simulate different degrees of hematoma, in this study, the dural position *μ*~a~ was set at a range of 0.5--1.5 cm^−1^ and the adjustment step length was 0.1 cm^−1^.

2.3. Data processing and VIP screening of detector position {#s0025}
-----------------------------------------------------------

### 2.3.1. Data processing {#s0030}

According to Monte Carlo simulation, the luminous fluxes from 30 detectors placed at different radial positions from the light source were obtained. The ΔOD values of the detectors located within a 2.0--5.0 cm range with a 0.1 cm interval were calculated using Eq. [(4)](#e0020){ref-type="disp-formula"}. White noise is the random noise evenly distributed throughout the frequency spectrum with similar energy density and it exists in every detector. In this study, white noise was added to the ΔOD value in each detector according to the signal-to-noise ratio of the spectrometer in order to verify the robustness of the model. The USB2000+ fiber optic spectrometer from Ocean Optics is the most classic and popular spectrometer. It is suitable for many types of scientific studies and industrial applications. Its signal-to-noise ratio is 250:1. This ratio was used as a basis for the added white noise in this study. The PLS model was built using ΔOD values with added white noise from the detectors at each position as the input and the dural position *μ*~a~ as the output.

### 2.3.2. VIP screening of detector position {#s0035}

The VIP values of the difference in luminous flux ΔOD of detectors at different position versus the hematoma degree y were calculated based on the definition of the importance of variable projection of Eq. [(1)](#e0005){ref-type="disp-formula"}. The result is shown in [Fig. 2](#f0010){ref-type="fig"}.Fig. 2Radial distribution map of VIP.

The importance of each detector position on modeling can be viewed clearly using the VIP radial distribution map between 2.0 cm and 5.0 cm. Based on the characteristics of the brain structure as well as the VIP values, the chosen screening criterion for the detector position was that the normalized VIP value should be greater than 0.14. As shown in the radial distribution map of VIP in [Fig. 4](#f0020){ref-type="fig"}, the detector positions at 2.1 cm, 2.4 cm, 3.4 cm and 4.2 cm met the screening criteria. After VIP screening, the ΔOD values with the added white noise from the 4 detector positions were used as the input and the dural position *μ*~a~ as the output to establish the PLS model. The model verification method used was Leave-One-Out Cross Validation (LOOCV) and the model extraction fraction was 2.

3. Results and discussions {#s0040}
==========================

One experimental group and two control groups were used for this study. The experimental group was modeled using the VIP screening result. Control group 1 underwent PLS modeling using the ΔOD values from 30 detectors positioned at 2.0--5.0 cm with a 0.1 cm division as the input. Control group 2 underwent PLS modeling using the ΔOD values from 4 detectors located at 2.0 cm, 3.0 cm, 4.0 cm and 5.0 cm from the light source as the input and the dural position *μ*~a~ as the output. White noise was added to all models. The prediction results are shown in [Fig. 3](#f0015){ref-type="fig"}. [Fig. 3](#f0015){ref-type="fig"}(a) shows the results of the dural position *μ*~a~ prediction PLS model using the 4 selected VIP detectors. The model correlation was 99.48%; the mean error was 0.0200 cm^−1^; and the maximum error was 0.0414 cm^−1^. [Fig. 3](#f0015){ref-type="fig"}(b) shows the results of the dural position *μ*~a~ prediction PLS model of control group 1 using 30 detectors. The model correlation was 98.54%; the mean error was 0.0339 cm^−1^; and the maximum error was 0.0684 cm^−1^. [Fig. 3](#f0015){ref-type="fig"}(c) shows the results of the dural position *μ*~a~ prediction PLS model of control group 2 using 4 equally positioned detectors. The model correlation was 97.26%; the mean error was 0.0474 cm^−1^; and the maximum error was 0.0814 cm^−1^.Fig. 3*μ*~a~ predicated result with PLS. (a) *μ*~a~ predicated result with PLS model of experiment group. (b) *μ*~a~ predicated result with PLS model of control group 1. (c) *μ*~a~ predicated result with PLS model of control group 2.

The comparison of the relative errors (RE) among the three different model groups is shown in [Fig. 4](#f0020){ref-type="fig"}. The experiment group modeled via the PLS model used the ΔOD from 4 detectors with distances of 2.1 cm, 2.4 cm, 3.4 cm and 4.2 cm from the light source respectively to predict the dural position *μ*~a~. The RE was 2.06%. The control group 1 PLS model used the ΔOD from 30 detectors to predict the dural position *μ*~a~. The RE was 4.08%. Control group 2 PLS model used the ΔOD from 4 equally distanced detectors 2.0 cm, 3.0 cm, 4.0 cm and 5.0 cm from the light source to predict the dural position *μ*~a~. The RE was 5.68%. These results show that the model with the 4 detectors selected through the VIP method has the smallest RE. Therefore, satisfactory results of dural position *μ*~a~ prediction using VIP for the selection of detector position and number can be obtained.Fig. 4Map of RE.

4. Conclusion {#s0045}
=============

In this study, variable importance in the projection (VIP) method was used to simplify the model that uses differential near-infrared luminous flux for the detection of the position *μ*~a~. A highly accurate model with smaller relative error can be built using a smaller number of detectors using VIP screening for detectors with different radial distances from the light source. This study is important for the miniaturization and portability of brain hematoma detection device research. It plays an important role in the promotion of portable near-infrared brain hematoma detector applications in complex environments and provides an important reference and a new direction for the research on relationship between the location of detectors and source and the effective depth of photons.

The work described in this paper was fully supported by a grant from the Tianjin Science and Technology Project of China (No. 13ZCZDSY20100) and the National Science Foundation of Tianjin (No. 15JCQNJC14500), besides, the paper was supported by the Key Laboratory of Biomedical Effects of Nanomaterials and Nano-safety, CAS (No: NSKF201602).

Peer review under responsibility of King Saud University.
