Requirements for Coregistration Accuracy in On-Scalp MEG

Non-rigid-registration-based positioning and labelling of triaxial OPMs on a flexible cap for wearable magnetoencephalography

BACKGROUND

Recent advances in highly sensitive miniaturized optically pumped magnetometers (OPMs) have enabled the development of wearable magnetoencephalography (MEG) offering great flexibility in experimental setting. The OPM array for wearable MEG is typically attached to a flexible cap and exhibits a variable spatial layout across different subjects, which imposes challenges concerning the efficient positioning and labelling of OPMs.

NEW METHOD

A pair of reflective markers are affixed to each triaxial OPM sensor above its cable to determine its location and sensitive axes. A non-rigid registration of optically digitized marker locations with a pre-labelled template of marker locations is performed to map newly digitized markers to OPMs.

RESULTS

The positioning and labelling of 66 OPM sensors could be completed within 35 s. Across ten experiments, all OPMs were accurately labelled, and the mean test-retest errors were 0.48 mm for sensor locations and 0.20 degree for sensitive axes. By combining six OPMs' positions with their respective recordings, magnetic dipoles inside a phantom were located with a mean error of 5.5 mm, and the best fitted dipole for MEG with auditory stimuli presented was located on a subject's primary auditory cortex.

COMPARISON WITH EXISTING METHODS

The proposed method reduces the reliance on error-prone and laborious manual operations inherent in existing methods, therefore significantly improving the efficiency of OPM positioning and labelling on a flexible cap.

CONCLUSION

We developed a method for the precise and rapid positioning and labelling triaxial OPMs on a flexible cap, thereby facilitating the practical implementation of wearable OPM-MEG.

Co-registration of OPM-MCG signals with CT using optical scanning

Magnetocardiography (MCG) can be used to noninvasively measure the electrophysiological activity of myocardial cells. The high spatial resolution of magnetic source localization can precisely determine the location of cardiomyopathy, which is of great significance for the diagnosis and treatment of cardiovascular disease. To perform magnetic source localization, MCG data must be co-registered with anatomical images. We propose a co-registration method that can be applied to OPM-MCG systems. In this method, the sensor array and the trunk of the subject are scanned using structured light-scanning technology, and the scan results are registered with the reconstructed structure using computed tomography (CT). This can increase the number of effective cloud points acquired and reduce the interference from respiratory motion. The scanning bed of the OPM-MCG system was modified to be consistent with the CT device, ensuring that the state of the body remains consistent between the cardiac magnetometry measurements and CT scans.

Optical Co-Registration Method of Triaxial OPM-MEG and MRI

The advent of optically pumped magnetometers (OPMs) facilitates the development of on-scalp magnetoencephalography (MEG). In particular, the triaxial OPM emerged recently, making simultaneous measurements of all three orthogonal components of vector fields possible. The detection of triaxial magnetic fields improves the interference suppression capability and achieves higher source localization accuracy using fewer sensors. The source localization accuracy of MEG is based on the accurate co-registration of MEG and MRI. In this study, we proposed a triaxial co-registration method according to combined principal component analysis and iterative closest point algorithms for use of a flexible cap. A reference phantom with known sensor positions and orientations was designed and constructed to evaluate the accuracy of the proposed method. Experiments showed that the average co-registered position errors of all sensors were approximately 1 mm and average orientation errors were less than 2.5° in the X -and Y orientations and less than 1.6° in the Z orientation. Furthermore, we assessed the influence of co-registration errors on the source localization using simulations. The average source localization error of approximately 1 mm reflects the effectiveness of the co-registration method. The proposed co-registration method facilitates future applications of triaxial sensors on flexible caps.

Optimising the sensing volume of OPM sensors for MEG source reconstruction

Magnetoencephalography (MEG) based on optically pumped magnetometers (OPMs) has been hailed as the future of electrophysiological recordings from the human brain. In this work, we investigate how the dimensions of the sensing volume (the vapour cell) affect the performance of both a single OPM-MEG sensor and a multi-sensor OPM-MEG system. We consider a realistic noise model that accounts for background brain activity and residual noise. By using source reconstruction metrics such as localization accuracy and time-course reconstruction accuracy, we demonstrate that the best overall sensitivity and reconstruction accuracy are achieved with cells that are significantly longer and wider that those of the majority of current commercial OPM sensors. Our work provides useful tools to optimise the cell dimensions of OPM sensors in a wide range of environments.

OMMR: Co-registration toolbox of OPM-MEG and MRI

Magnetoencephalography (MEG) based on optically pumped magnetometers (OPM-MEG) has shown better flexibility in sensor configuration compared with the conventional superconducting quantum interference devices-based MEG system while being better suited for all-age groups. However, this flexibility presents challenges for the co-registration of MEG and magnetic resonance imaging (MRI), hindering adoption. This study presents a toolbox called OMMR, developed in Matlab, that facilitates the co-registration step for researchers and clinicians. OMMR integrates the co-registration methods of using the electromagnetic digitization system and two types of optical scanners (the structural-light and laser scanner). As the first open-source co-registration toolbox specifically for OPM-MEG, the toolbox aims to standardize the co-registration process and set the ground for future applications of OPM-MEG.

On-scalp magnetocorticography with optically pumped magnetometers: Simulated performance in resolving simultaneous sources

Currently, the gold standard for high-resolution mapping of cortical electrophysiological activity is invasive electrocorticography (ECoG), a procedure that carries with it the risk of serious morbidity and mortality. Due to these risks, the use of ECoG is largely limited to pre-surgical mapping in intractable epilepsy. Nevertheless, many seminal studies in neuroscience have utilized ECoG to explore domains such as visual perception, attention, auditory processing, and sensorimotor behavior. Studies such as these, occurring in patients with epilepsy rather than healthy controls, may lack generalizability, and are limited by the placement of the electrode arrays over the presumed seizure focus. This manuscript explores the use of optically pumped magnetometers (OPMs) to create a non-invasive alternative to ECoG, which we refer to as magnetocorticography. Because prior ECoG studies reveal that most cognitive processes are driven by multiple, simultaneous independent neuronal assemblies, we characterize the ability of a theoretical 56-channel dense OPM array to resolve simultaneous independent sources, and compare it to currently available SQUID devices, as well as OPM arrays with inter-sensor spacings more typical of other systems in development. Our evaluation of this theoretical system assesses many potential sources of error, including errors of sensor calibration and position. In addition, we investigate the influence of geometrical and anatomical factors on array performance. Our simulations reveal the potential of high-density, on-scalp OPM MEG devices to localize electrophysiological brain responses at unprecedented resolution for a non-invasive device.

Calibration and Localization of Optically Pumped Magnetometers Using Electromagnetic Coils

In this paper, we propose a method to estimate the position, orientation, and gain of a magnetic field sensor using a set of (large) electromagnetic coils. We apply the method for calibrating an array of optically pumped magnetometers (OPMs) for magnetoencephalography (MEG). We first measure the magnetic fields of the coils at multiple known positions using a well-calibrated triaxial magnetometer, and model these discreetly sampled fields using vector spherical harmonics (VSH) functions. We then localize and calibrate an OPM by minimizing the sum of squared errors between the model signals and the OPM responses to the coil fields. We show that by using homogeneous and first-order gradient fields, the OPM sensor parameters (gain, position, and orientation) can be obtained from a set of linear equations with pseudo-inverses of two matrices. The currents that should be applied to the coils for approximating these low-order field components can be determined based on the VSH models. Computationally simple initial estimates of the OPM sensor parameters follow. As a first test of the method, we placed a fluxgate magnetometer at multiple positions and estimated the RMS position, orientation, and gain errors of the method to be 1.0 mm, 0.2°, and 0.8%, respectively. Lastly, we calibrated a 48-channel OPM array. The accuracy of the OPM calibration was tested by using the OPM array to localize magnetic dipoles in a phantom, which resulted in an average dipole position error of 3.3 mm. The results demonstrate the feasibility of using electromagnetic coils to calibrate and localize OPMs for MEG.

Imaging somatosensory cortex responses measured by OPM-MEG: Variational free energy-based spatial smoothing estimation approach

In recent years, optically pumped magnetometer (OPM)-based magnetoencephalography (MEG) has shown potential for analyzing brain activity. It has a flexible sensor configuration and comparable sensitivity to conventional SQUID-MEG. We constructed a 32-channel OPM-MEG system and used it to measure cortical responses to median and ulnar nerve stimulations. Traditional magnetic source imaging methods tend to blur the spatial extent of sources. Accurate estimation of the spatial size of the source is important for studying the organization of brain somatotopy and for pre-surgical functional mapping. We proposed a new method called variational free energy-based spatial smoothing estimation (FESSE) to enhance the accuracy of mapping somatosensory cortex responses. A series of computer simulations based on the OPM-MEG showed better performance than the three types of competing methods under different levels of signal-to-noise ratios, source patch sizes, and co-registration errors. FESSE was then applied to the source imaging of the OPM-MEG experimental data.

Co-registration Comparison of On-Scalp Magnetoencephalography and Magnetic Resonance Imaging

Magnetoencephalography (MEG) can non-invasively measure the electromagnetic activity of the brain. A new type of MEG, on-scalp MEG, has attracted the attention of researchers recently. Compared to the conventional SQUID-MEG, on-scalp MEG constructed with optically pumped magnetometers is wearable and has a high signal-to-noise ratio. While the co-registration between MEG and magnetic resonance imaging (MRI) significantly influences the source localization accuracy, co-registration error requires assessment, and quantification. Recent studies have evaluated the co-registration error of on-scalp MEG mainly based on the surface fit error or the repeatability error of different measurements, which do not reflect the true co-registration error. In this study, a three-dimensional-printed reference phantom was constructed to provide the ground truth of MEG sensor locations and orientations relative to MRI. The co-registration performances of commonly used three devices—electromagnetic digitization system, structured-light scanner, and laser scanner—were compared and quantified by the indices of final co-registration errors in the reference phantom and human experiments. Furthermore, the influence of the co-registration error on the performance of source localization was analyzed via simulations. The laser scanner had the best co-registration accuracy (rotation error of 0.23° and translation error of 0.76 mm based on the phantom experiment), whereas the structured-light scanner had the best cost performance. The results of this study provide recommendations and precautions for researchers regarding selecting and using an appropriate device for the co-registration of on-scalp MEG and MRI.

Contactless measurements of retinal activity using optically pumped magnetometers

Optically pumped magnetometers (OPMs) have been adopted for the measurement of brain activity. Without the need to be cooled to cryogenic temperatures, an array of these sensors can be placed more flexibly, which allows for the recording of neuronal structures other than neocortex. Here we use eight OPM sensors to record human retinal activity following flash stimulation. We compare this magnetoretinographic (MRG) activity to the simultaneously recorded electroretinogram of the eight participants. The MRG shows the familiar flash-evoked potentials (a-wave and b-wave) and shares a highly significant amount of information with the electroretinogram (both in a simultaneous and separate measurement). We conclude that OPM sensors have the potential to become a contactless alternative to fiber electrodes for the measurement of retinal activity. Such a contactless solution can benefit both clinical and neuroscientific settings.

Automatic coregistration of MRI and on-scalp MEG

BACKGROUND

Recent progress in optically pumped magnetometers (OPMs) and high-temperature superconducting quantum interference devices (SQUIDs) has facilitated the development of an on-scalp magnetoencephalography (MEG) system that offers high signal intensity and flexibility at a lower cost. While the on-scalp sensor array has high flexibility, it brings new challenges to accurate sensor-to-brain coregistration, which is essential for MEG source localization.

NEW METHOD

A novel automatic filtering algorithm based on plane segmentation was proposed to locate on-scalp MEG sensors in 3D images reconstructed from optical scanning. Global image registration was employed for the automatic alignment of anatomical images and sensor positions.

RESULTS

Seventy-one sensor dummies on the scalp were located and registered to brain anatomical images. The deviations of the sensor location and orientation from the averaged result of 10 measurements were less than 1 mm and 0.6°, respectively. The entire process could be completed in less than 4 min.

COMPARISON WITH EXISTING METHODS

Compared with existing methods that involve various manual procedures, such as moving digitizers to fiducials and repeatedly pulling out sensors, our proposed coregistration method is more efficient and accurate.

CONCLUSION

An automatic method for the coregistration of anatomical structure and on-scalp sensors that will have a large impact on the practical use of on-scalp MEG is developed.

Pragmatic spatial sampling for wearable MEG arrays

Several new technologies have emerged promising new Magnetoencephalography (MEG) systems in which the sensors can be placed close to the scalp. One such technology, Optically Pumped MEG (OP-MEG) allows for a scalp mounted system that provides measurements within millimetres of the scalp surface. A question that arises in developing on-scalp systems is: how many sensors are necessary to achieve adequate performance/spatial discrimination? There are many factors to consider in answering this question such as the signal to noise ratio (SNR), the locations and depths of the sources, density of spatial sampling, sensor gain errors (due to interference, subject movement, cross-talk, etc.) and, of course, the desired spatial discrimination. In this paper, we provide simulations which show the impact these factors have on designing sensor arrays for wearable MEG. While OP-MEG has the potential to provide high information content at dense spatial samplings, we find that adequate spatial discrimination of sources (< 1 cm) can be achieved with relatively few sensors (< 100) at coarse spatial samplings (~ 30 mm) at high SNR. After this point approximately 50 more sensors are required for every 1 mm improvement in spatial discrimination. Comparable discrimination for traditional cryogenic systems require more channels by these same metrics. We also show that sensor gain errors have the greatest impact on discrimination between deep sources at high SNR. Finally, we also examine the limitation that aliasing due to undersampling has on the effective SNR of on-scalp sensors.

Detection of human auditory evoked brain signals with a resilient nonlinear optically pumped magnetometer

•

We have built and tested an Optically Pumped Magnetometer based on non-linear magnetooptical rotation.

•

We present a new approach to OPM- based MEG using a modular NOPM sensor.

•

Our sensor is resilient to non-zero and non-uniform magnetic field environments and crosstalk free.

•

We demonstrate the operation of the NOPM sensor by measuring auditory response and calculating time-frequency representation of power.

Optically Pumped Magnetometers (OPMs) have been hailed as the future of human magnetoencephalography, as they enable a level of flexibility and adaptability that cannot be obtained with systems based on superconductors. While OPM sensors are already commercially available, there is plenty of room for further improvements and customization. In this work, we detected auditory evoked brain fields using an OPM based on the nonlinear magneto-optical rotation (NMOR) technique. Our sensor head, containing only optical and non-magnetizable elements, is connected to an external module including all the electronic components, placed outside the magnetically shielded room. The use of the NMOR allowed us to detect the brain signals in non-zero magnetic field environments. In particular, we were able to detect auditory evoked fields in a background field of 70 nT. We benchmarked our sensor with conventional SQUID sensors, showing comparable performance. We further demonstrated that our sensor can be employed to detect modulations of brain oscillations in the alpha band. Our results are a promising stepping-stone towards the realization of resilient OPM-based magnetoencephalography systems that do not require active compensation.

Multi-channel whole-head OPM-MEG: Helmet design and a comparison with a conventional system

Magnetoencephalography (MEG) is a powerful technique for functional neuroimaging, offering a non-invasive window on brain electrophysiology. MEG systems have traditionally been based on cryogenic sensors which detect the small extracranial magnetic fields generated by synchronised current in neuronal assemblies, however, such systems have fundamental limitations. In recent years, non-cryogenic quantum-enabled sensors, called optically-pumped magnetometers (OPMs), in combination with novel techniques for accurate background magnetic field control, have promised to lift those restrictions offering an adaptable, motion-robust MEG system, with improved data quality, at reduced cost. However, OPM-MEG remains a nascent technology, and whilst viable systems exist, most employ small numbers of sensors sited above targeted brain regions. Here, building on previous work, we construct a wearable OPM-MEG system with ‘whole-head’ coverage based upon commercially available OPMs, and test its capabilities to measure alpha, beta and gamma oscillations. We design two methods for OPM mounting; a flexible (EEG-like) cap and rigid (additively-manufactured) helmet. Whilst both designs allow for high quality data to be collected, we argue that the rigid helmet offers a more robust option with significant advantages for reconstruction of field data into 3D images of changes in neuronal current. Using repeat measurements in two participants, we show signal detection for our device to be highly robust. Moreover, via application of source-space modelling, we show that, despite having 5 times fewer sensors, our system exhibits comparable performance to an established cryogenic MEG device. While significant challenges still remain, these developments provide further evidence that OPM-MEG is likely to facilitate a step change for functional neuroimaging.

Ultrasensitive Magnetic Field Sensors for Biomedical Applications

The development of magnetic field sensors for biomedical applications primarily focuses on equivalent magnetic noise reduction or overall design improvement in order to make them smaller and cheaper while keeping the required values of a limit of detection. One of the cutting-edge topics today is the use of magnetic field sensors for applications such as magnetocardiography, magnetotomography, magnetomyography, magnetoneurography, or their application in point-of-care devices. This introductory review focuses on modern magnetic field sensors suitable for biomedicine applications from a physical point of view and provides an overview of recent studies in this field. Types of magnetic field sensors include direct current superconducting quantum interference devices, search coil, fluxgate, magnetoelectric, giant magneto-impedance, anisotropic/giant/tunneling magnetoresistance, optically pumped, cavity optomechanical, Hall effect, magnetoelastic, spin wave interferometry, and those based on the behavior of nitrogen-vacancy centers in the atomic lattice of diamond.

On-scalp MEG sensor localization using magnetic dipole-like coils: A method for highly accurate co-registration

Source modelling in magnetoencephalography (MEG) requires precise co-registration of the sensor array and the anatomical structure of the measured individual's head. In conventional MEG, the positions and orientations of the sensors relative to each other are fixed and known beforehand, requiring only localization of the head relative to the sensor array. Since the sensors in on-scalp MEG are positioned on the scalp, locations of the individual sensors depend on the subject's head shape and size. The positions and orientations of on-scalp sensors must therefore be measured at every recording. This can be achieved by inverting conventional head localization, localizing the sensors relative to the head - rather than the other way around. In this study we present a practical method for localizing sensors using magnetic dipole-like coils attached to the subject's head. We implement and evaluate the method in a set of on-scalp MEG recordings using a 7-channel on-scalp MEG system based on high critical temperature superconducting quantum interference devices (high-T_c SQUIDs). The method allows individually localizing the sensor positions, orientations, and responsivities with high accuracy using only a short averaging time (≤ 2 ​mm, < 3° and < 3%, respectively, with 1-s averaging), enabling continuous sensor localization. Calibrating and jointly localizing the sensor array can further improve the accuracy of position and orientation (< 1 ​mm and < 1°, respectively, with 1-s coil recordings). We demonstrate source localization of on-scalp recorded somatosensory evoked activity based on co-registration with our method. Equivalent current dipole fits of the evoked responses corresponded well (within 4.2 ​mm) with those based on a commercial, whole-head MEG system.

Non-Invasive Functional-Brain-Imaging with an OPM-based Magnetoencephalography System

A non-invasive functional-brain-imaging system based on optically-pumped-magnetometers (OPM) is presented. The OPM-based magnetoencephalography (MEG) system features 20 OPM channels conforming to the subject's scalp. We have conducted two MEG experiments on three subjects: assessment of somatosensory evoked magnetic field (SEF) and auditory evoked magnetic field (AEF) using our OPM-based MEG system and a commercial MEG system based on superconducting quantum interference devices (SQUIDs). We cross validated the robustness of our system by calculating the distance between the location of the equivalent current dipole (ECD) yielded by our OPM-based MEG system and the ECD location calculated by the commercial SQUID-based MEG system. We achieved sub-centimeter accuracy for both SEF and AEF responses in all three subjects. Due to the proximity (12 mm) of the OPM channels to the scalp, it is anticipated that future OPM-based MEG systems will offer enhanced spatial resolution as they will capture finer spatial features compared to traditional MEG systems employing SQUIDs.

Optically pumped magnetometers: From quantum origins to multi-channel magnetoencephalography

Optically Pumped Magnetometers (OPMs) have emerged as a viable and wearable alternative to cryogenic, superconducting MEG systems. This new generation of sensors has the advantage of not requiring cryogenic cooling and as a result can be flexibly placed on any part of the body. The purpose of this review is to provide a neuroscience audience with the theoretical background needed to understand the physical basis for the signal observed by OPMs. Those already familiar with the physics of MRI and NMR should note that OPMs share much of the same theory as the operation of OPMs rely on magnetic resonance. This review establishes the physical basis for the signal equation for OPMs. We re-derive the equations defining the bounds on OPM performance and highlight the important trade-offs between quantities such as bandwidth, sensor size and sensitivity. These equations lead to a direct upper bound on the gain change due to cross-talk for a multi-channel OPM system.

Potential of on-scalp MEG: Robust detection of human visual gamma-band responses

Electrophysiological signals recorded intracranially show rich frequency content spanning from near‐DC to hundreds of hertz. Noninvasive electromagnetic signals measured with electroencephalography (EEG) or magnetoencephalography (MEG) typically contain less signal power in high frequencies than invasive recordings. Particularly, noninvasive detection of gamma‐band activity (>30 Hz) is challenging since coherently active source areas are small at such frequencies and the available imaging methods have limited spatial resolution. Compared to EEG and conventional SQUID‐based MEG, on‐scalp MEG should provide substantially improved spatial resolution, making it an attractive method for detecting gamma‐band activity. Using an on‐scalp array comprised of eight optically pumped magnetometers (OPMs) and a conventional whole‐head SQUID array, we measured responses to a dynamic visual stimulus known to elicit strong gamma‐band responses. OPMs had substantially higher signal power than SQUIDs, and had a slightly larger relative gamma‐power increase over the baseline. With only eight OPMs, we could obtain gamma‐activity source estimates comparable to those of SQUIDs at the group level. Our results show the feasibility of OPMs to measure gamma‐band activity. To further facilitate the noninvasive detection of gamma‐band activity, the on‐scalp OPM arrays should be optimized with respect to sensor noise, the number of sensors and intersensor spacing.

Optical Co-registration of MRI and On-scalp MEG

To estimate the neural generators of magnetoencephalographic (MEG) signals, MEG data have to be co-registered with an anatomical image, typically an MR image. Optically-pumped magnetometers (OPMs) enable the construction of on-scalp MEG systems providing higher sensitivity and spatial resolution than conventional SQUID-based MEG systems. We present a co-registration method that can be applied to on-scalp MEG systems, regardless of the number of sensors. We apply a structured-light scanner to create a surface mesh of the subject’s head and the sensor array, which we fit to the MR image. We quantified the reproducibility of the mesh and localised current dipoles with a phantom. Additionally, we measured somatosensory evoked fields (SEFs) to median nerve stimulation and compared the dipole positions between on-scalp and SQUID-based systems. The scanner reproduced the head surface with <1 mm error. Phantom dipoles were localised with 2.1 mm mean error. SEF dipoles corresponding to the P35m response for OPMs were well localised to the somatosensory cortex, while SQUID dipoles for two subjects were erroneously localised to the motor cortex. The developed co-registration method is inexpensive, fast and can easily be applied to on-scalp MEG. It is more convenient than traditional co-registration methods while also being more accurate.

On-scalp MEG system utilizing an actively shielded array of optically-pumped magnetometers

The spatial resolution of magnetoencephalography (MEG) can be increased from that of conventional SQUID-based systems by employing on-scalp sensor arrays of e.g. optically-pumped magnetometers (OPMs). However, OPMs reach sufficient sensitivity for neuromagnetic measurements only when operated in a very low absolute magnetic field of few nanoteslas or less, usually not reached in a typical magnetically shielded room constructed for SQUID-based MEG. Moreover, field drifts affect the calibration of OPMs. Static and dynamic suppression of interfering fields is thus necessary for good-quality neuromagnetic measurements with OPMs. Here, we describe an on-scalp MEG system that utilizes OPMs and external compensation coils that provide static and dynamic shielding against ambient fields.

In a conventional two-layer magnetically shielded room, our coil system reduced the maximum remanent DC-field component within an 8-channel OPM array from 70 to less than 1 nT, enabling the sensors to operate in the sensitive spin exchange relaxation-free regime. When compensating field drifts below 4 Hz, a low-frequency shielding factor of 22 dB was achieved, which reduced the peak-to-peak drift from 1.3 to 0.4 nT and thereby the standard deviation of the sensor calibration from 1.7% to 0.5%. Without band-limiting the field that was compensated, a low-frequency shielding factor of 43 dB was achieved.

We validated the system by measuring brain responses to electric stimulation of the median nerve. With dynamic shielding and digital interference suppression methods, single-trial somatosensory evoked responses could be detected. Our results advance the deployment of OPM-based on-scalp MEG in lighter magnetic shields.