Next Article in Journal
GridHTM: Grid-Based Hierarchical Temporal Memory for Anomaly Detection in Videos
Next Article in Special Issue
Spatial Smoothing Effect on Group-Level Functional Connectivity during Resting and Task-Based fMRI
Previous Article in Journal
An IoT-Based GeoData Production System Deployed in a Hospital
Previous Article in Special Issue
Lightweight Compound Scaling Network for Nasopharyngeal Carcinoma Segmentation from MR Images
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Impact of Anatomical Variability on Sensitivity Profile in fNIRS–MRI Integration

by
Augusto Bonilauri
1,*,
Francesca Sangiuliano Intra
2,
Francesca Baglio
3 and
Giuseppe Baselli
1
1
Department of Electronics, Information and Bioengineering, Politecnico di Milano, 20133 Milan, Italy
2
Faculty of Education, Free University of Bolzano-Bozen, 39042 Brixen, Italy
3
IRCCS Fondazione Don Carlo Gnocchi ONLUS, CADITER, 20148 Milan, Italy
*
Author to whom correspondence should be addressed.
Sensors 2023, 23(4), 2089; https://doi.org/10.3390/s23042089
Submission received: 23 December 2022 / Revised: 6 February 2023 / Accepted: 7 February 2023 / Published: 13 February 2023

Abstract

:
Functional near-infrared spectroscopy (fNIRS) is an important non-invasive technique used to monitor cortical activity. However, a varying sensitivity of surface channels vs. cortical structures may suggest integrating the fNIRS with the subject-specific anatomy (SSA) obtained from routine MRI. Actual processing tools permit the computation of the SSA forward problem (i.e., cortex to channel sensitivity) and next, a regularized solution of the inverse problem to map the fNIRS signals onto the cortex. The focus of this study is on the analysis of the forward problem to quantify the effect of inter-subject variability. Thirteen young adults (six males, seven females, age 29.3 ± 4.3) underwent both an MRI scan and a motor grasping task with a continuous wave fNIRS system of 102 measurement channels with optodes placed according to a 10/5 system. The fNIRS sensitivity profile was estimated using Monte Carlo simulations on each SSA and on three major atlases (i.e., Colin27, ICBM152 and FSAverage) for comparison. In each SSA, the average sensitivity curves were obtained by aligning the 102 channels and segmenting them by depth quartiles. The first quartile (depth < 11.8 (0.7) mm, median (IQR)) covered 0.391 (0.087)% of the total sensitivity profile, while the second one (depth < 13.6 (0.7) mm) covered 0.292 (0.009)%, hence indicating that about 70% of the signal was from the gyri. The sensitivity bell-shape was broad in the source–detector direction (20.953 (5.379) mm FWHM, first depth quartile) and steeper in the transversal one (6.082 (2.086) mm). The sensitivity of channels vs. different cortical areas based on SSA were analyzed finding high dispersions among subjects and large differences with atlas-based evaluations. Moreover, the inverse cortical mapping for the grasping task showed differences between SSA and atlas based solutions. In conclusion, integration with MRI SSA can significantly improve fNIRS interpretation.

1. Introduction

Functional near-infrared spectroscopy (fNIRS) is gaining relevance in the imaging of cortical activity in a variety of tasks [1,2,3,4]. Compared to functional MRI (fMRI), it allows for the performance of low-cost and ecological measurements of brain activity in an open environment, also minimizing discomfort to participants. Furthermore, fNIRS provides a direct estimation of both the concentration changes in oxygenated HbO 2 and deoxygenated hemoglobin HbR , hence providing deeper insight to the hemodynamic response to neural activation compared to the blood oxygen level-dependent signal of fMRI. Moreover, fNIRS is well-suited to multimodal integration with other neuroimaging techniques such as EEG [5] and MRI [6,7].
Acquisitions are based on a grid of optodes alternating near-infrared light sources and detectors, whose neighboring pairs form a measurement channel. However, fNIRS’ practical advantages and noninvasiveness requirements for clinical applications must face two major technological issues: (i) the need of accurate signal processing to reduce artifacts and other confounding factors of physiological origin [8,9,10,11]; and (ii) the fairly large source–detector distance needed to reach the cortex; hence limiting the spatial resolution of fNIRS. The present study addresses the latter issue by focusing on an fNIRS system with optodes placed according to the international 10/5 system [12]. In this case, the large source–detector distance (about 3–4 cm) has a lower spatial resolution compared to the structure of gyri, sulci, and functional areas to be investigated, thus supporting a wider application of methods integrating the fNIRS measures with the subject specific anatomy (SSA) of the cortical surface from an anatomical MRI scan.
An analysis of the potential impact of SSA integration is based on two working hypotheses: (i) the average sensitivity of channels may show sharp drops with cortical depth and lateral displacements [13,14]; (ii) the gyri and sulci pattern over cortical surface may cause a large subject-specific variability [15]. As result, the coupling between a cortical area and a channel may experience large variations, which should be considered when referring channel-wise fNIRS data to cortical areas.
The specific aim of this study is to quantify in a group of young healthy adults the impact of the above channel sensitivity drop and SSA influence. More generally, the driving reason is to highlight the potential impact of fNIRS integration with the SSA. Indeed, in clinical neurological studies, fNIRS might show easier applicability, particularly in longitudinal studies, compared to fMRI, which remains the actual gold-standard in functional evaluations [16]. Moreover, in neurological patients, an anatomical MRI is routinely available, thus permitting the SSA extraction and the fNIRS improved interpretation.
Several studies already address the problem of estimating fNIRS channels vs. cortical sensitivity (i.e., forward problem) using analytical approaches [17,18], or (as performed in the present study) using Monte Carlo simulations [19,20,21], which require the computation of the subject-specific optical profile. Additionally, the inverse problem of mapping the surface measures onto the cortical surface has intensively been studied [20]. Within this line of research, other studies have investigated methods for the optimal placement of optodes [22,23,24,25], and the performance of high density diffuse optical tomography (HD-DOT) systems to obtain reconstructions of cortical activation with spatial resolution comparable to fMRI [26,27,28,29,30,31,32].
Nevertheless, fNIRS–MRI integration is not systematically employed in clinical research and practice. Hence, in this study, we employ a standard 10/5 optode placement focusing on the forward problem addressing the factors influencing fNIRS sensitivity profile with respect to individual anatomical variability. Namely, a set of Monte Carlo simulations of light propagation according to MRI SSA is performed in a group of healthy young adults. We propose a numerical method to quantify the sensitivity profile with respect to channel displacements and scalp cortex distance (Section 2.3) and the coupling between channels and cortical parcellations of interest (Section 2.4). We also compared SSA results to the atlas-based anatomy (ABA) of three major atlases, which are often used in cases where the SSA was not available [19,33].

2. Materials and Methods

2.1. Dataset and Experimental Set-Up

At IRCCS Fondazione Don Carlo Gnocchi, 13 healthy young adults (6 males, 7 females, age 29.3 ± 4.3) underwent both fNIRS acquisitions and anatomical MRI scans. The experiment was approved by the IRCCS Fondazione Don Carlo Gnocchi ethical committee, and all volunteers provided their informed consent. An expert neuroradiologist examined each MRI to exclude the presence of any pathology.
Anatomical MRI scans were acquired using a 3T Siemens Prisma scanner and included a high-resolution, T1-weighted 3D image (MPRAGE, resolution = 0.8 × 0.8 × 0.8 mm3), which was processed using the Freesurfer package to extract scalp, skull, cerebrospinal fluid (CSF), grey (GM), and white matter (WM) surfaces.
The focus of this study is on the analysis of the forward problem and the impact of SSA variability, ahead of the regularization implied by the inverse solution. Nonetheless, application examples of fNIRS cortical mapping are shown on representative subjects performing a motor task. Data were acquired using a continuous-wave system at λ = 760 nm and λ = 850 nm wavelengths (NIRScoutX 32 × 32, NIRx Medizintechnik, Berlin, Germany) that employs 32 LED sources and 32 avalanche photodiode detectors placed according to the international 10/5 EEG system (Figure 1). Source and detector pairs were combined into a total of 102 measurement channels.
The functional task consisted of a block-designed motor-grasping task, where participants had to alternatively move their left or right hand [8]. The fNIRS signals were preprocessed using the Nirstorm package (https://github.com/Nirstorm/nirstorm (accessed on 5 February 2023)) in Brainstorm [34] to extract the relative changes in oxygenated Δ HbO 2   and deoxygenated Δ HbR hemoglobin concentrations, denoising, and artifact removal. Namely, fNIRS data were firstly preprocessed to remove channels with a coefficient of variation greater than 10%; raw intensity signals were converted into optical density variations, corrected for motion artifacts using the temporal derivative distribution repair (TDDR) algorithm [35], and bandpass-filtered at (0.01–0.08) Hz.
The fNIRS protocol required 40 to 45 min on average, including the placement of the optodes, quality control of channels, and signal acquisition during the considered task, while the MRI protocol average 15 min required on. The pre-processing of fNIRS and MRI data were not time consuming. Conversely, a significant computational burden is referred to—the sensitivity matrix computation over SSA (Section 2.2). Nonetheless, the overall experimental procedure is suitable for clinical applications since the computation of the sensitivity matrix can be performed offline after a routine anatomical MRI is acquired. Moreover, the sensitivity matrix can be reused for several fNIRS trials in a follow-up protocol in most cases where a stationary anatomical condition can be assumed.

2.2. Sensitivity Matrix Computation

Anatomical MRI data were processed in Nirstorm to compute the forward model at both wavelengths λ :
Δ O D λ = A λ Δ μ a λ
Matrix A λ elements represent the cortical surface-to-channel sensitivity coefficients at both wavelengths. The channel vector of optical density data is Δ O D λ = y i , with i = 1 ,   2 ,   ,   I and I = 102 number of channels in our configuration. The vector Δ μ a λ = Δ μ α , j λ , j = 1 , 2 , , J represents the absorption coefficient changes associated to the j t h cortical element CE j activity, which in turn is related to hemoglobin concentration changes as Δ μ α , j λ = α HbO 2 λ Δ HbO 2 + α HbR λ Δ HbR , where the α are the molar absorption coefficients of the respective chromophores at wavelength λ . The solution of this linear algebraic system is hence transferred from the surface channel level to the level of cortical elements CE j . Importantly, the CE j are automatically defined in Freesurfer by a vectorized mesh representation of about 15,000 elements for both subject-specific and atlas-based anatomies.
The optodes were placed on the scalp at the subject’s MRI-native coordinates and a virtual 10/5 cap model was co-registered to the specific anatomy, assuring the correspondence of the relevant landmarks. The resulting grid of sources (red) and detectors (green) is shown in Figure 1.
For the sole purpose of a quality check, all scalp locations and optode landmarks were registered to the Montreal Neurological Institute (MNI) space, thus observing the low dispersion of the clusters of corresponding optodes through subjects and atlases (see Supplementary Material for details).
Matrix A λ was estimated using Monte Carlo (MC) simulations in the native MRI space of SSAs, or in the atlas coordinates for ABAs. The forward problem was approached with GPU-accelerated MC simulations [36] with 5 × 10 7 photons per optode. A 5-layer head model was considered: scalp, skull, CSF, GM, and WM. The optical parameters shown in Table 1 were set from the work of Tak et al. [37] and Eggbrecht et al. [27]. Next, the channel-wise sensitivity profile was obtained from the respective source and detector fluencies [38], hence interpolated onto the pial surface (i.e., interface between GM and CSF, from now on simply referred to as cortical surface) using a Voronoi-based method [39] and smoothed with a 2 mm FWHM gaussian kernel.
The same analysis was applied to three atlases. This ABA approach had the dual aim of assessing the validity of applications in which the absence of an individual SSA is surrogated by an ABA [19,33] and also to provide three benchmarks to our SSA group dispersion analysis. The selected atlases were Colin27 [40] due to its wide adoption in fNIRS applications [14] and software [7,41,42], ICBM152 [43] as a standard MNI template, and FSAverage as the default brain template implemented in the Freesurfer package [44].
A qualitative view was provided by mapping the total sensitivity (TS) of all channels vs. cortical elements CE j . This was readily computed by summing the values in each column of matrix A :
T S λ = T S j λ = i = 1 I A i , j λ
The inverse problem (i.e., projecting surface measures onto the cortex) is ill-conditioned and undetermined, thus requiring heavy regularization. For the given application examples (see Section 3.4), we employed the solution based on the minimum norm estimate approach (for both λ , dropped for clarity), according to the work of Machado et al. [24]:
x ^ = A T h ^ 1 C 1 A + h ^ 2 C 2 1 A T h ^ 1 C 1 1 y
where h ^ 1 , h ^ 2 are regularization hyperparameters estimated through the restricted maximum likelihood method, while C 1 , C 2 are the covariance matrices, respectively, employed to model measurement noise and the a priori distribution of absorption changes.
Notably, Equation (1) indicates the size of regularization to solve for the ill-conditioned and undetermined problem presented in Equation (3), since estimating an order of 10.000 unknown cortical elements’ values from an order of 100 surface measures. Conversely, the present study focuses the forward sensitivity matrix A itself, based on the sole SSA (or atlas) and the simulated optical properties. Namely, two aspects are considered: (i) the drop in sensitivity with lateral displacement, statistically analyzed across channels and subjects by sensitivity displacement surfaces (SDS, Section 2.3); (ii) the variability of sensitivity coefficients grouped by cortical areas into area to channel (A2Ch) sensitivities (Section 2.4).

2.3. Sensitivity Displacement Surfaces

This section describes the proposed method to statistically assess the mean 3D profile of sensitivity with respect to displacements from the channel center (CC) at increasing cortical depths. For this purpose, the longitudinal γ (i.e., in the source to detector direction) and transverse τ (i.e., orthogonal direction to γ ) curvilinear distances from the CC were defined for each i t h channel considering the scalp curvature specific to each subject. The depth d was defined as the distance from a given γ ,   τ scalp location along the local normal to the scalp down to a cortex element C E j , for which sensitivity A i , j was also recorded. Since the cortical element index j is univocally determined by the channel i and the displacement γ ,   τ , the sampled A i , j can be written as A i γ ,   τ . Hence, in a given subject (sbj) and in the frame of the i t h -channel, a depth map d i s b j γ ,   τ and a sensitivity map A i sbj γ ,   τ were derived. Coordinates γ ,   τ were sampled over a 40 × 40 mm2 field by 2 × 2 mm2 bins (with some approximation at borders, due to the curvilinear head shape).
To gain a statistical representation of the sensitivity profile A i sbj γ ,   τ , d , depth and sensitivity values from all 102 channels were superimposed by aligning the respective frames. Thus, for each γ ,   τ coordinate the relationship between depth and sensitivity could be extracted from 102 samples. A segmentation by depth was fixed according to the quartiles of the overall sampled depths in the subject: d Q s b j , with the quartile index Q = 1 ; 2 ;   3 ;   4 ranging from the least cortical depths (i.e., gyri) to the deepest ones (i.e., sulci). Thus, it was possible to map the average sensitivity profiles at progressive depths according to four sensitivity displacement Surfaces (SDS) averaged across all channels for a single subject:
S D S Q s b j γ ,   τ = mean d d Q s b j A i s b j γ ,   τ , d
To quantify the expected drop in sensitivity with depth, the integral value of the four SDS, named integral under the surface (IUS), were computed and expressed in percentage values:
I U S Q s b j % = I U S Q s b j I U S 1 s b j + I U S 2 s b j + I U S 3 s b j + I U S 4 s b j · 100
where Q = 1 ; 2 ;   3 ;   4 indicates the quartile of cortical depth.
The shape of the SDS curves was assessed by the full width at half-maximum (FWHM) along the longitudinal FWHM λ and transverse FWHM τ directions. Namely, FWHM λ and FWHM τ were computed by considering the profile of SDS curves along these directions and calculating the distance between points at half of maximum value. Values for the first and second quartiles of cortical depth are reported, since almost flat SDSs were found at further depths.

2.4. Area to Channel Sensitivity

The interpretation of surface fNIRS data are mainly addressed to detect the activations of cortical areas expected to be involved during the execution of a task or in response to specific stimuli. In this perspective, the cortical elements belonging to each area may be grouped as subset of cortical elements   C E j composing an a r e a a (i.e., cortical parcellation), which are identified by the subset of indexes J a j   |   C E j a r e a a , with a = 1 ,   2 N a r e a s . Accordingly, the columns of A λ relevant to the same a r e a a are summed up into the a t h -column of a reduced sensitivity matrix B λ :
B λ = B i , a λ = j J a A i , j λ
The coefficients B i , a λ express the value of area to channel (A2Ch) sensitivities. This reduced model is adopted here to evaluate the sensitivity profile of a parcellation a r e a a due to the specific i t h -channel and across SSA vs. ABA cases.
The final A2Ch values were normalized using the sum of channel-wise sensitivity coefficients across cortical elements and expressed as percentage ratio.
A 2 C h % i , a = B i , a A i = j J a A i , j j = 1 J A i , j
As results, A 2 C h % i , a values span in the [0, 1] range and express the relative influence of an area over a measurement channel (i.e., area to channel coupling) and its variability through subjects.
The A 2 Ch % values were analyzed on a representative set of Brodmann areas (BA) (Figure 2A–D panels): bilateral BA1–2–3–4 as primary sensorimotor cortex (SensoriMotor Network, SMN), BA17 as primary (Visual Network 1, VIS1), and BA18-19 as secondary visual areas of occipital cortex (Visual Network 2, VIS2), BA46-9 as high-level, and BA45-47 as low-level sites of executive functions of dorsolateral– (PFC1) and ventrolateral–prefrontal cortex (PFC2), respectively. The correspondence between BAs with anatomical areas, along with the notation employed in the results section, is presented in Table 2.
We also considered the bilateral precentral and postcentral gyri as additional anatomical parcellations on the Destrieux atlas [45] (Figure 2E–G panels). Separate analyses of the selected BAs and Destrieux areas were performed. In both cases, A 2 Ch % values above a threshold of 0.2 are presented and considered as significant channel-wise values (i.e., at least of 20% of the sensitivity profile of the considered area is explained by that channel), while empty channels (i.e., those distant from any considered areas) are omitted.

3. Results

3.1. Sensitivity Maps

Figure 3 provides an example of sensitivity profile of channel S9D7 (i.e., the representation of A i , j i = S 9 D 7 at source 9 and detector 7), while the red dot represents the S9D7 midpoint as channel position. Considering the log10 scale down to 1/100 (i.e., blue in the color scale), it is evident that the sensitivity profile drops in a short range of depth from the top of gyri to the bottom of sulci. Moreover, there is a noticeable change in the sensed pattern using SSAs (A–C panels) and ABAs (D–F panels), whose variations can be better-quantified with A 2 Ch % indexes (Section 3.3).
The total sensitivity T S j maps on the cortical surface are shown in Figure 4 for three subjects and the three atlases. Values are normalized to the maximum sensitivity (i.e., red color scale) and in log 10 · scale down to 100-fold less (i.e., blue color scale). The qualitative inspection of these maps anticipates the quantitative results of the following analyses. The most evident feature is that the drop in sensitivity with depth shown above for a single channel (S9D7, as example) is maintained, also summing up all channel sensitivities. This confirms that the signal recorded on the scalp surface is almost entirely determined by the gyri. Secondly, along the sagittal plane, the inflation shows a part of the mesial cortex, which is a gray color scale since it is below the considered sensitivity range, hence negligible. Finally, similarities and differences among SSAs (A–C panels) and among atlases (D–panels) are highlighted. Overall, the maps presented in Figure 3 and Figure 4 attain to MC simulation at λ = 760 nm; however, λ = 850 nm provided similar qualitative information.

3.2. Sensitivity Displacement Surfaces

Table 3 provides a summary of the percentage IUS values (median and IQR) at increasing quartiles of cortical depth. Results quantify how the sensitivity profile drops with depth. Additionally, Table 4 presents the depth separation values between each subsequent quartile: d Q 1 / Q 2 , d Q 2 / Q 3 , and d Q 3 / Q 4 , respectively. The IUS 1 reflects the integral of the sensitivity profile to cortical anatomy at the top of gyri and shows that the sensitivity profile is mostly confined to this cortical range. Values for SSAs had a median of IUS 1   = 0.39, while lower values are found for most ABAs: 0.3 for Colin27, 0.36 for ICBM152. Conversely, the IUS 1 over the sensitivity profile of FSAverage is comparable to the one of SSAs.
Considering the sum of the two upper quartiles IUS 1 + IUS 2 values of 0.683 were found for SSAs, 0.596 for Colin27, 0.656 for ICBM152 and 0.681 for FSAaverage, which further confirms that the top of the sulcogyral pattern should be considered in the interpretation of fNIRS (i.e., almost of 70% of the sensitivity profile is entirely confined within this cortical range). Moreover, Table 4 highlights the steep drop of sensitivity with cortical depth. The quartile separation values in SSAs present a median decrease of about 2 mm per quartile (i.e., 11.8 mm at d Q 1 / Q 2 , 13.6 mm at d Q 2 / Q 3 , 15.7 mm at d Q 3 / Q 4 ).
Considering the shape of the SDS plots, it is possible to evaluate the longitudinal and transversal sensitivity decay. Figure 5 displays this analysis relevant to two SSAs. The drop in amplitude with the depth quartile (from SDS 1 to SDS 4 ), which is quantified by the IUS% indexes, is visually confirmed. The sensitivity profile with displacement reveals two major features: (i) high sensitivity is kept for a wide portion of the longitudinal γ displacement, yet with a drop approaching the source and the detector point; (ii) a sensibly steeper drop is seen in the transverse direction τ .
These features were common to all SSAs and are quantified in Table 5 by the full width at half-maximum (FWHM) along the longitudinal FWHM γ and transverse FWHM τ directions, respectively. Values for S D S 1 and S D S 2 are reported since they justify most of the sensitivity. The remarks advanced about Figure 5 are confirmed: (i) FWHM γ covers about 21 mm of the longitudinal axis, which is about 60% of the source-detector average distance; (ii) the sensitivity profile of SSAs is sensibly confined along the longitudinal direction by the steeper transverse drop with FWHM τ about threefold narrower than FWHM y .
Conversely, this result is not valid for the ABAs case, FWHM τ , 1 and FWHM τ , 2 being almost comparable or even higher than FWHM γ , 1 and FWHM γ , 2 , respectively. Therefore, this result suggests that ABAs provide a more uniform drop in the sensitivity profile along longitudinal and transverse direction, possibly due to the smoother anatomy resulting from merging many subjects.

3.3. Area to Channel Sensitivity

The A 2 Ch % coefficients are evaluated here, considering the sample areas defined in Section 2.4. Variability is presented by the SSA group statistics (median (IQR) and compared to the values of the three ABAs. Since the result shown in Section 3.2 displayed negligible differences in sensitivities at λ = 760 nm and 850 nm, the following analyses will be shown only for the former wavelength. Table 6 and Table 7 report the results of A 2 Ch % -coefficients for the selected Destrieux areas, while Table 8, Table 9, Table 10, Table 11 and Table 12 report functional ROIs associated with selected BAs. In both cases, results for SSAs are reported as median (IQR) values across subjects.
To limit the size of these tables, channels displaying an A 2 Ch % < 0.2 in all columns (i.e., SSA median and ABA values) were omitted. Despite this cutoff, the tables show that the considered areas are sensed by a fairly long list of channels, generally around 10. This indicates that the contribution of many channels enters with varying weights the quantification of the activation of a specific cortical area provided by the inverse problem.
Remarkably, in several areas, one or more channels were found with a median A 2 Ch % close to or above 0.5, hence indicating that those channels prevalently sensed that area. However, in many of these cases, a wide IQR was found. Hence, interpreting that channel as primarily indicative of that area is only worth it for some of the subjects.
SSA results were also compared to ABAs. The tables highlight A 2 Ch % values for which the SSA median differed by more than 0.2 from any one of the ABA cases in bold. Such differences, when largely present, indicate that surrogating the SSA with an ABA might provide misleading results when mapping surface data to cortical anatomy. Namely, the highest rate of SSA vs. ABA differences (i.e., highest number of highlighted table rows) is found around the VIS regions (Table 8 and Table 9), while the lowest rate around the PFC regions (Table 11 and Table 12). An intermediate case was represented by motor and sensorimotor regions, as highlighted in Table 6 and Table 7 for Destrieux areas and Table 10 for the SMN region. Although the analysis was limited to few representative ROIs, it is worth noting that the IQR range of A 2 Ch % coefficients across SSAs is often comparable or higher than the median value, hence indicating that there is an intrinsic variability regarding the identification of cortical anatomy.

3.4. Example of Cortical Image Reconstruction

Figure 6 provides an example of cortical image reconstruction of the block-averaged precentral and postcentral gyri response to the right-hand motor task in a representative subject (i.e., subject #9). Surface measures, averaged from 10 repetitions, where projected according to the own SSA (A–C panels) and compared with those resulting from employing the three considered atlases: Colin27 (D–F panels), ICBM152 (G–I panels), and FSAaverage (J–L panels). The color maps represent the peak Δ HbO 2 and Δ HbR values (left and mid column, respectively). The averaged responses are shown in the right column relevant to the contralateral precentral and postcentral gyri (integral in the areas). As expected, the time course shapes were negligibly affected by the adopted anatomy, since they are bound to the task paradigm and to the hemodynamic response function. Conversely, large differences were seen in their amplitudes, in keeping with the differences in the cortical maps.
Namely, the peak of Δ HbO 2 response (650 μ mol / L ) was found in the precentral gyrus if the correct SSA was applied. Conversely, the main activation was shifted on the postcentral gyrus for the ABAs, with large amplitude differences: 450 μ mol / L for Colin27, 980 μ mol / L for ICBM152, and 700 μ mol / L for FSAverage. Importantly, the ratio between the postcentral and the precentral gyrus shows sensible differences, being about 2:1 Colin27 and FSAverage, while up to 3:1 for ICBM152.
This result confirms that the anatomy considered in the inverse problem solution strongly influences the localization and quantification of functional responses, which is explained by the forward problem analyses. Namely, the above results relevant to channel S9D7 can be revisited under this perspective. Qualitatively, in Figure 3, by the sensitivity patterns on different anatomies. Quantitively, by the variable A2Ch% values for the Destrieux left precentral and postcentral gyri in Table 6 and for the L-SMN in Table 10.

4. Discussion

In this work, we showed an approach to (i) statistically assess the distribution of fNIRS sensitivity profile with respect to channel displacements and scalp to cortex distance; (ii) quantify the coupling between channels and selected functional/anatomical cortical areas (i.e., A2Ch% sensitivity). The proposed method addressed the multimodal integration of fNIRS and MRI techniques to further support the applicability of fNIRS in clinical research by alleviating the problems related to its limited spatial resolution and indirect sensing of cortical areas. Indeed, clinical applications would require a specific identification of cortical activations in patients with impaired neurovascular activity. The pairing of fNIRS with other recognized clinical neuroimaging techniques, such as anatomical MRI, is essential for progress in this direction while simultaneously taking advantage of its portability and potential applicability in various clinical contexts. For these reasons, our analysis was limited to multichannel fNIRS equipment with optodes placed at 10/5 locations, while the ongoing progresses relevant to more sophisticated instrumentation, optimized optode placement, and/or high-density probes were out of our scope.
Results for SDSs and IUS values provided a numerical and visual outcome regarding the spatial sensitivity distribution of the employed fNIRS probes along longitudinal and transverse direction at varying cortical depths. Results showed the characteristic features of the sensing field in SSAs: (i) along the longitudinal axis (i.e., source–detector direction) the sensitivity profile preserved meaningful values along a high percentage of the entire source–detector distance (i.e., the FWHM γ is about 21 mm across both the first and second quartiles of depth), while along the transverse axis we found a steeper decay in sensitivity (i.e., FWHM τ is about 6 mm across both the first and second quartiles of depth); (ii) a very wide drop in sensitivity is seen with depth, such that the bottom of sulci was attenuated about 100fold compared to the top of the gyri. This last result is quantified by the integral values of SDS plots according to quartiles of cortical depth. Approximately, the first depth quartile (i.e., top of gyri) provided about 40% of the sensitivity, the second about 30%, the third about 20%, and the fourth (i.e., bottom of sulci) about 10%. This result showed limited group dispersion across the SSAs, as quantified by low IQRs, since the SDSs were obtained by combining all channels, thus compensating the individual differences in the patterns of gyri and sulci. However, it should be stressed that this result refers to a group of healthy young adults with similar age span. Individual variability of fNIRS sensing with depth may conversely show even wider dispersions incase significant age spans were considered, or even neuropathologies implying brain atrophy [46,47,48].
The comparison of IUS values with the ABAs cases delivered similar results except for the Colin27 atlas, which provided an underestimation of the first quartile. In general, the proposed analysis is in line with recent works. Namely, Tian and Liu [49] proposed ICBM152 as reference for a depth compensation algorithm to improve the localization of functional activation. Liu et al. [50] adopted SSAs, indicating that fNIRS surface measurements can infer deep-brain activity comparable to simultaneous fMRI acquisitions. Strangman et al. [14,51] mapped the depth sensitivity profile over a Colin27 function of source–detector separation using MC simulations and quantified the influence of scalp and skull thickness. Moreover, the interpretation of fNIRS findings can be highly affected by anatomical variability and optode placement [15].
Regarding the A 2 Ch % analysis, in the fNIRS literature, it is a common practice to constrain scalp measurements to the cortical surface [21,33]. However, to the best of our knowledge, a numerical assessment of the forward problem and individual anatomic variability is usually overlooked, while directly assessing the inverse problem solution and regularization. The analysis was necessarily limited to a few Brodmann areas (BAs), yet representative of frontal, sensorimotor, and occipital regions. Separately, the precentral and postcentral gyri from the Destrieux atlas were also considered, given the different logic of this atlas based on anatomical features. To summarize the huge and sparse intersection between areas and channels, we provided a list of channels for each area, showing A 2 Ch % 20 % . The major result was the complex portrait of channels sensing each of the considered areas that provided information about the area activation while solving the inverse problem. Importantly, the SSA group analysis revealed high IQRs through the subjects. Finally, for many A 2 Ch % , the result provided by one or more of the three atlases sensibly differed from the SSA group median. The major deviation of the ABAs from the SSA median were found in the occipital visual regions, decreased in the sensorimotor region, and even more in the frontal region.
So far, fNIRS–MRI fusion is still seldom employed in clinical research and practice. Notwithstanding, the ecological measurements provided by fNIRS could support the evaluation of brain disorders [52,53], functional recovery in stroke [54], epilepsy [55], transcranial electrical stimulation [56], implementations of brain–computer interface and neurofeedback [57,58], monitoring of chronic neurological diseases [59], and the mapping of brain plasticity in rehabilitation [60,61,62].
The results of this study show that the further computational effort needed for cortical fNIRS mapping may correct the nonnegligible individual variability of the sensitivity of channels vs. cortical areas, thus improving interpretation via the inclusion of subject-specific anatomy. Indeed, future developments of this work will addressed the inverse problem of providing a more stable solution by reducing the number of sensitivity matrix coefficients. Namely, the proposed method will allow for the definition of a spatial prior on cortical depth to avoid the inversion of sensitivity matrix elements not significantly contributing to the sensitivity profile (i.e., the SDS analysis of Section 2.3) and/or that do not properly map a cortical area of interest (i.e., the A 2 Ch % analysis of Section 2.4).

5. Conclusions

We reconsidered the current problem of fNIRS image reconstruction by introducing methods for the numerical quantification of the sensitivity profile with respect to cortical depth and coupling with cortical areas of interest. This step requires the estimation of a sensitivity matrix of surface channels vs. cortical elements, whose dimension challenges qualitative and quantitative analyses. We also limited the analysis to a multichannel fNIRS equipment according to a 10/5 system configuration of optodes. The choice was derived from the need for providing concrete indications for the clinical research use of fNIRS. An indication was obtained using an analysis of the cortical surface sensed on average by each scalp channel, which is limited in the transverse direction (orthogonal to the source–detector line) and drops dramatically with depth. The sharp sensing volume and the subjective variability of gyri and sulci, even in a group of young healthy adults, provided high dispersion of the coupling between channels and cortical areas of major interest, as indicated by the area of channel sensitivity coefficients. The analysis performed on three major atlases showed that such standard anatomies may provide a limited surrogate for subject-specific anatomies according to different regions of interest. In conclusion, moving towards the application of fNIRS in clinical contexts would greatly benefit from the integration of fNIRS data with the subject-specific MRI anatomy as routine practice.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/s23042089/s1.

Author Contributions

Conceptualization, A.B., F.S.I., F.B. and G.B.; methodology, F.B. and G.B; software, A.B.; validation, A.B., F.S.I. and G.B.; formal analysis, A.B. and G.B.; investigation, A.B. and F.S.I.; resources, F.B.; data curation, A.B., F.S.I. and F.B.; writing—original draft preparation, A.B. and F.S.I.; writing—review and editing, F.B. and G.B.; visualization, F.B. and G.B.; supervision, F.B. and G.B.; project administration, F.B. and G.B; funding acquisition G.B. All authors have read and agreed to the published version of the manuscript.

Funding

This work was supported by the Italian Ministry of Health (“Ricerca Corrente” program 2022–2024).

Institutional Review Board Statement

This study was conducted according to the guidelines of the Declaration of Helsinki and approved by the Ethics Committee of Fondazione Don Carlo Gnocchi ONLUS (protocol code 26/2020/CE_FdG/FC/SA; date of approval 20 October 2020).

Informed Consent Statement

Informed consent was obtained from all subjects involved in the study.

Data Availability Statement

The data are not publicly available due to privacy issues.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Yücel, M.A.; Selb, J.J.; Huppert, T.J.; Franceschini, M.A.; Boas, D.A. Functional Near Infrared Spectroscopy: Enabling routine functional brain imaging. Curr. Opin. Biomed. Eng. 2017, 4, 78–86. [Google Scholar] [CrossRef] [PubMed]
  2. Boas, D.A.; Elwell, C.E.; Ferrari, M.; Taga, G. Twenty years of functional near-infrared spectroscopy: Introduction for the special issue. Neuroimage 2014, 85, 1–5. [Google Scholar] [CrossRef] [PubMed]
  3. Gramigna, V.; Pellegrino, G.; Cerasa, A.; Cutini, S.; Vasta, R.; Olivadese, G.; Martino, I.; Quattrone, A. Near-Infrared Spectroscopy in Gait Disorders: Is It Time to Begin? Neurorehabilit. Neural Repair 2017, 31, 402–412. [Google Scholar] [CrossRef]
  4. Cutini, S.; Moro, S.B.; Bisconti, S. Functional near Infrared Optical Imaging in Cognitive Neuroscience: An Introductory Review. J. Near Infrared Spectrosc. 2012, 20, 75–92. [Google Scholar] [CrossRef]
  5. Chiarelli, A.M.; Zappasodi, F.; Di Pompeo, F.; Merla, A. Simultaneous functional near-infrared spectroscopy and electroencephalography for monitoring of human brain activity and oxygenation: A review. Neurophotonics 2017, 4, 041411. [Google Scholar] [CrossRef] [PubMed]
  6. Forbes, S.H.; Wijeakumar, S.; Eggebrecht, A.T.; Magnotta, V.A.; Spencer, J.P. Processing pipeline for image reconstructed fNIRS analysis using both MRI templates and individual anatomy. Neurophotonics 2021, 8, 025010. [Google Scholar] [CrossRef]
  7. Aasted, C.M.; Yücel, M.; Cooper, R.; Dubb, J.; Tsuzuki, D.; Becerra, L.; Petkov, M.; Borsook, D.; Dan, I.; Boas, D. Anatomical guidance for functional near-infrared spectroscopy: AtlasViewer tutorial. Neurophotonics 2015, 2, 020801. [Google Scholar] [CrossRef]
  8. Bonilauri, A.; Intra, F.S.; Baselli, G.; Baglio, F. Assessment of fNIRS Signal Processing Pipelines: Towards Clinical Applications. Appl. Sci. 2021, 12, 316. [Google Scholar] [CrossRef]
  9. Brigadoi, S.; Ceccherini, L.; Cutini, S.; Scarpa, F.; Scatturin, P.; Selb, J.; Gagnon, L.; Boas, D.; Cooper, R. Motion artifacts in functional near-infrared spectroscopy: A comparison of motion correction techniques applied to real cognitive data. Neuroimage 2014, 85, 181–191. [Google Scholar] [CrossRef]
  10. Strangman, G.; Franceschini, M.A.; Boas, D.A. Factors affecting the accuracy of near-infrared spectroscopy concentration calculations for focal changes in oxygenation parameters. Neuroimage 2003, 18, 865–879. [Google Scholar] [CrossRef]
  11. Tachtsidis, I.; Scholkmann, F. False positives and false negatives in functional near-infrared spectroscopy: Issues, challenges, and the way forward. Neurophotonics 2016, 3, 031405. [Google Scholar] [CrossRef]
  12. Jurcak, V.; Tsuzuki, D.; Dan, I. 10/20, 10/10, and 10/5 systems revisited: Their validity as relative head-surface-based positioning systems. Neuroimage 2007, 34, 1600–1611. [Google Scholar] [CrossRef]
  13. Tsuzuki, D.; Dan, I. Spatial registration for functional near-infrared spectroscopy: From channel position on the scalp to cortical location in individual and group analyses. Neuroimage 2014, 85, 92–103. [Google Scholar] [CrossRef]
  14. Strangman, G.E.; Li, Z.; Zhang, Q. Depth Sensitivity and Source-Detector Separations for Near Infrared Spectroscopy Based on the Colin27 Brain Template. PLoS ONE 2013, 8, e66319. [Google Scholar] [CrossRef] [PubMed]
  15. Zhai, X.; Santosa, H.; Huppert, T.J. Using anatomically defined regions-of-interest to adjust for head-size and probe alignment in functional near-infrared spectroscopy. Neurophotonics 2020, 7, 035008. [Google Scholar] [CrossRef]
  16. Scarapicchia, V.; Brown, C.; Mayo, C.; Gawryluk, J.R. Functional magnetic resonance imaging and functional near-infrared spectroscopy: Insights from combined recording studies. Front. Hum. Neurosci. 2017, 11, 419. [Google Scholar] [CrossRef]
  17. Schweiger, M.; Arridge, S. The Toast++ software suite for forward and inverse modeling in optical tomography. J. Biomed. Opt. 2014, 19, 040801. [Google Scholar] [CrossRef]
  18. Dehghani, H.; Eames, M.E.; Yalavarthy, P.K.; Davis, S.C.; Srinivasan, S.; Carpenter, C.M.; Pogue, B.W.; Paulsen, K.D. Near infrared optical tomography using NIRFAST: Algorithm for numerical model and image reconstruction. Commun. Numer. Method. Eng. 2009, 25, 711–732. [Google Scholar] [CrossRef]
  19. Cooper, R.J.; Caffini, M.; Dubb, J.; Fang, Q.; Custo, A.; Tsuzuki, D.; Fischl, B.; Wells, W.; Dan, I.; Boas, D. Validating atlas-guided DOT: A comparison of diffuse optical tomography informed by atlas and subject-specific anatomies. Neuroimage 2012, 62, 1999–2006. [Google Scholar] [CrossRef] [PubMed]
  20. Tremblay, J.; Martínez-Montes, E.; Vannasing, P.; Nguyen, D.K.; Sawan, M.; Lepore, F.; Gallagher, A. Comparison of source localization techniques in diffuse optical tomography for fNIRS application using a realistic head model. Biomed. Opt. Express 2018, 9, 2994. [Google Scholar] [CrossRef]
  21. Boas, D.A.; Dale, A.M. Simulation study of magnetic resonance imaging–guided cortically constrained diffuse optical tomography of human brain function. Appl. Opt. 2005, 44, 1957. [Google Scholar] [CrossRef]
  22. Cutini, S.; Scatturin, P.; Zorzi, M. A new method based on ICBM152 head surface for probe placement in multichannel fNIRS. Neuroimage 2011, 54, 919–927. [Google Scholar] [CrossRef] [PubMed]
  23. Morais, G.A.Z.; Balardin, J.B.; Sato, J.R. fNIRS Optodes’ Location Decider (fOLD): A toolbox for probe arrangement guided by brain regions-of-interest. Sci. Rep. 2018, 8, 3341. [Google Scholar] [CrossRef] [PubMed]
  24. Machado, A.; Cai, Z.; Pellegrino, G.; Marcotte, O.; Vincent, T.; Lina, J.-M.; Kobayashi, E.; Grova, C. Optimal positioning of optodes on the scalp for personalized functional near-infrared spectroscopy investigations. J. Neurosci. Methods 2018, 309, 91–108. [Google Scholar] [CrossRef] [PubMed]
  25. Machado, A.; Marcotte, O.; Lina, J.M.; Kobayashi, E.; Grova, C. Optimal optode montage on electroencephalography/functional near-infrared spectroscopy caps dedicated to study epileptic discharges. J. Biomed. Opt. 2014, 19, 026010. [Google Scholar] [CrossRef]
  26. Wheelock, M.D.; Culver, J.P.; Eggebrecht, A.T. High-density diffuse optical tomography for imaging human brain function. Rev. Sci. Instrum. 2019, 90, 051101. [Google Scholar] [CrossRef]
  27. Eggebrecht, A.T.; White, B.R.; Culver, J.P. A quantitative spatial comparison of high-density diffuse optical tomography and fMRI cortical mapping. Neuroimage 2012, 61, 1120–1128. [Google Scholar] [CrossRef] [PubMed]
  28. Eggebrecht, A.T.; Ferradal, S.L.; Robichaux-Viehoever, A.; Hassanpour, M.S.; Dehghani, H.; Snyder, A.Z.; Hershey, T.; Culver, J.P. Mapping distributed brain function and networks with diffuse optical tomography. Nat. Photonics 2014, 8, 448–454. [Google Scholar] [CrossRef]
  29. Cui, X.; Bray, S.; Bryant, D.M.; Glover, G.H.; Reiss, A.L. A quantitative comparison of NIRS and fMRI across multiple cognitive tasks. Neuroimage 2011, 54, 2808–2821. [Google Scholar] [CrossRef]
  30. Zhang, X.; Toronov, V.Y.; Webb, A.G. Simultaneous integrated diffuse optical tomography and functional magnetic resonance imaging of the human brain. Opt. Express 2005, 13, 5513. [Google Scholar] [CrossRef]
  31. Cai, Z.; Uji, M.; Aydin, Ü.; Pellegrino, G.; Spilkin, A.; Delaire, É.; Abdallah, C.; Lina, J.-M.; Grova, C. Evaluation of a personalized functional near-infra-red optical tomography workflow using maximum entropy on the mean. Hum. Brain Mapp. 2021, 42, 4823–4843. [Google Scholar] [CrossRef] [PubMed]
  32. Hoshi, Y.; Yamada, Y. Overview of diffuse optical tomography and its clinical applications. J. Biomed. Opt. 2016, 21, 091312. [Google Scholar] [CrossRef] [PubMed]
  33. Custo, A.; Boas, D.; Tsuzuki, D.; Dan, I.; Mesquita, R.; Fischl, B.; Grimson, E.; Wells, W. Anatomical atlas-guided diffuse optical tomography of brain activation. Neuroimage 2010, 49, 561–567. [Google Scholar] [CrossRef] [PubMed]
  34. Tadel, F.; Baillet, S.; Mosher, J.C.; Pantazis, D.; Leahy, R.M. Brainstorm: A User-Friendly Application for MEG/EEG Analysis. Comput. Intell. Neurosci. 2011, 2011, 879716. [Google Scholar] [CrossRef]
  35. Fishburn, F.A.; Ludlum, R.S.; Vaidya, C.J.; Medvedev, A.V. Temporal Derivative Distribution Repair (TDDR): A motion correction method for fNIRS. Neuroimage 2019, 184, 171–179. [Google Scholar] [CrossRef]
  36. Fang, Q.; Boas, D.A. Monte Carlo Simulation of Photon Migration in 3D Turbid Media Accelerated by Graphics Processing Units. Opt. Express 2009, 17, 20178. [Google Scholar] [CrossRef] [PubMed]
  37. Tak, S.; Kempny, A.M.; Friston, K.J.; Leff, A.P.; Penny, W.D. Dynamic causal modelling for functional near-infrared spectroscopy. Neuroimage 2015, 111, 338–349. [Google Scholar] [CrossRef] [PubMed]
  38. Boas, D.A.; Culver, J.P.; Stott, J.J.; Dunn, A.K. Three dimensional Monte Carlo code for photon migration through complex heterogeneous media including the adult human head. Opt. Express 2002, 10, 159. [Google Scholar] [CrossRef]
  39. Grova, C.; Makni, S.; Flandin, G.; Ciuciu, P.; Gotman, J.; Poline, J.B. Anatomically informed interpolation of fMRI data on the cortical surface. Neuroimage 2006, 31, 1475–1486. [Google Scholar] [CrossRef]
  40. Collins, D.L.; Zijdenbos, A.P.; Kollokian, V.; Sled, J.G.; Kabani, N.J.; Holmes, C.J.; Evans, A.C. Design and construction of a realistic digital brain phantom. IEEE Trans. Med. Imaging 1998, 17, 463–468. [Google Scholar] [CrossRef]
  41. Huppert, T.J.; Diamond, S.G.; Franceschini, M.A.; Boas, D.A. HomER: A review of time-series analysis methods for near-infrared spectroscopy of the brain. Appl. Opt. 2009, 48, D280. [Google Scholar] [CrossRef] [PubMed]
  42. Santosa, H.; Zhai, X.; Fishburn, F.; Huppert, T. The NIRS Brain AnalyzIR Toolbox. Algorithms 2018, 11, 73. [Google Scholar] [CrossRef]
  43. Fonov, V.; Evans, A.C.; Botteron, K.; Almli, C.R.; McKinstry, R.C.; Collins, D.L. Unbiased average age-appropriate atlases for pediatric studies. Neuroimage 2011, 54, 313–327. [Google Scholar] [CrossRef]
  44. Fischl, B.; Sereno, M.I.; Tootell, R.B.H.; Dale, A.M. High-resolution intersubject averaging and a coordinate system for the cortical surface. Hum. Brain Mapp. 1999, 8, 272–284. [Google Scholar] [CrossRef]
  45. Destrieux, C.; Fischl, B.; Dale, A.; Halgren, E. Automatic parcellation of human cortical gyri and sulci using standard anatomical nomenclature. Neuroimage 2010, 53, 1–15. [Google Scholar] [CrossRef]
  46. Whiteman, A.C.; Santosa, H.; Chen, D.F.; Perlman, S.; Huppert, T. Investigation of the sensitivity of functional near-infrared spectroscopy brain imaging to anatomical variations in 5- to 11-year-old children. Neurophotonics 2017, 5, 011009. [Google Scholar] [CrossRef]
  47. Scholkmann, F.; Wolf, M. General equation for the differential pathlength factor of the frontal human head depending on wavelength and age. J. Biomed. Opt. 2013, 18, 105004. [Google Scholar] [CrossRef]
  48. Bonnéry, C.; Desjardins, M.; Pouliot, P.; Lesage, F.; Leclerc, P.-O.; Hoge, R.D.; Bherer, L. Changes in diffusion path length with old age in diffuse optical tomography. J. Biomed. Opt. 2012, 17, 056002. [Google Scholar] [CrossRef]
  49. Tian, F.; Liu, H. Depth-compensated diffuse optical tomography enhanced by general linear model analysis and an anatomical atlas of human head. Neuroimage 2014, 85, 166–180. [Google Scholar] [CrossRef]
  50. Liu, N.; Cui, X.; Bryant, D.M.; Glover, G.H.; Reiss, A.L. Inferring deep-brain activity from cortical activity using functional near-infrared spectroscopy. Biomed. Opt. Express 2015, 6, 1074. [Google Scholar] [CrossRef] [Green Version]
  51. Strangman, G.E.; Zhang, Q.; Li, Z. Scalp and skull influence on near infrared photon propagation in the Colin27 brain template. Neuroimage 2014, 85, 136–149. [Google Scholar] [CrossRef] [PubMed]
  52. Chen, W.-L.; Wagner, J.; Heugel, N.; Sugar, J.; Lee, Y.-W.; Conant, L.; Malloy, M.; Heffernan, J.; Quirk, B.; Zinos, A.; et al. Functional Near-Infrared Spectroscopy and Its Clinical Application in the Field of Neuroscience: Advances and Future Directions. Front. Neurosci. 2020, 14, 724. [Google Scholar] [CrossRef] [PubMed]
  53. Irani, F.; Platek, S.M.; Bunce, S.; Ruocco, A.C.; Chute, D. Functional Near Infrared Spectroscopy (fNIRS): An Emerging Neuroimaging Technology with Important Applications for the Study of Brain Disorders. Clin. Neuropsychol. 2007, 21, 9–37. [Google Scholar] [CrossRef]
  54. Yang, M.; Yang, Z.; Yuan, T.; Feng, W.; Wang, P. A Systemic Review of Functional Near-Infrared Spectroscopy for Stroke: Current Application and Future Directions. Front. Neurol. 2019, 10, 58. [Google Scholar] [CrossRef]
  55. Peng, K.; Pouliot, P.; Lesage, F.; Nguyen, D.K. Multichannel continuous electroencephalography-functional near-infrared spectroscopy recording of focal seizures and interictal epileptiform discharges in human epilepsy: A review. Neurophotonics 2016, 3, 031402. [Google Scholar] [CrossRef]
  56. Yang, D.; Shin, Y.-I.; Hong, K.-S. Systemic Review on Transcranial Electrical Stimulation Parameters and EEG/fNIRS Features for Brain Diseases. Front. Neurosci. 2021, 15, 629323. [Google Scholar] [CrossRef]
  57. Obrig, H NIRS in clinical neurology—A ‘promising’ tool? Neuroimage 2014, 85, 535–546. [CrossRef]
  58. Mihara, M.; Miyai, I. Review of functional near-infrared spectroscopy in neurorehabilitation. Neurophotonics 2016, 3, 031414. [Google Scholar] [CrossRef]
  59. Bonilauri, A.; Intra, F.S.; Pugnetti, L.; Baselli, G.; Baglio, F. A Systematic Review of Cerebral Functional Near-Infrared Spectroscopy in Chronic Neurological Diseases—Actual Applications and Future Perspectives. Diagnostics 2020, 10, 581. [Google Scholar] [CrossRef]
  60. Baglio, F.; Pirastru, A.; Bergsland, N.; Cazzoli, M.; Tavazzi, E. Neuroplasticity mediated by motor rehabilitation in Parkinson’s disease: A systematic review on structural and functional MRI markers. Rev. Neurosci. 2022, 33, 213–226. [Google Scholar] [CrossRef]
  61. Tavazzi, E.; Bergsland, N.; Cattaneo, D.; Gervasoni, E.; Laganà, M.M.; Dipasquale, O.; Grosso, C.; Saibene, F.L.; Baglio, F.; Rovaris, M. Effects of motor rehabilitation on mobility and brain plasticity in multiple sclerosis: A structural and functional MRI study. J. Neurol. 2018, 265, 1393–1401. [Google Scholar] [CrossRef] [PubMed]
  62. Allali, G.; Blumen, H.M.; Devanne, H.; Pirondini, E.; Delval, A.; Van De Ville, D. Brain imaging of locomotion in neurological conditions. Neurophysiol. Clin. 2018, 48, 337–359. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (AD) Graphical representation of NIRS source (S, red dots) and detector (D, green dots) positioning over scalp surface ((A)—right hemisphere view; (B)—left hemisphere view; (C)—superior view; (D)—occipital view).
Figure 1. (AD) Graphical representation of NIRS source (S, red dots) and detector (D, green dots) positioning over scalp surface ((A)—right hemisphere view; (B)—left hemisphere view; (C)—superior view; (D)—occipital view).
Sensors 23 02089 g001
Figure 2. (AD) Cortical parcellations according to Brodmann areas on Colin27 atlas: BA1-2-3-4 (green), BA17 (light violet), BA18-19 (orange), BA45-47 (red), and BA46-9 (dark violet). (EG) Cortical parcellations according to Destrieux on Colin27 atlas: precentral gyrus (dark blue) and postcentral gyrus (light blue). White dots indicate channel positions over scalp surface.
Figure 2. (AD) Cortical parcellations according to Brodmann areas on Colin27 atlas: BA1-2-3-4 (green), BA17 (light violet), BA18-19 (orange), BA45-47 (red), and BA46-9 (dark violet). (EG) Cortical parcellations according to Destrieux on Colin27 atlas: precentral gyrus (dark blue) and postcentral gyrus (light blue). White dots indicate channel positions over scalp surface.
Sensors 23 02089 g002
Figure 3. Graphical representation of channel-wise sensitivity profile placed over motor areas (i.e., S9D7) across three SSAs (AC) and ABAs ((D)—Colin27, (E)—ICBM152, (F)—FSAverage). The respective cortical parcellations of precentral (dark blue) and postcentral (light blue) gyri are also indicated. Sensitivity values are expressed as log 10 · adimensional units. Cortical surfaces were inflated at 40% for visualization purposes.
Figure 3. Graphical representation of channel-wise sensitivity profile placed over motor areas (i.e., S9D7) across three SSAs (AC) and ABAs ((D)—Colin27, (E)—ICBM152, (F)—FSAverage). The respective cortical parcellations of precentral (dark blue) and postcentral (light blue) gyri are also indicated. Sensitivity values are expressed as log 10 · adimensional units. Cortical surfaces were inflated at 40% for visualization purposes.
Sensors 23 02089 g003
Figure 4. Maps of total sensitivity in three SSAs (AC) and ABAs ((D)—Colin27, (E)—ICBM152, (F)—FSAverage) according to anteroposterior direction. Sensitivity values are expressed in log 10 · scale. Cortical surfaces are inflated at 40%. All maps refer to λ = 760 nm, while maps for λ = 850 nm are not shown since they display similar patterns.
Figure 4. Maps of total sensitivity in three SSAs (AC) and ABAs ((D)—Colin27, (E)—ICBM152, (F)—FSAverage) according to anteroposterior direction. Sensitivity values are expressed in log 10 · scale. Cortical surfaces are inflated at 40%. All maps refer to λ = 760 nm, while maps for λ = 850 nm are not shown since they display similar patterns.
Sensors 23 02089 g004
Figure 5. Graphical representation of SDS plots regarding two SSAs at increasing quartiles of cortical depth. (A,B) Lateral and top views of SDS plots regarding the first subject. (C,D) Lateral and top views of SDS plots regarding the second subject.
Figure 5. Graphical representation of SDS plots regarding two SSAs at increasing quartiles of cortical depth. (A,B) Lateral and top views of SDS plots regarding the first subject. (C,D) Lateral and top views of SDS plots regarding the second subject.
Sensors 23 02089 g005
Figure 6. Example of block-averaged image reconstruction in a representative subject (i.e., subject #9) associated to right-hand grasping. From top to bottom: SSA case (AC); Colin27 atlas (DF); ICBM152 atlas (GI); FSAverage atlas (JL). From left to right: image reconstruction of Δ HbO 2   μ mol / L concentrations at peak response (A,D,G,J); image reconstruction of Δ HbR   μ mol / L concentrations at peak response (B,E,H,K); mean profile of Δ HbO 2 and Δ HbR of cortical vertices at precentral and postcentral gyri according to Destrieux (C,F,I,L).
Figure 6. Example of block-averaged image reconstruction in a representative subject (i.e., subject #9) associated to right-hand grasping. From top to bottom: SSA case (AC); Colin27 atlas (DF); ICBM152 atlas (GI); FSAverage atlas (JL). From left to right: image reconstruction of Δ HbO 2   μ mol / L concentrations at peak response (A,D,G,J); image reconstruction of Δ HbR   μ mol / L concentrations at peak response (B,E,H,K); mean profile of Δ HbO 2 and Δ HbR of cortical vertices at precentral and postcentral gyri according to Destrieux (C,F,I,L).
Sensors 23 02089 g006
Table 1. Optical properties of the five-layered model employed in Monte Carlo forward problem (absorption coefficient μ α scattering coefficient μ s , anisotropy factor g, refraction index n).
Table 1. Optical properties of the five-layered model employed in Monte Carlo forward problem (absorption coefficient μ α scattering coefficient μ s , anisotropy factor g, refraction index n).
760 nm850 nm
μ a μ s μ a μ s gn
Scalp0.0176.7270.0195.8180.891.37
Skull0.01168.5450.01397.6360.891.37
CSF0.0042.7270.0042.7270.891.37
GM0.0187.5990.01926.1650.891.37
WM0.016710.8250.02089.1880.891.37
Table 2. Summary of employed left (L-) and right (R-) hemisphere areas and their correspondence between Brodmann areas and anatomical areas.
Table 2. Summary of employed left (L-) and right (R-) hemisphere areas and their correspondence between Brodmann areas and anatomical areas.
Functional LabelFunctional AreasAnatomical Areas
L-SMN/R-SMNBA1-2-3-4SensoriMotor Network
L-VIS1/R-VIS1BA17Visual Network 1
(Occipital cortex)
L-VIS2/R-VIS2BA18-19Visual Network 2
(Occipital cortex)
L-PFC1/R-PFC1BA46-9Prefrontal cortex 1
(Dorsolateral)
L-PFC2/R-PFC2BA45-47Prefrontal cortex 2
(Ventrolateral)
Table 3. Summary of percentage Integral Under the Surface (IUS) of SDSs at increasing quartile of cortical depth across all channels. Reported values are expressed as median (IQR) across all SSAs.
Table 3. Summary of percentage Integral Under the Surface (IUS) of SDSs at increasing quartile of cortical depth across all channels. Reported values are expressed as median (IQR) across all SSAs.
λ (nm) IU S 1 (%) IU S 2 (%) IU S 3 (%) IU S 4 (%)
SSA7600.391 (0.087)0.292 (0.009)0.215 (0.035)0.108 (0.044)
Colin270.2990.2700.2610.170
ICBM1520.3660.2900.2280.115
FSAverage0.3960.2850.2100.110
SSA8500.395 (0.089)0.293 (0.01)0.214 (0.033)0.106 (0.043)
Colin270.3000.2700.2600.169
ICBM1520.3670.2900.2270.115
FSAverage0.4000.2840.2080.108
Table 4. Separation values between depth quartiles. Results for SSA case are expressed as median (IQR) across all subjects.
Table 4. Separation values between depth quartiles. Results for SSA case are expressed as median (IQR) across all subjects.
d Q 1 / Q 2 (mm) d Q 2 / Q 3 (mm) d Q 3 / Q 4 (mm)
SSA11.8 (0.7)13.6 (0.7)15.7 (1.2)
Colin2714.10016.10018.000
ICBM15212.90014.90017.100
FSAverage17.20019.80022.795
Table 5. Full width at half-maximum of SDS plots along the longitudinal ( FWHM λ ) and transverse direction ( FWHM τ ) at the first and second quartiles of cortical depth. Results are reported as median (IQR) for SSAs case.
Table 5. Full width at half-maximum of SDS plots along the longitudinal ( FWHM λ ) and transverse direction ( FWHM τ ) at the first and second quartiles of cortical depth. Results are reported as median (IQR) for SSAs case.
λ (nm) FWH M γ , 1 (mm) FWH M τ , 1 (mm) FWH M γ , 2 (mm) FWH M τ , 2 (mm)
SSA76020.953 (5.379)6.082 (2.086)21.147 (5.091)6.187 (4.01)
Colin2729.18918.33121.05824.856
ICBM15228.23625.03424.54320.815
FSAverage18.73814.10820.74621.326
SSA85020.067 (5.4)6.082 (1.634)21.312 (4.718)6.332 (4.01)
Colin2728.77218.33121.05824.856
ICBM15228.23625.03424.54320.815
FSAverage18.29214.10820.52321.326
Table 6. Summary of A2Ch coefficients associated with left hemisphere Destrieux postcentral and precentral gyri. Results are reported for wavelength λ   = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 6. Summary of A2Ch coefficients associated with left hemisphere Destrieux postcentral and precentral gyri. Results are reported for wavelength λ   = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelG_postcentralLG_precentralL
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S11D310.766 (0.4)0.6570.5590.953
S11D70.291 (0.417)0.6140.5160.4850.431 (0.199)0.1070.1010.462
S11D90.682 (0.167)0.3150.0730.475
S7D7 0.042 (0.121)0.3090.4310.318
S8D5 0.196 (0.218)0.2570.2720.393
S8D80.157 (0.133)0.2580.2530.1850.429 (0.159)0.1440.1610.402
S9D50.022 (0.068)0.2030.2670.0570.343 (0.232)0.4870.5420.740
S9D70.121 (0.143)0.3000.6230.2440.576 (0.187)0.5670.3150.708
S9D80.477 (0.191)0.4710.2970.5520.357 (0.203)0.0090.0000.312
S9D90.689 (0.065)0.5910.1980.672
Table 7. Summary of A2Ch coefficients associated to right hemisphere Destrieux postcentral and precentral gyri. Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 7. Summary of A2Ch coefficients associated to right hemisphere Destrieux postcentral and precentral gyri. Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelG_postcentralRPrecentral
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S21D21 0.04 (0.159)0.3070.5310.201
S22D19 0.144 (0.154)0.2660.4450.438
S22D22 0.396 (0.145)0.2330.3000.341
S23D190.026 (0.05)0.0690.2620.0260.329 (0.135)0.6410.5400.655
S23D210.067 (0.103)0.2710.4800.0820.523 (0.237)0.5910.4790.771
S23D220.461 (0.16)0.5080.4760.5040.358 (0.129)0.2000.0370.332
S23D230.599 (0.257)0.3890.3490.6790.143 (0.086)0.0650.0040.248
S25D210.259 (0.23)0.5620.3810.3550.431 (0.184)0.2010.0900.507
S25D230.513 (0.356)0.3220.0590.673
S25D320.546 (0.333)0.7300.1900.7830.378 (0.457)0.0230.0000.217
Table 8. Summary of A2Ch coefficients associated to selected left hemisphere functional ROIs: primary (VIS1) and secondary visual network (VIS2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 8. Summary of A2Ch coefficients associated to selected left hemisphere functional ROIs: primary (VIS1) and secondary visual network (VIS2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelL-VIS1L-VIS2
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S13D12 0.033 (0.069)0.4060.3050.068
S13D13 0.298 (0.206)0.7560.8220.528
S14D14 0.422 (0.202)0.8200.7930.819
S15D130.100 (0.162)0.4560.4500.1680.873 (0.195)0.5440.5500.819
S15D140.087 (0.112)0.3630.2880.1650.883 (0.189)0.6370.7120.835
S15D150.527 (0.121)0.5700.4360.4170.478 (0.142)0.4300.5640.583
S15D300.589 (0.311)0.9610.9710.8000.353 (0.348)0.0350.0290.165
S16D14 0.988 (0.032)1.0001.0000.997
S16D15 0.973 (0.069)0.0000.0000.992
S31D150.413 (0.229)0.2690.0000.1130.56 (0.276)0.7311.0000.887
S31D300.729 (0.379)0.7690.6140.480
S32D13 0.731 (0.128)0.9630.9720.687
S32D300.092 (0.116)0.2960.2970.1280.581 (0.37)0.4660.3920.372
Table 9. Summary of A2Ch coefficients associated to selected right hemisphere functional ROIs: primary (VIS1) and secondary visual network (VIS2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 9. Summary of A2Ch coefficients associated to selected right hemisphere functional ROIs: primary (VIS1) and secondary visual network (VIS2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelR-VIS1R-VIS2
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S30D29 0.982 (0.034)1.0000.0000.992
S31D290.278 (0.122)0.3590.0000.2260.522 (0.251)0.6410.0000.769
S32D27 0.611 (0.189)0.8450.8550.798
S27D27 0.251 (0.178)0.8920.8630.545
S28D26 0.022 (0.022)0.3830.2510.297
S28D28 0.416 (0.269)0.8130.7220.911
S29D270.156 (0.220)0.6140.6470.2830.825 (0.192)0.3860.3530.714
S29D280.131 (0.148)0.3720.2640.2040.846 (0.167)0.6280.7360.796
S29D290.564 (0.173)0.6420.5090.3920.429 (0.215)0.3580.4780.608
S29D300.456 (0.153)0.7920.8610.7680.414 (0.374)0.0460.0670.172
S30D28 0.978 (0.118)1.0001.0001.000
Table 10. Summary of A2Ch coefficients associated to selected functional ROIs associated to sensorimotor network (SMN). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 10. Summary of A2Ch coefficients associated to selected functional ROIs associated to sensorimotor network (SMN). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelL-SMNR-SMN
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S11D310.956 (0.256)0.8670.8910.992
S11D70.581 (0.360)0.8840.8580.787
S11D90.922 (0.189)0.5860.3480.641
S12D80.532 (0.253)0.2580.3450.427
S12D90.482 (0.205)0.1770.0370.202
S8D80.210 (0.199)0.1860.3900.245
S9D50.104 (0.188)0.2920.5030.241
S9D70.266 (0.276)0.7260.8690.628
S9D80.756 (0.182)0.8770.9040.727
S9D90.913 (0.116)0.8620.4860.834
S22D22 0.227 (0.137)0.2110.3410.156
S23D19 0.081 (0.080)0.2310.3800.098
S23D21 0.270 (0.200)0.7630.8540.340
S23D22 0.743 (0.167)0.8610.9370.611
S23D23 0.904 (0.091)0.8590.7620.877
S25D21 0.555 (0.227)0.8870.8760.605
S25D23 0.933 (0.106)0.5890.4640.907
S25D25 0.446 (0.238)0.0990.0670.430
S25D32 0.863 (0.214)0.8600.8000.870
S26D22 0.701 (0.114)0.6230.6290.259
S26D23 0.474 (0.238)0.2240.0830.350
S27D23 0.280 (0.174)0.0640.0230.334
Table 11. Summary of A2Ch coefficients associated to selected left hemisphere functional ROIs: dorsolateral (PFC1) and ventrolateral prefrontal cortex (PFC2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 11. Summary of A2Ch coefficients associated to selected left hemisphere functional ROIs: dorsolateral (PFC1) and ventrolateral prefrontal cortex (PFC2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelL-PFC1L-PFC2
SSAColin27ICBM152FSAverage SSAColin27ICBM152FSAverage
S4D20.521 (0.156)0.5760.5130.375
S2D20.341 (0.25)0.2910.3630.211
S2D30.892 (0.108)0.8920.9060.796
S3D1 0.402 (0.171)0.4160.2680.355
S3D40.257 (0.196)0.2830.1970.2390.745 (0.186)0.7100.7970.729
S4D30.891 (0.174)0.8470.7741.000
S5D40.308 (0.091)0.2900.3190.2110.692 (0.102)0.6940.6650.782
S5D50.638 (0.087)0.7720.7850.4520.252 (0.077)0.1190.0780.405
S6D4 0.976 (0.096)0.9520.9851.000
S6D6 0.678 (0.177)0.4960.6341.000
S7D30.619 (0.382)0.4120.3600.612
S7D50.742 (0.194)0.5400.5820.621
S8D4 0.863 (0.2)0.7670.7290.763
Table 12. Summary of A2Ch coefficients associated to selected right-hemisphere functional ROIs: dorsolateral (PFC1) and ventrolateral prefrontal cortex (PFC2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
Table 12. Summary of A2Ch coefficients associated to selected right-hemisphere functional ROIs: dorsolateral (PFC1) and ventrolateral prefrontal cortex (PFC2). Results are reported for wavelength λ = 760 nm. Bold ABA values indicate channels presenting a 20% difference in A2Ch coefficients from the SSA median.
ChannelR-PFC1R-PFC2
SSAColin27ICBM152FSAverageSSAColin27ICBM152FSAverage
S17D170.87 (0.166)0.8620.9270.781
S17D20.284 (0.251)0.2110.3510.204
S18D16 0.253 (0.176)0.4040.2970.394
S18D180.383 (0.163)0.2900.1850.3910.615 (0.174)0.6980.8080.575
S19D170.983 (0.051)0.9460.9470.993
S19D180.356 (0.193)0.2060.2600.3300.637 (0.184)0.7820.7270.669
S19D190.730 (0.072)0.5720.6990.6280.202 (0.145)0.2500.1840.326
S20D18 0.929 (0.055)0.9350.9840.939
S20D20 0.761 (0.257)0.4790.4971.000
S21D170.602 (0.308)0.4530.3630.620
S21D190.619 (0.194)0.5400.4950.711
S22D18 0.864 (0.145)0.8410.8080.809
S22D19 0.202 (0.164)0.1600.1050.232
S22D20 0.283 (0.196)0.1410.0450.312
S23D190.416 (0.188)0.1900.1300.239
S4D170.783 (0.242)0.7540.7650.785
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Bonilauri, A.; Sangiuliano Intra, F.; Baglio, F.; Baselli, G. Impact of Anatomical Variability on Sensitivity Profile in fNIRS–MRI Integration. Sensors 2023, 23, 2089. https://doi.org/10.3390/s23042089

AMA Style

Bonilauri A, Sangiuliano Intra F, Baglio F, Baselli G. Impact of Anatomical Variability on Sensitivity Profile in fNIRS–MRI Integration. Sensors. 2023; 23(4):2089. https://doi.org/10.3390/s23042089

Chicago/Turabian Style

Bonilauri, Augusto, Francesca Sangiuliano Intra, Francesca Baglio, and Giuseppe Baselli. 2023. "Impact of Anatomical Variability on Sensitivity Profile in fNIRS–MRI Integration" Sensors 23, no. 4: 2089. https://doi.org/10.3390/s23042089

APA Style

Bonilauri, A., Sangiuliano Intra, F., Baglio, F., & Baselli, G. (2023). Impact of Anatomical Variability on Sensitivity Profile in fNIRS–MRI Integration. Sensors, 23(4), 2089. https://doi.org/10.3390/s23042089

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop