Web Portal: A Database of Realistic Simulated and Empirical MEG Data for Testing Algorithms

 

VI
V2/V3
I. LOG
IPS
S.LOG
DLPFC
AC
RHC
Set 1.A
80 ms
90 ms
 
100 ms
     
30 nAm
30 nAm
 
30 nAm
    
1.B
80 ms
90 ms
 
100 ms
     
50 nAm
50 nAm
 
50 nAm
    
Set 2
 
90 ms
90 ms
100 ms
      
15 nAm
30 nAm
30 nAm
    
Set 3.A
80 ms
90 ms
 
100 ms
     
30 nAm
30 nAm
 
30 nAm
    
3.B
80 ms
90 ms
 
100 ms
     
50 nAm
50 nAm
 
50 nAm
    
Set 4
 
90 ms
90 ms
100 ms
100 ms
300 ms*
400 ms*
   
15 nAm
30 nAm
20 nAm
30 nAm
20 nAm
30 nAm
 
Set 5
 
90 ms
90 ms
100 ms
100 ms
300 ms*
400 ms*
80 ms
  
15 nAm
30 nAm
20 nAm
30 nAm
20 nAm
30 nAm
51 nAm
*DLPFC and AC were treated as ramping activity peaking later in time
In the simulated example shown in Fig. 1, a Freesurfer-segmented gray matter/white matter boundary for the simulations was imported into MRIVIEW (Fig. 1a), although the segmentation may also be accomplished within MRIVIEW. The simulated activation timecourses (signal) are shown Fig. 1b. In each case, 100 single trials of real spontaneous background activity were averaged together as the noise trial for each of the 5 participants and for each of the MEG systems (Fig. 1c). Then the signal was embedded within the averaged noise file (Fig. 1d). For all simulated datasets on the web portal, a spherical head model was used for the simulations and modeled data; however, a boundary element model (BEM) is also available in MRIVIEW.
A272390_1_En_14_Fig1_HTML.gif
Fig. 1
A Freesurfer-segmented gray matter/white matter boundary for the simulations (shown in red) was imported into MRIVIEW from which patches (a) of simulated activity (b) were generated. 100 passes of spontaneous activity or noise (c) were identified using CTF software (Data Editor) and averaged together using MEGAN. The simulated activity was embedded within the averaged noise file (d) and saved in netCDF format (i.e., a netMEG file in MEGAN). Reprinted from Aine et al. (2012) with permission from Springer
Table 2 shows actual source locations, CSST estimated source locations, and errors when either noise was absent (no-noise) or empirical noise was present for visual simulated data Set 4. CTF head-centered coordinate system is used, where -x points out the back of the head, +y points out the left ear, and +z points out the top of the head. Average error across the 6 sources was 0.1 mm for the no-noise condition and 6.8 mm for the real noise condition. Standard deviation (SDev) is shown for estimated solutions for real-noise simulated data. This table demonstrates that the presence of real noise significantly affects source localization accuracy; however, our CSST solution for the real noise condition was still good for this complicated dataset, and inconsistent with previous critiques of dipole-modeling approaches that state dipole methods cannot accurately localize more than a few point sources of activity. Further, Table 3 lists CSST output when varying the model order (i.e. number of fitted dipoles) for a 3-dipole simulated dataset. The solutions (1–4 Dipoles) shown are for real spontaneous noise. Timecourses (shown as absolute values, bottom) are from the 4-dipole fit to 3-source data. In Table 3 the entries P1, P2, and P3 correspond to the Pk 1, Pk 2, Pk 3 timecourses. Notice that the noise timecourse is low-amplitude and without structure. As this table shows, under-modeling (1- and 2-dipoles) results in large localization errors. In contrast, localization errors are often reduced when over-modeling by 1 dipole (i.e., 4-dipoles for this 3-source dataset). Fortunately, noise sources are often easy to identify by a lack of timecourse structure and low amplitude (lower right panel).
Table 2
Actual and CSST estimated (“no-noise” and “real-noise”) locations for a 6-source, realistic simulation
SET 4 6 sources
Source V3
Error (mm)
Source I. LOG
Error (mm)
Source IPS
Error (nm)
X
Y
Z
 
X
Y
Z
 
X
Y
Z
 
Actual
−70.0
5.9
75.8
 
−59.7
33.2
42.9
 
−22.1
38.3
82.6
 
No noise
−69.7
6.0
75.9
0.3
−59.8
33.3
42.9
0.1
−22.1
38.2
82.7
0.1
Real noise
−61.3
4.3
74.1
9.0
−55.6
31.7
44.5
4.6
−18.7
28.7
71.8
14.8
SDev (Real)
0.3
2.3
1.6
 
1.6
0.3
0.4
 
1.5
1.7
0.9
 
 
Source R. frontal
 
Source AC
 
Source S. LOG
 
 
X
Y
Z
 
X
Y
Z
 
X
Y
Z
 
Actual
58.1
−41.5
46.2
 
74.1
−7.0
47.8
 
−31.3
−40.7
60.3
 
No noise
58.1
−41.7
46.2
0.1
74.0
−7.1
47.6
0.2
−31.4
−40.8
60.3
0.1
Real noise
58.5
−43.2
44.0
2.8
72.6
−9.9
46.7
3.4
−27.5
−36.1
59.1
6.1
SDev (Real)
0.1
0.1
0.4
 
0.5
0.3
0.2
 
0.5
0.6
0.6
 
Table 3
Sample output from an automated routine for determining best-fits to 3-source simulated data
 
Source location
Loc error mm (STD)
Peak amplitude error nAm
Peak latency error ms
Avg Loc error mm
   
Pk 1
Pk 2
Pk 3
Pk 1
Pk 2
Pk 3
 
Real spontaneous noise
1 Dip-V3
17.5 (0.15)
12.4
42.5
19.6
3.0
1.0
18.0
17.5
2 Dip-Vl
9.25 (0.11)
1.6
4.7
1.0
4.0
2.0
2.0
 
2 Dip-IPS
7.22 (0.08)
5.1
21.5
11.0
7.0
4.0
18.0
8.23
3 Dip-V1
4.93 (0.13)
1.1
0.45
2.2
5.0
2.0
1.0
 
3 Dip-V3
4.98 (0.12)
1.9
7.1
4.1
13.0
0.0
25.0
 
3 Dip-IPS
2.32 (0.05)
1.3
4.6
2.8
6.0
2.0
15.0
4.08
4 Dip-V1
3.11 (0.14)
0.03
3.6
0.80
4.0
2.0
6.0
 
4 Dip-V3
3.51 (0.14)
3.2
5.1
2.5
1.0
1.0
29.0
 
4 Dip-IPS
1.56 (0.05)
1.4
5.6
3.6
5.0
2.0
18.0
2.73
4-Dip-N
Noise
 
A272390_1_En_14_Figa_HTML.gif
Set 6 (remaining sets are not shown in Table 1) includes late activity (e.g., 400–600 ms) that was synchronous across four cortical sites (V1, I. LOG, IPS, and DLPFC), as is seen in working memory studies (Aine et al. 2006). The upper left panel of Fig. 2 displays the locations of the cortical patches (cortical patches are located at the cross-hairs) while the timecourses assigned to the cortical patches are shown beneath the MRIs. The averaged waveforms (128 trials with signals embedded in real spontaneous noise) seen across the 275 channels of the CTF MEG system are shown in the middle left column. CSST source locations are shown in the upper right panel (see tabled values). The table shows the coordinates of the actual sources, the estimated source locations, and the errors using Euclidean distance. Net source orientation errors were 42.0° for V1, 58.2° for I. LOG, 20.9° for IPS and 48.0° for the DLPFC sources. However, summarizing absolute orientation error is challenging since the original sources consisted of patches of cortex with the orientation of the patch activity conforming to the cortical folds. The middle right panel shows the estimated timecourses and source locations. The average localization error across all 4 sources was 6.7 mm with the greatest error for the I. LOG source. The cross-correlations between timecourses are shown in the bottom row of Fig. 2. We examined early activity first (200–350 ms–bottom left panel) which shows that V1 activity correlated highly with I. LOG, regions showing the initial spike-like activity (~280 ms). IPS and DLPF cross-correlations were also highly correlated with near zero-lag. The maximal correlation coefficients of the other pairs of sources were lower in value and were not near zero-lag. In contrast, the late activity (350–600 ms—bottom right panel) shows higher zero-lag correlation coefficients for activity between the 4 brain regions (i.e., late activity was synchronous across brain regions) with IPS and DLPFC revealing the highest correlation coefficient. This dataset is also suitable for examining coherence either between sensors or between reconstructed sources.
Apr 2, 2016 | Posted by in MAGNETIC RESONANCE IMAGING | Comments Off on Web Portal: A Database of Realistic Simulated and Empirical MEG Data for Testing Algorithms

Full access? Get Clinical Tree

Get Clinical Tree app for offline access