皮层源分析已成为脑活动分析的重要工具。源分析的目的是重建头皮上脑电图信号的皮层发生器(源)。源重建的质量取决于正问题的精度和反问题的精度。通过头皮上传感器位置的3，使用适当的成像模式来描述头部的几何形状D当确定头部模型的每种组织类型的真实导电性值时，地图可以确定准确的电极位置。这些参数有助于定义真实的头部模型。在这里，我们描述了重建记录在头皮上的脑电图信号皮层发生器的必要步骤。我们提供与事件相关的电位(ERPs)在一个6个月大的婴儿身上进行面部治疗的例子。我们讨论了不同的用途ERP调整源分析所需的措施。该方法可用于研究不同年龄受试者的认知任务。
脑电图（EEG）皮层源分析已成为定位整个生命周期脑功能反应的重要神经影像学工具。功能磁共振成像（fMRI）它是最广泛应用于认知任务中定位大脑功能反应的技术，最近被用于儿科人群的研究。最近，研究人员利用功能磁共振成像（fMRI）研究清醒婴儿和儿童在执行视觉任务、注意力和统计学习任务时的神经反应。这些有前途的首次应用fMRI该协议仍受高磨损率和小样本量的影响。另一方面，晚期儿童和老年参与者的扫描成功率较高。
脑电图是婴儿研究中最重要的神经成像方法之一，允许其高时间分辨率EEG捕捉动态神经反应。此外，EEG毫秒量级事件可以在不受影响的情况下捕获高时间分辨率fMRI依赖于缓慢的血流动力学响应(~秒)限制。通过应用源分析分析计算来识别记录在头皮上的神经发生器EEG空间分辨率差的特点。使用头部头皮记录进行皮层源分析，并通过计算推断出在皮层中产生信号的来源。因此，根据脑电图信号提供的功能活动(如振幅和潜伏期的变化)，fMRI研究大脑发育的源分析方法可以通过结构解剖和先进的数据处理和分析来确定。源定位分析可分为偶极子法或分布式源法。基于脑电位可以模拟为单个或多个离散偶极子的假设，其位置由先验定义。该方法已应用于认知过程的研究，包括儿童和成人晚期听觉诱发电位的发育变化。在另一项研究中，偶极子源模型用于研究6个月大婴儿的听觉变化检测机制。在颞上和额叶区域定位声音识别过程中产生的负和正峰值。在一项类似的研究中，研究人员评估了6个月大的婴儿在处理语音信息时大脑活动的时间过程。结果表明，婴儿的反应与成人相似，但激活时间不同。与颞区活动相比，前扣带皮层的早期激活可能取决于边缘和皮层区域的成熟差异，这对语音功能反应的分布有重要影响。
分布式源建模在定义源空间上的空间分布，分布式源建模方法通常局限于皮层表面。该分析计算了分布在头皮活动模式下的大脑区域的电流密度。一些研究实施了儿童脑电图记录的分布式源方法，旨在调查广泛的认知任务。源分析可作为评估脑电功能连通性的一种方法。脑功能连接通常在头皮-电极水平上完成，并通过相位耦合显示不同电极之间的相关性。源级可以完成与脑电的功能连接。这揭示了皮层源而不是头皮位置之间的连接。
总之，这些结果显示了源分析如何广泛用于认知任务中，并提供了关于大脑功能组织的进一步信息。儿科使用源分析的一个提醒是，需要真正的头部模型来反映个人头部的准确表征。这方面的黄金标准是利用个体参与者的结构核磁共振成像构建源分析的基本属性。这种方法可以解释头部几何形状的个体差异，这可能比成年人对儿童样本的影响更大，因为头部和大脑在年轻时的生长速度不同。比如大脑和头骨的关系在成人和婴儿之间有很大的不同。足月出生后，头盖骨被囟门和缝合线分开几个月。人颅骨的这一特征有助于大脑的生长，影响隔间的传导特性。皮质脊髓液在婴儿期和学步期(CSF)增加有助于改变大脑-头皮的距离，对源信号强度有重要影响。此外，如果使用真实的个体头部模型，反解会变得更加准确，尤其是在不能用球形模型描述的区域(如颞叶)和临床样本(如病变、肿瘤)的研究中。研究人员应对电极位置程序采取额外的预防措施。脑电图网络放置在不同的年轻参与者之间，因此应该避免使用标准的电极蒙太奇。
在本文中，我们总结了高密度脑电图皮层发生器的重建过程，并讨论了一些开发研究人员可供选择的解决方案。我们利用6个月大的参与者和12个月大的代表性婴儿的数据进行分布式源分析。将采用该源分析方法P400 ERP研究组件的振幅值。当单个结构磁共振成像或数字电极位置不可用时，我们将讨论源分析的替代方法。最后，我们将简要讨论非功能活动ERP源分析所需的调整措施(如频带、连通性)。
头皮上记录的脑电图电位变化是由皮层活动中突触后电位引起的。头皮上的电位变化是各种电源活动的总和，因此头皮上的电极位置与底层皮层源之间没有真正的关系。皮质源分析方法已发展成为确定皮质活动的定量工具(s)。一旦计算出正演模型成分，就可以构建逆空间滤波器来解决反问题。逆模型使用逆空间滤波器记录头皮的脑电图/ERP乘法运算，产生脑内皮层源的模式。这一步从最小的电极组中产生许多元素。这一过程被定义为逆问题，其解是待定的。将空间逆计算限制在理论分布范围内，通过选择参数来解决这个问题。
我们使用分布式源分析方法，不需要假设偶极子的数量；源空间是通过使用逼真的头部模型创建的，来自MRI ；反约束采用eLORETA方法。正反计算矩阵示意图如图1所示。图2显示源分析方法的可视化表示，可按以下步骤总结：
1）EEG和MRI预处理。如果不能准确进行。MRI可选择单体积的替代方案。
2)创建电极位置3D坐标，并与MRI头部体积共配准。
3)采用有限元法或边界元法MRI头部体积分割。
4)在源空间局限于灰质的情况下计算导频场矩阵。
5)精确低分辨率脑电磁成像(eLORETA)解决反问题的算法。
6)源估计和体素功率计算。
7）ROI源值统计分析中空间滤波器的应用。
图1 正(上)和反(下)模型示意图。这些元素对应于上一行的左右两侧EEG/ERP头皮记录、电极位置、参与者头部MRI以电流密度图表示的衍生几何形状、电导率、源空间和源发生器。
图2 源分析方法的示意图。 准确重构源的关键是从原始脑电信号中获取每个参与者的最小伪影数据。分割后的脑电图记录和平均值ERP数据集都可以作为源分析过程的输入。脑电预处理过程超出了当前工作范围。其它细节可以在其他地方找到。可靠的神经发生器可以通过脑电图信号重建，其中来自非脑源的噪声(如眼影)最小化。此外，考虑到头皮分布的个体差异，高密度脑电图系统开发更广泛的源发生器。在之前的研究中，我们提供了蒙太奇无法捕捉到的头皮电位表征，表明高密度脑电图记录。因此，建议使用高密度(如64通道或更高)和分布式蒙太奇来实现源定位技术。建立真实的头部模型是准确定位头皮脑电分布神经发生器的重要步骤。儿科研究参与者不能直接应用成人解剖学和参数。单个MRI使用扫描克服了使用成人头部模型分析儿童脑电图数据来源造成的定位错误。
头部模型源分析的金标准是单个参与者MRI。然而，fMRI获取具有挑战性，特别是对于非常年轻的参与者，因此限制了获取高质量图像的可能性。当个体核磁共振不能使用时，可以使用两种方法来获得真正的头部模型。首先，可以免费提供MRI数据库(例如，神经发育MRI数据库）。该解决方案需要收集参与者头部的精确尺寸，并在一组核磁共振成像中找到理想的匹配。为了尽可能准确地描述复杂的头部形状，我们建议获得多个头部测量值。头部测量的示意图表示在图3中。头部周长、前后距离、左右距离可用柔性不可拉伸卷尺和头部卡尺精确测量。第二，个体MRI替代选择是使用适合年龄的平均模板。从多个受试者那里获数据构建磁共振成像模板。对于发展研究，选择MRI为了捕捉发育中大脑的适当特征，模板包括被调查年龄范围内的参与者。因此，即使没有参与者的头部测量数据，年龄合适的平均模板也可以作为头部模型构建的替代品MRI。
图3 为了选择接近大小的头部测量示意图MRI。头部的宽度和长度可以用卡尺(顶板)获得。柔性胶带可以测量头围(黄线)、矢状(橙线)、冠状(红线)、前(蓝线)和后(绿线)的半围。
图4 与单个MRI电极位置的3D渲染视图(右下角的面板)。参与者的特定电极位置可以从放置在网络上的图片(红框和电极)或放置在所有电极上的数字表示(绿框和电极)来确定。 电极位置与MRI体积的准确配置是逼真头部模型的重要要求。由于头部形状的特殊差异，特别是在婴儿和儿童时期，电极的平均位置不能反映实际的网状位置。因此，应实施定位系统来确定头皮上电极的分析位置。许多设备在参与者头部使用不同的方法3D在空间中定位电极，如电磁场或图片。在没有定位系统的情况下，参与者的网状图片可以用来在三维空间中定位电极。头部周围几个电极的位置可以作为重建所有电极位置的基准标记。我们建议128通道Hydrocel测量地线传感网(HGSN)顶点、中耳、小翼、左右耳前分别使用5个基准标记Cz、17、75、57和100个电极。可以拍摄放置网的照片，定位受试者头部的基准标记，手动放置MRI上。基准点的坐标用于重建剩余电极的位置，并共同注册MRI头部。我们为单个人提供了一个MRI标记基准的图片示例(图4)。
通过引入场矩阵得到源分析的正解。铅场矩阵(图1)根据源模型、头部模型和电极放置位置计算。头部模型描述了考虑的隔间的头部几何形状和电导率值。个性化的头部模型来自参与者MRI适合年龄的电导率值的精确分割和使用。进行了一些模拟研究，以确定头部模型对源定位精度的影响。主要结论是，与简化的头部模型(如三腔室同心圆)相比，应优先考虑用真实的头部模型，并且至少应建模5个腔室。
在确定震源定位所模拟的不同隔室的电导率值时，还应考虑其他因素。不同头部介质的组成不同，会导致组织电导率的年龄变化。成人和幼儿的脑颅骨电导率的估计有很高的变异性，而其余组织的电导率则相对一致。颅骨室因其渐进性骨化而发生重要变化。此外，婴儿时期的颅室传导率在其表面并不均匀，直到2岁左右才完全闭合。已经实施了几种方法来估计幼儿的最佳导电性值。
当发生重要的发育结构变化时，头部组织间室电导率值建模的不准确性对电活动定位的影响需要进行进一步的研究。有限元方法模型中的每个介质都可以用六面体或四面体网格表示。在对成人参与者的研究中，六面体和四面体网格都表现良好，并显示出最小的源分析精度差异。与直接由MRI分割过程生成的六面体模型相比，四面体化过程更加复杂和耗时。在这里，我们在头部模型中区分了10个不同的区隔（即头皮、颅骨、硬脑膜、肌肉、眼睛、鼻腔、脑脊液、GM、WM和无髓鞘轴突），并报告了两种网格类型获得的来源结果。头部模型的六面体和四面体网格表示如图5所示。
图5 单个MRI分割输出用六面体(面板a)和四面体(面板b)网格表示。 电极的电位是根据头部模型配置计算的。可以解释测量信号在源分析中被重建。反问题的产生是由于电极与反源模型之间的病态关系，即很少有电极记录来自多个源的活动。有许多建议的逆方法程序（例如，最小范数估计）。我们使用分布式源分析方法，不需要假设偶极子的数量。源空间是通过使用逼真的头部模型来创建的，该模型来源于MRI 。
源激活在源空间的每个体素中用一个由三个参数组成的向量表示，空间中的每个维度都有一个参数。然后，通过计算每个体素位置的三维向量(这里称为CDR方法)的大小，将当前密度重建减少到单个值。采用这种方法，无论脑电信号的极性如何，原始脑电数据中的高值都会导致CDR高值。或者，奇异值分解(SVD)可以用作投影方法。与CDR方法不同的是，重构源保留了脑电活动的极性。SVD方法被认为可以在ERP源分析中提供稳健的结果，并在进行时频分解和连通性分析的研究中具有优势。
收集1例6月龄婴儿(M, 186天)的数据。报告了5个具有代表性的12个月婴儿(1名男性，年龄为376.4天)的聚合结果。所有婴儿都是足月(至少妊娠38周，出生体重至少2500克)，出生时健康，没有已知的发育异常。参与者主要是白人和中等社会经济地位。根据南卡罗莱纳大学机构审查委员会的伦理批准，获得了知情的父母同意。所有婴儿都被动地观看垂直和倒置的脸和房子，500毫秒内随机呈现，没有替换。对刺激序列使用了可变的试验间隔，当婴儿将目光从屏幕上移开时，芝麻街视频被作为吸引注意力的元素呈现。在脑电任务开始前进行头部测量。在EEG任务结束时，通过GPS拍摄图片和电极位置(见图4)。
受试者的MRI体积在Siemens Magnetom Prisma 3.0 T扫描仪上单独采集，采用20通道头线圈。采用MP-RAGE协议，使用以下参数获取全头t1加权图像:重复时间(TR) = 3000 ms，回波时间(TE) = 4.44 ms，翻转角度(FA) = 9◦，视野(FoV) = 192 × 192像素，体素尺寸= 1.0 × 1.0 × 1.0 mm。总采集时间为9 min 36 s。使用持续4 min和32 s的SPC序列获得t2加权图像(TR = 4500 ms, TE = 13 ms, FoV = 192 ×192 pixels，体素大小= 1.0 ×1.0 ×2.5 mm)。将MRI分割成10种不同的组织类型，获得了一个逼真的头部模型(图5)。每种组织类型的电导率值如下:头皮0.35 S/m，颅骨0.0132 S/m，脑脊液1.79 S/m, WM 0.2 S/m, GM 0.33 S/m，硬脑膜0.33 S/m，肌肉0.35 S/m，眼睛0.5 S/m，鼻腔0.0048 S/m。
根据解剖立体定位图谱定义解剖ROI。这些图谱被用来定义17个与婴儿面部处理的神经来源假说有关的解剖区域。所选ROI的3D渲染表示如图6所示。
图6 不同图谱。每个图谱都登记在6个月大的MRI中。选择枕颞和顶叶感兴趣区域。
图7展示了典型的6个月婴儿顶叶和顶枕通道ERP线图。P400在刺激发作后约400-600毫秒之间，在右侧电极（例如P10和PO10）上的振幅较大。、类似的头皮图对面部的响应比房屋显示出更大和右侧的活性。
在源网格的每个点上计算估计的源，这里限制在GM和眼睛。然后将反解应用于数据，并在网格的每个点计算CDR和SVD值。生成的三维矩阵的维度由网格中的点数、条件数和时间点数定义。
图7 典型6个月婴儿顶叶和顶枕通道ERP线图。整个片段的ERP活动被绘制成刺激类型的函数(顶部)。头皮图谱显示了所有刺激类型的P400成分的峰值活动(即刺激开始后432 ms)。 可以将一组特定的ROI应用于重建的源活动，并用于进一步的分析。ROI集可以从现有的立体定向地图集或fMRI任务中受试者的功能反应中创建。图8以3D绘制图的形式报告了P400峰值对直立面响应的分布活性结果。线图报告了在特定ROI中四种实验条件下的源活动。颞上沟和后扣带回的活动与ERP反应模式相似。重构的P400 的ERP在峰值时似乎存在一种特异的反转效应。通过使用接近参与者头部大小的MRI和适合年龄的MRI模板，我们得到了类似的解决方案。六面体网格计算结果见图9。
在图10中，我们比较了三种不同头部模型获得的后扣带回的源活动。三种解决方案的反应模式相似。自我核磁共振和近距离核磁共振的结果比年龄相适应的核磁共振平均数据的结果更相似。然而，只有自核磁共振和近核磁共振解决方案似乎是敏感的反转效应。这些结果需要在更大样本量的研究中进一步探索。总之，使用根据参与者头部特征和电极位置量身定制的头部模型进行源分析程序将提供更准确的解决方案。尽管如此，用适合年龄的头部模型获得的源解决方案被认为比使用成人头部模型更合适。最后，我们使用适合年龄的MRI模板和头部模型的六面体网格，将源分析程序应用于一小群12个月大的受测者中。在受测者层面应用正解和逆解，并在参与者之间平均得出结果。图11描述了ERP和CDR值的聚合结果。
图8 在单个t1加权MRI上绘制P400 ERP峰值对直立面孔的响应活动。(a)使用报告的两种网格类型(六面体，四面体)和源测量(CDR, SVD)的单个头部模型重建源活动。(b)线图描绘了四种实验条件下整个试验期间的源活动。
图9 (a)近磁共振头部模型(左面板)和平均模板头部模型(右面板)源分析结果的线图和3D渲染图。(b) 3D渲染图显示重构源在接近或年龄合适的t1加权MRI上绘制的P400峰值面响应。©线图显示了四种实验条件下的源活动。
图10 受测者后扣带回的CDR活动与MRI类型(图a)和面部取向(图b)有关。注意，y轴刻度被设置为最大化每个分布。400 - 600 ms之间的平均振幅作为面向的函数绘制。在P400时间窗内，无论是自我MRI还是近距离MRI，人脸反转效应似乎更明显。
图11 从12个月大的小样本受测者中收集的数据。顶部面板报告了四种实验条件下的顶叶和顶叶枕通道的ERP反应，以及直立面孔P400峰值的头皮图谱。底部面板显示了四个ROI中的CDR活动作为刺激类型的函数，以及P400峰值处垂直面的体素源图的3D呈现。
本文概述了脑电图源分析的过程，并特别参考了使用儿科人群获得的数据进行此类分析所需的调整。脑电源重建的关键是获得精确的电极位置图和逼真的头部模型。所有电极位置的数字地图(即GPS解决方案)是获得所有电极的3D坐标的最精确的方法，这是非常小的婴儿可以接受的程序。几个基准电极的数字坐标可以用来重建剩余位置的位置。当没有数字化电极或基准位置时，训练有素的研究人员可以使用参与者的脑电图网络放置图片，手动将基准电极放置在MRI上。我们建议在脑电图记录结束时从头部的各个侧面拍摄照片。根据成人研究，个性化电极放置要优于标准电极放置，以减少源定位误差。
头部模型的建立是源分析研究的另一个关键步骤。使用逼真的头部模型有利于最小化源定位误差，而不考虑所调查的年龄组。然而，儿童MRI和MRI模板的缺乏导致发育研究人员使用成人头部模型。这种方法增加了定位误差的解剖差异，不同年龄的参与者。为此，至关重要的是，在实验期间采取准确的头部测量。当无法获得头部测量值时，可使用适合年龄的MRI模板。所有参与者使用相同的头部模型和源空间来计算头皮脑电活动的神经发生器。因此，潜在的个体差异可能会被忽视，导致定位误差的增加。然而，使用平均模板构建头部模型所产生的小定位误差可以通过在预先定义的解剖ROI内平均源活动来补偿。我们与一位有代表性的参与者的比较结果表明，与年龄相适应的头部模型通过自身MRI、近距离MRI或MRI模板的形式在后扣带回中产生了相似的源激活模式。进一步的研究是必要的，以更好地量化由三种替代头部模型解决方案产生的定位误差。
总之，源分析的金标准方法利用了参与者MRI中电极位置的数字地图和个性化头部模型。可以采用其他方法来克服每个参与者获得电极位置和MRI的限制。使用神经导航系统或网状布局的图片，可以从几个电极标记的坐标重建电极坐标地图。可以使用接近参与者头部大小的MRI或与年龄相适应的MRI模板来创建头部和源模型。
参考文献：Cortical Source Analysis of Event-Related Potentials: A Developmental Approach