CN109758132A - 一种血管狭窄的检测方法及设备 - Google Patents
一种血管狭窄的检测方法及设备 Download PDFInfo
- Publication number
- CN109758132A CN109758132A CN201910142018.8A CN201910142018A CN109758132A CN 109758132 A CN109758132 A CN 109758132A CN 201910142018 A CN201910142018 A CN 201910142018A CN 109758132 A CN109758132 A CN 109758132A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- curved section
- diameter
- edge curved
- point
- Prior art date
- Legal status (The legal status is an assumption and is not a legal conclusion. Google has not performed a legal analysis and makes no representation as to the accuracy of the status listed.)
- Pending
Links
- 238000001514 detection method Methods 0.000 title claims abstract description 24
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 308
- 230000000630 rising effect Effects 0.000 claims abstract description 183
- 238000000034 method Methods 0.000 claims abstract description 26
- 239000008280 blood Substances 0.000 claims description 14
- 210000004369 blood Anatomy 0.000 claims description 14
- 230000007423 decrease Effects 0.000 claims description 9
- 230000002792 vascular Effects 0.000 claims description 7
- 239000000203 mixture Substances 0.000 description 6
- 208000031481 Pathologic Constriction Diseases 0.000 description 4
- 238000010586 diagram Methods 0.000 description 4
- 208000037804 stenosis Diseases 0.000 description 4
- 230000036262 stenosis Effects 0.000 description 4
- 238000005516 engineering process Methods 0.000 description 2
- 239000000463 material Substances 0.000 description 2
- 210000004556 brain Anatomy 0.000 description 1
- 201000010099 disease Diseases 0.000 description 1
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 1
- 230000002526 effect on cardiovascular system Effects 0.000 description 1
- 238000007689 inspection Methods 0.000 description 1
- 238000005259 measurement Methods 0.000 description 1
- 230000000250 revascularization Effects 0.000 description 1
- 230000011218 segmentation Effects 0.000 description 1
- 208000024891 symptom Diseases 0.000 description 1
Landscapes
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本发明公开了一种血管狭窄的检测方法及设备,所述方法包括:获取血管直径曲线图;计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值;根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段;将所述血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为所述血管的一处狭窄区域。本发明通过对血管直径进行梯度计算,再根据血管直径梯度值查找上升沿曲线段和对应下降沿曲线段的方法来确定血管狭窄区域,实现了血管狭窄检测的自动化,不需人工辨别血管狭窄情况,省时省力,而且精准度高。
Description
技术领域
本发明涉及血管图像技术领域,尤其涉及一种血管狭窄的检测方法及设备。
背景技术
在现代医学技术领域中,血管图像技术能够帮助医生了解病人的血管状态,有助于及时发现和诊断各种疾病,因此对医生有着重要的临床价值和实际意义。随着年龄的上升,血管出现问题的概率会越来越高,而血管狭窄是其中一个非常普遍的问题,包括如脑血管狭窄和心血管狭窄等都是中老年人高发病状。
现代医疗技术通过采用血管重建技术等能够自动建立血管模型并得到各种易于观察的血管拉直分割图,甚至自动测量出血管的直径曲线图像,便于医生检查血管情况。但是,在血管的直径曲线图建立完成后通常还需要人工辨别血管病况,查找狭窄区域,不仅耗时耗力,而且很容易出现遗漏或误差。
发明内容
本发明实施例为了有效克服现有技术所存在的上述缺陷,创造性地提供一种血管狭窄的检测方法及设备。
本发明一方面提供一种血管狭窄的检测方法及设备,所述方法包括:获取血管直径曲线图;计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值;根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段;将所述血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为所述血管的一处狭窄区域。
在一个实施例中,根据本发明的上述方法,在确定所述狭窄区域之后,所述方法还包括:针对每一处所述狭窄区域,根据所述狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;将所述血管最小直径除以所述血管理论直径,得到对应所述狭窄区域的血管狭窄率。
在一个实施例中,根据本发明的上述方法,所述计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值包括:遍历所述血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到所述血管直径梯度值。
在一个实施例中,根据本发明的上述方法,所述根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段包括:沿血管的根部向末端遍历所述血管直径梯度值,找到血管直径梯度值大于0的第一上升点;针对每一个所述第一上升点,分别从所述第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将所述第一上升点和所述第一下降点之间的血管直径曲线段确定为所述上升沿曲线段;从所述第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将所述第二上升点和所述第一上升点之间的血管直径曲线段确定为与所述上升沿曲线段对应的下降沿曲线段;其中,所述第一方向为自血管根部向末端的延伸方向,所述第二方向为自血管末端向根部的延伸方向。
在一个实施例中,根据本发明的上述方法,所述处理所述上升沿曲线段和所述下降沿曲线段数据,得到血管理论直径包括:针对每一个上升沿曲线段及对应的下降沿曲线段,计算所述上升沿曲线段的最大血管直径,作为第一血管直径;计算所述下降沿曲线段的最大血管直径,作为第二血管直径;对所述第一血管直径和所述第二血管直径进行加权平均计算,得到所述血管理论直径。
在一个实施例中,根据本发明的上述方法,所述处理所述上升沿曲线段和所述下降沿曲线段数据,得到血管最小直径包括:针对每一个上升沿曲线段及对应的下降沿曲线段,找到所述上升沿曲线段和所述下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。
本发明另一方面提供一种血管狭窄的检测设备,所述设备包括:曲线获取模块,用于获取血管直径曲线图;梯度计算模块,用于计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值;曲线查找模块,用于根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段;狭窄确定模块,用于将所述血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为所述血管的一处狭窄区域。
在一个实施例中,根据本发明的上述设备,所述设备还包括:狭窄计算模块,用于针对每一处所述狭窄区域,根据所述狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;将所述血管最小直径除以所述血管理论直径,得到对应所述狭窄区域的血管狭窄率。
在一个实施例中,根据本发明的上述设备,所述梯度计算模块进一步用于,遍历所述血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到所述血管直径梯度值。
在一个实施例中,根据本发明的上述设备,所述曲线查找模块包括:梯度查找子模块,用于沿血管的根部向末端遍历所述血管直径梯度值,找到血管直径梯度值大于0的第一上升点;上升查找子模块,用于针对每一个所述第一上升点,分别从所述第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将所述第一上升点和所述第一下降点之间的血管直径曲线段确定为所述上升沿曲线段;下降查找子模块,用于从所述第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将所述第二上升点和所述第一上升点之间的血管直径曲线段确定为与所述上升沿曲线段对应的下降沿曲线段;其中,所述第一方向为自血管根部向末端的延伸方向,所述第二方向为自血管末端向根部的延伸方向。
在一个实施例中,根据本发明的上述设备,所述狭窄计算模块包括:上升计算子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,计算所述上升沿曲线段的最大血管直径,作为第一血管直径;下降计算子模块,用于计算所述下降沿曲线段的最大血管直径,作为第二血管直径;加权计算子模块,用于对所述第一血管直径和所述第二血管直径进行加权平均计算,得到所述血管理论直径。
在一个实施例中,根据本发明的上述设备,所述狭窄计算模块包括:最小直径查找子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,找到所述上升沿曲线段和所述下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。
本发明通过获取血管直径曲线图,再针对获得的血管直径曲线图计算血管的直径梯度,得到血管直径梯度值。由于健康的血管从根部到末端总是越来越细的,其直径曲线是均匀变化的函数,当血管出现狭窄区域时,其直径曲线会出现上升点,因此本发明根据血管直径梯度值查找出上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段,再将血管中的每一个上升沿曲线段和对应于每一个上升沿曲线段组成的血管区域确定为血管的一处狭窄区域。通过对血管直径进行梯度计算,再根据血管直径梯度值查找上升沿曲线段和对应下降沿曲线段的方法来确定血管狭窄区域,实现了血管狭窄检测的自动化,不需人工辨别血管狭窄情况,省时省力,而且精准度高。
附图说明
通过参考附图阅读下文的详细描述,本发明示例性实施方式的上述以及其他目的、特征和优点将变得易于理解。在附图中,以示例性而非限制性的方式示出了本发明的若干实施方式,其中:
在附图中,相同或对应的标号表示相同或对应的部分。
图1为本发明实施例一种血管狭窄的检测方法的实现流程示意图;
图2为本发明实施例一种血管狭窄的检测设备的组成结构示意图。
具体实施方式
为使本发明的目的、特征、优点能够更加的明显和易懂,下面将结合本发明实施例中的附图,对本发明实施例中的技术方案进行清楚、完整地描述,显然,所描述的实施例仅仅是本发明一部分实施例,而非全部实施例。基于本发明中的实施例,本领域技术人员在没有做出创造性劳动前提下所获得的所有其他实施例,都属于本发明保护的范围。
图1为本发明实施例一种血管狭窄的检测方法的实现流程示意图,请参考图1。本发明实施例提供的一种血管狭窄的检测方法,方法包括:
步骤101,获取血管直径曲线图;
步骤102,计算血管直径曲线图中的血管直径梯度,得到血管直径梯度值;
步骤103,根据血管直径梯度值查找上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段;
步骤104,将血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为血管的一处狭窄区域。
本发明实施例通过步骤101获取血管直径曲线图,再针对获得的血管直径曲线图通过步骤102计算血管的直径梯度,得到血管直径梯度值。由于健康的血管从根部到末端总是越来越细的,其直径曲线是均匀变化的函数,当血管出现狭窄区域时,其直径曲线会出现上升点,因此本发明实施例通过步骤103根据血管直径梯度值查找出上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段,再通过步骤104将血管中的每一个上升沿曲线段和对应于每一个上升沿曲线段组成的血管区域确定为血管的一处狭窄区域。通过对血管直径进行梯度计算,再根据血管直径梯度值查找上升沿曲线段和对应下降沿曲线段的方法来确定血管狭窄区域,实现了血管狭窄检测的自动化,不需人工辨别血管狭窄情况,省时省力,而且精准度高。
在一可实施方式中,在确定狭窄区域之后,方法还包括:针对每一处狭窄区域,根据狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;将血管最小直径除以血管理论直径,得到对应狭窄区域的血管狭窄率。
本发明实施例通过对每一处狭窄区域进行计算,得到每一处狭窄区域的血管理论直径和血管最小直径,并进一步计算得到每一处狭窄区域的血管狭窄率,从而精准得到每一处狭窄区域的血管狭窄情况,并可以根据实际情况设置血管狭窄率阈值,将超过血管狭窄率阈值的数据进行输出或处理。使用本方案进行狭窄检测速度快,省时省力,同时,只要获取的血管直径数据测量准确,本发明可检测出所有狭窄,准确率较高。
在一可实施方式中,计算血管直径曲线图中的血管直径梯度,得到血管直径梯度值包括:遍历血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到血管直径梯度值。
由于健康的血管是从根部到末端越来越细的,因此,本方案中采用沿血管的根部向末端方向来计算每一个像素点相对于前一个像素点的血管直径差值,这样能够检测到图像上的每一个像素点,而且得到的血管直径梯度值较便于计算。也可根据实际需要选择其他的方向进行计算,只要能够使得血管直径曲线图像中的每一个像素点都能被计算到,从而保证血管狭窄检测的完整性和准确性即可。
在一可实施方式中,根据血管直径梯度值查找上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段包括:沿血管的根部向末端遍历血管直径梯度值,找到血管直径梯度值大于0的第一上升点;针对每一个第一上升点,分别从第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将第一上升点和第一下降点之间的血管直径曲线段确定为上升沿曲线段;从第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将第二上升点和第一上升点之间的血管直径曲线段确定为与上升沿曲线段对应的下降沿曲线段;其中,第一方向为自血管根部向末端的延伸方向,第二方向为自血管末端向根部的延伸方向。
本发明实施例通过沿血管的根部向末端方向先查找到每一个血管直径梯度值大于0的第一上升点,再将每一个上升点作为起始点,沿第一方向逐点进行查找,所遇到的血管直径梯度值大于0的上升点和血管直径梯度值等于0的平台点都与第一上升点位于同一个上升沿中,直至找到血管直径梯度值小于0的第一下降点,并将第一上升点和第一下降点之间的这整段血管直径曲线段确定为上升沿曲线段,这样确保了每一个上升点和上升沿曲线段都能够被检测计算到,有效保证检测结果的完整性。由于每一个狭窄都由一个上升沿曲线段和一个下降沿曲线段组成,这里将第一上升点作为出发点沿血管末端向根部反向查找,所遇到的血管直径梯度值小于0的下降点和血管直径梯度值等于0的平台点都位于同一个下降沿中,直至找到一个血管直径梯度值大于0的第二上升点,说明下降沿结束,将第二上升点和第一上升点之间的这整段血管直径曲线段确定为与上升沿曲线段对应的下降沿曲线段。这样,能够准确查找出每一个狭窄的上升沿曲线段和对应的下降沿曲线段,有效保证检测结果的完整性和准确性。
在一可实施方式中,处理上升沿曲线段和下降沿曲线段数据,得到血管理论直径包括:针对每一个上升沿曲线段及对应的下降沿曲线段,计算上升沿曲线段的最大血管直径,作为第一血管直径;计算下降沿曲线段的最大血管直径,作为第二血管直径;对第一血管直径和第二血管直径进行加权平均计算,得到血管理论直径。通过对上升沿曲线段和对应的下降沿曲线段分别计算最大血管直径,然后进行加权平均计算获得血管理论直径数据的方式,能够自动计算出每一处狭窄的血管理论直径,并且有效提高血管理论直径数据的准确性。
在一可实施方式中,处理上升沿曲线段和下降沿曲线段数据,得到血管最小直径包括:针对每一个上升沿曲线段及对应的下降沿曲线段,找到上升沿曲线段和下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。由于血管的狭窄由一个上升沿曲线段和一个下降沿曲线段组成,下降沿曲线段的最低点即上升沿的起始点,也就是上升沿曲线段和下降沿曲线段的连接点,通过自动查找每一个上升沿曲线段及对应的下降沿曲线段之间的连接点,并将此连接点数据作为每一处狭窄的血管最小直径的方法,简单快捷,检测精度较高。
图2为本发明实施例一种血管狭窄的检测设备的组成结构示意图,请参考图2。本发明实施例提供的一种血管狭窄的检测设备包括:
曲线获取模块201,用于获取血管直径曲线图;
梯度计算模块202,用于计算血管直径曲线图中的血管直径梯度,得到血管直径梯度值;
曲线查找模块203,用于根据血管直径梯度值查找上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段;
狭窄确定模块204,用于将血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为血管的一处狭窄区域。
本发明实施例通过曲线获取模块201获取血管直径曲线图,再针对获得的血管直径曲线图通过梯度计算模块202计算血管的直径梯度,得到血管直径梯度值。由于健康的血管从根部到末端总是越来越细的,其直径曲线是均匀变化的函数,当血管出现狭窄区域时,其直径曲线会出现上升点,因此本发明实施例通过曲线查找模块203根据血管直径梯度值查找出上升沿曲线段和对应于每一个上升沿曲线段的下降沿曲线段,再通过狭窄确定模块204将血管中的每一个上升沿曲线段和对应于每一个上升沿曲线段组成的血管区域确定为血管的一处狭窄区域。通过对血管直径进行梯度计算,再根据血管直径梯度值查找上升沿曲线段和对应下降沿曲线段的方法来确定血管狭窄区域,实现了血管狭窄检测的自动化,不需人工辨别血管狭窄情况,省时省力,而且精准度高。
在一可实施方式中,设备还包括:狭窄计算模块205,用于针对每一处狭窄区域,根据狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;将血管最小直径除以血管理论直径,得到对应狭窄区域的血管狭窄率。
本发明实施例通过狭窄计算模块205对每一处狭窄区域的曲线段进行计算,得到每一处狭窄区域的血管理论直径和血管最小直径,并进一步计算得到每一处狭窄区域的血管狭窄率,从而精准得到每一处狭窄区域的血管狭窄情况,并可以根据实际情况设置血管狭窄率阈值,将超过血管狭窄率阈值的数据进行输出或处理。使用本方案进行狭窄检测速度快,省时省力,同时,只要获取的血管直径数据测量准确,本发明可检测出所有狭窄,准确率较高。
在一可实施方式中,梯度计算模块202进一步用于,遍历血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到血管直径梯度值。
由于健康的血管是从根部到末端越来越细的,因此,本方案中梯度计算子模块采用沿血管的根部向末端方向来计算每一个像素点相对于前一个像素点的血管直径差值,这样能够检测到图像上的每一个像素点,而且得到的血管直径梯度值较便于计算。也可根据实际需要选择其他的方向进行计算,只要能够使得血管直径曲线图像中的每一个像素点都能被计算到,从而保证血管狭窄检测的完整性和准确性即可。
在一可实施方式中,曲线查找模块203包括:梯度查找子模块,用于沿血管的根部向末端遍历血管直径梯度值,找到血管直径梯度值大于0的第一上升点;上升查找子模块,用于针对每一个第一上升点,分别从第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将第一上升点和第一下降点之间的血管直径曲线段确定为上升沿曲线段;下降查找子模块,用于从第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将第二上升点和第一上升点之间的血管直径曲线段确定为与上升沿曲线段对应的下降沿曲线段;其中,第一方向为自血管根部向末端的延伸方向,第二方向为自血管末端向根部的延伸方向。
本发明实施例通过梯度查找子模块沿血管的根部向末端方向先查找到每一个血管直径梯度值大于0的第一上升点,再将每一个上升点作为起始点,通过上升查找子模块沿第一方向逐点进行查找,所遇到的血管直径梯度值大于0的上升点和血管直径梯度值等于0的平台点都与第一上升点位于同一个上升沿中,直至找到血管直径梯度值小于0的第一下降点,并将第一上升点和第一下降点之间的这整段血管直径曲线段确定为上升沿曲线段,这样确保了每一个上升点和上升沿曲线段都能够被检测计算到,有效保证检测结果的完整性。由于每一个狭窄都由一个上升沿曲线段和一个下降沿曲线段组成,这里通过下降查找子模块将第一上升点作为出发点沿血管末端向根部反向查找,所遇到的血管直径梯度值小于0的下降点和血管直径梯度值等于0的平台点都位于同一个下降沿中,直至找到一个血管直径梯度值大于0的第二上升点,说明下降沿结束,将第二上升点和第一上升点之间的这整段血管直径曲线段确定为与上升沿曲线段对应的下降沿曲线段。这样,能够准确查找出每一个狭窄的上升沿曲线段和对应的下降沿曲线段,有效保证检测结果的完整性和准确性。
在一可实施方式中,狭窄计算模块205包括:上升计算子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,计算上升沿曲线段的最大血管直径,作为第一血管直径;下降计算子模块,用于计算下降沿曲线段的最大血管直径,作为第二血管直径;加权计算子模块,用于对第一血管直径和第二血管直径进行加权平均计算,得到血管理论直径。通过上升计算子模块和下降计算子模块对上升沿曲线段和对应的下降沿曲线段分别计算最大血管直径,然后通过加权计算子模块进行加权平均计算获得血管理论直径数据的方式,能够自动计算出每一处狭窄的血管理论直径,并且有效提高血管理论直径数据的准确性。
在一可实施方式中,狭窄计算模块205包括:最小直径查找子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,找到上升沿曲线段和下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。由于血管的狭窄由一个上升沿曲线段和一个下降沿曲线段组成,下降沿曲线段的最低点即上升沿的起始点,也就是上升沿曲线段和下降沿曲线段的连接点,通过最小直径查找子模块自动查找每一个上升沿曲线段及对应的下降沿曲线段之间的连接点,并将此连接点数据作为每一处狭窄的血管最小直径的方法,简单快捷,检测精度较高。
在本说明书的描述中,参考术语“一个实施例”、“一些实施例”、“示例”、“具体示例”、或“一些示例”等的描述意指结合该实施例或示例描述的具体特征、结构、材料或者特点包含于本发明的至少一个实施例或示例中。而且,描述的具体特征、结构、材料或者特点可以在任一个或多个实施例或示例中以合适的方式结合。此外,在不相互矛盾的情况下,本领域的技术人员可以将本说明书中描述的不同实施例或示例以及不同实施例或示例的特征进行结合和组合。
此外,术语“第一”、“第二”仅用于描述目的,而不能理解为指示或暗示相对重要性或者隐含指明所指示的技术特征的数量。由此,限定有“第一”、“第二”的特征可以明示或隐含地包括至少一个该特征。在本发明的描述中,“多个”的含义是两个或两个以上,除非另有明确具体的限定。
以上,仅为本发明的具体实施方式,但本发明的保护范围并不局限于此,任何熟悉本技术领域的技术人员在本发明揭露的技术范围内,可轻易想到变化或替换,都应涵盖在本发明的保护范围之内。因此,本发明的保护范围应以权利要求的保护范围为准。
Claims (12)
1.一种血管狭窄的检测方法,其特征在于,所述方法包括:
获取血管直径曲线图;
计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值;
根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段;
将所述血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为所述血管的一处狭窄区域。
2.根据权利要求1所述的方法,其特征在于,在确定所述狭窄区域之后,所述方法还包括:
针对每一处所述狭窄区域,根据所述狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;
将所述血管最小直径除以所述血管理论直径,得到对应所述狭窄区域的血管狭窄率。
3.根据权利要求1或2所述的方法,其特征在于,所述计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值包括:
遍历所述血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到所述血管直径梯度值。
4.根据权利要求1或2所述的方法,其特征在于,所述根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段包括:
沿血管的根部向末端遍历所述血管直径梯度值,找到血管直径梯度值大于0的第一上升点;
针对每一个所述第一上升点,分别从所述第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将所述第一上升点和所述第一下降点之间的血管直径曲线段确定为所述上升沿曲线段;
从所述第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将所述第二上升点和所述第一上升点之间的血管直径曲线段确定为与所述上升沿曲线段对应的下降沿曲线段;
其中,所述第一方向为自血管根部向末端的延伸方向,所述第二方向为自血管末端向根部的延伸方向。
5.根据权利要求2所述的方法,其特征在于,所述处理所述上升沿曲线段和所述下降沿曲线段数据,得到血管理论直径包括:
针对每一个上升沿曲线段及对应的下降沿曲线段,计算所述上升沿曲线段的最大血管直径,作为第一血管直径;
计算所述下降沿曲线段的最大血管直径,作为第二血管直径;
对所述第一血管直径和所述第二血管直径进行加权平均计算,得到所述血管理论直径。
6.根据权利要求2或5所述的方法,其特征在于,所述处理所述上升沿曲线段和所述下降沿曲线段数据,得到血管最小直径包括:
针对每一个上升沿曲线段及对应的下降沿曲线段,找到所述上升沿曲线段和所述下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。
7.一种血管狭窄的检测设备,其特征在于,所述设备包括:
曲线获取模块,用于获取血管直径曲线图;
梯度计算模块,用于计算所述血管直径曲线图中的血管直径梯度,得到血管直径梯度值;
曲线查找模块,用于根据所述血管直径梯度值查找上升沿曲线段和对应于每一个所述上升沿曲线段的下降沿曲线段;
狭窄确定模块,用于将所述血管中每一个上升沿曲线段及对应的下降沿曲线段组成的血管区域确定为所述血管的一处狭窄区域。
8.根据权利要求7所述的设备,其特征在于,所述设备还包括:
狭窄计算模块,用于针对每一处所述狭窄区域,根据所述狭窄区域内的上升沿曲线段和下降沿曲线段的血管直径,得到血管理论直径和血管最小直径;将所述血管最小直径除以所述血管理论直径,得到对应所述狭窄区域的血管狭窄率。
9.根据权利要求7或8所述的设备,其特征在于,所述梯度计算模块进一步用于,遍历所述血管直径曲线图像中沿血管的根部向末端的每一个像素点,计算每一个像素点相对于前一个像素点的血管直径差值,得到所述血管直径梯度值。
10.根据权利要求7或8所述的设备,其特征在于,所述曲线查找模块包括:
梯度查找子模块,用于沿血管的根部向末端遍历所述血管直径梯度值,找到血管直径梯度值大于0的第一上升点;
上升查找子模块,用于针对每一个所述第一上升点,分别从所述第一上升点出发,沿第一方向逐点查找,直至找到血管直径梯度值小于0的第一下降点,将所述第一上升点和所述第一下降点之间的血管直径曲线段确定为所述上升沿曲线段;
下降查找子模块,用于从所述第一上升点出发沿第二方向逐点查找,直至找到血管直径梯度值大于0的第二上升点,将所述第二上升点和所述第一上升点之间的血管直径曲线段确定为与所述上升沿曲线段对应的下降沿曲线段;
其中,所述第一方向为自血管根部向末端的延伸方向,所述第二方向为自血管末端向根部的延伸方向。
11.根据权利要求8所述的设备,其特征在于,所述狭窄计算模块包括:
上升计算子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,计算所述上升沿曲线段的最大血管直径,作为第一血管直径;
下降计算子模块,用于计算所述下降沿曲线段的最大血管直径,作为第二血管直径;
加权计算子模块,用于对所述第一血管直径和所述第二血管直径进行加权平均计算,得到所述血管理论直径。
12.根据权利要求8或11所述的设备,其特征在于,所述狭窄计算模块包括:
最小直径查找子模块,用于针对每一个上升沿曲线段及对应的下降沿曲线段,找到所述上升沿曲线段和所述下降沿曲线段的连接点,得到连接点处的血管直径为血管最小直径。
Priority Applications (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910142018.8A CN109758132A (zh) | 2019-02-26 | 2019-02-26 | 一种血管狭窄的检测方法及设备 |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
CN201910142018.8A CN109758132A (zh) | 2019-02-26 | 2019-02-26 | 一种血管狭窄的检测方法及设备 |
Publications (1)
Publication Number | Publication Date |
---|---|
CN109758132A true CN109758132A (zh) | 2019-05-17 |
Family
ID=66457130
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
CN201910142018.8A Pending CN109758132A (zh) | 2019-02-26 | 2019-02-26 | 一种血管狭窄的检测方法及设备 |
Country Status (1)
Country | Link |
---|---|
CN (1) | CN109758132A (zh) |
Cited By (2)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112487342A (zh) * | 2020-11-25 | 2021-03-12 | 苏州润迈德医疗科技有限公司 | 精确获取狭窄病变区间的方法、系统及存储介质 |
CN114359128A (zh) * | 2021-09-10 | 2022-04-15 | 数坤(北京)网络科技股份有限公司 | 一种血管狭窄的检测方法、装置及计算机可读介质 |
Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101028198A (zh) * | 2006-01-17 | 2007-09-05 | 西门子公司 | 在检查区域内的图像数据的基础上检查患者的血管的方法 |
CN101991420A (zh) * | 2009-08-27 | 2011-03-30 | 上海西门子医疗器械有限公司 | 一种疑似血管狭窄区域的定位方法和系统 |
US20150268039A1 (en) * | 2014-03-18 | 2015-09-24 | Medis Associated B.V. | Method and device for determining deviation in pressure in a blood vessel |
US20160166209A1 (en) * | 2014-12-16 | 2016-06-16 | Siemens Healthcare Gmbh | Method and System for Personalized Non-Invasive Hemodynamic Assessment of Renal Artery Stenosis from Medical Images |
CN107613874A (zh) * | 2015-03-31 | 2018-01-19 | 新加坡科技研究局 | 用于评估血管狭窄的方法和装置 |
-
2019
- 2019-02-26 CN CN201910142018.8A patent/CN109758132A/zh active Pending
Patent Citations (5)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN101028198A (zh) * | 2006-01-17 | 2007-09-05 | 西门子公司 | 在检查区域内的图像数据的基础上检查患者的血管的方法 |
CN101991420A (zh) * | 2009-08-27 | 2011-03-30 | 上海西门子医疗器械有限公司 | 一种疑似血管狭窄区域的定位方法和系统 |
US20150268039A1 (en) * | 2014-03-18 | 2015-09-24 | Medis Associated B.V. | Method and device for determining deviation in pressure in a blood vessel |
US20160166209A1 (en) * | 2014-12-16 | 2016-06-16 | Siemens Healthcare Gmbh | Method and System for Personalized Non-Invasive Hemodynamic Assessment of Renal Artery Stenosis from Medical Images |
CN107613874A (zh) * | 2015-03-31 | 2018-01-19 | 新加坡科技研究局 | 用于评估血管狭窄的方法和装置 |
Non-Patent Citations (2)
Title |
---|
徐智: "心血管造影图像的二维信息处理及其三维重建研究", 《中国优秀博硕士学位论文全文数据库 (博士)信息科技辑》 * |
王一丁等: "《数字图像处理》", 31 August 2015 * |
Cited By (3)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
CN112487342A (zh) * | 2020-11-25 | 2021-03-12 | 苏州润迈德医疗科技有限公司 | 精确获取狭窄病变区间的方法、系统及存储介质 |
WO2022109909A1 (zh) * | 2020-11-25 | 2022-06-02 | 苏州润迈德医疗科技有限公司 | 精确获取狭窄病变区间的方法、系统及存储介质 |
CN114359128A (zh) * | 2021-09-10 | 2022-04-15 | 数坤(北京)网络科技股份有限公司 | 一种血管狭窄的检测方法、装置及计算机可读介质 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
EP3473177A1 (en) | Baseline impedance maps for tissue proximity indications | |
US10275881B2 (en) | Semi-automated image segmentation system and method | |
US20150297373A1 (en) | Method and apparatus for automated determination of a lumen contour of a stented blood vessel | |
RU2466679C2 (ru) | Визуализация сосудистой структуры | |
CN109758132A (zh) | 一种血管狭窄的检测方法及设备 | |
CN109872321A (zh) | 一种血管狭窄检测方法及设备 | |
CN109886982B (zh) | 一种血管图像分割方法、装置及计算机存储设备 | |
CN109288537B (zh) | 评估血流储备分数的系统、方法、设备及存储介质 | |
CN111932554A (zh) | 一种肺部血管分割方法、设备及存储介质 | |
EP4159131A1 (en) | Method and device for automatically processing blood vessel image | |
CN112215060A (zh) | 一种基于霍夫变换的高精度机械式仪表示数识别方法 | |
CN109345546A (zh) | 一种冠脉体数据模型分割方法及设备 | |
CN107392891B (zh) | 血管树提取方法、装置、设备及存储介质 | |
CN104688190A (zh) | 检测冠状动脉内支架贴壁情况的装置 | |
Tahoces et al. | Automatic detection of anatomical landmarks of the aorta in CTA images | |
CN115049807A (zh) | 肺部血管模型的建立方法、装置及服务器 | |
US11741574B2 (en) | Method and device for processing blood vessel image on basis of user input | |
EP2152151A2 (en) | Method, apparatus and system for detection of arterial stiffness and artery tonus by pulse curve geometry analysis | |
Cao et al. | Automatic identification of side branch and main vascular measurements in intravascular optical coherence tomography images | |
CN113724315A (zh) | 眼底视网膜血管宽度测量方法、电子设备及计算机可读存储介质 | |
Lu et al. | Automatic detection of the coronary arterial contours with sub-branches from an x-ray angiogram | |
CN110349203B (zh) | 一种超声设备及其超声图像中血管直径的测量方法 | |
CN106796725A (zh) | 一种血管脊线追踪方法及装置 | |
CN109875595B (zh) | 一种颅内血管状态检测方法及装置 | |
CN109919913B (zh) | 一种冠状动脉的半径计算方法、终端及存储介质 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
PB01 | Publication | ||
PB01 | Publication | ||
SE01 | Entry into force of request for substantive examination | ||
SE01 | Entry into force of request for substantive examination | ||
CB02 | Change of applicant information |
Address after: 100102 No. 501 No. 12, 5th floor, No. 6, Wangjing Dongyuan District 4, Chaoyang District, Beijing Applicant after: Shukun (Beijing) Network Technology Co.,Ltd. Address before: 100102 No. 501 No. 12, 5th floor, No. 6, Wangjing Dongyuan District 4, Chaoyang District, Beijing Applicant before: SHUKUN (BEIJING) NETWORK TECHNOLOGY Co.,Ltd. |
|
CB02 | Change of applicant information | ||
RJ01 | Rejection of invention patent application after publication |
Application publication date: 20190517 |
|
RJ01 | Rejection of invention patent application after publication |