血管造影中介下先心病术中血管辨识的装置及方法制造方法及图纸

一种血管造影中介下先心病术中血管辨识的装置及方法，其装置包括术前的CT图像预处理模块、分割模块、三维重建模块、二维平面化模块；术中的造影成像模块、配准初始化模块、量子更新模块、个体变异模块、变异后处理模块以及图像配准模块、和血管辨识结果展示模块。本发明专利技术方法的步骤为：使用分割术前CT图像，对分割完的CT图像进行三维可视化重建，得到分割的术前心脏和血管的二维图像，将分割的术前心脏和血管的二维图像与术中造影像配准，并展示配准后的结果。本发明专利技术配准的结果具有较高的精度和丰富的细节信息，可为进行先心病手术中的医生提供患者的血管的准确位置，降低手术的难度。

【技术实现步骤摘要】

本专利技术属于医疗
，更进一步涉及医学影像处理与应用
中的一种血管造影中介下的先心病手术中血管辨识的装置及方法。本专利技术利用血管造影中介下的先心病术中血管辨识的装置中的图像提取模块，对先心病患者的术前CT图像进行分割，然后将分割后得到的心脏及周围血管的三维模型进行平面化处理，得到若干张CT二维图像。再利用本专利技术的装置术中造影成像模块获取先心病患者的心脏及周围血管的术中造影图像，再利用本专利技术的装置中的图像配准模块将CT二维图像和术中造影图像进行配准，所获得的配准图像可以为手术中的医生提供血管辨识参考。
技术介绍
在先天性心脏病手术中，能为手术中的医生提供血管辨识辅助的有术前CT图像、MR图像等，但是由于患儿个体较小，血管与健康人差异较大，并且开胸之后血肉模糊，医生难以单凭一种图像完成血管辨识。对于在先心病手术中的血管辨识这一问题，并没有相应的专利或文献予以说明。而与之相近的一些配准技术，也存在诸多缺陷，因此并不适用于术中这一场景。王婕妤等人在其发表的论文“基于改进光流场和尺度不变特征变换的非刚性医学图像配准”(电子与信息学报，vol.35,No.5,May.2013)中提出了一种加入尺度不变特征变换特征的光流场模型的配准方法。该方法首先使用尺度不变特征变换方法对源图像和目标图像进行特征点的提取及匹配，并且使用随机抽样一致方法去除部分误匹配点，接着为特征点赋权值，然后在确定特征点的速度场后建立特征匹配项，极小化得到的能量函数完成图像配准。该方法的不足之处在于，两幅图像的特征点不易找到，更不易准确匹配，且内部纹理细节不明显。周其慧等人在其发表的论文“基于多分辨率和Powell算法的医学图像配准”(计算机应用研究，vol.30,No.4,Apr.2013)中提出了一种改进的Powell算法与多分辨策略相结合的配准方法。该方法先将待配准的两幅图像进行两层小波分解，得到图像的近似分量，然后配准最低分辨率的近似分量，再以上一层的结果作为粗略位置，进行下一轮循环，直到最底层的源图像，最终得到搜索结果完成图像配准。该方法存在的不足之处是，尽管该算法配准所需的时间相对于基本的Powell算法已有所减少，但仍不够快，而且仍然容易陷入局部最优结果，无法保证配准的精度。
技术实现思路
本专利技术的目的在于克服上述现有技术的不足，提出了一种血管造影中介下的先心病手术中血管辨识的装置及方法。实现本专利技术目的的思路是，分割出先心病患者术前CT图像中的血管和心脏，并将其进行三维可视化重建，对重建得到的模型进行平面化处理，得到术前二维的心脏和血管图像，在得到术中心脏和血管的造影像后，将两种图像配准，并展示配准的结果。为实现上述目的，本专利技术的具体步骤如下：本专利技术的装置包括术前CT图像预处理模块、术前CT图像分割模块、术前CT图像三维重建模块、术前CT图像二维平面化模块、术中造影成像模块、术中配准初始化模块、术中量子更新模块、术中个体变异模块、术中变异后处理模块、图像配准模块、血管辨识结果展示模块。所述的术前CT图像预处理模块，用于从术前心脏造影CT图像中选择一幅主动脉瓣CT图像作为初始图像；标记初始图像中心脏和血管所在的位置，并将其作为图割方法中无向图的权值；所述的术前CT图像分割模块，是利用图割方法，对术前CT图像预处理模块得到的初始图像进行分割，得到分割后的参考模板；用参考模板中心脏和血管所在的位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的前一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最前一幅图像，若是，则执行下一步，否则，继续分割图像；用参考模板中的心脏和血管位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的后一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最后一幅图像，若是，则将所有的参考模板传入三维重建模块，否则，继续分割图像；所述的术前CT图像三维重建模块，是将所有的参考模板加载到具有三维可视化功能的软件中；利用软件的体绘制选项，对所有的参考模板中的心脏和血管，进行三维可视化重建；所述的术前CT图像二维平面化模块，是将先心病患者手术时拟平躺于水平手术台且面向医生时的心脏位置作为基准，以先心病患者心脏的主动脉窦的中心为坐标原点；以过坐标原点且垂直于手术台平面的直线为Z轴，将先心病患者面部所在的方向确定为Z轴正方向；以过原点且垂直于脊柱中心线所在平面的直线为X轴，将指向先心病患者左臂的方向确定为X轴正方向；以过原点且垂直于X轴和Z轴所在平面的直线为Y轴，将指向先心病患者头部的方向确定为Y轴正方向，完成术前心脏和血管标准坐标系的构建；构建术前心脏和血管标准坐标系完成后，保持建立标准坐标系时分割的术前心脏和血管的位置不变，从Z轴正半轴看向负半轴的方向上，得到一张分割的术前心脏和血管的二维图像；以Y轴为旋转轴，将分割的术前心脏和血管分别向X轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；以X轴为旋转轴，将分割的术前心脏和血管分别向Y轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；所述的术中造影成像模块，用于获取先心病患者心脏和血管的术中造影像；所述的术中配准初始化模块，是将从术前CT图像二维平面化模块中获得的五张分割的术前心脏和血管的二维图像中任选一张图像；分别将配准所需要的水平位移种群、垂直位移种群、旋转角度种群中的每个种群初始化为5个个体，每个个体在其取值范围内取一个随机值；利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像与术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；将水平位移种群、垂直位移种群的复制规模设定为50，旋转角度种群的复制规模设定为20，根据每种个体组合的互信息占所有个体组合的互信息总和的比例，对个体组合中的个体按比例进行复制，生成复制后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中量子更新模块，是利用量子更新方法，对术中配准初始化模块复制后的水平位移种群、垂直位移种群、旋转角度种群中的个体进行量子更新，得到更新后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中个体变异模块，是将每个种群中个体发生变异的概率设定为20％，对术中量子更新模块更新后的水平位移种群、垂直位移种群、旋转角度种群中变异的个体在其取值范围内取随机值，得到变异后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中变异后处理模块，是对术中个体变异模块获得的变异后的水平位移种群、垂直位移种群、旋转角度种群，利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像和术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；判断当前的极大互信息是否小于配准初始化模块的极大互信息，若是，则获得极大互信息中的最大值，并执行下一步，否则，执行水平位移种群、垂直位移种群、旋转角度种群的复制步骤；判断是否选完五张分割的术前心脏和血管的二维图像，若是，则执行造影像配准步骤，否则，执行得到图像的最大互信息的步骤本文档来自技高网...

【技术保护点】
一种血管造影中介下先心病术中血管辨识的装置，包括术前CT图像预处理模块、术前CT图像分割模块、术前CT图像三维重建模块、术前CT图像二维平面化模块、术中造影成像模块、术中配准初始化模块、术中量子更新模块、术中个体变异模块、术中变异后处理模块、图像配准模块、血管辨识结果展示模块，其中：所述的术前CT图像预处理模块，用于从术前心脏造影CT图像中选择一幅主动脉瓣CT图像作为初始图像；标记初始图像中心脏和血管所在的位置，并将其作为图割方法中无向图的权值；所述的术前CT图像分割模块，是利用图割方法，对术前CT图像预处理模块得到的初始图像进行分割，得到分割后的参考模板；用参考模板中心脏和血管所在的位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的前一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最前一幅图像，若是，则执行下一步，否则，继续分割图像；用参考模板中的心脏和血管位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的后一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最后一幅图像，若是，则将所有的参考模板传入三维重建模块，否则，继续分割图像；所述的术前CT图像三维重建模块，是将所有的参考模板加载到具有三维可视化功能的软件中；利用软件的体绘制选项，对所有的参考模板中的心脏和血管，进行三维可视化重建；所述的术前CT图像二维平面化模块，是将先心病患者手术时拟平躺于水平手术台且面向医生时的心脏位置作为基准，以先心病患者心脏的主动脉窦的中心为坐标原点；以过坐标原点且垂直于手术台平面的直线为Z轴，将先心病患者面部所在的方向确定为Z轴正方向；以过原点且垂直于脊柱中心线所在平面的直线为X轴，将指向先心病患者左臂的方向确定为X轴正方向；以过原点且垂直于X轴和Z轴所在平面的直线为Y轴，将指向先心病患者头部的方向确定为Y轴正方向，完成术前心脏和血管标准坐标系的构建；构建术前心脏和血管标准坐标系完成后，保持建立标准坐标系时分割的术前心脏和血管的位置不变，从Z轴正半轴看向负半轴的方向上，得到一张分割的术前心脏和血管的二维图像；以Y轴为旋转轴，将分割的术前心脏和血管分别向X轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；以X轴为旋转轴，将分割的术前心脏和血管分别向Y轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；所述的术中造影成像模块，用于获取先心病患者心脏和血管的术中造影像；所述的术中配准初始化模块，是将从术前CT图像二维平面化模块中获得的五张分割的术前心脏和血管的二维图像中任选一张图像；分别将配准所需要的水平位移种群、垂直位移种群、旋转角度种群中的每个种群初始化为5个个体，每个个体在其取值范围内取一个随机值；利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像与术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；将水平位移种群、垂直位移种群的复制规模设定为50，旋转角度种群的复制规模设定为20，根据每种个体组合的互信息占所有个体组合的互信息总和的比例，对个体组合中的个体按比例进行复制，生成复制后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中量子更新模块，是利用量子更新方法，对术中配准初始化模块复制后的水平位移种群、垂直位移种群、旋转角度种群中的个体进行量子更新，得到更新后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中个体变异模块，是将每个种群中个体发生变异的概率设定为20％，对术中量子更新模块更新后的水平位移种群、垂直位移种群、旋转角度种群中变异的个体在其取值范围内取随机值，得到变异后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中变异后处理模块，是对术中个体变异模块获得的变异后的水平位移种群、垂直位移种群、旋转角度种群，利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像和术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；判断当前的极大互信息是否小于配准初始化模块的极大互信息，若是，则获得极大互信息中的最大值，并执行下一步，否则，执行水平位移种群、垂直位移种群、旋转角度种群的复制步骤；判断是否选完五张分割的术前心脏和血管的二维图像，若是，则执行造影像配准步骤，否则，执行得到图像的最大互信息的步骤；所述的图像配准模块，是从五张分割的术前心脏和血管的二维图像中找出与术中造影像的最大互信息；从五张分割的术前心脏和血管的二维图像中找出得到最大互信息的分割的术前心脏和血管的二维图像；从五张分割的术前心脏和血管的二维图像...

【技术特征摘要】
1.一种血管造影中介下先心病术中血管辨识的装置，包括术前CT图像预处理模块、术前CT图像分割模块、术前CT图像三维重建模块、术前CT图像二维平面化模块、术中造影成像模块、术中配准初始化模块、术中量子更新模块、术中个体变异模块、术中变异后处理模块、图像配准模块、血管辨识结果展示模块，其中：所述的术前CT图像预处理模块，用于从术前心脏造影CT图像中选择一幅主动脉瓣CT图像作为初始图像；标记初始图像中心脏和血管所在的位置，并将其作为图割方法中无向图的权值；所述的术前CT图像分割模块，是利用图割方法，对术前CT图像预处理模块得到的初始图像进行分割，得到分割后的参考模板；用参考模板中心脏和血管所在的位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的前一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最前一幅图像，若是，则执行下一步，否则，继续分割图像；用参考模板中的心脏和血管位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的后一幅图像进行分割，得到并保存更新后的参考模板；判断当前参考模板更新时所选取的图像是否为所有CT图像的最后一幅图像，若是，则将所有的参考模板传入三维重建模块，否则，继续分割图像；所述的术前CT图像三维重建模块，是将所有的参考模板加载到具有三维可视化功能的软件中；利用软件的体绘制选项，对所有的参考模板中的心脏和血管，进行三维可视化重建；所述的术前CT图像二维平面化模块，是将先心病患者手术时拟平躺于水平手术台且面向医生时的心脏位置作为基准，以先心病患者心脏的主动脉窦的中心为坐标原点；以过坐标原点且垂直于手术台平面的直线为Z轴，将先心病患者面部所在的方向确定为Z轴正方向；以过原点且垂直于脊柱中心线所在平面的直线为X轴，将指向先心病患者左臂的方向确定为X轴正方向；以过原点且垂直于X轴和Z轴所在平面的直线为Y轴，将指向先心病患者头部的方向确定为Y轴正方向，完成术前心脏和血管标准坐标系的构建；构建术前心脏和血管标准坐标系完成后，保持建立标准坐标系时分割的术前心脏和血管的位置不变，从Z轴正半轴看向负半轴的方向上，得到一张分割的术前心脏和血管的二维图像；以Y轴为旋转轴，将分割的术前心脏和血管分别向X轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；以X轴为旋转轴，将分割的术前心脏和血管分别向Y轴的正方向和反方向偏移五度，从Z轴正半轴看向负半轴的方向上，分别得到两张分割的术前心脏和血管的二维图像；所述的术中造影成像模块，用于获取先心病患者心脏和血管的术中造影像；所述的术中配准初始化模块，是将从术前CT图像二维平面化模块中获得的五张分割的术前心脏和血管的二维图像中任选一张图像；分别将配准所需要的水平位移种群、垂直位移种群、旋转角度种群中的每个种群初始化为5个个体，每个个体在其取值范围内取一个随机值；利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像与术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；将水平位移种群、垂直位移种群的复制规模设定为50，旋转角度种群的复制规模设定为20，根据每种个体组合的互信息占所有个体组合的互信息总和的比例，对个体组合中的个体按比例进行复制，生成复制后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中量子更新模块，是利用量子更新方法，对术中配准初始化模块复制后的水平位移种群、垂直位移种群、旋转角度种群中的个体进行量子更新，得到更新后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中个体变异模块，是将每个种群中个体发生变异的概率设定为20％，对术中量子更新模块更新后的水平位移种群、垂直位移种群、旋转角度种群中变异的个体在其取值范围内取随机值，得到变异后的水平位移种群、垂直位移种群、旋转角度种群；所述的术中变异后处理模块，是对术中个体变异模块获得的变异后的水平位移种群、垂直位移种群、旋转角度种群，利用个体组合的互信息计算方法，计算所选的分割的术前心脏和血管的二维图像和术中造影像配准的互信息，得到当前的极大互信息和当前的极大互信息的个体组合；判断当前的极大互信息是否小于配准初始化模块的极大互信息，若是，则获得极大互信息中的最大值，并执行下一步，否则，执行水平位移种群、垂直位移种群、旋转角度种群的复制步骤；判断是否选完五张分割的术前心脏和血管的二维图像，若是，则执行造影像配准步骤，否则，执行得到图像的最大互信息的步骤；所述的图像配准模块，是从五张分割的术前心脏和血管的二维图像中找出与术中造影像的最大互信息；从五张分割的术前心脏和血管的二维图像中找出得到最大互信息的分割的术前心脏和血管的二维图像；从五张分割的术前心脏和血管的二维图像中找出得到最大互信息的个体组合；用最大互信息的个体组合中的水平位移、垂直位移、旋转角度，对最大互信息的分割的术前心脏和血管的二维图像与术中造影像进行配准；所述的血管辨识结果展示模块，是将术前三维参考模板的角度调整到与配准所用的分割的术前心脏和血管的二维图像的角度一致后，将先心病患者术前血管的位置输出到显示器展示。2.一种血管造影中介下先心病术中血管辨识的方法，包括如下步骤：(1)分割术前心脏造影CT图像：(1a)术前CT图像预处理模块，从术前心脏造影CT图像中选择一幅主动脉瓣CT图像作为初始图像；(1b)术前CT图像预处理模块标记初始图像中心脏和血管所在的位置，并将其作为图割方法中无向图的权值；(1c)术前CT图像分割模块利用图割方法，对初始图像进行分割，得到分割后的参考模板；(1d)术前CT图像分割模块用参考模板中心脏和血管所在的位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的前一幅图像进行分割，得到并保存更新后的参考模板；(1e)术前CT图像分割模块判断当前参考模板更新时所选取的图像是否为所有术前心脏造影CT图像的最前一幅图像，若是，则执行步骤(1f)，否则，执行步骤(1d)；(1f)术前CT图像分割模块用参考模板中心脏和血管所在的位置作为图割方法的新的权值，采用图割方法，对当前参考模板分割图像时的后一幅图像进行分割，得到并保存更新后的参考模板；(1g)判断当前参考模板更新时所选取的图像是否为所有术前心脏造影CT图像的最后一幅图像，若是，则执行步骤(2)，否则，执行步骤(1f)；(2)重建术前三维参考模板：(2a)术前CT图像三维重建模块将所有参考模板加载到具有三维可视化功能的软件中；(2b)术前CT图像三维重建模块利用软件的体绘制选项，对参考模板中的心脏和血管，进行三维可视化重建；(3)构建术前心脏和血管标准坐标系：(3a)术前CT图像二维平面化模块将先心病患者手术时拟平躺于水平手术台且面向医生时的心脏位置作为基准，以先心病患者心脏的主动脉窦的中心为坐标原点；(3b)术前CT图像二维平面化模块，以过坐标原点且垂直于手术台平面的直线为Z轴，将先心病患者面部所在的方向确定为Z轴正方向；(3c)术前CT图像二维平面化模块，以过原点且垂直于脊柱中心线所在平面的直线为X轴，将指向先心病患者左臂的方向确定为X轴正方向；(3d)术前CT图像二维平面化模块，以过原点且垂直于X轴和Z轴所在平面的直线为Y轴，将指向先心病患者头部的方向确定为Y轴正方向，完成术前心脏和血管标准坐标系的构建；(4)将分割的术前心脏和血管平面化：(4a)保持建立标准坐标系时分割的术前心脏和血管的位置不变，术前CT图像二维平面化模块，从Z轴正半轴看向负半轴的方向上，得到一张分割的术前心脏和血管的二维图像；(4b)以Y轴为旋转轴，术前CT图像二维平面化模...

【专利技术属性】
技术研发人员：黄力宇，王芮东，陈煜，黄美萍，李军，庄建，
申请(专利权)人：西安电子科技大学，
类型：发明
国别省市：陕西;61