[发明专利]三维大脑磁共振图像的大脑皮层表面上脑沟盆地分割方法

权利要求书

1.三维大脑磁共振图像的大脑皮层表面上脑沟盆地分割方法，其特征在于：

步骤1对三维大脑核磁共振图像进行预处理：利用可变性模型方法去除脑壳，利用配准方法去除非大脑组织，利用高斯混合模型方法对大脑图像进行组织分割，得到白质，灰质和脑脊髓液三种组织类型表示的图像；

步骤2大脑皮层表面重建：利用Marching Cubes方法从组织分割后的大脑图像中重构三角化的大脑皮层表面；

步骤3利用有限差分方法估计三角化的大脑皮层表面上每个顶点的最大主曲率和最大主方向；

步骤4：在满足v(x)·n(x)＝0的条件下，最小化能量函数

得到最大主方向流场v(x)＝(u(x)，v(x)，w(x))：其中λ是权重参数为0.05～0.15，是梯度算子，p(x)是顶点x的最大主方向，f(x)是顶点x的最大主曲率强度的递减函数，f(x)＝|c(x)|，其中c(x)是顶点x的最大主曲率；n(x)是顶点x的法向量，S表示大脑皮层表面上所有顶点的集合；所述的v(x)·n(x)＝0限定最大主方向流场在大脑皮层表面上的切平面内；

步骤5利用大脑皮层表面上最大主方向流场跟踪方法，得到脑沟盆地分割：具体步骤如下：在最大主方向流场中，计算大脑皮层表面上任意一个顶点x所流经的下一个顶点其中x_i是顶点x的第一圈邻域；如果两个连续流向量的夹角小于π/2，利用以上公式计算下一个流经的顶点，流场跟踪步骤继续执行；否则，流场跟踪步骤就停止，到达了一个脑沟底部区域，将流场跟踪停止顶点标记为脑沟底部顶点，同时将该流场跟踪停止顶点记录为流路径所经过的所有顶点的对应脑沟底部顶点；利用连通成分分析方法，将所有同一连通的脑沟底部区域标记为同一颜色，将所有对应于同一个脑沟底部顶点的顶点标记为与其对应脑沟底部顶点相同的颜色，得到大脑皮层表面上脑沟盆地的分割。

2.根据权利要求1所述的三维大脑磁共振图像的大脑皮层表面上脑沟盆地分割方法，其特征在于：所述最小化能量函数的计算通过将v看作是时间的方程，计算得到最大主方向流场，具体步骤如下：

步骤1：根据v(x，t)(t＝0，1，2...)，将方程

得到弥散的最大主方向场v(x，t+1)，上述公式中，v_t(x，t)表示v(x，t)对于时间t的偏导数，拉普拉斯项其中x_i是顶点x周围一圈的邻域，是一个权重系数计算为顶点x和x_i之间距离的倒数；同理估计和

步骤2：投影弥散的最大主方向场v(x，t+1)到切平面内；

步骤3：归一化投影的弥散的最大主方向场v(x，t+1)；

步骤4：将v(x，t)替换为投影的弥散的归一化的v(x，t+1)，重复步骤1至3共30～50次迭代停止。