一种立体脑电图电极定位方法

摘要

本发明提出了一种立体脑电图电极定位方法，采用CT影像进行定位：S1：获取电极植入后的CT影像进行配准；S2：将配准得到的多层图像进行二值化处理；S3：将包含目标电极的所有图像叠加为一张图片；S4：根据步骤S3得出的图片确定所述目标电极的起始坐标与终点坐标，计算电极位点坐标；S5：将步骤S4得出的电极位点坐标进行空间坐标变换，根据变换后的坐标定位脑区位置。定位过程仅用到手术后的CT影像即可工作，无需用到病人的多组影像数据；而且，定位的进行不受得到的CT影像分辨率影响，且计算量小，工作效率高。进一步，计算步骤的应用能够适用于不同系统上的开发场景，可移植性强，定位方法实用且有效。

说明书

技术领域

本发明涉及医学图像处理技术领域，尤其涉及一种通过CT影像数据进行的立体脑电图电极定位方法。

背景技术

立体定向脑电图(SEEG)是一种微创外科手术，用于识别大脑中癫痫发作的起始区域。医生在手术时将含有多个位点的电极放置在目标脑区，然后对电极采集到的脑电信号进行分析，以精确定位癫痫疾病的致痫区域。立体定向脑电图可以做到精确可靠、有针对性的治疗，并能够发现常规脑电图(EEG)测试无法发现的大脑深处的癫痫发作部位。但是，在分辨率不高的MRI和CT影像中电极和电极上每个位点的位置难以定位，对后期对影像的处理和分析上有较大的影响。所以，如何有效定位电极位点所在脑区的位置，是一项亟待解决的问题。

现有技术中采用先以计算机进行识别初步筛选电极位点，然后通过人工确认及筛选得到待选电极信号的定位方式，确保对电极位点的定位准确性。这种方式相比于传统的采用人工标注的方式提升了定位准确性，但是，定位过程中需要运用到病人的多种影像数据、对影像数据的分辨率要求高，且需要对电极图像信号和体素进行处理；而且，由于不同人体的结构存在微小差异，定位过程中需要进行大量的运算，且针对不同的情况需要重新设置大量的参数，过程复杂。

发明内容

有鉴于此，为了克服上述现有技术的缺陷，本发明提出了一种能够在分辨率不高的影像中进行立体脑电图电极定位的方法。

具体地，所述立体脑电图电极定位方法包括：

S1：获取电极植入后的CT影像进行配准；

S2：将配准得到的多层图像进行二值化处理；

S3：将包含目标电极的所有图像叠加为一张图片；

S4：根据步骤S3得出的图片确定所述目标电极的起始坐标与终点坐标，计算电极位点坐标；

S5：将步骤S4得出的电极位点坐标进行空间坐标变换，根据变换后的坐标定位脑区位置。优选地，使用标准模板对所述CT影像进行配准。

进一步，步骤S2还包括，对二值化处理得到的图像进行去噪处理，提升所得到的图像质量。

步骤S4的具体执行过程包括：

S41：确定步骤S3得到的图片中，所述目标电极的起始点与终点的二维坐标；

S42：根据所述起始点与所述终的点二维坐标及所述目标电极的起始点和终点的所在的图像编号，得出所述起始点与所述终点的三维坐标；

S43：根据所述起始点与所述终点的三维坐标计算所述电极位点的坐标。

步骤S43的具体算法为：

根据起始坐标(X

signX＝sign(X

oignY-oign(Y

signZ＝sign(Z

具体地，所述电极中位点数量N为包含所述目标电极的冠状面图像数量。进一步，对步骤S43计算得出的三维坐标的数值进行取整。

步骤S5的具体执行过程包括：

S51：根据公式(X,Y,Z)＝T*(Xk,Yk,Zk)+b将步骤S4得到的空间坐标变换以转换为标准模板中的坐标，其中，T为线性映射矩阵，b为平移表示需要的向量；

S52：将配准后的CT影像与标准模板叠加；

S53：根据变换后的坐标定位脑区位置。

所述二值化的具体过程包括：设定阈值，将步骤S1所获得的图像像素与所述阈值比较，将将大于所述阈值的部分的像素调整为255，小于所述阈值的部分的像素调整为0。

基于上述立体脑电图电极定位方法，本发明还提供了一种立体脑电图电极定位系统，包括：

配准模块，用于将所获取的电极植入后的CT影像进行配准；

二值化模块，用于将配准得到的多层图像进行二值化处理；

图像叠加模块，用于将包含所述目标电极的所有图像叠加为一张图片；

坐标计算模块，用于根据步骤S3得出的图片确定所述目标电极的起始点坐标与终点坐标，计算电极位点坐标；

坐标变换模块，用于将步骤S4得出的电极位点坐标进行空间坐标变换，根据变换后的坐标定位脑区位置。

综上所述，本发明的立体脑电图定位方法只用到手术后CT影像即可工作，不需要用到病人的多种影像数据，且不需要设置大量的参数，泛化能力强；在CT影像分辨率较低的情况下依然可以使用，不受得到的CT影像分辨率的影响，提高实验的工作效率。进一步，计算步骤的应用能够适用于不同系统上的开发场景，可移植性强，定位方法实用且有效。

附图说明

为了更清楚地说明本申请实施例中的技术方案，下面将对实施例描述中所需要使用的附图作简单地介绍，显而易见地，下面描述中的附图仅仅是本申请的一些实施例，对于本领域普通技术人员来讲，在不付出创造性劳动的前提下，还可以根据这些附图获得其他的附图。

图1为本发明立体脑电图电极定位方法的工作流程示意图；

图2为二值化后某一电极CT冠状面图像；

图3为包含图2所示电极的所有CT冠状面图像的拼接图

图4为实施例2中的目标电极的部分冠状面图像；

图5为本发明立体脑电图电极定位系统的结构示意图。

具体实施方式

下面将结合本发明实施例中的附图，对本发明实施例中的技术方案进行清楚、完整地描述，显然，所描述的实施例仅是本发明一部分实施例，而不是全部的实施例。基于本发明中的实施例，本领域普通技术人员在没有做出创造性劳动前提下所获得的所有其他实施例，都属于本发明保护的范围。

本发明提供了一种立体脑电图电极定位方法，采用参考影像将电极植入后获得的CT影像进行配准；对配准得到的每一层图像都进行二值化处理；确定需要定位的目标电极跨越的每层图像，对其进行拼接；根据目标电极的起始点坐标和终点坐标计算得到CT影像中的电极位点的坐标；将CT影像中的电极位点坐标进行空间坐标变换转换为标准模板中的坐标，根据标准模板中的坐标定位找到电极对应的脑区位置。

实施例1

参见说明书附图1，本发明的立体脑电图电极定位方法的具体工作流程如下：

S1：将电极植入后获得的CT影像配准到参考影像。其中，参考影像采用标准模板，配准过程可以使用开源工具fs l或spm。在本实施例中，采用fs l工具中的线性配准功能将病人手术后的CT影像配准到CT的模板上。具体地，进行配准时对部分参数的选择为：-in:需要配准的影像，-ref:参考影像，-omat:配准过程保存的仿射矩阵，-out:配准结果；其中，仿射变换作用为使得配准过程中目标影像与参考影像在体素层面对齐。其他如代价函数加权(cost function weighting)使得不同体素有不同权重，自由度(degrees of freedom)，插值方法，成本函数等参数选用默认设置。

在一些实施例中，使用MNI152模板作为标准模板进行配准操作。上述CT模板和MNI152模板具有相同结构。可选地，还可以采用ICBM152模板作为标准模板。

S2：配准之后的CT影像形成多层的图像，先将每一层图像进行二值化处理。二值化的具体过程包括：设定阈值，将步骤S1所获得的图像像素与阈值比较，将大于阈值的部分的像素调整为255，小于阈值的部分的像素调整为0。二值化之后的图像中颅骨和电极的像素值为255，显示为白色，背景像素值为0，显示为黑色。在本实施例中，采用0TSU二值化方法处理。然后对二值化之后的图像进行去噪，提升图像质量。具体地，可以使用数字图像处理中的形态学操作开操作。步骤S2执行完成得到的图像如图2所示。

S3：确定目标电极跨越的每层图像，将电极的起始和终点位置之间所有的二值化图像叠加为一张图片，获得可以清晰显示此电极形状的图片。可选地，上述图像叠加可以采用0penCV中的add()，或是其他图像叠加方法进行，例如，resize()和addWeighted()。拼接得到的图片如图3所示，白色电极清晰可见。

S4：将步骤S3得到的拼接后的图片视为电极从三维位置向二维的投影，在拼接后图片中计算出电极的起始点和终点的二维坐标，则根据投影还原出原始三维的坐标，即为电极在配准后的CT影像中的坐标位置。例如，计算名称为A，电极位点个数为N的电极位点坐标的过程如下：将电极上所有的位点标记为A

则运用以下公式计算所有电极位点的X方向坐标：

oignX-oign(x

同理，电极位点的Y方向坐标的计算过程为：

signY＝sign(Y

电极位点的Z方向坐标的计算过程为：

signZ＝sign(Z

S5：得到CT影像中的电极位点坐标之后，由于CT/MNI体素位置坐标(voxellocation)和MNI空间坐标(MNI coordinates)的偏差方向不同，需要将坐标进行空间坐标变换：(X，Y，Z)＝T*(Xk，Yk，Zk)+b，其中，T为线性映射矩阵，b为平移表示需要的向量。在具体使用公式时，需要根据实际情况选择不同的变换矩阵。之后，将配准后的CT影像与MNI152标准模板叠加，根据变换后的坐标定位到脑区所在位置即可得到此位点的脑区位置。重复此流程即可得到所有电极位点的脑区坐标。

实施例2

本实施例提供了对一个目标电极进行定位的具体计算过程。

具体地，本实施例对一个编号为86-99的目标电极进行定位，上述电极具有共14个冠状面图像，参见说明书附图4，为从中提取的7张图像，可以看到电极在不同图像中的变化。由于CT影像是由若干张图像形成的影像集，每张图片只能反应出电极的一部分，所以需要将图像叠加以展示出完整电极的形状。

电极跨越编号86-99，可得到MNI体素空间位置Z坐标为86-99，则Z

signX＝sign(X

同理可得：

然后，进行空间仿射变换

如将Talairach坐标转换为MNI坐标，作如下变换：

如将体素位置坐标转换为MNI坐标，由于体素坐标只能为正整数和0，需计算round(X

则得到

其中，X，Y，Z为电极位点在MNI体素位置空间坐标，X′，Y′，Z′为得到的MNI空间坐标。

实施例3

本实施例为一种运用实施例1所提供的立体脑电图电极定位方法的立体脑电图电极定位系统。

参见说明书附图5，立体脑电图电极定位系统包括：

配准模块，用于将所获取的电极植入后的CT影像进行配准；

二值化模块，用于将配准得到的多层图像进行二值化处理，随后对二值化之后的图像进行去噪；

图像叠加模块，用于将包含目标电极的所有图像叠加为一张图片；

坐标计算模块，用于根据步骤S3得出的图片确定目标电极的起始坐标与终点坐标，计算电极位点坐标；

坐标变换模块，用于将步骤S4得出的电极位点坐标进行空间坐标变换，根据变换后的坐标定位脑区位置。

综上所述，本发明提供了一种能够在分辨率不高的影像中进行立体脑电图电极定位的方法，只用到手术后CT影像即可工作，不需要用到病人的多种影像数据，且不需要设置大量的参数，泛化能力强；不受得到的CT影像分辨率的影响，在CT影像分辨率较低的情况下依然可以使用，提高实验的工作效率；可以根据计算得到的电极位点的坐标找到位点所在的脑区位置，是之后在脑区层面的处理分析的前提条件；而且，本方法除配准步骤外其他计算步骤在所有平台都可以使用，可移植性强，适应不同系统上的开发场景，因而更加实用有效。

以上所述仅为本发明的较佳实施例，并不用以限制本发明，除了以上实施例以外，还可以具有不同的变形例，以上实施例的技术特征可以相互组合，凡在本发明的精神和原则之内，所作的任何修改、等同替换、改进等，均应包含在本发明的保护范围之内。