The development of constrained source localisation algorithms for human brain imaging

This work sets out to evaluate the potential benefits and pit-falls in using a priori information to help solve the Magnetoencephalographic (MEG) inverse problem. In chapter one the forward problem in MEG is introduced, together with a scheme that demonstrates how a priori information can be incorporated into the inverse problem. Chapter two contains a literature review of techniques currently used to solve the inverse problem. Emphasis is put on the kind of a priori information that is used by each of these techniques and the ease with which additional constraints can be applied. The formalism of the FOCUSS algorithm is shown to allow for the incorporation of a priori information in an insightful and straightforward manner. In chapter three it is described how anatomical constraints, in the form of a realistically shaped source space, can be extracted from a subject’s Magnetic Resonance Image (MRI). The use of such constraints relies on accurate co-registration of the MEG and MRI co-ordinate systems. Variations of the two main co-registration approaches, based on fiducial markers or on surface matching, are described and the accuracy and robustness of a surface matching algorithm is evaluated. Figures of merit introduced in chapter four are shown to given insight into the limitations of a typical measurement set-up and potential value of a priori information. It is shown in chapter five that constrained dipole fitting and FOCUSS outperform unconstrained dipole fitting when data with low SNR is used. However, the effect of errors in the constraints can reduce this advantage. Finally, it is demonstrated in chapter six that the results of different localisation techniques give corroborative evidence about the location and activation sequence of the human visual cortical areas underlying the first 125ms of the visual magnetic evoked response recorded with a whole head neuromagnetometer.

Topographic mapping and source localization of the pattern onset visual evoked magnetic response

The topography of the visual evoked magnetic response (VEMR) to a pattern onset stimulus was investigated using 4 check sizes and 3 contrast levels. The pattern onset response consists of three early components within the first 200ms, CIm, CIIm and CIIIm. The CIIm is usually of high amplitude and is very consistent in latency within a subject. Half field (HF) stimuli produce their strongest response over the contralateral hemisphere; the RHF stimulus exhibiting a lower positivity (outgoing field) and an upper negativity (ingoing field), rotated towards the midline. LHF stimulation produced the opposite response, a lower negative and an upper positive. Larger check sizes produce a single area of ingoing and outgoing field while smaller checks produce on area of ingoing and outgoing field over each hemisphere. Latency did not appear to vary with change in contrast but amplitudes increased with increasing contrast. A more detailed topographic study incorporating source localisation procedures suggested a source for CIIm - 4cm below the scalp, close to the midline with current flowing towards the lateral surface. Similar depth and position estimates but with opposite polarity were obtained for the pattern shift P100m previously. Hence, the P100m and the CIIm may originate in similar areas of visual cortex but reveal different aspects of visual processing. © 1992 Human Sciences Press, Inc.

Functional localisation of human sensory-motor cortex using magnetoencephalography

The 19 channel Neuromagnetometer system in the Clinical Neurophysiology Unit at Aston University is a multi-channel system, unique in the United Kingdom. A bite bar head localisation and MRI co-registration strategy which enabled accurate and reproducible localisation of MEG data into cortical space was developed. This afforded the opportunity to study magnetic fields of the human cortex generated by stimulation of peripheral nerve, by stimulation of visceral sensory receptors and by those evoked through voluntary finger movement. Initially, a study of sensory-motor evoked data was performed in a healthy control population. The techniques developed were then applied to patients who were to undergo neurosurgical intervention for the treatment of epilepsy and I or space occupying lesions. This enabled both validation of the effective accuracy of source localisation using MEG as well as to determine the clinical value of MEG in presurgical assessment of functional localisation in human cortex. The studies in this thesis have demonstrated that MEG can repeatedly and reliably locate sources contained within a single gyrus and thus potentially differentiate between disparate gyral activation. This ability is critical in the clinical application of any functional imaging technique; which is yet to be fully validated by any other 'non-invasive' functional imaging methodology. The technique was also applied to the study of visceral sensory representation in the cortex which yielded important data about the multiple cortical representation of visceral sensory function.

Can you tell your clunis from your cubitus? A benchmark for functional imaging

Advances in functional brain imaging have allowed the development of new investigative techniques with clinical application—ranging from presurgical mapping of eloquent cortex to identifying cortical regions involved in religious experiences. Similarly a variety of methods are available to referring physicians, ranging from metabolic measures such as functional magnetic resonance imaging and positron emission tomography to measurements based on electrical activity such as electroencephalography and magnetoencephalography. However, there are no universal benchmarks by which to judge between these methods. In this study we attempt to develop a standard for functional localisation, based on the known functional organisation of somatosensory cortex. Studies have shown spatially distinct sites of brain activity in response to stimulation of various body parts. Generally these studies have focused on areas with large cortical representations, such as the index finger and face. We tested the limits of magnetoencephalography source localisation by stimulation of body parts, namely the clunis and the cubitus, that map to proximal and relatively poorly represented regions of somatosensory cortex.

Visual evoked magnetic responses (VEMR) to flash and pattern reversal stimulation

The problems of using a single channel magnetometer (BTi, Model 601) in an unshielded clinical environment to measure visual evoked magnetic responses (VEMR) were studied. VEMR to flash and pattern reversal stimuli were measured in 100 normal subjects. Two components, the P100M to pattern reversal and P2M to flash, were measured successfully in the majority of patients. The mean latencies of these components in different decades of life were more variable than the visual evoked potentials (VEP) that have been recorded to these stimuli. The latency of the P100M appeared to increase significantly after about 55 years of age whereas little change occurred for the flash P2M. The effects of blur, check size, stimulus size and luminance intensity on the latency and amplitude of the VEMR were studied. Blurring a small (32') check significantly increased latency whereas blurring a large (70') check had little effect on latency. Increasing check size significantly reduced latency of the P100M but had little effect on amplitude. Increasing the field size decreases the latency and increases the amplitude of the P100M. Within a normal subject, most of the temporal variability of the P100M appeared to be associated with run to run variation rather than between recording sessions on the same day or between days. Reproducibility of the P100M was improved to a degree by employing a magnetically shielded room. Increasing flash intensity decreases the latency and increases the amplitude of the P2M component. The magnitude of the effects of varying stimulus parameters on the VEMR were frequently greater than is normally seen in the VEP. The topography of the P100M and P2M varied over the scalp in normal subjects. Full field responses to a large check could be explained as approximately the sum of the half field responses and were consistent with the cruciform model of the visual cortex. Preliminary source localisation data suggested a shallower source in the visual cortex for the flash P2M compared with the P100M. The data suggest that suitable protocols could be devised to obtain normative data of sufficient quality to use the VEMR to flash and pattern clinically.

Effective electromagnetic noise cancellation with beamformers and synthetic gradiometry in shielded and partly shielded environments

The major challenge of MEG, the inverse problem, is to estimate the very weak primary neuronal currents from the measurements of extracranial magnetic fields. The non-uniqueness of this inverse solution is compounded by the fact that MEG signals contain large environmental and physiological noise that further complicates the problem. In this paper, we evaluate the effectiveness of magnetic noise cancellation by synthetic gradiometers and the beamformer analysis method of synthetic aperture magnetometry (SAM) for source localisation in the presence of large stimulus-generated noise. We demonstrate that activation of primary somatosensory cortex can be accurately identified using SAM despite the presence of significant stimulus-related magnetic interference. This interference was generated by a contact heat evoked potential stimulator (CHEPS), recently developed for thermal pain research, but which to date has not been used in a MEG environment. We also show that in a reduced shielding environment the use of higher order synthetic gradiometry is sufficient to obtain signal-to-noise ratios (SNRs) that allow for accurate localisation of cortical sensory function.

Localising the auditory N1m with event-related beamformers:localisation accuracy following bilateral and unilateral stimulation

The auditory evoked N1m-P2m response complex presents a challenging case for MEG source-modelling, because symmetrical, phase-locked activity occurs in the hemispheres both contralateral and ipsilateral to stimulation. Beamformer methods, in particular, can be susceptible to localisation bias and spurious sources under these conditions. This study explored the accuracy and efficiency of event-related beamformer source models for auditory MEG data under typical experimental conditions: monaural and diotic stimulation; and whole-head beamformer analysis compared to a half-head analysis using only sensors from the hemisphere contralateral to stimulation. Event-related beamformer localisations were also compared with more traditional single-dipole models. At the group level, the event-related beamformer performed equally well as the single-dipole models in terms of accuracy for both the N1m and the P2m, and in terms of efficiency (number of successful source models) for the N1m. The results yielded by the half-head analysis did not differ significantly from those produced by the traditional whole-head analysis. Any localisation bias caused by the presence of correlated sources is minimal in the context of the inter-individual variability in source localisations. In conclusion, event-related beamformers provide a useful alternative to equivalent-current dipole models in localisation of auditory evoked responses.