Open Access
5 July 2024 Flexible circuit-based spatially aware modular optical brain imaging system for high-density measurements in natural settings
Edward Xu, Morris Vanegas, Miguel Mireles, Artem Dementyev, Ashlyn McCann, Meryem Yücel, Stefan A. Carp, Qianqian Fang
Author Affiliations +
Abstract

Significance

Functional near-infrared spectroscopy (fNIRS) presents an opportunity to study human brains in everyday activities and environments. However, achieving robust measurements under such dynamic conditions remains a significant challenge.

Aim

The modular optical brain imaging (MOBI) system is designed to enhance optode-to-scalp coupling and provide a real-time probe three-dimensional (3D) shape estimation to improve the use of fNIRS in everyday conditions.

Approach

The MOBI system utilizes a bendable and lightweight modular circuit-board design to enhance probe conformity to head surfaces and comfort for long-term wearability. Combined with automatic module connection recognition, the built-in orientation sensors on each module can be used to estimate optode 3D positions in real time to enable advanced tomographic data analysis and motion tracking.

Results

Optical characterization of the MOBI detector reports a noise equivalence power of 8.9 and 7.3 pW/Hz at 735 and 850 nm, respectively, with a dynamic range of 88 dB. The 3D optode shape acquisition yields an average error of 4.2 mm across 25 optodes in a phantom test compared with positions acquired from a digitizer. Results for initial in vivo validations, including a cuff occlusion and a finger-tapping test, are also provided.

Conclusions

To the best of our knowledge, the MOBI system is the first modular fNIRS system featuring fully flexible circuit boards. The self-organizing module sensor network and automatic 3D optode position acquisition, combined with lightweight modules (18 g/module) and ergonomic designs, would greatly aid emerging explorations of brain function in naturalistic settings.

1.

Introduction

The widespread adoption of functional neuroimaging techniques has fundamentally advanced our understanding of how the human brain perceives, interprets, and responds to stimuli.1 However, brain activation exhibits complex patterns and dynamics, many of which are only apparent when studied in their native environments.2 This has posed a challenge for many established neuromonitoring techniques, such as functional magnetic resonance imaging (fMRI)3 and magnetoencephalography (MEG),4 which, despite tremendous technological innovation over recent years, still suffer from a key limitation—the need for subjects to remain relatively motionless within a highly confined space. Such limitations have greatly restricted our ability to understand the human brain in everyday situations and interactions.5,6 Although another commonly used neuroimaging technique, electroencephalography (EEG), complements fMRI/MEG with excellent portability, it is limited by relatively poor spatial resolution.5,7 These inherent disadvantages of gold-standard neuromonitoring modalities have motivated active exploration and development of new techniques that are more suited to measuring brain activities in naturalistic and dynamic environments.

In recent years, researchers have turned to functional near-infrared spectroscopy (fNIRS) to address this technological gap.8 fNIRS is an emerging functional neuroimaging technique that measures the hemodynamic changes associated with brain activity9 using safe, non-ionizing near-infrared light.8 With its relatively low cost and ability to probe a rich set of physiological parameters, fNIRS has seen rapid adoption in recent years6,10 and has been used in studies ranging from psychiatric conditions,11,12 language,1315 cognitive neurodevelopment,1619 stroke recovery,20 education,21 pain detection,22 and even brain–computer interfaces.2325 However, many early fNIRS studies employed mobile, cart-sized instruments that housed the source and detector units and optical fibers that transmit light to and from the subject’s scalp.2629 Although such fiber-based fNIRS systems are significantly more portable than conventional fMRI and MEG, the fragility of optical fibers and physical constraints upon subject mobility still preclude them from investigating brain activity in unrestricted environments.1

In response to this, fNIRS devices featuring improved wearability and portability have been actively developed over the past decade and deployed in both research and clinical settings.30 In particular, the rapid rise of modular fNIRS architectures,31,32 alongside a diverse range of lightweight and wearable fNIRS instruments,3339 represents a promising step toward extending fNIRS-based brain monitoring to address the aforementioned challenges. Modular fNIRS systems are wearable probes comprised of repeating source and detector circuits called modules.40 The use of repeating modules not only greatly simplifies system fabrication and lowers overall production costs but also offers improved flexibility in probe design to readily adapt to the diverse coverage requirements and regions of interest (ROIs) of different fNIRS experiments.41 In addition, wearable fNIRS devices have increasingly replaced cumbersome fiber optics with compact sources and detectors placed at the scalp surface to allow for greater optode density and more robust system designs. Modular fNIRS architectures with the ability to measure from both intra- and inter-module (optodes on different modules) source-detector (SD) pairs further improve spatial sampling density and overlapping channels, which have been shown to enhance both spatial resolution and depth specificity.38,42 However, compared with traditional “monolithic” fNIRS systems, modular headgear requires additional considerations, including optode layout design and modular connectivity, to optimize measurement coverage. Our group has recently released an open-source toolbox specifically addressing this need.43

While the transition toward modular and wearable fNIRS systems represents an important initial step toward achieving robust measurements of brain activity in naturalistic settings, additional challenges remain to be addressed. For example, the increased subject and headgear movement under dynamic experimental conditions can result in strong motion-related artifacts.5,44 While some of these artifacts can be detected and corrected using signal processing techniques45,46 or alleviated through mechanical approaches such as glue47 or spring-loaded optode designs,48 effectively and efficiently coupling sources and detectors to the head remains a key obstacle.46 In addition, model-based tomographic image reconstruction techniques using data from multiple channels overlapping the same area have been shown to improve the spatial localization accuracy, signal contrast, and measurement robustness of recovered brain activations.49,50 However, this approach requires three-dimensional (3D) optode positions51 at the very least and, in some cases, 3D anatomical scans. The typical approach for acquiring 3D optode positions involves manual measurements using 3D digitization systems, providing a snapshot of optode locations without accounting for movement or positioning changes over the course of an experiment.37,52 More recent solutions leveraging photogrammetry offer greater portability and constant sampling over time but still require a set of external cameras that restrict mobility to areas within the cameras’ field of view.38,51 Integrating automated, real-time measurements of 3D optode positions into a wearable fNIRS system would further assist motion correction and advance tomographic fNIRS to mainstream use.

To address the aforementioned challenges, we report here a lightweight, re-configurable, high-density fNIRS modular optical brain imaging (MOBI) system using head-conforming flexible/bendable circuits with integrated 3D optode position tracking capabilities.53 This ultra-compact and fiberless system adopts a diamond-shaped flexible printed circuit board (fPCB) module that improves the robustness of optode-scalp coupling.54,55 In addition, a dense peer-to-peer (P2P) network automatically determines module-to-module connectivity and topology, facilitating the implementation of spatial multiplexing strategies that increase data acquisition frame rates. Finally, each flexible MOBI module contains an inertial measurement unit (IMU), which can be used in combination with the automatically determined connection topology to produce real-time, 3D tracking of optode positions. This built-in position-sensing capability greatly simplifies tomographic fNIRS data acquisition, shortens setup times, and provides real-time metrics for assessing optode motion.

In Secs. 2.12.3, we provide an overview of the MOBI system, including technical 5 specifications for the MOBI modules and supporting hardware components as well as implementation strategies for novel features such as connection topology recognition, IMU-based optode positioning, and data acquisition speed improvement through spatial multiplexing. We further outline the characterization and validation protocols (Sec. 2.4) and results of the MOBI system (Sec. 3.1), head conformability of flexible- versus rigid-circuit probes (Sec. 3.3), accuracy of automatically determined locations (Sec. 3.4), and preliminary phantom and in vivo experiments (Secs. 3.1, 3.5, and 3.6).

2.

Methods

2.1.

Module Design

A MOBI module consists of a double-sided fPCB fabricated on polyimide-based films [Fig. 1(c)], with each module housing two detectors (OPT101, Texas Instruments, Dallas, Texas, United States) with built-in trans-impedance amplifiers and three dual-wavelength light-emitting diode (LED) sources at 735 and 850 nm (Marubeni, Chiyoda, Japan). The LEDs are driven by a single programmable constant current driver (LT3092, Analog Devices Inc., Wilmington, Massachusetts, United States) capable of delivering 0 to 100 mA of current coupled to a digital-to-analog converter (DAC, Microchip Technology Inc., Chandler, Arizona, United States) and controlled by a digital multiplexer (NX3L4051, NXP Semiconductors, Eindhoven, Netherlands) using a spatial multiplexing encoding strategy, as detailed in Sec. 2.4.2 (Fig. 2).

Fig. 1

Photos of the MOBI module showing the (a) top view without a silicone cover, (b) bottom view showing light guides attached to optodes and a black silicone enclosure, and (c) MOBI flexible circuit board before adding components.

NPh_11_3_035002_f001.png

Fig. 2

Schematic diagram of a single MOBI module. The microcontroller uses an internal I2C protocol to communicate with components on a single board. A P2P network allows communication between neighboring modules. Associated components are color-coded: red for sources, blue for detectors, green for the IMU, and yellow for the microcontroller.

NPh_11_3_035002_f002.png

The positions of the LEDs and detectors are labeled in Fig. 1(a). This arrangement produces one 8 mm, one 30 mm, and four 24.5 mm dual-wavelength channels within a single module (referred to as intra-module channels hereafter). In addition, a 3 mm diameter acrylic optical fiber (Edmund Optics, Barrington, New Jersey, United States) is attached to each source and detector, serving as a light guide to couple light to and from the scalp [Fig. 1(b)]. At the center of each module, we also mount a nine-axis IMU that tracks the module’s orientation (Bosch Sensortec, Kusterdingen, Germany). In addition, three low-dropout voltage regulators are used to control the power source voltage for the source, detector, and auxiliary measurement components. All optodes and sensors on the MOBI module are controlled by an onboard system-on-a-chip (BC832, Fanstel Corp., Scottsdale, Arizona, United States) with an integrated microcontroller (nRF52832, Nordic Semiconductor, Trondheim, Norway). The microcontroller also incorporates a 12-bit, 2-channel analog-to-digital converter (ADC) with eight built-in gain settings, which change the effective input range of the ADC, that samples from the detectors and IMU. Each optical measurement averages 175 samples, providing a single-channel, single-wavelength acquisition time of 7.5 ms and a modular frame rate, including all six dual-wavelength channels, dark measurements, and IMU measurements, of 11 Hz. An inter-integrated circuit (I2C) protocol is used to communicate between various components (Fig. 2).

The five optodes (3× LEDs and 2× detectors) are located on the side of the module facing the scalp, while all other electronic components are placed on the opposite side to assist with signal isolation, heat dissipation, and comfort. MOBI modules have a channel density of 0.56 dual-wavelength channels per square centimeter. We designed our module using a diamond or dual-equilateral-triangle shape with an edge length of 5 cm, a choice motivated by the superior tesselation of curved head surfaces with triangular meshes. Each module is encapsulated in a cast flexible black silicone enclosure to provide a smooth touch, blockage of stray light, and protection for electronic components during use. Each module has four flexible printed circuit (FPC) connectors that allow multiple modules to be connected using FPC cables [Fig. 1(a)]. A P2P serial network protocol allows each module to communicate with up to four connected neighboring modules. Five 8 mm diameter holes next to each optode are used for securing modules to a headgear and potentially hair adjustment. Each module, including all components, light guides, and the silicone cover, has a total weight of around 18 g.

2.2.

System Architecture

A compact master node, measuring 45×30  mm2, is used to connect to and control an arbitrary number of daisy-chained MOBI modules (Fig. 3). When the master node is powered by a laptop using a single universal serial bus (USB) cable, it can support up to 6× MOBI modules. An external power supply can be used instead when using more modules is desired. In addition, the master node utilizes I2C communication to supply power to and acquire data from each module [Fig. 3(a)]. It incorporates a microcontroller development board (Teensy 4.0, PJRC, Sherwood, Oregon, United States), a voltage regulator, an FPC connector, a 2-pin Japan Solderless Terminal (JST, Jihlava, Czech Republic) connector, and two switches. A micro-USB cable connects the master board to a computer for serial communication. The switches are used to manually toggle the power supply (between micro-USB and an external battery) to the modules, and the JST connector allows an external trigger board [Fig. 3(c)] to be connected for synchronizing auxiliary signals such as experimental triggers. This trigger board is based on a simple microcontroller (ATmega328P, Microchip Technology, Chandler, Arizona, United States) and transmits transistor–transistor logic (TTL) signals to the master node. The master node and trigger boards are individually encased inside 3D-printed covers.

Fig. 3

Schematic diagram of (a) a multi-module MOBI probe. Optional external power source and trigger board not shown. We also include (b) a photo of the master node with its circuit board exposed and (c) a photo of the hardware trigger board.

NPh_11_3_035002_f003.png

Data acquisition of a MOBI probe is controlled through a MATLAB (Mathworks, Natick, Massachusetts, United States) based graphical user interface (GUI) on the computer, which allows users to perform pre-experimental calibration and setup, acquire and visualize measurements, and save data into the community-driven shared near-infrared spectroscopy format (SNIRF).56 Optical signal quality can be manually adjusted from the GUI by setting the source current and detector gain settings (ranging from 1/6 to 4× gain), with real-time signal display providing instantaneous feedback, including visual indicators for channels near saturation or noise floor. In addition, the GUI accepts software triggers transmitted through LabStreamingLayer (LSL),57 a middleware for synchronizing experimental data streams, and provides visual cues of trigger arrival within the data visualization windows. MOBI also supports wireless communication of data and user-defined commands between the master node and control computer through a WiFi co-processor breakout board (ESP32, Adafruit Industries, New York, New York, United States) for the master node and transmission control protocol/Internet protocol (TCP/IP) communication. During WiFi-based data acquisition, a MOBI probe can be entirely powered by a battery without needing to physically connect to a computer.

2.3.

Automatic Features

2.3.1.

Automatic connection topology recognition

Our MOBI system can automatically recognize the connection topology among all modules to determine the number of modules and the orientation of each module in an arbitrary probe. An animation demonstrating this feature is shown in Fig. 4. Modules can be dynamically added to an existing probe by connecting one of its four FPC connections with another module’s FPC connector using an FPC cable. Upon start-up, each module samples all four of its P2P communication channels to determine if one or multiple modules are connected to its FPC connectors. This sampling of inputs and outputs allows the P2P serial network to automatically determine how each module is positioned relative to others. Each module is also programmed with a global module ID. With this information, the computer can determine the spatial orientation of each module in the probe based on the connection topology, and the master node can implement spatial multiplexing patterns generated using the modular optode configuration analyzer (MOCA) toolbox.43

Fig. 4

Demonstration of automatic determination of module-to-module connectivity, shown in the top-left overlay, using two five-module configurations (Video 1, MP4, 5.18 MB [URL: https://doi.org/10.1117/1.NPh.11.3.035002.s1]).

NPh_11_3_035002_f004.png

2.3.2.

Automatic 3D optode position estimation using IMU

In addition, our MOBI system can automatically determine the 3D position of each source and detector in the probe without the use of an external digitizer. A demonstration of this feature is also provided in Fig. 4. The previously determined connection topology, known FPC cable lengths, and orientation measurements acquired from each onboard IMU are used together to generate a 3D piece-wise spherical model to estimate the surface upon which the probe is placed. The detailed methods are described in Sec. 2.4.4. Module geometry information and optode layouts within a module can then be superimposed upon this surface to derive the 3D location of each optode. IMU measurements are also recorded over the duration of an experiment and can be used for robust temporal signal rejection through real-time monitoring of optode movements during use.

2.3.3.

Spatial multiplexing groupings

Our MOBI modules utilize a spatial multiplexing encoding strategy. Source LEDs are grouped into spatial multiplexing group (SMG) based on their spatial distributions during the probe design phase using our MOCA toolbox.43 LEDs within the same SMG can be illuminated simultaneously due to sufficient spatial separation, provided by the user when creating the probe, and negligible cross-talk while LEDs across SMGs are illuminated sequentially. This approach allows us to accelerate data acquisition by a factor referred to as spatial multiplexing ratio, or SMR, defined as the total number of source LEDs in a probe divided by the number of SMGs. For example, an SMR of 2 suggests that the data acquisition is twofold faster than turning on every source sequentially.

2.4.

System Characterization Protocols

2.4.1.

Optical performance characterization

The MOBI system’s signal quality is evaluated using several performance metrics. Detector sensitivity is characterized using a dual-axis filter wheel setup, shining light from a source LED through 36 levels of attenuation ranging from optical density (OD) ratings between 0 and 4.8 and sampling with a MOBI detector. From these measurements, we can calculate noise equivalence power (NEP), defined as the value at which the optical signal and noise are of equal magnitude. This also allows us to estimate the detector’s dynamic range as the ratio of maximal detectable power over the minimal detectable signal (i.e., NEP). In addition, we characterize MOBI’s performance on the probe level by measuring the signal-to-noise ratio (SNR) at various SD separations (SDSs). A linear three-module probe was created, providing SDSs ranging between 24 and 85 mm. Measurements were performed on an optical phantom (μs=4.7  cm1 and μa=0.063  cm1 at 830 nm), acquiring 100 samples for each channel and wavelength to compute the respective SNRs. Both the source LED current and detector gain settings are tuned to maximize the number of usable channels.

2.4.2.

Sampling rate enhancement using spatial multiplexing

The ability to improve probe sampling rates utilizing SMGs is conceptually demonstrated via our MOCA toolbox using three distinct probe configurations as examples. All tested configurations contain identical numbers of modules (N=5) and only differ in their spatial connections. Using MOCA, we calculate the number of valid intra- and inter-module channels that fall within an SDS threshold, set to 35 mm for this analysis, and determine the lowest number of SMGs needed. We compare the channel numbers, SMGs, and probe data frame rate with these five-module probes to exemplify the importance of spatial multiplexing in modular probe design.

2.4.3.

Optode-scalp conformity improvement

We previously fabricated rigid PCB versions of our flexible modules with the same dimensions and diamond shape, providing an opportunity to quantitatively assess the enhancement in optode-scalp conformity within the constraints of our module design. We first 3D-print a hemisphere phantom of radius 100 mm to simulate a simplified head surface. Identical V-shaped five-module probes are then created using either rigid or flexible modules and placed atop the hemisphere and held in place using a mesh hair net. The 3D positions of all 25 optodes are manually digitized for both rigid and flexible probes using a 3D digitizer (Fastrak, Polhemus, Colchester, Vermont, United States) and then used to compute the distance of each optode to the center of the sphere. Four landmarks along the bottom circumference of the hemisphere are also digitized to estimate the position of the sphere center. The distances between the digitized optode positions from either probe are compared with the expected radial distance (105 mm, consisting of the sphere radius and the module thickness) when the modules are perfectly conforming to the surface.

2.4.4.

Automated optode localization algorithm and accuracy assessment

A fast algorithm was developed to recover 3D optode/probe positions using the acquired IMU readings and the known probe connection topology, explained in a graphical illustration shown in Fig. 8(a). The 3D optode positions are estimated in two steps. First, we select any of the modules, preferably the middle module of a series of inter-connected modules, as the reference (R) with a specified origin (pR), and estimate the center position of the module (P) immediately adjacent to the reference module using the below algorithm: we first use the quaternion readings from the IMUs in modules R and P to compute the angle of rotation between them, θ, and their respective normal vectors, (nR) and (nP), which are the quaternion rotation applied to the (0, 0, 1) unit vector. By assuming that the surface between every adjacent module can be approximated by a spherical surface, we can use θ and the known cable length (LP,R) between the two modules (55 mm in this experiment) to determine the radius of the sphere, r=LP,R/θ. The module center of P can be then computed by

Eq. (1)

pP=pRr(nRnP),
where pP and pR are the respective module center positions, nP and nR are the respective normal vectors, and r is the radius of the sphere. Once we determine the center positions of P using R as reference, we continue applying the above algorithm to all subsequently connected modules and propagate the positions to every connected module. To enhance the accuracy of our estimation, we constrain modules that are connected in a straight line topologically to the same plane by averaging the pair-wised planes between each estimated module and the reference module and projecting the quaternion-derived normal vectors onto this plane. Using the updated normal vectors, we recompute the module center positions starting again from the reference module.

In the second step, we estimate the 3D positions of the optodes from the center of each module based on the two-dimensional (2D) module layout and radii used to estimate module positions in the previous step. This is achieved using the module’s local coordinate system, with the origin located at the module center with the short axis serving as the x-axis. The 2D distance between each optode to the module center is known, defining the arc length between the origin and the optode over a local spherical surface with a radius calculated as the average of the pair-wised spherical radii estimated in the previous step. Once the 3D positions are computed in the local coordinate system, they are then translated to the probe’s world coordinate using the module center’s position and normal vector.

The accuracy of IMU-derived 3D optode positions is assessed using the same flexible V-shaped five-module probe and hemisphere phantom used in the previous section. Manual digitization of optode positions is again performed using the Polhemus digitizer and repeated five times to produce the “ground truth” from the averaged positions and the standard error of multiple digitizations from a single operator. The central module lying at the tip of the V-shaped probe serves as the reference module for IMU-derived optode locations and for registering IMU-derived and manually digitized positions.

2.4.5.

Initial in vivo validation

The MOBI system’s ability to measure hemodynamic changes is evaluated in human subjects using two experimental protocols. First, we perform a blood pressure cuff occlusion experiment on a 31-year-old adult male subject while collecting simultaneous measurements using MOBI and a commercial fNIRS system (Brite23, Artinis, Elst, Netherlands) as a reference to assess the trend and magnitude of recovered physiological changes. A single MOBI module and single-channel Artinis probe are placed on the medial forearm. A blood pressure cuff is positioned on the upper arm, which is rested at the same level as the heart. The cuff is first inflated to 100 mmHg and then to 220 mmHg for 75 s each to achieve venous and arterial occlusions, respectively. Measurements from the 30 mm channels of both the MOBI and Artinis systems are acquired simultaneously at 10 Hz in a dark room.

In addition, the evoked response for a finger-tapping experiment is measured in a 31-year-old, right-handed adult male. The task consists of 10 repetitions of 20-s task periods, during which the non-dominant hand sequentially taps the thumb to the index, middle, ring, and little fingers of the same hand, and 25-s rest periods, with 60-s rest periods added before and after the entire task. The subject sits in a comfortable position with their eyes closed during the entire experiment and is verbally instructed to start and stop tapping. During rest, the subject is instructed to place both hands on their laps. A probe consisting of two MOBI modules is created and symmetrically placed over the ipsilateral and contralateral primary motor cortices,58 positioned halfway between C1 and C3 for the left hemisphere and between C2 and C4 for the right hemisphere using the EEG 10-20 system.59 The modules are mounted and secured using a flexible, 3D-printed head cap developed in-house [Fig. 10(a)].60 Our cap contains embedded 10-10 position markers computed from an adult atlas to guide module placement and features a wire-frame mesh to securely and conveniently anchor MOBI modules in any position using 3D-printed locking pins and clips. This combination of wire-frame and locking pins also allows for the rapid re-configuration of MOBI modules between experimental probe designs. A black light-blocking fabric is placed over the head cap to further block ambient light.

Acquired data are converted into the SNIRF format and loaded into Homer361 for processing. A simple processing stream is employed—the raw data are converted to OD, a band-pass filter with pass-band frequencies between 0.01 and 0.4 Hz is applied, and the OD is converted to blood concentration using the modified Beer–Lambert law with the partial path-length factor set to 1.62 The hemodynamic response function (HRF) is then computed with the general linear model (GLM) using the short separation most closely correlated with each measurement channel to regress physiological signals.63

3.

Results

3.1.

System Characterization

Characterizing the optical sensitivity of the MOBI modules, as shown in Fig. 5(a), demonstrates a noise equivalent power of 8.90  pW/Hz at 735 nm and 7.28  pW/Hz at 850 nm. In addition, the detector saturates at 250 and 220 nW for 735 and 850 nm, respectively, providing a dynamic range of around 88 dB for both wavelengths. Figure 5(b) shows the SNR versus SDS over an optical phantom. For this measurement, the source current is attenuated to 50%, and no detector gain is applied to maximize the number of non-saturated channels. Under these conditions, intra-module channels (24.8 and 30 mm) exhibit an SNR of 65 to 72 dB, and SNRs of greater than 40 dB are measured in SDSs of up to 50 mm.

Fig. 5

System characterization showing (a) measured signal versus incident optical power, providing noise equivalent power at both 735 and 850 nm and (b) SNR of versus SDS from a three-module probe, with each point representing a unique channel.

NPh_11_3_035002_f005.png

Although we do not have enough completed modules to physically create and test a full-head MOBI probe, we use our MOCA toolbox and MOBI module specifications to predict some of the key metrics for a full-head probe configuration. Based on the head surface area computed from an adult atlas, we estimate that a total of 20 MOBI modules are needed. Using the MOCA toolbox, this full-head probe provides 372 dual-wavelength channels, including both intra- and inter-module channels. This increases the achievable channel density to 1.72 channels per square centimeter. In addition, with MOBI’s spatial multiplexing encoding strategy and an SDS cutoff of 50 mm, this full-head probe results in 14 spatial multiplexing groups and a full-frame rate of 2.4 Hz. In terms of power, assuming each source is operated at the maximum 100 mA to produce an illumination of 12.3 and 9.6 mW at 735 and 850 nm, respectively, at the output of the coupler,64 the total power draw of the full-head probe is estimated at 2.31 W, resulting in a 2.85-h battery life when a 3.7 V 2000 mAh battery is used.

3.2.

Spatial Multiplexing Groups and Data Acquisition Speed

As described in Sec. 2.4.2, three example five-module probe layouts are tested and shown in Fig. 6. The predicted performance using MOCA for each probe is compared in Table 1. When using a sequential acquisition scheme, the full-frame rate of the probe is inversely proportional to the number of sources in the probe, resulting in 2.2 Hz for all three configurations. Using MOCA, we can divide the sources into SMGs so that LEDs belonging to different SMGs can be turned on simultaneously. In Table 1, we report different probe-level metrics reported by MOCA. It is evident that SMGs can significantly enhance the full-frame rate of probes, while module layouts that are more spatially packed require higher numbers of SMGs.

Fig. 6

Three example probe layouts composed of five identical MOBI modules. Optodes are represented by small red circles (sources) and blue crosses (detectors). Each layout has multiple spatial multiplexing groups determined based on the global proximity of sources to each other. Red dashed circles show which sources are simultaneously on for each layout’s first spatial multiplexing group.

NPh_11_3_035002_f006.png

Table 1

Full-frame rate of a probe depends on the layout of the modules within a probe. The three layouts in Fig. 6 result in the channels, groupings, and full-frame rates below.

Layout 1Layout 2Layout 3
Number of modules555
Frame rate for sequential acquisition (Hz)2.22.22.2
Number of channels385355
Number of spatial multiplexing groups577
Improvement ratio3.02.22.2
Spatial encoding full-frame rate (Hz)6.64.84.8

3.3.

Quantifying Optode Conformity to Head Surface

Figure 8(b) shows a photo of the five-module probe, constructed using either rigid or flexible modules [see Figs. 7(a) and 7(b)], pressed upon a hemispherical phantom with a mesh hair net. If the module is perfectly conforming to the phantom surface, the distance between the top of the module and the center of the sphere is expected to be 105 mm, consisting of the 100 mm dome radius and 5 mm thickness of the module. For the rigid-circuit probe, the average measured distance is 107.0±3.7mm, with nine optodes exhibiting errors of more than 5 mm from the expected distance. The use of flexible probes lowers the average measured distance to 104.8±2.8mm. In addition, measured distances for the flexible modules are more centrally distributed [Fig. 7(f)] when compared with those of the rigid modules [Fig. 7(c)].

Fig. 7

Comparison between (a) a rigid and (d) flexible module implementation in terms of conformity to a hemisphere phantom. Specifically, we report the distances to the sphere center measured for the (b) rigid and (e) flexible five-module probes using a digitizer. The histogram of the resulting distances for a total of 25 optodes is shown for the (c) rigid and (f) flexible probes. With a 5-mm module thickness, a reading of 105 mm is expected for a perfectly conforming probe.

NPh_11_3_035002_f007.png

3.4.

Evaluating the Accuracy of Automatic Optode 3D Positioning

In Fig. 8(c), we plot and compare the 3D optode positions acquired using the Polhemus digitizer with those computed using the MOBI module’s IMU readings and a piece-wise spherical model. The IMU-derived sources (red cross) and detector (blue cross) positions are overlaid on the manually digitized positions (cyan circles). Figure 8(d) plots the error between manually and automatically determined optode positions across all 25 optodes of the probe, indicating an average error of 4.2 mm, and Fig. 8(e) shows the breakdown by module. In comparison, the standard deviation across the five repetitions of digitizer readings is 1.8 mm.

Fig. 8

Schematic and results for automated optode position estimation, showing (a) a diagram of the piece-wise spherical algorithm to sequentially estimate module positions, pi (i=1,2,), from a reference module location, pR, using their quaternion-derived unitary normal vectors, ni, and the computed radius ri between adjacent modules; (b) a photo of a V-shaped MOBI probe secured by a hair net over a hemisphere phantom; and (c) estimated module positions (green) and optode positions (red—source, blue—detector) using IMU data versus manually digitized positions (purple circles). We also report the error of IMU-derived optode positions (d) on the probe level and (e) per module.

NPh_11_3_035002_f008.png

3.5.

Physiological Measurements of Venous and Arterial Occlusion

We plot the estimated changes in hemoglobin concentrations during the dual-pressure blood cuff occlusion experiment in Fig. 9. During venous occlusion, oxygenated (HbO) and deoxygenated (HbR) hemoglobin show increasing trends in the results from both systems [Figs. 9(a) and 9(b)]. The arterial occlusion resulted in a negative correlation between HbO and HbR, as demonstrated by the horizontal total hemoglobin concentration (HbT) line. Upon releasing pressure, there is a brief increase in HbO and HbT accompanied by decreasing HbR.

Fig. 9

Results from a dual-pressure blood occlusion experiment using (a) a single MOBI module and (b) a single Artinis channel placed on the forearm. Venous (100 mmHg) and arterial (220 mmHg) occlusions lasted 75 s each prior to release.

NPh_11_3_035002_f009.png

3.6.

Evoked Response in Motor Execution Task

In Fig. 10(a), we include a photo showing two modules affixed to the MOBI 3D-printed head cap for the finger-tapping experiment, while Fig. 10(b) provides the HRFs corresponding with each channel. A notable hemodynamic response is shown between detector D1 and source S2 in the contralateral (right) hemisphere, with a peak increase in HbO of 40  μMmm 8  s after stimulus onset.

Fig. 10

Validation of MOBI system in a finger-tapping experiment, showing (a) a photo of a two-module probe mounted over the wire-frame head cap and (b) resulting hemodynamic responses for each channel, with red, green, and blue lines, indicating HbO, HbR, and HbT changes, respectively. Source LEDs (red) and detectors (blue) are numbered in both plots.

NPh_11_3_035002_f010.png

4.

Discussion

While modular fNIRS devices convey critical advantages in portability, wearability, and re-configurability, the MOBI system presented here is, to the best of our knowledge, the first 3D-aware and fully fPCB-based system. These features are expected to further facilitate the translation of fNIRS from laboratory and research settings to real-world applications and provide an expanded toolkit that empowers fNIRS studies in naturalistic environments. In combination with our 3D-printed, wire-frame head cap, the flexible MOBI modules freely conform to any head shape, helping provide more robust optode-scalp contact even in environments with increased motion. Sampling from the IMUs embedded in each module also provides real-time, module-specific information on subject and probe movement while enabling the automatic determination of optode 3D positioning. With further improvements, this 3D optode localization function could readily enable tomographic fNIRS data analysis for more accurate localization of cerebrovascular activation while dramatically shortening probe setup time without the need for manual digitization. Additional support for wireless data collection over WiFi and easy integration with data synchronization protocols such as LSL uniquely position MOBI to aid the exploration of the human brain’s day-to-day functions.

Characterization of the MOBI system’s optical sensitivity, as shown in Fig. 5(a), reports an NEP of 7.3 and 8.9  pW/Hz at 735 and 850 nm illumination, respectively, which is capable of measuring SDSs of up to 45 to 50 mm in human head tissue based on previously reported studies.40 A dynamic range of 88  dB, in conjunction with programmable source currents and detector gain settings that expand the effective range of measurable signals, also allows MOBI to simultaneously sample from long, inter-module channels and short separation channels for physiological noise regression65 without saturating the detectors.

The MOBI probe design pays specific attention to enhanced portability, long-term wearability, and accelerated experimental setup. The combination of the 3D printed wire-frame head cap, rapid clip-on module mounting, and automatic optode 3D position sensing makes it possible to potentially shorten the experimental setup time by half or even more. With a measured 9.8   input resistance for each module, MOBI probes can be completely powered by a lightweight battery pack, theoretically capable of accommodating up to 75 modules on a 5 V power supply before dropping below the 3.3 V necessary to drive the microcontrollers. Each module weighs only 18 g, meaning that a 10-module probe, including a head cap, modules, master node, and a battery pack, would total less than 300 g, making MOBI well-suited for long-term wear.

Full-head modular probes require not only lightweight modules and large dynamic ranges but also an encoding strategy to ensure fast data frame rates. The results reported in Fig. 6 and Table 1 show how a spatial multiplexing encoding strategy can improve a full probe’s sampling rate. As long as sources are adequately spatially dispersed, their contributions to a detector’s readings avoid cross-talk. Table 1 also shows that tightly packed modules require more SMGs, thus lowering the frame rate, compared with layouts that are more spatially separate even though the total number of modules remains constant. In this work, the SMG assignment is computed using a 2D probe representation. It is also possible to perform this on a 3D curved head surface, although we anticipate a potential slight increase in required SMG numbers due to shortened distances. Although spatial multiplexing is a software feature and does not require extensive hardware, it does require knowledge of the location of all sources and detectors to identify the SMGs for a particular probe.

In addition, based on the comparisons between manual digitization shown in Fig. 7, it is not surprising that rigid PCB modules yield a higher average optode position than the expected value when compared with its flexible counterpart, as a result of poorer conformity to a curved surface. We want to note that previously reported rigid or flex-rigid modular fNIRS systems usually address such issues using a smaller module size36,66 or spring-loaded couplers48 in conjunction with flexible cables between rigid modules.41 MOBI’s fully flexible modules offer broader intra-module SDS ranges compared with modules with a smaller area, and its optode-to-scalp coupling is expected to be further improved when combined with spring-loaded couplers similar to those used among other systems.

Figure 8 provides initial validation for the accuracy of our internal IMU-based optode positioning system. It demonstrates that the average error for MOBI’s IMU-based optode position estimates is around 4.2 mm, and the calculation requires less than a second to complete. As expected, optode localization exhibits lower errors near the reference position (module 3) and increases when moving away from the reference. While this error is higher than that of repeated manual digitizations, which shows an average standard deviation of 1.8 mm over five repetitions, there are a few points of note. First, inaccuracies in manually digitized optode positions may compound the error of the IMU-derived positions reported here. As shown in Fig. 7(f), digitized positions for a flexible probe exhibit deviations of up to 5 mm from the expected distance of 105 mm between optodes and the center of the hemispherical phantom. The error of IMU-derived positions is computed with respect to each optode’s digitized positions, potentially resulting in an overestimation of said error. The higher localization error for module 5 in Fig. 8(e) may be attributed to sub-optimal coupling that results from the non-uniform pressure applied by the elastic hair net when securing modules onto the large hemispherical dome phantom. The development of our 3D-printed headgear, as shown in Fig. 10(a), is specifically intended to address such issues in human studies. Finally, the piece-wise spherical approximation approach used to compute the position of each module with respect to other modules is a naive approach toward position estimation with respect to an unknown extrinsic reference using intrinsic quaternions. To further improve 3D optode position estimation accuracy, we would consider strategies developed for the study of human kinemetics67 or optimization-based methods.68

We would like to note here as well that this method requires knowledge of the probe’s connection topology, inter-module distances, module geometry, and optode layout. While we use fixed FPC cable lengths to connect adjacent modules in this experiment, it is perceivable that experiments targeting spatially distant ROIs may require varying lengths of inter-module cables, which must be provided to obtain accurate 3D optode positions. In addition, anatomically labeled headgear, such as our 3D-printed cap, serves as the bridge for registering the computed relative optode positions onto a head surface by allowing to affix the modular probe between desired head landmark positions. The use of emerging augmented reality and real-time head landmark tracking workflows could further shorten fNIRS probe setup time and improve probe placement accuracy.69 Although the use of orientation sensors in modular fNIRS systems has been previously reported for motion artifact correction,7072 our reported results are among the first for extending IMU readings for 3D probe shape acquisition.

Preliminary experiments in human subjects demonstrate that the MOBI system is capable of recovering changes in human physiology. Figure 9 shows that simultaneous measurements using MOBI and a commercial fNIRS system during an arm cuff occlusion experiment recover similar hemodynamic changes compared with previous research results,39 such as the continuous increase of HbO and HbT during venous occlusion at 100 mmHg and a flattened HbT during arterial occlusion at 220 mmHg. In addition, both systems capture the transient increase in HbO following the release of arterial occlusion. The results from a single subject during a finger-tapping task also produce a distinctive and expected hemodynamic response73 in the contralateral hemisphere following physiological signal regression using GLM.

Future work on the MOBI platform includes applying MOBI in experiments with more complex tasks, particularly in settings involving greater subject mobility. This would also allow further evaluation of the accuracy of IMU-derived versus manually digitized optode positions and assess MOBI’s ability to capture changes in probe/optode movement over the course of an experiment. In addition, the data acquisition frame rate is currently limited by the relatively low sampling rate and bit-depth of the ADC (12-bit at 32 ksps) included in the microcontroller. This limitation can be addressed by adopting higher-performance sources and detectors. Increasing the source LED power, while still keeping it under safety limits, and adopting a high-speed ADC in the detection hardware would lead to significant improvements in data frame rate. We are also considering using an analog front-end (AFE) to enable signal amplification, improve sampling rates, and perform onboard signal filtering. Any resulting improvements to SNR and sensitivity would also enable measurements across longer SDSs and thereby increase sampling density.

5.

Conclusion

We have designed and validated a lightweight, fiber-less, diamond-shaped modular fNIRS system with features tailored toward its use in natural environments. Its lightweight and flexible circuit-based module design allows the system to conform to the scalp for improved optode-scalp coupling during use. In addition, a P2P communication network allows for the automatic determination of modular connection topology for the implementation of spatial multiplexing groupings to increase the frame rate of any probe configuration relative to sequential encoding methods. This, in conjunction with the 3D orientation sensors embedded in each module, also allows for the determination of optode positions without the need for external hardware. Our MOBI modules were validated against a commercial system in blood pressure cuff occlusion and finger-tapping tests. The MOBI fNIRS system directly addresses ergonomic considerations, enables fast experimental setup, and achieves the robust measurements necessary for advancing fNIRS toward further tackling challenges in naturalistic neuromonitoring tasks.

Disclosures

The 3D optode position estimation using IMU sensor networks was disclosed previously in a US patent (US 11,246,529).

Code and Data Availability

Our open-source MATLAB toolbox for designing and optimizing modular fNIRS probes,43 as used in this work, can be freely accessed at http://github.com/COTILab/MOCA.

Acknowledgments

The authors would like to thank Dr. David Boas for his insightful comments and recommendations regarding improving signal quality. This research is supported by the National Institutes of Health (NIH) (Grant Nos. R01-EB026998 and R01-GM114365).

References

1. 

J. E. Crum, “Future applications of real-world neuroimaging to clinical psychology,” Psychol. Rep., 124 (6), 2403 –2426 https://doi.org/10.1177/0033294120926669 PYRTAZ (2021). Google Scholar

2. 

S. Ladouce et al., “Understanding minds in real-world environments: toward a mobile cognition approach,” Front. Hum. Neurosci., 10 694 https://doi.org/10.3389/fnhum.2016.00694 (2017). Google Scholar

3. 

S. Heinzel et al., “Variability of (functional) hemodynamics as measured with simultaneous fNIRS and fMRI during intertemporal choice,” NeuroImage, 71 125 –134 https://doi.org/10.1016/j.neuroimage.2012.12.074 NEIMEF 1053-8119 (2013). Google Scholar

4. 

J. Gross, “Magnetoencephalography in cognitive neuroscience: a primer,” Neuron, 104 (2), 189 –204 https://doi.org/10.1016/j.neuron.2019.07.001 NERNET 0896-6273 (2019). Google Scholar

5. 

F. Herold et al., “Functional near-infrared spectroscopy in movement science: a systematic review on cortical activity in postural and walking tasks,” Neurophotonics, 4 (4), 041403 https://doi.org/10.1117/1.NPh.4.4.041403 (2017). Google Scholar

6. 

M. A. Yücel et al., “Functional near infrared spectroscopy: enabling routine functional brain imaging,” Curr. Opin. Biomed. Eng., 4 78 –86 https://doi.org/10.1016/j.cobme.2017.09.011 (2017). Google Scholar

7. 

P. Sirpal et al., “fNIRS improves seizure detection in multimodal EEG-fNIRS recordings,” J. Biomed. Opt., 24 (5), 051408 https://doi.org/10.1117/1.JBO.24.5.051408 JBOPFO 1083-3668 (2019). Google Scholar

8. 

D. A. Boas et al., “Twenty years of functional near-infrared spectroscopy: introduction for the special issue,” NeuroImage, 85 1 –5 https://doi.org/10.1016/j.neuroimage.2013.11.033 NEIMEF 1053-8119 (2014). Google Scholar

9. 

M. Ferrari and V. Quaresima, “A brief review on the history of human functional near-infrared spectroscopy (fNIRS) development and fields of application,” NeuroImage, 63 (2), 921 –935 https://doi.org/10.1016/j.neuroimage.2012.03.049 NEIMEF 1053-8119 (2012). Google Scholar

10. 

V. Quaresima and M. Ferrari, “A mini-review on functional near-infrared spectroscopy (fNIRS): where do we stand, and where should we go?,” Photonics, 6 (3), 87 https://doi.org/10.3390/photonics6030087 (2019). Google Scholar

11. 

A. C. Ehlis et al., “Application of functional near-infrared spectroscopy in psychiatry,” NeuroImage, 85 478 –488 https://doi.org/10.1016/j.neuroimage.2013.03.067 NEIMEF 1053-8119 (2014). Google Scholar

12. 

V. Kumar et al., “Functional near infra-red spectroscopy (fNIRS) in schizophrenia: a review,” Asian J. Psychiatry, 27 (2017), 18 –31 https://doi.org/10.1016/j.ajp.2017.02.009 (2017). Google Scholar

13. 

V. Quaresima, S. Bisconti and M. Ferrari, “A brief review on the use of functional near-infrared spectroscopy (fNIRS) for language imaging studies in human newborns and adults,” Brain Lang., 121 (2), 79 –89 https://doi.org/10.1016/j.bandl.2011.03.009 (2012). Google Scholar

14. 

S. Rossi et al., “Shedding light on words and sentences: near-infrared spectroscopy in language research,” Brain Lang., 121 (2), 152 –163 https://doi.org/10.1016/j.bandl.2011.03.008 (2012). Google Scholar

15. 

E. S. Jackson et al., “A fNIRS investigation of speech planning and execution in adults who stutter,” Neuroscience, 406 73 –85 https://doi.org/10.1016/j.neuroscience.2019.02.032 (2019). Google Scholar

16. 

R. N. Aslin, M. Shukla and L. L. Emberson, “Hemodynamic correlates of cognition in human infants,” Annu. Rev. Psychol., 66 (1), 349 –379 https://doi.org/10.1146/annurev-psych-010213-115108 ARPSAC 0066-4308 (2015). Google Scholar

17. 

R. E. Vanderwert and C. A. Nelson, “The use of near-infrared spectroscopy in the study of typical and atypical development,” NeuroImage, 85 264 –271 https://doi.org/10.1016/j.neuroimage.2013.10.009 NEIMEF 1053-8119 (2014). Google Scholar

18. 

T. Wilcox and M. Biondi, “fNIRS in the developmental sciences,” Wiley Interdiscipl. Rev. Cogn. Sci., 6 (3), 263 –283 https://doi.org/10.1002/wcs.1343 (2015). Google Scholar

19. 

M. Soltanlou et al., “Applications of functional near-infrared spectroscopy (fNIRS) in studying cognitive development: the case of mathematics and language,” Front. Psychol., 9 277 https://doi.org/10.3389/fpsyg.2018.00277 1664-1078 (2018). Google Scholar

20. 

M. Yang et al., “A systemic review of functional near-infrared spectroscopy for stroke: current application and future directions,” Front. Neurol., 10 58 https://doi.org/10.3389/fneur.2019.00058 (2019). Google Scholar

21. 

G. Brockington et al., “From the laboratory to the classroom: the potential of functional near-infrared spectroscopy in educational neuroscience,” Front. Psychol., 8 1840 https://doi.org/10.3389/fpsyg.2018.01840 1664-1078 (2018). Google Scholar

22. 

D. Lopez-Martinez et al., “Pain detection with fNIRS-measured brain signals: a personalized machine learning approach using the wavelet transform and Bayesian hierarchical modeling with dirichlet process priors,” in 8th Int. Conf. Affect. Comput. and Intell. Interaction Workshops and Demos (ACIIW), 304 –309 (2019). https://doi.org/10.1109/ACIIW.2019.8925076 Google Scholar

23. 

N. Naseer and K. S. Hong, “fNIRS-based brain-computer interfaces: a review,” Front. Hum. Neurosci., 9 3 https://doi.org/10.3389/fnhum.2015.00003 (2015). Google Scholar

24. 

S. Ahn and S. C. Jun, “Multi-modal integration of EEG-fNIRS for brain-computer interfaces—current limitations and future directions,” Front. Hum. Neurosci., 11 1 –6 https://doi.org/10.3389/fnhum.2017.00503 (2017). Google Scholar

25. 

K. S. Hong, M. J. Khan and M. J. Hong, “Feature extraction and classification methods for hybrid fNIRS-EEG brain-computer interfaces,” Front. Hum. Neurosci., 12 246 https://doi.org/10.3389/fnhum.2018.00246 (2018). Google Scholar

26. 

F. Scholkmann et al., “A review on continuous wave functional near-infrared spectroscopy and imaging instrumentation and methodology,” NeuroImage, 85 6 –27 https://doi.org/10.1016/j.neuroimage.2013.05.004 NEIMEF 1053-8119 (2014). Google Scholar

27. 

A. T. Eggebrecht et al., “Mapping distributed brain function and networks with diffuse optical tomography,” Nat. Photonics, 8 (6), 448 –454 https://doi.org/10.1038/nphoton.2014.107 NPAHBY 1749-4885 (2014). Google Scholar

28. 

M. D. Wheelock, J. P. Culver and A. T. Eggebrecht, “High-density diffuse optical tomography for imaging human brain function,” Rev. Sci. Instrum., 90 (5), 051101 https://doi.org/10.1063/1.5086809 RSINAK 0034-6748 (2019). Google Scholar

29. 

N. L. Everdell et al., “A frequency multiplexed near-infrared topography system for imaging functional activation in the brain,” Rev. Sci. Instrum., 76 (9), 093705 https://doi.org/10.1063/1.2038567 RSINAK 0034-6748 (2005). Google Scholar

30. 

A. von Lühmann et al., “Toward neuroscience of the everyday world (new) using functional near-infrared spectroscopy,” Curr. Opin. Biomed. Eng., 18 100272 https://doi.org/10.1016/j.cobme.2021.100272 (2021). Google Scholar

31. 

A. Curtin and H. Ayaz, “The age of neuroergonomics: towards ubiquitous and continuous measurement of brain function with fNIRS,” Jpn. Psychol. Res., 60 (4), 374 –386 https://doi.org/10.1111/jpr.12227 (2018). Google Scholar

32. 

J. L. Park, P. A. Dudchenko and D. I. Donaldson, “Navigation in real-world environments: new opportunities afforded by advances in mobile brain imaging,” Front. Hum. Neurosci., 12 361 https://doi.org/10.3389/fnhum.2018.00361 (2018). Google Scholar

33. 

A. von Luhmann et al., “M3BA: a mobile, modular, multimodal biosignal acquisition architecture for miniaturized EEG-NIRS-based hybrid BCI and monitoring,” IEEE Trans. Biomed. Eng., 64 (6), 1199 –1210 https://doi.org/10.1109/TBME.2016.2594127 IEBEAX 0018-9294 (2017). Google Scholar

34. 

R. Zimmermann et al., “Silicon photomultipliers for improved detection of low light levels in miniature near-infrared spectroscopy instruments,” Biomed. Opt. Express, 4 (5), 659 –666 https://doi.org/10.1364/BOE.4.000659 BOEICL 2156-7085 (2013). Google Scholar

35. 

T. Funane et al., “Rearrangeable and exchangeable optical module with system-on-chip for wearable functional near-infrared spectroscopy system,” Neurophotonics, 5 (1), 011007 https://doi.org/10.1117/1.NPh.5.1.011007 (2017). Google Scholar

36. 

D. G. Wyser et al., “Wearable and modular functional near-infrared spectroscopy instrument with multidistance measurements at four wavelengths,” Neurophotonics, 4 (4), 041413 https://doi.org/10.1117/1.NPh.4.4.041413 (2017). Google Scholar

37. 

H. Zhao et al., “A wide field-of-view, modular, high-density diffuse optical tomography system for minimally constrained three-dimensional functional neuroimaging,” Biomed. Opt. Express, 11 (8), 4110 https://doi.org/10.1364/BOE.394914 BOEICL 2156-7085 (2020). Google Scholar

38. 

E. E. Vidal-Rosas et al., “Evaluating a new generation of wearable high-density diffuse optical tomography (HD-DOT) technology via retinotopic mapping in the adult brain,” in Opt. InfoBase Conf. Pap., 025002 (2021). https://doi.org/10.1117/1.NPh.8.2.025002 Google Scholar

39. 

G. Liu et al., “Development of a miniaturized and modular probe for fNIRS instrument,” Lasers Med. Sci., 37 (4), 2269 –2277 https://doi.org/10.1007/s10103-021-03493-w (2022). Google Scholar

40. 

D. Chitnis et al., “Functional imaging of the human brain using a modular, fibre-less, high-density diffuse optical tomography system,” Biomed. Opt. Express, 7 (10), 4275 –4288 https://doi.org/10.1364/BOE.7.004275 BOEICL 2156-7085 (2016). Google Scholar

41. 

H. Zhao and R. J. Cooper, “Review of recent progress toward a fiberless, whole-scalp diffuse optical tomography system,” Neurophotonics, 5 (1), 011012 https://doi.org/10.1117/1.NPh.5.1.011012 (2017). Google Scholar

42. 

N. M. Gregg et al., “Brain specificity of diffuse optical imaging: improvements from superficial signal regression and tomography,” Front. Neuroenerg., 2 14 https://doi.org/10.3389/fnene.2010.00014 (2010). Google Scholar

43. 

M. Vanegas, M. Mireles and Q. Fang, “MOCA: a systematic toolbox for designing and assessing modular functional near-infrared brain imaging probes,” Neurophotonics, 9 (1), 017801 https://doi.org/10.1117/1.NPh.9.1.017801 (2022). Google Scholar

44. 

S. Perrey, “Possibilities for examining the neural control of gait in humans with fNIRS,” Front. Physiol., 5 10 –13 https://doi.org/10.3389/fphys.2014.00204 FROPBK 0301-536X (2014). Google Scholar

45. 

R. Cooper et al., “A systematic comparison of motion artifact correction techniques for functional near-infrared spectroscopy,” Front. Neurosci., 6 147 https://doi.org/10.3389/fnins.2012.00147 1662-453X (2012). Google Scholar

46. 

S. Brigadoi et al., “Motion artifacts in functional near-infrared spectroscopy: a comparison of motion correction techniques applied to real cognitive data,” NeuroImage, 85 181 –191 https://doi.org/10.1016/j.neuroimage.2013.04.082 NEIMEF 1053-8119 (2014). Google Scholar

47. 

M. A. Yücel et al., “Reducing motion artifacts for long-term clinical nirs monitoring using collodion-fixed prism-based optical fibers,” NeuroImage, 85 192 –201 https://doi.org/10.1016/j.neuroimage.2013.06.054 NEIMEF 1053-8119 (2014). Google Scholar

48. 

A. von Lühmann et al., “Toward a wireless open source instrument: functional near-infrared spectroscopy in mobile neuroergonomics and BCI applications,” Front. Hum. Neurosci., 9 617 https://doi.org/10.3389/fnhum.2015.00617 (2015). Google Scholar

49. 

D. K. Joseph et al., “Diffuse optical tomography system to image brain activation with improved spatial resolution and validation with functional magnetic resonance imaging,” Appl. Opt., 45 8142 –8151 https://doi.org/10.1364/AO.45.008142 APOPAI 0003-6935 (2006). Google Scholar

50. 

B. R. White and J. P. Culver, “Quantitative evaluation of high-density diffuse optical tomography: in vivo resolution and mapping performance,” J. Biomed. Opt., 15 (2), 026006 https://doi.org/10.1117/1.3368999 JBOPFO 1083-3668 (2010). Google Scholar

51. 

I. Mazzonetto et al., “Smartphone-based photogrammetry provides improved localization and registration of scalp-mounted neuroimaging sensors,” Sci. Rep., 12 (1), 10862 https://doi.org/10.1038/s41598-022-14458-6 SRCEC3 2045-2322 (2022). Google Scholar

52. 

P. Pinti et al., “A review on the use of wearable functional near-infrared spectroscopy in naturalistic environments,” Jpn. Psychol. Res., 60 (4), 347 –373 https://doi.org/10.1111/jpr.12206 (2018). Google Scholar

53. 

M. Vanegas et al., “A modular, fiberless, 3-D aware, flexible-circuit-based wearable fNIRS system,” in Biophotonics Congr.: Biomed. Opt. 2020 (Transl., Microsc., OCT, OTS, BRAIN), BM3C.3 (2020). https://doi.org/10.1364/BRAIN.2020.BM3C.3 Google Scholar

54. 

C. H. Bartkowski et al., “Variable source emitter smart textile headgear design for functional near infrared spectroscopy,” in IEEE Eng. in Med. and Biol., (July 2019). https://doi.org/10.13140/RG.2.2.14834.58562 Google Scholar

55. 

T. Muehlemann, D. Haensse and M. Wolf, “Wireless miniaturized in-vivo near infrared imaging,” Opt. Express, 16 (14), 10323 https://doi.org/10.1364/OE.16.010323 OPEXFF 1094-4087 (2008). Google Scholar

56. 

S. Tucker et al., “Introduction to the shared near infrared spectroscopy format,” Neurophotonics, 10 (1), 013507 https://doi.org/10.1117/1.NPh.10.1.013507 (2022). Google Scholar

57. 

C. Kothe, “Lab streaming layer (LSL),” https://github.com/sccn/labstreaminglayer (). Google Scholar

58. 

L. M. Silva et al., “Localizing the primary motor cortex of the hand by the 10-5 and 10-20 systems for neurostimulation: an MRI study,” Clin. EEG Neurosci., 52 (6), 427 –435 https://doi.org/10.1177/1550059420934590 (2021). Google Scholar

59. 

M. R. Nuwer et al., “IFCN standards for digital recording of clinical EEG,” Electroencephalogr. Clin. Neurophysiol., 106 (3), 259 –261 https://doi.org/10.1016/S0013-4694(97)00106-5 ECNEAZ 0013-4694 (1998). Google Scholar

60. 

A. McCann et al., “Designing anatomically derived, 3-D printable head caps for functional neuroimaging,” in Optica Biophotonics Congr.: Biomed. Opt. Congr., (2024). Google Scholar

61. 

T. J. Huppert et al., “HomER: a review of time-series analysis methods for near-infrared spectroscopy of the brain,” Appl. Opt., 48 (10), D280 –D298 https://doi.org/10.1364/AO.48.00D280 APOPAI 0003-6935 (2009). Google Scholar

62. 

D. A. Boas, A. M. Dale and M. A. Franceschini, “Diffuse optical imaging of brain activation: approaches to optimizing image sensitivity, resolution, and accuracy,” NeuroImage, 23 S275 –S288 https://doi.org/10.1016/j.neuroimage.2004.07.011 NEIMEF 1053-8119 (2004). Google Scholar

63. 

M. A. Yücel et al., “Short separation regression improves statistical significance and better localizes the hemodynamic response obtained by near-infrared spectroscopy for tasks with differing autonomic responses,” Neurophotonics, 2 (3), 035005 https://doi.org/10.1117/1.NPh.2.3.035005 (2015). Google Scholar

64. 

E. Xu et al., “Evaluating the effect of optical couplers on fNIRS light delivery,” in fNIRS 2022, The Society for functional Near-Infrard Spectroscopy, (2022). Google Scholar

65. 

L. Gagnon et al., “Short separation channel location impacts the performance of short channel regression in NIRS,” NeuroImage, 59 (3), 2518 –2528 https://doi.org/10.1016/j.neuroimage.2011.08.095 NEIMEF 1053-8119 (2012). Google Scholar

66. 

H. Zhao et al., “Design and validation of a mechanically flexible and ultra-lightweight high-density diffuse optical tomography system for functional neuroimaging of newborns,” Neurophotonics, 8 (1), 015011 https://doi.org/10.1117/1.NPh.8.1.015011 (2021). Google Scholar

67. 

H. G. Kortier et al., “Assessment of hand kinematics using inertial and magnetic sensors,” J Neuroeng. Rehabil., 11 70 https://doi.org/10.1186/1743-0003-11-70 (2014). Google Scholar

68. 

T. McGrath and L. Stirling, “Body-worn IMU human skeletal pose estimation using a factor graph-based optimization framework,” Sensors, 20 (23), 6887 https://doi.org/10.3390/s20236887 SNSRES 0746-9462 (2020). Google Scholar

69. 

F.-Y. Yen, Y.-A. Lin and Q. Fang, “Real-time guidance for fNIRS headgear placement using augmented reality,” in Optica Biophotonics Congr.: Biomed. Opt. 2024, BW1C.6 (2024). Google Scholar

70. 

A. von Lühmann et al., “Improved physiological noise regression in fNIRS: a multimodal extension of the general linear model using temporally embedded canonical correlation analysis,” NeuroImage, 208 116472 https://doi.org/10.1016/j.neuroimage.2019.116472 NEIMEF 1053-8119 (2020). Google Scholar

71. 

M. R. Siddiquee et al., “Movement artefact removal from NIRS signal using multi-channel IMU data,” BioMed. Eng. OnLine, 17 120 https://doi.org/10.1186/s12938-018-0554-9 (2018). Google Scholar

72. 

J. Virtanen et al., “Accelerometer-based method for correcting signal baseline changes caused by motion artifacts in medical near-infrared spectroscopy,” J. Biomed. Opt., 16 (8), 087005 https://doi.org/10.1117/1.3606576 JBOPFO 1083-3668 (2011). Google Scholar

73. 

M. Plichta et al., “Event-related functional near-infrared spectroscopy (fNIRS) based on craniocerebral correlations: reproducibility of activation?,” Hum. Brain Mapp., 28 (8), 733 –741 https://doi.org/10.1002/hbm.20303 HBRME7 1065-9471 (2007). Google Scholar

Biographies of the authors are not available.

CC BY: © The Authors. Published by SPIE under a Creative Commons Attribution 4.0 International License. Distribution or reproduction of this work in whole or in part requires full attribution of the original publication, including its DOI.
Edward Xu, Morris Vanegas, Miguel Mireles, Artem Dementyev, Ashlyn McCann, Meryem Yücel, Stefan A. Carp, and Qianqian Fang "Flexible circuit-based spatially aware modular optical brain imaging system for high-density measurements in natural settings," Neurophotonics 11(3), 035002 (5 July 2024). https://doi.org/10.1117/1.NPh.11.3.035002
Received: 12 February 2024; Accepted: 7 June 2024; Published: 5 July 2024
Advertisement
Advertisement
KEYWORDS
Sensors

Brain

Flexible circuits

Head

Optical circuits

Multiplexing

3D metrology

RELATED CONTENT

MEG beamforming: magnetic source imaging
Proceedings of SPIE (March 13 2009)

Back to Top