Evaluation of cortical current density imaging methods using intracranial electrocorticograms and functional MRI

OBJECTIVE EEG source imaging provides important information regarding the underlying neural activity from noninvasive electrophysiological measurements. The aim of the present study was to evaluate source reconstruction techniques by means of the intracranial electrocorticograms (ECoGs) and functional MRI. METHODS Five source imaging algorithms, including the minimum norm least square (MNLS), LORETA with L(p)-norm (p equal to 1, 1.5 and 2), sLORETA, the minimum L(p)-norm (p equal to 1 and 1.5; when p=2, the MNLS method is mathematically equivalent to the minimum L(p)-norm) and L(1)-norm (the linear programming) methods, were evaluated in a group of 10 human subjects, in a paradigm with somatosensory stimulation. Cortical current density (CCD) distributions were estimated from the scalp somatosensory evoked potentials (SEPs), at approximately 30 ms following electrical stimulation of median nerve at the wrist. Realistic geometry boundary element head models were constructed from the MRIs of each subject and used in the CCD analysis. Functional MRI results obtained from a motor task and sensory stimulation in all subjects were used to identify the central sulcus, motor and sensory areas. In three patients undergoing neurosurgical evaluation, ECoGs were recorded in response to the somatosensory stimulation, and were used to help determine the central sulcus and the sensory cortex. RESULTS The CCD distributions estimated by the L(p)-norm and LORETA-L(p) methods were smoother when the p values were high. The LORETA based on the L(1)-norm performed better than the LORETA-L(2) method for imaging well localized sources such as the P30 component of the SEP. The mean and standard deviation of the distance between the location of maximum CCD value and the central sulcus, estimated by the minimum L(p)-norm (with p equal to 1), L(1)-norm (the Linear programming) and LORETA-L(p) (with p equal to 1) methods, were 4, 7, 7 mm and 3, 4, 2 mm, respectively (after converting into Talairach coordinates). The mean and standard deviation of the aforementioned distance, estimated by the MNLS, LORETA with L(p)-norm (p equal to 1.5 and 2.0), sLORETA and the minimum L(p)-norm (p equal to 1.5) methods, were over 11 mm and 6 mm, respectively. CONCLUSIONS The present experimental study suggests that L(1)-norm-based algorithms provide better performance than L(2) and L(1.5)-norm-based algorithms, in the context of CCD imaging of well localized sources induced by somatosensory electrical stimulation of median nerve at the wrist.