一种脑血管造影中伪影去除方法及系统技术方案

本发明专利技术提供了一种脑血管造影中伪影去除方法及系统，按照C型臂拍摄的时间，将脑部X光图像分为三类，分别为造影剂注入前脑部X光图像序列，注入造影剂并经过预设时间后的脑部X光图像序列，以及二者之间的脑部X光图像序列，对二者之间的脑部X光图像序列进行处理，得到只包括伪影的图像，再对通过传统方式获取的脑血管造影图像进行减影，进而实现对伪影的消除。本发明专利技术相较于传统的脑血管造影伪影去除方法，具有更高的去除伪影准确度，能够得到更为清晰的脑血管造影图像。清晰的脑血管造影图像。清晰的脑血管造影图像。

【技术实现步骤摘要】
一种脑血管造影中伪影去除方法及系统


[0001]本申请涉及医疗领域，尤其涉及脑血管造影中伪影去除方法及系统。

技术介绍

[0002]脑血管疾病是一种发病率高、容易复发的疾病，根据《2020中国卫生统计年鉴》的统计，2018年脑血管疾病的患病率为22.9
‰
，脑血管疾病已经成为严重危害人体健康的疾病。脑血管造影是诊断脑血管疾病的重要方式，尤其是数字减影(Digital Subtraction Angiography，DAS)技术的应用，可以获得高清晰度的脑血管图像，DSA被称为脑血管疾病诊断的金标准。通过脑血管造影可以直接观察脑血管是否存在畸形、狭窄、扩张等情况，DSA已经是脑血管疾病诊断中重要的方式。
[0003]脑血管造影是将造影剂注入到脑血管中，将造影剂注入前和注入后的两张图片相减，即可得到只包括脑血管的影像图片，其中将注入造影剂后拍摄的图像称为盈片或者活片，将未注入造影剂拍摄的图像称为蒙片(mask image)，减影也即将盈片和蒙片进行相减。但是由于设备本身、患者的呼吸等的影响，生成的脑血管图像会存在伪影，脑血管图像中的伪影会影响医生的判断，但是伪影是不可消除的，只能最大限度的减少。除了减少血管造影过程中周围环境的干扰降低伪影外，还可以通过后期处理的方式减少伪影，常用的伪影消除方式通过图像配准的方法尽量消除伪影。但是通过图像配准的方式消除伪影很多时候并不理想，有时候还会存在伪影的情况，影响医生的判断。

技术实现思路

[0004]为了尽可能的消除脑血管造影中的伪影，一方面，本专利技术提出了一种脑血管造影中伪影去除方法，所述方法包括以下步骤：
[0005]S1，获取脑部X光图像序列，将所述脑部X光图像分为三类，分别记为A
i
、B
j
、C
k
，其中A
i
表示造影剂注入前脑部X光图像序列，C
k
表示注入造影剂并经过预设时间后的脑部X光图像序列，B
j
表示A
i
和C
k
之间的脑部X光图像序列；其中i、j、k为正整数；
[0006]S2，根据脑部X光图像拍摄间隔从B
j
中选择n个第一盈片、第一蒙片对，对每个所述第一盈片、第一蒙片对分别在配准后进行减影处理，得到n个减影图像，对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，根据第一集合得到第一图像；其中，n为大于2的正整数；
[0007]S3，从A
i
、C
k
中分别选择第二蒙片和第二盈片，将第二盈片和第二蒙片在配准后进行减影处理，得到待处理减影图像，根据所述待处理减影图像和所述第一图像输出脑血管造影图像。
[0008]优选地，所述对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，具体包括以下过程：
[0009]S21，判断减影图像中是否存在像素点的灰度值都小于第二阈值的区域，若存在，将所述区域内所有像素点的灰度值调整为255，得到调整后的减影图像；
[0010]S22，对所述调整后的减影图像进行边缘检测，若边缘检测得到边缘长度之和大于第三阈值，将所述调整后的减影图像放入第一集合中。
[0011]优选地，所述根据第一集合得到第一图像，具体为：
[0012]以所述第一集合中第一图像为基准图像，将第一集合中其他图像分别和所述基准图像配准，累加所述基准图像和配准后的其他图像对应位置的像素点的灰度值得到累加灰度值H，计算平均灰度值h＝H/s，其中s表示第一集合中图像个数，根据每个像素点的平均灰度值得到第一图像。
[0013]优选地，所述根据所述待处理减影图像和所述第一图像输出脑血管造影图像，具体为：
[0014]对所述待处理减影图像进行特征检测，判断所述待处理减影图像是否存在伪影，若存在，则将所述待处理减影图像和所述第一图像配准后减影，输出脑血管造影图像，否则，直接将所述待处理减影图像输出。
[0015]优选地，所述配准包括第一配准方式和第二配准方式，所述第一配准方式仅包括粗配准过程，所述第二配准方式是先进行粗配准、再进行细配准。
[0016]优选地，在步骤S1前，还包括判断系统空闲资源情况，根据系统空闲资源情况确定后续处理过程涉及到配准到过程是采用第一配准方法，还是第二配准方法。
[0017]另一方面，本专利技术还提供了一种脑血管造影中伪影去除系统，所述系统包括以下模块：
[0018]脑血管图像分类模块，用于获取脑部X光图像序列，将所述脑部X光图像分为三类，分别记为A
i
、B
j
、C
k
，其中A
i
表示造影剂注入前脑部X光图像序列，C
k
表示注入造影剂并经过预设时间后的脑部X光图像序列，B
j
表示A
i
和C
k
之间的脑部X光图像序列；其中i、j、k为正整数；
[0019]第一图像获取模块，用于根据脑部X光图像拍摄间隔从B
j
中选择n个第一盈片、第一蒙片对，对每个所述第一盈片、第一蒙片对分别在配准后进行减影处理，得到n个减影图像，对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，根据第一集合得到第一图像；其中，n为大于2的正整数；
[0020]伪影去除模块，用于从A
i
、C
k
中分别选择第二蒙片和第二盈片，将第二盈片和第二蒙片在配准后进行减影处理，得到待处理减影图像，根据所述待处理减影图像和所述第一图像输出脑血管造影图像。
[0021]优选地，所述对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，具体为：
[0022]判断减影图像中是否存在像素点的灰度值都小于第二阈值的区域，若存在，将所述区域内所有像素点的灰度值调整为255，得到调整后的减影图像；
[0023]对所述调整后的减影图像进行边缘检测，若边缘检测得到边缘长度之和大于第三阈值，将所述调整后的减影图像放入第一集合中。
[0024]优选地，所述根据第一集合得到第一图像，具体为：
[0025]以所述第一集合中第一图像为基准图像，将第一集合中其他图像分别和所述基准图像配准，累加所述基准图像和配准后的其他图像对应位置的像素点的灰度值得到累加灰度值H，计算平均灰度值h＝H/s，其中s表示第一集合中图像个数，根据每个像素点的平均灰度值得到第一图像。
[0026]优选地，所述根据所述待处理减影图像和所述第一图像输出脑血管造影图像，具体为：
[0027]对所述待处理减影图像进行特征检测，判断所述待处理减影图像是否存在伪影，若存在，则将所述待处理减影图像和所述第一图像配准后减影，输出脑血管造影图像，否则，直接将所述待处理减影图像输出。
[0028]在造影剂注入过程中，C型臂拍摄的图像会先显示主血管，主血管相比微血管更本文档来自技高网...

【技术保护点】

【技术特征摘要】
1.一种脑血管造影中伪影去除方法，其特征在于，所述方法包括以下步骤：S1，获取脑部X光图像序列，将所述脑部X光图像分为三类，分别记为A
i
、B
j
、C
k
，其中A
i
表示造影剂注入前脑部X光图像序列，C
k
表示注入造影剂并经过预设时间后的脑部X光图像序列，B
j
表示A
i
和C
k
之间的脑部X光图像序列；其中i、j、k为正整数；S2，根据脑部X光图像拍摄间隔从B
j
中选择n个第一盈片、第一蒙片对，对每个所述第一盈片、第一蒙片对分别在配准后进行减影处理，得到n个减影图像，对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，根据第一集合得到第一图像；其中，n为大于2的正整数；S3，从A
i
、C
k
中分别选择第二蒙片和第二盈片，将第二盈片和第二蒙片在配准后进行减影处理，得到待处理减影图像，根据所述待处理减影图像和所述第一图像输出脑血管造影图像。2.如权利要求1所述的方法，其特征在于，所述对所述n个减影图像预处理，将预处理后的减影图像放入第一集合中，具体包括以下过程：S21，判断减影图像中是否存在像素点的灰度值都小于第二阈值的区域，若存在，将所述区域内所有像素点的灰度值调整为固定值，得到调整后的减影图像；S22，对所述调整后的减影图像进行边缘检测，若边缘检测得到边缘长度之和大于第三阈值，将所述调整后的减影图像放入第一集合中。3.如权利要求2所述的方法，其特征在于，所述根据第一集合得到第一图像，具体为：以所述第一集合中第一图像为基准图像，将第一集合中其他图像分别和所述基准图像配准，累加所述基准图像和配准后的其他图像对应位置的像素点的灰度值得到累加灰度值H，计算平均灰度值h＝H/s，其中s表示第一集合中图像个数，根据每个像素点的平均灰度值得到第一图像。4.如权利要求1所述的方法，其特征在于，所述根据所述待处理减影图像和所述第一图像输出脑血管造影图像，具体为：对所述待处理减影图像进行特征检测，判断所述待处理减影图像是否存在伪影，若存在，则将所述待处理减影图像和所述第一图像配准后减影，输出脑血管造影图像，否则，直接将所述待处理减影图像输出。5.如权利要求1
‑
4任一项所述的方法，其特征在于，所述配准包括第一配准方式和第二配准方式，所述第一配准方式仅包括粗配准过程，所述第二配准方式是先进行粗配准、再进行细配准。6.如权利要...

【专利技术属性】
技术研发人员：毛立武，康玉琪，
申请(专利权)人：郑州市中心医院，
类型：发明
国别省市：