一种SPECT三维重建图像左心室自动分割的方法技术

本发明专利技术公开了一种SPECT三维重建图像左心室自动分割的方法，通过对原始SPECT胸部三维图像线性插值进行等比例缩小，利用特征提取网络对缩小后图像提取刚性配准参数特征，利用空间变换网络及参数特征对SPECT图像进行自动转向获得标准视图的预测图像，从预测图像中切割中心部位获取心脏图像，并通过U‑NET网络进行图像自动分割获得标准视图下的左心室结构分割结果。本发明专利技术使用多任务学习的深度学习网络同步提取图像的位置特征和语义特征，并利用双网络特征的互相监督达到网络一体化训练的效果实现不同角度到标准视图的一体化自动转向、心脏定位及左心室的结构分割，减少了手动转向、分割的复杂性和人为误差，实现了图像操作的全自动并提高了准确性。

【技术实现步骤摘要】
一种SPECT三维重建图像左心室自动分割的方法
本专利技术涉及医学影像领域和深度学习领域，尤其涉及一种基于深度学习网络的SPECT三维重建图像左心室自动分割的方法。
技术介绍
SPECT心脏成像是目前临床诊断冠心病、心肌缺血等心血管疾病和疗效评价以及预后判断的金标准，其可非介入式的提供心肌组织的功能性信息来检测到尚未导致结构上变化的潜在病变，提供心肌组织更详细的功能活性信息。临床进行SPECT检查时，需要对重建后的SPECT图像进行一系列的操作和分析，其中左心室射血系数的计算是评价心脏功能的一个重要指标，其需要对左心室心室腔及心室壁进行分割，以提取不同心跳周期时的左心室腔体积进行计算。临床标准的SPECT心脏视图为短轴图SA方向，且由于左心室长轴与人体的长轴不平行，所以通常需要手动对重建SPECT图像进行旋转以获得标准SA视图，并在此视图中进行图像分割用以计算左心室射血系数。同时，基于标准SA视图的左心室图像可用于制备心脏极坐标图进行左心室心肌的活性分析。将图像从常规RA视图转至临床分析的心脏标准SA视图临床上常需要医师手动操作，此主观性的操作容易引入随机误差而影响分析准确性，且需消耗较长的手动操作时间。而对于左心室图像分割，目前常见的临床核医学心脏图像分析软件多采用常规图像处理的分割方法，例如基于左心室壁中心线的分割方法、基于左心室模型的分割方法、基于心脏图谱的分割方法、基于阈值或k均值聚类的分割方法等，由于SPECT图像分辨率较低，且心脏图像易受到呼吸和心跳运动影响出现图像运动边界模糊，目前的分割方法在进行分割时常出现分割精度低、对分割边缘提取不准确的问题从而进一步影响量化分析的精度。因此，对于SPECT心脏图像的临床处理和分析时，如何实现常规重建视图到特定的临床分析所用的标准SA视图的稳定、准确的图像自动转向和定位，如何在较低图像分辨率的SPECT心脏图像中实现对左心室结构的精准分割提取，是临床SPECT心脏处理中所面临的技术难题。
技术实现思路
本专利技术的目的在于针对现有技术的不足，提供一种基于深度学习网络的SPECT三维重建图像左心室自动分割的方法。本专利技术的目的是通过以下技术方案来实现的：一种SPECT三维重建图像左心室自动分割的方法，包括以下步骤：步骤一：对SPECT三维重建图像的常规视图RA进行缩小重采样，并将重采样后的缩小常规视图RA-r作为特征提取网络的输入，所述特征提取网络由卷积模块和全连接层组成，利用卷积模块对缩小常规视图RA-r进行特征提取，并全连接展开后形成6维特征向量T-r，所述6维特征向量T-r包含3个方向的平移参量和3个角度的旋转参量；步骤二：利用常规视图RA与缩小常规视图RA-r的等比例关系将T-r调整为特征向量T，其中旋转参量不变而平移参量等比例放大；步骤三：利用空间变换网络将特征向量T-r应用于缩小常规视图RA-r中得到预测缩小图像SA-r’，同时将特征向量T应用于常规视图RA中得到预测图像SA’；步骤四：以预测图像SA’的中心为中心，截取图像中的提取心脏部分，形成预测心脏图像SA-H’，同时，对心脏图像SA-H’进行图像梯度计算获得对应梯度图SA-G；步骤五：将心脏图像SA-H’和梯度图SA-G融合为双通道图像并经过三维U-NET网络的下采样及上采样提取图像特征并经过softmax层进行分割处理即获得预测的左心室结构分割结果F；其中，所述特征提取网络、空间变换网络和U-NET网络联合采用多任务共同学习训练，训练的总损失函数L=δL-par+μL-img+λL-seg；其中δ、μ和λ为权重系数；L-img为预测缩小图像SA-r’与缩小标准视图SA-r之间的图像损失函数、L-par为特征向量T-r与刚性配准参数P-r之间的参数损失函数，L-seg为预测分割结果F与SA方向分割标签G之间的标签损失函数；所述缩小标准视图SA-r通过手动转向SPECT三维重建图像的常规视图RA并按等比例缩小获得；所述刚性配准参数P-r为刚性配准算法计算的缩小常规视图RA-r与缩小标准视图SA-r之间的配准参数，包含3个方向的平移参量和3个角度的旋转参量；所述SA方向分割标签G包含左心室腔、心室壁及背景的3个数值，通过以标准视图SA中心为中心截取的心脏图像SA-H并手动勾画心脏左心室的心室腔及心室壁获得。进一步地，所述缩小常规视图RA-r、缩小标准视图SA-r的尺寸优选为64*64*64体素大小；心脏图像SA-H和预测心脏图像SA-H’应以涵盖整个心脏图像，同时尽量少的包含其他高强度器官信息为最佳，优选为32*32*32体素大小。进一步地，所述步骤三中，通过构建变换矩阵P=[RMT]获取预测图像，其中，M=[tx、ty、tz]表示位移矩阵，tx、ty、tz分别为特征向量中对应的3个方向的平移参量，R为旋转矩阵，将R从欧拉角R=[β、α、γ]转换为世界坐标系参数：。进一步地，所述梯度图SA-G获得的计算公式为：其中，gx,gy,gz分别为x、y、z方向的梯度，其计算公式为：其中i、j、k为图像内x、y、z三个方向的坐标索引。由于梯度图像可以更好的提供图像边界信息，同步读入梯度图可以优化网络对图像边缘的识别从而进一步提高网络对图像分割的效力。进一步地，所述特征提取网络的卷积模块、全连接层的个数均为3个，每个卷积模块包含卷积层和池化层；所述三维U-NET网络的下采样及上采样均为4次，其中下采样包含卷积层和池化层，上采样包含卷积层和反卷积层。进一步地，所述SA方向分割标签G中，将背景设置为0,左心室腔设置为1，左心室壁设置为2。进一步地，所述δ、μ和λ分别取值为1,100和10。进一步地，所述图像损失函数L-img采用均方差函数。进一步地，所述参数损失函数L-par采用绝对值损失函数L1或范数损失函数L2。进一步地，所述标签损失函数L-seg采用Dice-loss损失函数，其中，分别对左心室腔及左心室壁计算标签损失函数L-seg-1及L-seg-2，最终L-seg=L-seg-2+L-seg-1。本专利技术的有益效果是：本专利技术使用多任务学习的深度学习网络同步提取图像的位置特征和语义特征，并利用双网络特征的互相监督达到网络一体化训练的效果以实现不同角度到标准视图的一体化自动转向、心脏定位及左心室的结构分割，减少了手动转向、分割的复杂性和人为误差，实现了图像操作的全自动并提高了准确性，并且本专利技术还进一步采用梯度图像增加对图像数据边缘分割的准确性。附图说明图1是本专利技术SPECT三维图像左心室自动分割的流程示意图。图2是SPECT三维图像左心室自动分割模型结构中的自动转向和定位模块示意图。图3是SPECT三维图像左心室自动分割模型结构中的自动分割模块示意图。具体实施方式下面结合附图详细说明本专利技术。本专利技术提出的一种基于深度学习网络的SPECT三维本文档来自技高网...

【技术保护点】
1.一种SPECT三维重建图像左心室自动分割的方法，其特征在于，包括以下步骤：/n步骤一：对SPECT三维重建图像的常规视图RA进行缩小重采样，并将重采样后的缩小常规视图RA-r作为特征提取网络的输入，所述特征提取网络由卷积模块和全连接层组成，利用卷积模块对缩小常规视图RA-r进行特征提取，并全连接展开后形成6维特征向量T-r，所述6维特征向量T-r包含3个方向的平移参量和3个角度的旋转参量；/n步骤二：利用常规视图RA与缩小常规视图RA-r的等比例关系将T-r调整为特征向量T，其中旋转参量不变而平移参量等比例放大；/n步骤三：利用空间变换网络将特征向量T-r应用于缩小常规视图RA-r中得到预测缩小图像SA-r’，同时将特征向量T应用于常规视图RA中得到预测图像SA’；/n步骤四：以预测图像SA’的中心为中心，截取图像中的心脏部分，形成心脏图像SA-H’，同时，对心脏图像SA-H’进行图像梯度计算获得对应梯度图SA-G；/n步骤五：将心脏图像SA-H’和梯度图SA-G融合为双通道图像并经过三维U-NET网络的下采样及上采样提取图像特征并经过softmax层进行分割处理即获得预测的左心室结构分割结果F；/n其中，所述特征提取网络、空间变换网络和三维U-NET网络联合采用多任务共同学习训练，训练的总损失函数L = δL-par + μL-img + λL-seg；/n其中δ、μ和λ为权重系数；L-img为预测缩小图像SA-r’与缩小标准视图SA-r之间的图像损失函数、L-par为特征向量T-r与刚性配准参数P-r之间的参数损失函数，L-seg为预测分割结果F与SA方向分割标签G之间的标签损失函数；所述缩小标准视图SA-r通过手动转向SPECT三维重建图像的常规视图RA并按等比例缩小获得；所述刚性配准参数P-r为刚性配准算法计算的缩小常规视图RA-r与缩小标准视图SA-r之间的配准参数，包含3个方向的平移参量和3个角度的旋转参量；所述SA方向分割标签G包含左心室腔、心室壁及背景的3个数值，通过以标准视图SA中心为中心截取的心脏图像SA-H并手动勾画心脏左心室的心室腔及心室壁获得。/n...

【技术特征摘要】
1.一种SPECT三维重建图像左心室自动分割的方法，其特征在于，包括以下步骤：
步骤一：对SPECT三维重建图像的常规视图RA进行缩小重采样，并将重采样后的缩小常规视图RA-r作为特征提取网络的输入，所述特征提取网络由卷积模块和全连接层组成，利用卷积模块对缩小常规视图RA-r进行特征提取，并全连接展开后形成6维特征向量T-r，所述6维特征向量T-r包含3个方向的平移参量和3个角度的旋转参量；
步骤二：利用常规视图RA与缩小常规视图RA-r的等比例关系将T-r调整为特征向量T，其中旋转参量不变而平移参量等比例放大；
步骤三：利用空间变换网络将特征向量T-r应用于缩小常规视图RA-r中得到预测缩小图像SA-r’，同时将特征向量T应用于常规视图RA中得到预测图像SA’；
步骤四：以预测图像SA’的中心为中心，截取图像中的心脏部分，形成心脏图像SA-H’，同时，对心脏图像SA-H’进行图像梯度计算获得对应梯度图SA-G；
步骤五：将心脏图像SA-H’和梯度图SA-G融合为双通道图像并经过三维U-NET网络的下采样及上采样提取图像特征并经过softmax层进行分割处理即获得预测的左心室结构分割结果F；
其中，所述特征提取网络、空间变换网络和三维U-NET网络联合采用多任务共同学习训练，训练的总损失函数L=δL-par+μL-img+λL-seg；
其中δ、μ和λ为权重系数；L-img为预测缩小图像SA-r’与缩小标准视图SA-r之间的图像损失函数、L-par为特征向量T-r与刚性配准参数P-r之间的参数损失函数，L-seg为预测分割结果F与SA方向分割标签G之间的标签损失函数；所述缩小标准视图SA-r通过手动转向SPECT三维重建图像的常规视图RA并按等比例缩小获得；所述刚性配准参数P-r为刚性配准算法计算的缩小常规视图RA-r与缩小标准视图SA-r之间的配准参数，包含3个方向的平移参量和3个角度的旋转参量；所述SA方向分割标签G包含左心室腔、心室壁及背景的3个数值，通过以标准视图SA中心为中心截取的心脏图像SA-H并手动勾画心脏左心室的心室腔及心室壁获得。


2.根据权利要求1所述SPECT三维重建图像左心室自动分割的方法，其特征在于，所述缩小常规视图RA-r、缩小标准视图SA-r的尺寸为64*64*64体素大小；心脏图像SA-H和预测...

【专利技术属性】
技术研发人员：张铎，朱闻韬，韩璐，黄海亮，祁二钊，
申请(专利权)人：之江实验室，
类型：发明
国别省市：浙江;33