一种基于三维影像的脊柱分割方法技术

本发明专利技术公开了一种基于三维影像的脊柱分割方法，包括步骤：获取包含脊柱的正侧位数字重建放射影像；对所述正侧位数字重建放射影像进行分割得到各节脊柱区域图像与非脊柱区域图像；获取所述脊柱区域图像中各节脊柱的最小外接矩形，并根据其中心点得到三维影像中各节脊柱的空间关键点；以所述各节脊柱的空间关键点作为中心，以设定柱状体形式对三维影像沿脊柱高度方向进行分段切割得到各节脊柱分段切片图像；将各节脊柱分段切片图像进行拼接得到三维影像中脊柱的分割结果。本发明专利技术结合2D数字重建放射影像为输入，最终结果准确度更高，且预测耗时大大减小，在相同硬件配置上，算法推理速度更快。理速度更快。理速度更快。

【技术实现步骤摘要】
一种基于三维影像的脊柱分割方法


[0001]本专利技术涉及图像处理
，尤其涉及一种基于三维影像的脊柱分割方法。

技术介绍

[0002]目前常见医疗影像分割主要采用传统算法来完成，例如阈值法、区域增长法等。传统算法对特征设计要求比较高，算法建模难度比较大，且算法模型鲁棒性差，算法耗时也较长。
[0003]在医疗脊柱手术机器人上，为了准确完成对脊柱术中穿刺导航，要求准确分割出各节脊柱所在区域，而且要求每节脊柱具备不同的实例信息，因此需要对人体脊柱进行3D分割，但是存在以下困难：
[0004](1)医疗影像较复杂，各节脊柱间边界不清晰，各节脊柱之间相对相似，对检测所受干扰性比较大；
[0005](2)各节脊柱3D分割时要求边界严格分开，不能存在重叠现象，对分割精度要求较高；
[0006](3)要求各节脊柱满足准确分类信息，由于各节脊柱相似，对分类难度较大。

技术实现思路

[0007]专利技术目的：本专利技术针对上述不足，提出一种基于三维影像的脊柱分割方法，可以大大提高分割效率和分割准确度。
[0008]技术方案：
[0009]一种基于三维影像的脊柱分割方法，包括步骤：
[0010]获取包含脊柱的正侧位数字重建放射影像；
[0011]对所述正侧位数字重建放射影像进行分割得到各节脊柱区域图像与非脊柱区域图像；
[0012]获取所述各节脊柱区域图像中脊柱的最小外接矩形，并根据其中心点得到三维影像中各节脊柱的空间关键点；
[0013]以所述各节脊柱的空间关键点作为中心，以设定柱状体形式对三维影像沿脊柱高度方向进行分段切割得到各节脊柱分段切片图像；
[0014]将各节脊柱分段切片图像进行拼接得到三维影像中脊柱的分割结果。
[0015]所述获取所述各节脊柱区域图像中脊柱的最小外接矩形具体为：
[0016]根据分割得到正侧位数字重建放射影像中所述各节脊柱区域图像的各像素点的坐标，以该节脊柱的像素点在对应影像的二维坐标系的两个坐标轴中的最大坐标差值作为最小外接矩形的长与宽，据此得到各节脊柱区域图像中脊柱的最小外接矩形。
[0017]根据所述最小外接矩形的中心点基于最小二乘法计算得到三维影像中各节脊柱的空间关键点。
[0018]根据分割得到正侧位数字重建放射影像中所述各节脊柱区域图像的各像素点的
坐标包括以下步骤：
[0019]根据各节脊柱区域图像与非脊柱区域图像的像素级别进行像素分割，并对像素分割得到的各节脊柱区域图像作二值化掩膜处理，得到各节脊柱区域图像的像素矩阵，得到各节脊柱区域图像的各像素点的坐标。
[0020]对所述正侧位数字重建放射影像进行分割采用的是2D神经网络模型，通过若干2D脊柱图像进行训练得到。
[0021]所述设定柱状体为长方体，所述长方体的长、宽、高方向分别对应人体左右、前后和高度方向，且均大于各节脊柱中在对应方向的最大尺寸。
[0022]在对三维影像沿脊柱高度方向进行分段切割后，对得到各节脊柱分段切片图像进行二值化处理得到单节脊柱切片图像的3D掩膜，以得到的单节脊柱切片图像的3D掩膜进行拼接得到三维影像中脊柱的分割结果。
[0023]对三维影像分割得到的各节脊柱分段切片图像采用3D分割模型，所述3D分割模型采用3D脊柱图像中标注各节脊柱作为样本进行训练。
[0024]所述获取包含脊柱的正侧位数字重建放射影像采用如下方式：
[0025]通过输入的包含脊柱的三维影像生成正侧位数字重建放射影像，或者直接输入包含脊柱的正侧位数字重建放射影像。
[0026]有益效果：
[0027](1)本专利技术采用多模态输入，包括2D数字重建放射影像(DRR影像)与3D影像的切片图像，使得最终结果准确度更高。其中，以2D数字重建放射影像为输入，通过检测各节脊柱正侧位数字影像重建影像的中心点以及各节脊柱的外轮廓，具备较高定位与分割准确度以及较强鲁棒性；以3D影像的切片图像为输入，通过对3D影像的切片图像上各节脊柱进行语义分割，具备较高定位与分割准确度以及较强鲁棒性；
[0028](2)本专利技术结合2D数字重建放射影像为输入，预测耗时大大减小，在相同硬件配置上，算法推理速度更快。
附图说明
[0029]图1为本专利技术的流程图；
[0030]图2为三维影像示意图；
[0031]图3为正位数字重建放射影像示意图；
[0032]图4为侧位数字重建放射影像示意图；
[0033]图5为对正侧位数字重建放射影像进行分割的示意图；
[0034]图6为得到各节脊柱空间关键点的示意图；
[0035]图7为三维影像各切片图像的示意图；
[0036]图8为三维影像各切片图像生成掩膜的示意图；
[0037]图9为最终得到的脊柱三维分割结果图。
具体实施方式
[0038]下面结合附图和具体实施例，进一步阐明本专利技术。
[0039]本专利技术基于三维影像的脊柱分割方法如图1所示，包括步骤：
[0040](1)获取包含脊柱的正侧位数字重建放射影像；
[0041]本专利技术中，可以通过输入的包含脊柱的三维影像生成正侧位数字重建放射影像，输入的三维影像为原始临床包含脊柱的CT切片dicom数据，如图2所示，那么根据CT切片dicom数据特征可以模拟DRR过程自动生成正侧位CT DRR影像；
[0042]在本专利技术中，还可以直接输入经过前述得到的正侧位数字重建放射影像；
[0043](2)对步骤(1)得到的正侧位数字重建放射影像进行分割得到各节脊柱区域图像与非脊柱区域图像；
[0044]本专利技术中，通过预先训练得到的2D分割模型对正侧位数字重建放射影像进行分割，预先训练得到的2D分割模型采用神经网络模型，通过若干2D脊柱图像中标注各节脊柱作为样本进行训练得到；采用的神经网络可以选择2D卷积神经网络；本专利技术中，2D分割模型可以很大程度上提高关键点检测准确性与预测速度，增加算法落地可行性；
[0045](3)利用softmax函数对各节脊柱区域图像与非脊柱区域图像的像素级别做二分类以进行各节脊柱区域图像与非脊柱区域图像的像素分割，并对分割得到的各节脊柱区域作二值化掩膜处理，得到各节脊柱区域的像素矩阵；
[0046](4)根据步骤(3)得到的各节脊柱区域的像素矩阵计算得到在正侧位数字重建放射影像中各节脊柱区域中各像素点的坐标，以其中各节脊柱区域的像素点在对应数字重建放射影像的二维坐标系的两个坐标轴中的最大坐标差值作为该节脊柱的最小外接矩形的长与宽；以该最小外接矩形对患者正侧数字重建放射影像中该节脊柱区域进行框选；如图3、4所示；
[0047](5)以各节脊柱区域的最小外接矩形的中心点作为正侧位数字重建放射影像中对应的各节脊柱区域的中心点，如图5所示，例如某一节脊柱的中心点A点所示，分别可以检测到其正位的中心点A1(x,y)与侧位的中心点A2(z,y)；...

【技术保护点】

【技术特征摘要】
1.一种基于三维影像的脊柱分割方法，其特征在于：包括步骤：获取包含脊柱的正侧位数字重建放射影像；对所述正侧位数字重建放射影像进行分割得到各节脊柱区域图像与非脊柱区域图像；获取所述各节脊柱区域图像中脊柱的最小外接矩形，并根据其中心点得到三维影像中各节脊柱的空间关键点；以所述各节脊柱的空间关键点作为中心，以设定柱状体形式对三维影像沿脊柱高度方向进行分段切割得到各节脊柱分段切片图像；将各节脊柱分段切片图像进行拼接得到三维影像中脊柱的分割结果。2.根据权利要求1所述的基于三维影像的脊柱分割方法，其特征在于：所述获取所述各节脊柱区域图像中脊柱的最小外接矩形具体为：根据分割得到正侧位数字重建放射影像中所述各节脊柱区域图像的各像素点的坐标，以该节脊柱的像素点在对应影像的二维坐标系的两个坐标轴中的最大坐标差值作为最小外接矩形的长与宽，据此得到各节脊柱区域图像中脊柱的最小外接矩形。3.根据权利要求1或2所述的基于三维影像的脊柱分割方法，其特征在于：根据所述最小外接矩形的中心点基于最小二乘法计算得到三维影像中各节脊柱的空间关键点。4.根据权利要求2所述的基于三维影像的脊柱分割方法，其特征在于：根据分割得到正侧位数字重建放射影像中所述各节脊柱区域图像的各像素点的坐标包括以下步骤：根据各节脊柱区域图像与非脊柱区域图像的像素级别进行像素分割，并对...

【专利技术属性】
技术研发人员：沈林鹏，
申请(专利权)人：南京佗道医疗科技有限公司，
类型：发明
国别省市：