[发明专利]基于特征对齐域的SP-CTA图像冠脉分割方法及装置

权利要求书

1.基于特征对齐域的SP-CTA图像冠脉分割方法，其特征在于：其包括以下步骤：

(1)图像预处理：对源域图像及其标注数据进行缩放，使同一病人的CTA图像和SP-CTA图像具有相同的大小；

(2)在特征提取阶段，采用特征共享的方式使SP-CTA图像能够有效利用CTA图像空间学习的知识，进而加速网络的收敛；在卷积网络中使用特定于域的BN层；

(3)在训练阶段，CTA图像与SP-CTA图像分别以batch-size的大小加载并同时输入到冠脉分割网络中，网络利用共享卷积核及不同的BN层，得到两个域输入图像的血管预测概率图；

(4)在血管预测空间引入对抗性损失函数，网络中加入两个判别器，分别对分割网络输出的预测概率图与预测熵图进行对应模态的判别；

(5)分割网络采用三次下采样和三次上采样；卷积块由两个卷积核大小为3×3×3的卷积层，BN层以及LeakyReLU激活函数构成，卷积层的通道数分别设置为16，32，64，128，提取特征后，利用最大池化进行下采样；上采样阶段采用3D反卷积操作，通过学习相关参数恢复与原始图像大小相同的特征图；为了恢复下采样丢失的部分信息，在解码阶段通过跳连接融合低维特征；在最后的预测阶段，采用两个不同的卷积操作以预测不同域的血管概率。

2.根据权利要求1所述的基于特征对齐域的SP-CTA图像冠脉分割方法，其特征在于：所述步骤(1)中，通过调整图像z轴方向的spacing对图像进行插值操作实现重采样，插值操作首先根据目标图像和原图像的大小求出两张图像像素点的对应关系，坐标变换计算公式为：

其中，(X_src，Y_src，Z_src)为原图像中某个像素点m的空间坐标，(X_dst，Y_dst，Z_dst)为插值后的图像中与原图像点m对应像素点的坐标，W_src，H_src，L_src为原图像的宽度，高度和长度，W_dst，H_dst，L_dst为变换后图像的宽度，高度和长度；

分割金标准使用最邻近插值采样方式，m(x，y)为原图像中第x行第y列的像素，假设重采样后的图像中的某一像素点通过公式(1)映射到原图的P点，由于P点的坐标值为小数，不同于最邻近插值法的四舍五入取整，通过计算P点与其邻近的四个已知像素值的点之间的距离，然后利用四个点的灰度值以距离为权重的方式叠加计算出P点的像素值，计算公式表示为：

P＝u*v*Src(x，y)+u*(1-v)*Src(x+1，y)+(1-u)*v*Src(x，y+1)+(1-u)*(1-v)*Src(x+1，y+1) (2)。

3.根据权利要求2所述的基于特征对齐域的SP-CTA图像冠脉分割方法，其特征在于：所述步骤(1)中，源域图像经过重采样后，对两个域的数据进行标准化操作，采用3D的卷积神经网络，采用沿着z轴剪切的方式，实验中切片的个数为8。

4.根据权利要求3所述的基于特征对齐域的SP-CTA图像冠脉分割方法，其特征在于：所述步骤(3)中，假设在某一特定的层具有k个通道，其中某个特定通道的一个mini-batch的数据表示为B＝{x₁，x₂，x₃，...，x_m}，首先通过公式(3)与(4)分别计算此通道数据的平均值和方差，然后根据平均值和方差并利用式(5)对所有数据进行标准化，最后对变换后的数据进行平移和缩放处理以恢复网络初始所要学习的特征分布，计算公式为式(6)所示

采用全部训练数据的均值和方差，均值和方差计算公式如(7)和(8)所示，在训练阶段将利用移动平均法来获取全局统计量

E[x]＝E_B[μ_B] (7)

在测试阶段经过BN层处理后的数据分布如式(9)所示