基于个体特征的经颅脑图谱生成方法、导航方法及其系统与流程

本发明涉及一种个体经颅脑图谱生成方法，同时涉及一种基于个体特征的经颅脑图谱导航方法，还涉及相应的个体经颅脑图谱导航系统，属于认知神经科学技术领域。

背景技术：

当前，认知神经科学所面临的一个重要课题是建立脑功能和脑结构之间的对应关系。以功能磁共振成像(fmri)、功能近红外光谱(fnirs)为代表的功能脑成像技术，可以使神经科学的研究者们以非侵入的方式对活体人脑的功能进行观测。

包括经颅脑治疗和经颅脑成像在内的经颅脑映射技术，一直在快速发展，并在脑机制评估和脑损伤治疗方面表现出巨大的潜力。例如，经颅磁刺激(transcranialmagneticstimulation，简称为tms)通过外部施加的磁场，实现在高时空分辨率下对人脑局部活动的抑制或激发。经颅磁刺激的无创特性使其有助于评估特定脑区域/脑回路和行为之间的因果关系。经颅磁刺激也广泛用于治疗各种神经精神障碍，例如帕金森病、疼痛、成瘾等。值得一提的是，经颅磁刺激作为治疗耐药性抑郁症的临床干预手段，已经获得美国食品药品管理局(fda)的认证。

经颅脑映射技术通常在可见的头皮表面(称为操作空间)上操作。然而，经颅脑治疗和经颅脑成像的目标区域位于大脑内部(称为效用空间)，从操作者的角度来说，这是从外部不可见的。因此，给定目标位置或效用空间中的某个脑区域标号，在个体头皮表面上定位最佳的成像或刺激位点是非常困难的。这可能导致次优化、不一致的实验结果，甚至冲突的结论。例如，实施经颅磁刺激的线圈在治疗抑郁症的目标区域(例如背侧前额叶皮层)上的放置不一致可能导致治疗结果不同。类似地，fnirs光极的放置也是至关重要的。不适当的放置可能导致记录皮质位置的偏差，甚至可能覆盖错误的皮质区域。此外，在神经影像研究中的群组分析中，需要将来自不同个体的成像皮质位置彼此对应以进行不同个体之间结果的比较。目前，将fnirs光极正确放置在头皮上以覆盖感兴趣的皮层区域并且同时保持成像皮层位置在不同受试个体之间的对应性仍然具有挑战性。其中，可见的操作空间(头皮表面)和不可见的效用空间(大脑内部)之间的分离仍是以有效和准确应用这些经颅脑映射技术的最大挑战之一。

技术实现要素：

本发明所要解决的首要技术问题在于提供一种个体经颅脑图谱生成方法。

本发明所要解决的另一个技术问题在于提供一种基于个体特征的经颅脑图谱导航方法。

本发明所要解决的又一个技术问题在于提供一种个体经颅脑图谱导航系统。

为了实现上述目的，本发明采用下述的技术方案：

根据本发明实施例的第一方面，提供一种基于个体特征的颅面坐标系生成方法，包括如下步骤：

(11)在个体的头皮表面识别五个颅骨标记nz、iz、al、ar和cz；

(12)将头皮表面和通过nz、cz和iz的平面之间的相交曲线定义为颅骨赤道；

(13)给出头皮表面上的点p，经度曲线可以唯一地确定为头皮表面和通过al、ar和p的平面之间的相交曲线，p'是颅骨赤道与经度曲线之间的交叉点；

(14)上头皮的任意点p由一对非负实数(pe，pl)唯一确定：

pe＝lnz-p’/le，pe∈[01]

pl＝lal-p/lal-p-ar，pl∈[01]

其中，lnz-p'是沿着颅骨赤道的nz到p'的曲线长度，le是颅骨赤道的全长；lal-p是沿着全长为lal-p-ar的经度曲线从al到p的曲线长度。

其中较优地，上述颅面坐标系生成方法还包括步骤(15)：在标准半球上建立cpc空间；使用hammer-aitoff投影对带有cpc坐标的半球进行平面化，生成在扁平椭圆上呈现的cpc坐标系的地图。

根据本发明实施例的第二方面，提供一种基于个体特征的经颅脑图谱生成方法，包括如下步骤：

(1)按照上述的步骤在个体层面创建颅面坐标系；

(2)建立用于连接颅骨位置与大脑位置的经颅映射系统；

(3)使用马尔可夫链中的两步随机过程构建经颅脑图谱。

其中较优地，所述步骤(2)包括如下子步骤：

使用气球膨胀模型在个体空间中确定对应于给定的头皮表面任意点p的底层皮层位置c；

在所有皮层位置被空间标准化为mni空间之后，聚合所有(p，c)对，生成确定性的个体经颅脑映射模型。

根据本发明实施例的第三方面，提供一种基于个体特征的经颅脑图谱导航方法，包括如下步骤：

(1)在被试者头皮表面上进行若干个点的稀疏均匀采样，重建出被试者头皮表面的近似曲面；

(2)在所述近似曲面上建立cpc坐标系；

(3)将被试者头皮表面上移动的探笔位置实时转化到上述经颅脑图谱所在的cpc坐标系中，利用经颅脑图谱的信息指导经颅设备在头皮表面的放置。

其中较优地，所述步骤(1)中，采样点包括两个部分，第一部分是4个颅骨标记nz，iz，al和ar，第二部分由操作者手动在被试者头皮表面采集得到。

其中较优地，将基于所述采样点生成的离散点云表示在鼻凹-耳前坐标系统中；其中，鼻凹-耳前坐标系统的原点o定义在al-ar连线的中点，o点到nz的方向为x轴正方向，o点到al的方向为y轴正方向，而z轴正方向被定义在xy轴的叉积方向。

其中较优地，采用crust算法重建所述离散点云之间的拓补连接，得到三角面片集合tp{(i,j,k)|i,j,k＝1,2,...m}；其中，i，j，k分别是p中点的序号，每个三元组(i，j，k)表示具有拓扑相邻的三个点。

其中较优地，对所述三角面片集合tp{(i,j,k)|i,j,k＝1,2,...m}中的所有三角面片进行非线性插值，得到细分的三角网格；保留所述三角网格上的顶点，记作致密点云dp＝{xi,yi,zi|i＝1,2,3,...,ndp}，作为所述近似曲面。

根据本发明实施例的第四方面，提供一种个体经颅脑图谱导航系统，包括处理组件、显示组件、存储器、三维定位组件以及输入输出接口；其中，

所述处理组件分别连接三维定位组件、显示组件、存储器和输入输出接口，用于执行上述的经颅脑图谱导航方法中的全部或部分步骤；

所述显示组件与三维定位组件连接，用于实时显示三维定位组件在被试者头皮表面的当前状态及经颅脑图谱中的状态；

所述存储器，用于存储各种类型的数据；

所述三维定位组件，用于利用从被试者头皮表面采集的若干个点，重建出被试者头皮表面的近似曲面；

所述输入输出接口，用于实现所述个体经颅脑图谱导航系统与外部经颅设备的连接。

其中较优地，所述三维定位组件为fastrak磁场三维定位仪，包括一个磁场发射器和两个磁场探测器；其中一个磁场探测器为探笔形状，另一个磁场探测器是立方体形状。

其中较优地，所述磁场三维定位仪在使用时，所述磁场发射器固定在被试者身后的支架上，立方体形状的磁场探测器被固定在被试者颧骨位置，探笔形状的磁场探测器由实验者持握并在被试者头皮表面上标定位置。

与现有技术相比较，本发明首先提供了一种个体经颅脑图谱生成方法，可以将不可见的大脑内的图谱信息投射到可见的头皮表面上，使原本分离的操作空间和效用空间融合在一起。针对将经颅脑图谱应用到指导经颅设备的个体化放置中的问题，本发明进一步提供了一种个体经颅脑图谱导航系统。实验结果表明，利用本发明可以有效提高近红外测量光极感兴趣脑区的覆盖准确性，以及在不同人上皮层测量位置的一致性，从而提高近红外技术检测任务诱发脑活动的敏感性。

附图说明

图1(a)～图1(e)是经颅脑图谱的系列实施例示意图；

图2是从磁共振图像中识别颅骨标记的示意图；

图3是cpc坐标系的示意图；

图4给出了单个cpc坐标点的概率经颅脑映射示意图；

图5为本发明所提供的个体经颅脑图谱导航系统的结构示例图；

图6为磁场发射器参照系下，磁场探测器的空间位置和姿态表示示意图

图7为物理坐标系到头坐标系的空间变换示意图；

图8为鼻凹-耳前坐标系统的示意图；

图9为头皮表面重建过程的示意图，其中a显示了头皮表面采集的21个无结构散点，b显示了crust算法计算得到采样点的拓补关系，c显示了球面化处理得到的重建点云；

图10为本发明所提供的个体经颅脑图谱导航系统的使用状态示意图，其中a为显示组件所显示的屏幕截图，b为实验者进行实际操作的示意图。

具体实施方式

下面结合附图和具体实施例对本发明的技术内容做进一步的详细说明。

前已述及，对经颅成像装置可见的颅骨表面空间和对经颅成像装置不可见的颅内脑空间的分离是有效应用经颅脑映射技术的最大挑战之一。为了解决两个空间的对应问题，本发明首先提出了经颅脑图谱(transcranialbrainatlas，简称为tba)的概念。经颅脑图谱是建立在头皮表面上的脑图谱，它将不可见的大脑内的图谱信息投射到可见的头皮表面(尤其是头皮上表面)上，使研究者或者医生可以直接使用这些涉及脑结构和脑功能的图谱信息。

具体地说，本发明首先显式地构建一个标准的颅骨坐标系统，用来实现对不同个体颅骨表面空间的定量描述。然后，依据颅-脑对应关系在人群上一致性的假设，建立一个标准颅骨表面空间到脑图谱所在的标准脑空间之间的对应关系。最后，本发明从颅骨表面空间-标准脑空间，以及脑图谱提供的标准脑空间-脑分区标号空间这两部分的对应关系中，求解出颅骨表面空间-脑分区标号空间之间的对应关系。作为结果，本发明将对应的脑标准空间和脑图谱中的信息反向呈现到颅骨坐标系统，从而形成一种新型的“经颅脑图谱”。经颅脑图谱的一个重要性质是能够仅仅通过颅骨位置信息，直接推断对应脑分区标号的信息，因此可以粗略地理解为一个建立在头皮表面上的脑图谱。经颅脑图谱本质上是一个建立在坐标化的头皮表面上的大脑功能地图，即在一个坐标化的大脑空间中，让传统脑图谱对应每个脑部位置及其功能或解剖学标号，从而描绘出经颅脑映射技术可及的各个皮层位置及其对应的图谱标号，并将它们明确呈现在作为可见操作空间的头皮表面上。

经颅脑图谱可以将传统脑图谱中的先验脑分区信息，在人群水平颅-脑对应关系意义下映射到用于经颅脑成像装置放置的颅骨空间中。因此，可以看作是传统脑图谱在经颅脑成像技术领域的扩展。在经颅脑图谱的框架下，经颅数据在脑空间的定位可以转化为对经颅成像装置在颅骨空间的定位，这使得经颅脑映射技术的实时定位成为可能。此外，由于经颅脑图谱的标号信息是在颅骨空间显示的，这个特性非常有利于经颅脑图谱叠加到个体头皮表面进行显示，从而以直观的方式指导经颅成像装置在被试颅骨表面的放置。因此，经颅脑图谱的建立，解决了经颅脑映射技术中操作空间与效用空间分离的矛盾。

图1(a)～图1(e)是经颅脑图谱的系列实施例示意图。如图1(a)所示，假设在头皮表面上定义了统一的颅面坐标系，它描述了可以设置经颅脑映射技术的探头的所有可能的位置，并且将探头设置在给定的颅骨位置(图中黑点所在位置)。对于特定的个体(例如图1(b)中的受试者sub1)，设置在其头皮表面的探头将能够探及特定皮层位置/脑区域(图中黄点所在位置)。但在群体层面，考虑到个体之间的解剖结构的差异，这种颅-脑对应关系可能不是确定的。如图1(c)所示，在标准化到标准脑空间之后，我们可以捕获这种概率对应关系的空间分布(即黑色圆圈内的彩色区域)。通过给定从大脑图谱所获得的解剖信息，这种概率对应关系能够提供如何从颅骨位置访问每个大脑区域的群体层面的概率，并且将其作为先验知识。经颅脑图谱本质上是为颅面坐标系定义的整个大脑空间映射该先验知识。具体在图1所示的系列实施例中，如果仅针对每个颅骨位置考虑最可能探及的脑分区标号及其关联的概率，则如图1(d)所示的最大似然标号图(mlm)和如图1(e)所示的最大概率图(mpm)可以作为经颅脑映射技术的探头排列的有用指导。

在本发明的实施例中，经颅脑图谱的生成过程主要包括三个步骤：1)在个体层面创建颅面坐标系；2)建立用于连接颅骨位置与大脑位置的经颅映射系统；3)使用马尔可夫链中的两步随机过程构建经颅脑图谱。

下面对经颅脑图谱的详细生成过程进行说明。

1)在个体层面创建颅面坐标系(cranialsurfacecoordinatesystem)

颅面坐标系需要满足两个基本要求：第一，它应为个体头皮表面提供一对一的映射；第二，为了便于进行群体层面的研究，对于颅面坐标系中的每个位置，应该使来自不同个体的下层皮层位置在神经解剖学上基本一致。

cpc坐标系的基本思想是在头皮表面上构造一个类似于经纬线的坐标系统。与地理经纬线系统不同的是，cpc坐标系对“经纬度”的确定是利用两次表面比例测量的方式确定的。在本发明的实施例中，通过如下三个步骤在个体的头皮表面建立被称为连续比例坐标空间(简称为cpc空间)的二维比例坐标系：

11)在个体空间(naviespace)的头皮表面识别至少五个来源于10-20系统的颅骨标记nz、iz、al、ar和cz(参见图3(a))。在磁共振图像中识别颅骨标记的示例参见图2，其中，iz是人类头骨的外枕骨隆起，斜方肌附着于其上；al和ar是左耳和右耳前点，被识别为耳屏的峰值区域；nz被确定为鼻梁上部的凹陷位置；cz被确定为颅表面测地线al-cz-ar和nz-cz-iz的交点，并分别平分两条颅表面测地线。

12)将颅骨赤道被定义为头皮表面和通过nz、cz和iz的平面之间的相交曲线(即图3(a)中的曲线1)；

13)给出头皮表面上的点p，经度曲线(即图3(a)中的曲线2)可以唯一地确定为头皮表面和通过al、ar和p的平面之间的相交曲线，p'是颅骨赤道与经度曲线之间的交叉点。

在以上三步定义的基础上，上头皮(高于由nz、iz、al、ar点指定的曲线)的任意点p都可以由一对非负实数(pe，pl)唯一确定：

pe＝lnz-p’/le，pe∈[01](1)

pl＝lal-p/lal-p-ar，pl∈[01](2)

其中，lnz-p'是沿着颅骨赤道的nz到p'的曲线长度，le是颅骨赤道的全长(从nz到iz)；lal-p是沿着全长为lal-p-ar的经度曲线从al到p的曲线长度。如图3(b)所示，作为任意点的p点的表面位置由p'和p分别占这两条曲线的比例唯一地表示，计算公式参见公式(1)和公式(2)。

图3(b)是在头皮表面上建立的二维比例坐标系(简称为cpc坐标系)的示意图。其中，该二维比例坐标系为头皮表面上的任意点p向cpc空间提供一对一的映射。从上面的描述可以看出，cpc坐标系从本质上是在头皮表面上构造一个类似于经纬线的坐标系统。但与地理经纬线不同的是，cpc坐标系对“经纬度”的确定是利用两次表面比例测量的方式确定的。

基于颅骨标记(nz，iz，al，ar和cz)和cpc坐标系(与尺度和形状成比例)定义的比例关系所确立的对象间对应关系，可以在个体层面的头皮表面上建立合理的解剖学对应关系。为了从单个视角可视化整个头皮表面，如图3(c)所示，在标准半球上建立了一个特殊的cpc空间。然后，使用现有的hammer-aitoff投影对带有cpc坐标系的半球进行平面化，生成在扁平椭圆上呈现的cpc坐标系的地图，申请人将其称为bnu地图(beijingnormaluniversitymap)(参见图3(d))，其实质为标准cpc坐标系的二维投影图像。在bnu地图的基础上，任何与头皮表面相关的脑功能数据都可以在该地图中呈现，从而能够在不同项目、人群、实验室甚至不同成像模式之间实现有效的比较。

2)建立用于连接颅骨位置与大脑位置的经颅脑映射模型(transcranialbrainmapping，简称为tbm)

一旦在个体层面的头皮表面建立了cpc空间，使用成熟的气球膨胀模型(okamoto&dan，2005)，可以在个体空间(例如，个体3dmri图像)中确定对应于给定的头皮表面任意点p的底层皮层位置c。在所有皮层位置被空间归一化为标准脑空间(即mni空间)之后，聚合所有(p，c)对，可以生成确定性的个体经颅脑映射模型。然后，通过集成所有的个体经颅脑映射模型来生成群体层面上的概率经颅脑映射模型：

p(c|p)(3)

p(pe，pl)∈cpc，c(x，y，z)∈c，c是标准脑空间的子集，并包含经颅脑映射技术相关的所有皮层位置。当从给定坐标的头皮表面上任意点p(pe，pl)开始刺激或记录时，概率经颅脑映射模型给出每个靶向皮层位置c(x，y，z)的概率。

图4给出了单个cpc坐标所对应的概率经颅脑映射模型。当给出一对cpc坐标，例如p＝(0.4，0.6)，在每个个体层面的头皮表面上可以确定对应的点b，其所对应的皮层投射点c的位置使用成熟的气球膨胀模型在个体层面上的磁共振图像空间中识别出来。在图4中，下半部分的表格显示了任意点p＝(0.4，0.6)对应在概率经颅映射模型的不同位置的相应概率。

从颅骨空间一点到标号空间一个标号的映射，可以看做是一个两步映射。首先，从颅骨空间s映射到脑空间b，再从脑空间b映射到标号空间l。由于两步映射都是概率映射，因此，这个过程也可以看做一个两步随机过程。由于脑空间和脑分区标号的对应关系是由人脑自身的结构规律决定的，我们假定对于脑坐标空间中任意一点，其对应到各个脑区标号的概率是确定的，而与之前颅骨坐标空间到脑坐标空间的对应路径无关。因此，这个两步随机过程存在马尔可夫(markov)性质。

3)使用马尔可夫链中的两步随机过程构建经颅脑图谱(transcranialbrainatlas)

本领域技术人员知晓，脑图谱是在概率框架中构建的。例如，传统的脑图谱(例如mni图谱)描述的基本关系是条件概率：

p(l|b)(4)

其中，b(x，y，z)∈b，b是标准脑空间的子集，并且包含图谱的脑模板中的所有可能的脑组织点；l∈l，l包含所有可能的图谱标号，每个图谱标号都表示脑图谱中的特定脑区。对于给定的一对l和b，p(l|b)表示人脑中的位置b处出现图谱标号l可能性。

相应地，经颅脑图谱描述的基本关系也是条件概率：

p(l|p)(5)

其中p(pe，pl)∈cpc和l∈l。

在本发明的实施例中，可以通过马尔可夫链中的两步随机过程构建经颅脑图谱。具体地说，第一步：作为输入的给定点p(pe，pl)通过概率经颅映射p(c|p)被映射到标准脑空间中的皮层位置c(x，y，z)。第二步：特定的c(x，y，z)被映射到标号空间l中的标号l。研究者或者医生可以单独地通过皮层位置c(x，y，z)来预测图谱标号l，而无需考虑头皮表面上的位置p(pe，pl)，这里提供了公式(6)：

p(l|c，p)＝p(l|c)(6)

因此，两步随机过程使用了本领域技术人员常用的马尔科夫链。其中，假设皮层(域c)上的点是脑中(域b)的点的子集，p(l|p)可以使用公式(3)中的p(c|p)和公式(4)中的p(l|b)来计算。具体而言，如果p和c离散化，那么chapman-kolmogorov方程式表示：

其中，pi是cpc空间中的离散位置(即任意点，下同)，i＝1，2，...，np；cj是c在标准脑空间中的离散位置，j＝1，2，...，nc；lk是特定脑图谱的图谱标号，k＝1，2，...，n1。

需要说明的是，上述步骤所构建的经颅脑图谱是一种概率图谱，即当探头在给定坐标的头皮表面的任意位置p刺激或记录时，通过经颅脑图谱可以给出由l标记的每个靶向脑区域被探及的概率。它将不可见的脑内的图谱标号信息投射到可见的头皮上来，使得研究者或者医生可以直接使用这些脑结构和功能图谱信息，大大提升了脑图谱在经颅脑映射技术中的作用。

综上所述，本发明首先提供了一个基于两步马尔可夫链模型的概率框架作为经颅脑图谱的理论基础。其中，第一步是从cpc空间中的头皮位置到mni空间中的下划线皮层位置的颅骨-皮层映射(crainio-corticalmapping)。第二步是使用传统的大脑图谱构建经颅脑图谱，实质是一个从mni空间的皮层位置到图谱标号空间的映射。需要说明的是，本发明利用概率框架提供了一个可扩展的经颅脑图谱模型，上述第二步所使用的大脑图谱可以由任何其他脑图谱替代。虽然，在本发明中仅提供了三种基于宏观解剖学的脑图谱(ba图谱，aal图谱和lpba图谱)用于构建经颅脑图谱，但类似的功能图谱、连接图谱以及其他图谱都可以被用来为特定应用提供功能性的经颅脑图谱。

在现有技术中，经颅脑成像技术所面临的一个首要问题，就是如何将头盖骨测量所得到的功能信息定位到脑图谱中，即经颅定位问题。另外，与传统脑成像技术不同的是，这类技术在实际操作中，还需要事先规划经颅设备(包括经颅脑成像装置或经颅脑治疗装置)在头盖骨上的放置位置。因此，经颅脑成像技术定位问题的另一方面，是如何利用脑图谱中的解剖信息指导经颅设备的个体化放置，即经颅导航问题。

针对上述两个问题，本发明的整体解决思路是这样的：首先在经颅设备可见的颅骨空间构建一个标准化的坐标系，用来定量刻画整个经颅技术所能覆盖的整个头盖骨空间。然后基于大样本结构像数据，求解头盖骨空间到标准脑空间以及头盖骨空间到脑图谱标号空间的在人群水平的对应关系。以此方式解决经颅定位和经颅导航两方面的问题。下面对此展开详细具体的说明。

在经颅脑图谱的框架下，经颅数据在脑空间的定位可以转化为对经颅设备在颅骨空间的定位，这使得经颅脑映射技术的实时定位成为可能。此外，由于经颅脑图谱的标号信息是在颅骨空间显示的，这个特性非常利于将经颅脑图谱叠加到个体头皮表面进行显示，从而以直观的方式指导经颅设备在被试者头皮表面的放置。

为了将经颅脑图谱应用到解决经颅脑成像技术的实际定位和导航问题中，本发明提供了一种基于稀疏采样点的表面重建技术，并在其基础上实现了面向颅骨坐标系的自动定位算法。该自动定位算法能够仅通过在被试者头皮表面上的数十个点的稀疏均匀采样，重建出被试者头皮表面的近似曲面，并在此曲面上实现被试者个体cpc坐标系的重建。在上述自动定位算法的基础上，本发明提出了一种基于颅骨坐标系和经颅脑图谱的个体经颅脑图谱导航系统。该系统可以借助磁场三维定位仪的实时定位功能，将被试者头皮表面移动的探笔位置实时转化到经颅脑图谱所在的cpc坐标系中，以此实现计算空间的经颅脑图谱与物理空间的被试者头皮表面的融合，从而利用经颅脑图谱的信息指导经颅设备在头皮表面的放置。

如图5所示，本发明所提供的个体经颅脑图谱导航系统至少包括处理组件、显示组件、存储器、三维定位组件以及输入输出接口。其中，处理组件分别连接三维定位组件、显示组件、存储器和输入输出接口，用于控制整个个体经颅脑图谱导航系统的整体操作，执行本发明所提供的经颅脑图谱导航方法中的全部或部分步骤。在示例性实施例中，处理组件可以被一个或多个应用专用集成电路(asic)、数字信号处理器(dsp)、数字信号处理设备(dspd)、可编程逻辑器件(pld)、现场可编程门阵列(fpga)、控制器、微控制器、微处理器或其他电子元件实现，用于执行上述经颅脑图谱导航方法。显示组件与三维定位组件连接，用于实时显示三维定位组件在被试者头皮表面的当前物理状态及经颅脑图谱中的状态，以便指导经颅设备调整在被试者头皮表面的放置位置。存储器被配置为存储各种类型的数据以支持在个体经颅脑图谱导航系统上的操作。这些数据的示例包括用于在个体经颅脑图谱导航系统上操作的任何应用程序或方法的指令等。存储器可以由任何类型的易失性或非易失性存储设备或者它们的组合实现，如静态随机存取存储器(sram)，电可擦除可编程只读存储器(eeprom)，可擦除可编程只读存储器(eprom)，可编程只读存储器(prom)，只读存储器(rom)，磁存储器，快闪存储器等。三维定位组件用于利用从被试者头皮表面采集的若干个离散点，重建出被试者头皮表面的几何形状，并以毫米级精度实现将被试者头皮表面任意一点在颅骨比例坐标系统的定位。输入输出接口用于实现本个体经颅脑图谱导航系统与外部经颅设备，例如经颅脑治疗设备和经颅脑成像设备等的连接。除了上述功能模块之外，本个体经颅脑图谱导航系统还可以包括常规的电源模块及人机接口模块(例如键盘、鼠标等)，在此就不具体说明了。

在本发明的一个实施例中，为了获取物理空间中的空间定位，三维定位组件采用polhemus公司生产的fastrak磁场三维定位仪。这种磁场三维定位仪包括一个磁场发射器(transmitter)和两个磁场探测器(sensor)。其中一个磁场探测器被制作成探笔(stylus)的形状，另一个磁场探测器是立方体形状，为描述方便，我们分别称之为探笔和探测器。在磁场三维定位仪工作时，磁场发射器会产生一个充满局部磁场，这个磁场定义了一个半径为305cm的半球空间。在这个半球空间中，能够生成一个以磁场发射器为原点的参照系。在这个参照系中，探笔和探测器的空间位置和姿态都可以被定量表示。上述磁场三维定位仪在使用时，磁场发射器被固定在被试者身后的支架上，探测器被用弹性细绳固定在被试者颧骨位置，探笔由实验者持握并在被试者头皮表面上标定位置。

对于上述参照系中的一个磁场探测器，首先磁场探测器的空间位置可以由参照系的o-xyz坐标给出。然后，用欧拉角(eularangle)表述磁场探测器的姿态。另外，磁场三维定位仪也指定了一个原点和三个正方向，即定义了一个探测器特异的参照系o’-x’y’z’。假定最初探测器位于发射器参照系x正半轴上，且x’y’z’轴与xyz轴方向相同。则探测器可以通过三次旋转到达当前姿态。第一，在o-xy平面上，探测器沿z轴顺时针转过方位角(azimuth)ψ使得x’与当前x’在o-xy平面上的投影同向。第二，在o-x’z平面上，向z轴正反向转过仰角(elevation)θ，使得x’与当前x’方向重合。此时，z’轴和y‘轴位于图6中的z”和y”处。最后，在o’-y”z”面上，探测器延x’轴转过滚动角(roll)φ，使得z”和y”与当前的z’和y’轴重合。因此在物理空间中，磁场探测器的空间位置和姿态，可以由x、y、z、ψ、θ、φ六个维度唯一描述。

尽管磁场三维定位仪能够定位出物理空间中的位置信息。然而，在实际操作中，被试者的头部无法保持静止，因此采集到的头皮表面位置信息会受到被试者头动的影响。理想情况下，我们希望头皮表面位置的描述是建立在一个与被试者头部相对静止的参照系中。为了解决这个问题，我们将立方体形状的探测器与被试者头皮表面固定，形成一个相对头皮表面静止的参照系统，并且融合探笔和探测器的位置和姿态信息。

如图7所示，假定磁场发射器确定的坐标系为o-xyz，探测器参照系为o’-x’y’z’，探笔笔尖位置为p；此外我们额外定义一个参照系o’-xyz，该坐标系原点位于o’处但xyz轴方向与磁场发射器坐标系同向。已知p点在o-xyz坐标系下的坐标(x1,y1,z1)，o’点在o-xyz坐标系下的坐标为(x2,y2,z2)，姿态参数为(ψ,θ,φ)，我们需要求解p点在o’-x’y’z’的坐标(x3,y3,z3)。

首先，根据op，oo’和o’p的向量关系，我们可以得到

即

由于由经过三次旋转变换得到，存在变换关系：

使得

即

带入式(9)，得

于是有

(x3,y3,z3)＝(x1-x2,y1-y2,z1-z2)·a^-1(14)

式(14)给出了借助探测器坐标和姿态角将发射器参照系中的一点转换到头参照系的变换关系。借助这个变换关系，探笔笔尖定位到的位置，都可转换到随被试者头移动的头参照系。

上述经过磁场三维定位仪定位的头皮表面点，仍然表示在一个三维坐标系中。然而，个体经颅脑图谱导航系统的一个根本需求是将表示在三维立体空间中的点，在被试者个体的cpc坐标系中定位出来。为此，本发明采用基于稀疏采样的表面重建算法，利用从颅骨上采集到的零星散点，重建出整个颅骨的几何信息。然后，利用基于重建头皮表面的比例测量，将表示在头参照系中的颅骨点定位到被试者个体的cpc坐标系中。具体说明如下：

首先，我们对颅骨表面进行了采样。采样点包括两个部分，第一部分是4个颅骨地标点，即nz，iz，al和ar。第二部分被称为重建点，由操作者手动在被试者头皮表面采集得到。这里，4个颅骨地标点一方面和重建点一起被用在颅骨重建中，并定义了重建颅骨的边界；另一方面，在cpc测量中，4个颅骨地标点作为cpc系统构建的4个测量水准点，用于颅骨点的cpc坐标定位。在颅骨表面采样后，我们得到n个点的离散点云p{xi,yi,zi,i＝1,2,...,n}，其中n是包括4个颅骨地标点在内的所有颅骨采样点。为了方便表面重建算法的实施，我们对离散点云进行了一个预处理，将其表示在鼻凹-耳前坐标系统中。如图8所示，其中，鼻凹-耳前坐标系统的原点o定义在al-ar连线的中点，o点到nz的方向为x轴正方向，o点到al的方向为y轴正方向，而z轴正方向被定义在xy轴的叉积方向。

稀疏采点得到的结果是一个无结构的点云(图9中的a)，我们引入crust算法(参见amenta，n.，m.bern，andm.kamvysselis，anewvoronoi-basedsurfacereconstructionalgorithm，inproceedingsofthe25thannualconferenceoncomputergraphicsandinteractivetechniques.1998，acm.p.415-421.)重建出这些无结构点云之间的拓补连接。crust算法得到的结果是一个三角面片集合(trianglepatch)，tp{(i,j,k)|i,j,k＝1,2,...m}(图9中的b)，其中，i，j，k分别是p中点的序号，每个三元组(i，j，k)表示具有拓扑相邻的三个点。tp中的所有三角面片形成了点云p的一个外壳(crust)，crust算法可以保证这个外壳的与真实颅骨表面具有一致的拓补结构。

在稀疏采样下，点集p重建的外壳形状上相对头皮表面形状会过于“干瘪”。考虑到人类头皮表面形状应该是近似于椭球光滑曲面。因此我们引入一个球面化的过程对重建得到的外壳进行进一步完善。在球面化的过程中，我们对每一个三角面片进行多次迭代的非线性插值。经过对集合tp中的所有三角面片进行非线性插值，我们可以得到一个细分的三角网格。我们保留了三角网格上的顶点，记作致密点云dp＝{xi,yi,zi|i＝1,2,3,...,ndp}(图9中的c)，其中ndp是致密点云中的点数。对于推荐的～10量级的采样，ndp的量级大约是～10k。

颅骨表面重建得到的致密点云dp是被试者真实物理头皮表面形状的良好近似。并且，此时的致密点云dp和磁场三维定位仪记录到的4个颅骨水准点同处于一个虚拟的坐标空间，因此我们可以在这个dp上构建cpc坐标系统，并认为这个cpc坐标系统与直接在被试者物理头皮表面上得到的cpc系统近似。这样，对于磁场三维定位仪的探笔新纪录被试者物理头皮表面上的任意一点p，则经过上述坐标变换操作，也能实时进入到这个虚拟坐标空间，根据前述的cpc定位算法，我们可以实时计算出p点的cpc坐标。

在导航操作开始时，操作者使用磁场三维定位仪标定被试者头上的4个颅骨水准点nz，iz，al和ar，并以均匀的方式在被试者头皮采集多个稀疏点。然后，按照采样到的多个稀疏点重建被试者头皮表面并在重建头皮表面上构建cpc坐标。此后，探笔在被试者头皮表面的位置，会实时定位被试者自己的cpc坐标中。

由于被试者头皮表面的cpc坐标和个体经颅脑图谱导航系统中预置的3d头模型的cpc坐标是一一对应的，因此，探笔在被试者头皮表面的位置，将实时显示在个体经颅脑图谱导航系统中显示的3d头模型的头皮表面。这样，随着磁场三维定位仪的探笔在被试者头皮表面上的移动，3d头模型上的标记点也会相应移动。实验者按照个体经颅脑图谱导航系统给出的视觉反馈，通过移动探笔将3d头模型上的标记点导航到实验的感兴趣区域，此时的探笔在被试者头皮表面上物理位置就是基于经颅脑图谱的经颅设备的放置位置，如图10b所示。

综上所述，针对将经颅脑图谱应用到指导经颅设备的个体化放置中的问题，本发明进一步提供了一种个体经颅脑图谱导航系统。在一项对10-20地标点的标定实验中，我们验证了在个体经颅脑图谱导航系统指导下的真实实验中，操作者可以实现毫米级精度下对颅骨目标点的导航。最后，我们对基于经颅脑图谱的个体经颅脑图谱导航系统在真实经颅脑成像实验中的可行性进行了验证。实验结果表明，基于经颅脑图谱的导航方法及导航系统，可以有效提高近红外测量光极感兴趣脑区的覆盖准确性，以及在不同人上皮层测量位置的一致性，从而提高近红外技术检测任务诱发脑活动的敏感性。

另一方面，目前主流的功能近红外光谱(fnirs)设备采用多通道测量方式，通常由一个固定间距的光极阵列实现对脑活动的观测。这种观测实际上是一种离散的“多点观测”，不同观测点之间会相互影响，这主要是由于光极的形状限制。经颅脑图谱能够提供完整的颅-脑对应信息，便于实验者对光极放置位置进行全局规划。在一项真实的fnirs实验中，我们将基于经颅脑图谱的个体经颅脑图谱导航系统应用到多通道fnirs的光极放置中。实验结果表明，相比传统10-20方法，经颅脑图谱指导下的光极放置能够有效提高多通道fnirs测量通道对感兴趣脑区的覆盖准确性，以及测量通道在不同被试者上的观测位置的一致性，从而提高fnirs探测任务诱发脑活动的敏感性。因此，本发明所提供的经颅脑图谱导航方法及其系统对于保证光极在区域水平的覆盖准确性，甚至对于多区域覆盖(用同一台功能近红外光谱设备对多个脑区进行观测)也能体现出独特的优势。

以上对本发明所提供的基于个体特征的经颅脑图谱生成方法、导航方法及其系统进行了详细的说明。对本领域的一般技术人员而言，在不背离本发明实质精神的前提下对它所做的任何显而易见的改动，都将构成对本发明专利权的侵犯，将承担相应的法律责任。