[发明专利]基于图像子块支持向量机的脑血流信号计算方法

权利要求书

1.基于图像子块支持向量机的脑血流信号计算方法，其特征在于包括如下步骤：

步骤1、首先采集得到n个控制图像和标记图像对，记为{C_i,L_i},i＝1,2,....,n，其中C_i表示第i个控制图像，L_i表示与C_i相对应的标记图像；

步骤2、根据上一步得到的n个控制图像和标记图像对，对图像中的每一个体素逐一计算脑血流值，从而得到最终的脑血流图像，具体包括如下子步骤：

步骤1)、对于待计算的目标体素，分别在每个控制图像和标记图像中取出以该目标体素位置为中心的一个图像子块，从控制图像中取出的图像子块的类别标记为+1，从标记图像中取出的图像子块的类别标记为-1，接着将取出的图像子块拉成一个列向量p_k作为特征向量，对应的类别记为l_k(l_k∈{-1,+1})，得到一个训练数据集{p_k,l_k},k＝1,2,...,2n，利用该训练数据集计算对应的脑血流；

步骤2)、对每个特征向量p_k进行标准化得到x_k，使得p_k均值为0，标准差为1：

步骤3)、利用非线性支持向量机对数据集{x_k,l_k}进行分类：

s.t.l_k(w·φ(x_k)+b)≥1-ξ_k,k＝1,2,...,2n (2)

ξ_k≥0,k＝1,2,...,2n (3)

其中，C是一个平衡参数，φ(·)是一个非线性变换，w,b分别是超平面的法向量和截距，ξ_k是一个松弛变量；

步骤4)、优化公式(1)—(3)的对偶问题：

0≤α_i≤C，i＝1,2,...,2n (6)

其中，α_k是拉格朗日乘子，K(x_i,x_j)是正定核函数，采用二次多项式核，求解(4)-(6)得到拉格朗日乘子α_k，k＝1,2,...,2n的值，于是得到图像子块对应的脑血流信号

2.如权利要求1所述的基于图像子块支持向量机的脑血流信号计算方法，其特征在于：对每个体素计算脑血流值时，相邻的体素计算出来的图像子块脑血流会有重合部分，对于重合部分取它们的平均值作为最终的脑血流值。

3.如权利要求1所述的基于图像子块支持向量机的脑血流信号计算方法，其特征在于：步骤1所述的n≥40。

4.如权利要求1所述的基于图像子块支持向量机的脑血流信号计算方法，其特征在于：所述图像子块为三维立方体结构，图像子块的半径r的范围为2≤r≤5。