CN114209344A - 一种侧枝循环状态评估方法、装置、存储介质及电子设备 - Google Patents
一种侧枝循环状态评估方法、装置、存储介质及电子设备 Download PDFInfo
- Publication number
- CN114209344A CN114209344A CN202111308504.6A CN202111308504A CN114209344A CN 114209344 A CN114209344 A CN 114209344A CN 202111308504 A CN202111308504 A CN 202111308504A CN 114209344 A CN114209344 A CN 114209344A
- Authority
- CN
- China
- Prior art keywords
- blood vessel
- target
- collateral
- image
- blood
- 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
- 206010069729 Collateral circulation Diseases 0.000 title claims abstract description 70
- 238000011156 evaluation Methods 0.000 title claims abstract description 67
- 210000004204 blood vessel Anatomy 0.000 claims abstract description 177
- 238000000034 method Methods 0.000 claims abstract description 43
- 230000036770 blood supply Effects 0.000 claims description 67
- 238000003384 imaging method Methods 0.000 claims description 18
- 230000010412 perfusion Effects 0.000 claims description 16
- 238000005192 partition Methods 0.000 claims description 14
- 238000012545 processing Methods 0.000 claims description 11
- 238000004422 calculation algorithm Methods 0.000 claims description 9
- 238000004590 computer program Methods 0.000 claims description 9
- 238000000605 extraction Methods 0.000 claims description 6
- 238000001914 filtration Methods 0.000 claims description 6
- 238000007781 pre-processing Methods 0.000 claims description 4
- 238000001514 detection method Methods 0.000 claims 1
- 230000002490 cerebral effect Effects 0.000 description 28
- 210000001367 artery Anatomy 0.000 description 25
- 238000010968 computed tomography angiography Methods 0.000 description 7
- 230000004087 circulation Effects 0.000 description 6
- 238000005516 engineering process Methods 0.000 description 4
- 239000011159 matrix material Substances 0.000 description 4
- 238000012986 modification Methods 0.000 description 4
- 230000004048 modification Effects 0.000 description 4
- 210000002551 anterior cerebral artery Anatomy 0.000 description 3
- 210000001841 basilar artery Anatomy 0.000 description 3
- 230000017531 blood circulation Effects 0.000 description 3
- 210000004556 brain Anatomy 0.000 description 3
- 238000011161 development Methods 0.000 description 3
- 201000010099 disease Diseases 0.000 description 3
- 208000037265 diseases, disorders, signs and symptoms Diseases 0.000 description 3
- 210000003657 middle cerebral artery Anatomy 0.000 description 3
- 210000003388 posterior cerebral artery Anatomy 0.000 description 3
- 210000002385 vertebral artery Anatomy 0.000 description 3
- 208000032382 Ischaemic stroke Diseases 0.000 description 2
- 230000001154 acute effect Effects 0.000 description 2
- 230000002308 calcification Effects 0.000 description 2
- 235000019994 cava Nutrition 0.000 description 2
- 210000001638 cerebellum Anatomy 0.000 description 2
- 239000002872 contrast media Substances 0.000 description 2
- 238000009499 grossing Methods 0.000 description 2
- 238000013507 mapping Methods 0.000 description 2
- 210000003625 skull Anatomy 0.000 description 2
- 238000002560 therapeutic procedure Methods 0.000 description 2
- 210000003462 vein Anatomy 0.000 description 2
- 206010061751 Cerebrovascular stenosis Diseases 0.000 description 1
- 206010061216 Infarction Diseases 0.000 description 1
- 208000006011 Stroke Diseases 0.000 description 1
- 208000007536 Thrombosis Diseases 0.000 description 1
- 238000002583 angiography Methods 0.000 description 1
- 230000008321 arterial blood flow Effects 0.000 description 1
- 238000004364 calculation method Methods 0.000 description 1
- 210000001627 cerebral artery Anatomy 0.000 description 1
- 230000003727 cerebral blood flow Effects 0.000 description 1
- 208000026106 cerebrovascular disease Diseases 0.000 description 1
- 230000001447 compensatory effect Effects 0.000 description 1
- 238000002591 computed tomography Methods 0.000 description 1
- 238000010276 construction Methods 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000003111 delayed effect Effects 0.000 description 1
- 230000006870 function Effects 0.000 description 1
- 230000000004 hemodynamic effect Effects 0.000 description 1
- 230000007574 infarction Effects 0.000 description 1
- 230000000302 ischemic effect Effects 0.000 description 1
- 230000007774 longterm Effects 0.000 description 1
- 238000004393 prognosis Methods 0.000 description 1
- 230000010410 reperfusion Effects 0.000 description 1
- 230000000472 traumatic effect Effects 0.000 description 1
Images
Classifications
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/50—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications
- A61B6/507—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment specially adapted for specific body parts; specially adapted for specific clinical applications for determination of haemodynamic parameters, e.g. perfusion CT
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/02—Arrangements for diagnosis sequentially in different planes; Stereoscopic radiation diagnosis
- A61B6/03—Computed tomography [CT]
- A61B6/032—Transmission computed tomography [CT]
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/48—Diagnostic techniques
- A61B6/481—Diagnostic techniques involving the use of contrast agents
-
- A—HUMAN NECESSITIES
- A61—MEDICAL OR VETERINARY SCIENCE; HYGIENE
- A61B—DIAGNOSIS; SURGERY; IDENTIFICATION
- A61B6/00—Apparatus or devices for radiation diagnosis; Apparatus or devices for radiation diagnosis combined with radiation therapy equipment
- A61B6/52—Devices using data or image processing specially adapted for radiation diagnosis
- A61B6/5211—Devices using data or image processing specially adapted for radiation diagnosis involving processing of medical diagnostic data
Landscapes
- Health & Medical Sciences (AREA)
- Life Sciences & Earth Sciences (AREA)
- Engineering & Computer Science (AREA)
- Medical Informatics (AREA)
- Biomedical Technology (AREA)
- Surgery (AREA)
- High Energy & Nuclear Physics (AREA)
- Physics & Mathematics (AREA)
- Nuclear Medicine, Radiotherapy & Molecular Imaging (AREA)
- Optics & Photonics (AREA)
- Pathology (AREA)
- Radiology & Medical Imaging (AREA)
- Veterinary Medicine (AREA)
- Heart & Thoracic Surgery (AREA)
- Molecular Biology (AREA)
- Biophysics (AREA)
- Animal Behavior & Ethology (AREA)
- General Health & Medical Sciences (AREA)
- Public Health (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Pulmonology (AREA)
- Theoretical Computer Science (AREA)
- Dentistry (AREA)
- Oral & Maxillofacial Surgery (AREA)
- Apparatus For Radiation Diagnosis (AREA)
Abstract
本申请公开了一种侧枝循环状态评估方法、装置、存储介质以及电子设备，其中方法包括：获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。本申请中的侧枝循环状态评估方法能够使得评估结果更加准确，评估时间较短、效率较高。
Description
技术领域
本申请涉及医疗技术领域，特别涉及一种侧枝循环状态评估方法、装置、存储介质及电子设备。
背景技术
急性缺血性脑卒中作为全球首位致残及第二位致死性疾病，已成为我国病死率第一的疾病，其中急性缺血性脑卒中约占全部卒中的80％。随着血管内治疗技术的不断发展，针对具体患者建立个体化评估治疗方案，从而改善患者临床结局和预后已经成为亟待解决的问题。
脑侧枝循环是指脑血管狭窄或闭塞导致的供血区域缺血时，血流通过其他血管进行代偿性供血。侧支循环良好与将血管内血栓取出的手术过程和预测疾病发展情况密切相关，有效地侧支循环状态可以稳定梗死区脑血流量，延长治疗时间窗，改善再灌注治疗的远期发展情况。
现有一些侧枝状态评估方法包括基于有创影像技术，例如数字减影血管造影(Digital Subtraction Angiography，简称：DSA)，和基于无创影像技术，例如计算机断层血管造影(Computed tomography angiography,CTA)和计算机扫描灌注(Computedtomography perfusion,CTP)等，无创影像技术对首次影像检查期间作为符合血管内治疗患者的IA类被推荐。但是，目前，基于无创影像技术的侧枝评估方法存在不准确或者评估时间较长的问题，不适用急性期患者。
由此亟需一种侧枝循环状态评估方法，使得大脑侧枝循环评估更加准确及时。
发明内容
有鉴于此，本发明提供了一种侧枝循环状态评估方法、装置、存储介质及电子设备，主要目的在于解决目前侧枝循环状态评估不准确的问题。
为解决上述问题，本申请提供一种侧枝循环状态评估方法，包括：
获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；
提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定相应血管掩膜图像中患侧的充盈度，以获得若干充盈度；
至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
可选的，所述获取包含感兴趣区域的多个目标图像，具体包括：
基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描图像；
对各所述扫描图像进行预处理，获得与各所述扫描图像对应的第一图像；
对各所述第一图像进行血管增强处理，以获得各侧枝血管特征增强的各所述目标图像。
可选的，所述提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，具体包括：
利用最大密度投影算法对各所述目标图像分别进行处理，获得处理后的各目标图像；
采用非局部均值滤波的方式对处理后的各所述目标图像进行处理，获得各所述血管掩膜图像。
可选的，所述分别基于同一个血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度，具有包括：
将各所述血管掩膜图像与标准供血区分区模板图像进行匹配，以确定各所述血管掩膜图像中正常侧的第一目标供血区以及与所述第一目标供血区对应的第二目标供血区；所述第二目标供血区为患侧的供血区；
基于各所述血管掩膜图像中正常侧的第一目标供血区的灰度值以及所述第二目标供血区的灰度值，确定所述第二目标供血区的充盈度，以获得与所述各所述血管掩膜图像对应的所述若干充盈度。
可选的，所述至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果，具体包括：
基于各所述充盈度确定目标充盈度；
至少基于所述目标充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
可选的，所述侧枝循环状态评估方法，还包括：
基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值以及患侧侧枝血管的灰度值，确定患侧的侧枝血管充盈的延迟时间；
基于所述延迟时间对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
可选的，基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值以及患侧侧枝血管的灰度值，确定患侧的侧枝血管充盈的延迟时间，具体包括：
基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值，确定充盈度符合第一预设条件的至少一个第一血管掩膜图像；
基于各所述血管掩膜图像中患侧侧枝血管的灰度值，确定充盈度符合第二预设条件的至少一个第二血管掩膜图像；
基于各所述第一血管掩膜图像对应的扫描时间以及各所述第二血管掩膜图像对应扫描时间确定延迟时间。
为解决上述技术问题，本申请提供一种侧枝循环状态评估装置，包括：
获取模块，获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征
提取模块，提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
获得模块，用于分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
评估模块，用于至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
为解决上述技术问题，本申请提供一种存储介质，所述存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述任一项所述侧枝循环状态评估方法的步骤。
为解决上述技术问题，本申请提供一种电子设备，至少包括存储器、处理器，所述存储器上存储有计算机程序，所述处理器在执行所述存储器上的计算机程序时实现上述任一项侧枝循环状态评估方法的步骤。
本申请基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描数据(各期CTP数据)；本申请中评估方法区别于以往的任何侧支评估方法，通过确定侧枝循环充盈度，基于充盈度来确定侧枝循环状态为良好状态或者非良好状态，能够使得评估结果更加准确，评估时间较短、效率较高。为后续将该评估结果提供给临床医生作为参考奠定了基础。
上述说明仅是本发明技术方案的概述，为了能够更清楚了解本发明的技术手段，而可依照说明书的内容予以实施，并且为了让本发明的上述和其它目的、特征和优点能够更明显易懂，以下特举本发明的具体实施方式。
附图说明
通过阅读下文优选实施方式的详细描述，各种其他的优点和益处对于本领域普通技术人员将变得清楚明了。附图仅用于示出优选实施方式的目的，而并不认为是对本发明的限制。而且在整个附图中，用相同的参考符号表示相同的部件。在附图中：
图1为本申请实施例一种侧枝循环状态评估方法的流程图；
图2为本申请又一实施例侧枝循环状态评估方法的流程图。
具体实施方式
此处参考附图描述本申请的各种方案以及特征。
应理解的是，可以对此处申请的实施例做出各种修改。因此，上述说明书不应该视为限制，而仅是作为实施例的范例。本领域的技术人员将想到在本申请的范围和精神内的其他修改。
包含在说明书中并构成说明书的一部分的附图示出了本申请的实施例，并且与上面给出的对本申请的大致描述以及下面给出的对实施例的详细描述一起用于解释本申请的原理。
通过下面参照附图对给定为非限制性实例的实施例的优选形式的描述，本申请的这些和其它特性将会变得显而易见。
还应当理解，尽管已经参照一些具体实例对本申请进行了描述，但本领域技术人员能够确定地实现本申请的很多其它等效形式。
当结合附图时，鉴于以下详细说明，本申请的上述和其他方面、特征和优势将变得更为显而易见。
此后参照附图描述本申请的具体实施例；然而，应当理解，所申请的实施例仅仅是本申请的实例，其可采用多种方式实施。熟知和/或重复的功能和结构并未详细描述以避免不必要或多余的细节使得本申请模糊不清。因此，本文所申请的具体的结构性和功能性细节并非意在限定，而是仅仅作为权利要求的基础和代表性基础用于教导本领域技术人员以实质上任意合适的详细结构多样地使用本申请。
本说明书可使用词组“在一种实施例中”、“在另一个实施例中”、“在又一实施例中”或“在其他实施例中”，其均可指代根据本申请的相同或不同实施例中的一个或多个。
本申请实施例提供一种侧枝循环状态评估方法，具体可以应用大脑侧枝循环状态的评估。如图1所示，本实施例中的评估方法包括如下步骤：
步骤S101，获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；
本步骤中在获取多个目标图像之前，可以先基于CT灌注成像进行多次扫描，以此来获得与各次扫描对应的扫描数据，即获得若干组/期数的CTP数据，由此获得若干组/期数扫描图像。在获取扫描图像后，可以对扫描数据进行处理来获得目标图像，例如可采用SPM12软件对大脑前联合、后联合进行矫正，头骨剥离，配准到MNI空间，多尺度平滑滤波对各扫描图像进行预处理获得与各所述扫描图像对应的第一图像；可采用海森矩阵相关的血管增强算法对各第一图像进行血管增强处理以获得各侧枝血管特征增强的各所述目标图像，在获得目标图像后还可以利用造影剂在血管中的流动特征，进一步采用矩阵补全算法分离钙化背景来进行图像提取，以获取最终的各目标图像，所述目标图像中包括正常侧和患侧的侧枝血管特征。
步骤S102，提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
本步骤中，利用最大密度投影算法沿着轴向对各所述目标图像进行处理，便于在图像上呈现血管形态和灰度分布，获得与之对应的血管掩膜图像；对预定数期CTP数据中的每一次扫描数据/各期扫描数据重复上述步骤，分别获得对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征。
步骤S103，分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定相应血管掩膜图像中患侧的充盈度，以获得若干充盈度；
本步骤中的正常侧和患侧是通过比较各动脉血管分区中血管灰度均值差异得到的，灰度均值低的一侧为患侧，另一侧为正常侧，本步骤中，充盈度具体是指患侧血管血流量灰度值占正常侧血管血流量灰度值的百分比。
步骤S104，至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
本步骤中，根据充盈度对侧枝循环状态进行评估，例如当充盈度大于预设值时说明侧枝循环状态良好，当充盈度小于预设值时则说明侧支循环状态不良。
本实施例中的方法，通过基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描数据(各期CTP数据)，然后通过确定侧枝循环的充盈度，基于充盈度来确定侧枝循环状态为良好状态或者非良好状态，能够使得评估结果更加准确，评估时间较短、效率较高。
本申请在上述实施例的基础上，为了使得对侧枝循环状态的评估更加准确，本申请又一实施例提供一种侧枝循环状态评估方法，如图2所示，包括如下步骤：
步骤S201，基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描图像；
在具体实施过程中，可以根据实际需要设定扫描次数。例如设定扫描次数为19次，以此来获得19期CTP扫描数据/扫描图像，在获取扫描图像时，还可以进一步结合不同医院不同设备的扫描方案，来获得与各扫描次数对应的扫描图像。预定的扫描次数可以覆盖CT灌注成像扫描的整个周期，CT灌注成像是对感兴趣区层面进行连续CT扫描,从而获得感兴趣区时间-密度曲线，扫描的整个周期可以根据感兴趣区时间-密度曲线确定，以脑侧枝血管评估为例，通过采集CT灌注成像完整周期的数据，通过重建图像获取不同期图像的正常侧和患侧的血管的灰度特征，以获取血管的充盈度，而不用人为的评经验去选择动脉期，此后每间隔8S分别选取静脉期、静脉晚期数据，提升了数据的准确性。
步骤S202，对各所述扫描图像进行预处理，获得预处理后的第一图像；
本步骤可以使用SPM12软件对各所述扫描图像进行大脑前联合和大脑后联合矫正，剥离头骨，配准到MNI空间后使用多尺度平滑滤波去除干扰信息获得预处理后的第一图像。
步骤S203，对各所述第一图像进行血管增强处理，以获得各侧枝血管特征增强的各所述目标图像。
本步骤中可采用海森矩阵相关的血管增强算法对各第一图像进行血管增强处理以获得各侧枝血管特征增强的各所述目标图像，在获得目标图像后还可以利用造影剂在血管中的流动特征，进一步采用矩阵补全算法分离钙化背景来进行图像提取，以获取最终的各目标图像，所述目标图像中包括正常侧和患侧的侧枝血管特征。
步骤S204，提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
本步骤在具体实施过程中，可以沿着轴向空间维度进行灰度的最大密度投影，采用非局部均值滤波的方式对投影后的各所述目标图像进行处理，获得各所述血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
步骤S205，分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
本步骤中将各所述血管掩膜图像与标准供血区分区模板图像进行匹配，以确定各所述血管掩膜图像中正常侧的第一目标供血区以及第二目标供血区；所述第二目标供血区为患侧的供血区；基于各所述血管掩膜图像中正常侧的第一目标供血区的灰度值以及与所述第一目标供血区对应的第二目标供血区的灰度值，确定所述第二目标供血区的充盈度，以获得所述若干充盈度。所述标准供血区分区模板具体可以为大脑的分区模板，该模板中是包含有小脑后下动脉、小脑上动脉、椎动脉和基底动脉分支、脉络膜前动脉、豆纹动脉、大脑前动脉、大脑中动脉、大脑后动脉的各个标准分区的模板，即该分区模板中包含有各供血区的位置信息与供血区名称的对应关系。当然该标准供血区分区模板还可以为其他部位对应的分区模板。当患侧出现在其他动脉供血区时，标准供血区分区模板为与之对应的供血区分区模板。
本步骤在实施过程中，供血区的确定方式具体可以为将各所述血管掩膜图像与预定的标准供血区分区模板图像进行匹配，以确定各所述血管掩膜图像中正常侧的第一目标供血区以及与正常侧对应的患侧的第二目标供血区；所述第二供血区与所述第一供血区一一对应；基于各所述血管掩膜图像中正常侧的第一目标供血区的灰度值以及与所述第一目标供血区对应的第二目标供血区的灰度值，确定所述第二目标供血区的充盈度，以获得所述若干期侧枝的充盈度。本步骤中以患侧为脑动脉为例进行说明，所述第一目标供血区可以为如下任意一种供血区：小脑后下动脉、小脑上动脉、椎动脉和基底动脉分支、脉络膜前动脉、豆纹动脉、大脑前动脉、大脑中动脉、大脑后动脉；第二目标供血区可以为如下任意一种供血区：小脑后下动脉、小脑上动脉、椎动脉和基底动脉分支、脉络膜前动脉、豆纹动脉、大脑前动脉、大脑中动脉、大脑后动脉。在确定充盈度时，可以针对某一组/对供血区来确定患侧的充盈度；也可以同时分别基于多组/对供血区来确定患侧的充盈度。例如当第一目标供血区为小脑后下动脉时，第二目标供血区为小脑后下动脉；则后将正常侧的小脑后下动脉的灰度值与患侧小脑后下动脉的灰度值进行比较，以确定各所述血管掩膜图像中患侧小脑后下动脉的充盈度，由此获得若干充盈度。同理，当第一目标供血区为小脑后下动脉和小脑上动脉时，第二目标供血区为小脑后下动脉和小脑上动脉，则后续将正常侧的小脑后下动脉的灰度值与患侧小脑后下动脉的灰度值进行比较，将正常侧的小脑后上动脉的灰度值与患侧小脑后上动脉的灰度值进行比较，以确定各所述血管掩膜图像中患侧小脑后下动脉的充盈度，以及各所述血管掩膜图像中患侧小脑后上动脉的充盈度。后续可以根据各小脑后下动脉的充盈度，确定第一目标充盈度，根据各小脑后上动脉的充盈度确定第二目标充盈度，最后根据第一目标充盈度和第二目标充盈度确定侧支循环状态。第一目标供血区的灰度值和第二目标供血区的灰度值可以是对应供血区域内的平均灰度值，也可以是对应供血区域内某一点的灰度值。优选地，第一目标供血区的灰度值和第二目标供血区的灰度值是对应供血区域内的平均灰度值。
由于不同的血流量在图像中灰度值也会不同，灰度值越高代表血流量越大，因此本步骤中可以通过确定患侧第二目标供血区灰度值占正常侧第一目标供血区灰度值的百分比，根据百分比来确定患侧相对于正常侧的供血区的充盈度，使得充盈度的确定更加准确。
步骤S206，基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值以及患侧侧枝血管的灰度值，确定患侧的侧枝血管充盈的延迟时间；
本步骤在具体实施过程中，具体基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值，确定充盈度符合第一预设条件的至少一个第一血管掩膜图像；基于各所述血管掩膜图像中患侧侧枝血管的灰度值，确定充盈度符合第二预设条件的至少一个第二血管掩膜图像；基于各第一血管掩膜图像对应的扫描时间以及各所述第二血管掩膜图像对应扫描时间确定延迟时间，延迟时间等于第二血管掩膜图像对应的扫描时间与第一血管掩膜图像对应扫描时间的差值。所述的延迟时间(Time delay)具体还可以表示为：Time delay＝|Plm-Prm|*Tx。Tx时间间隔，由数据扫描方案中规定的P_lm、P_rm之间时间差确定，Plm为充盈度符合第一预设条件的一个第一血管掩膜图像的期数，Prm为充盈度符合第二预设条件的一个第二血管掩膜图像的期数。需要说明的是，上述Tx为等时间间隔。在实际数据扫描过程中，相邻期的时间间隔也可以不相同，上述公式是Tx为等时间间隔时的示例性表达。在本实施例中，第一目标供血区的灰度值和第二目标供血区的灰度值可以是动脉血管的灰度平均值。
步骤S207，基于所述充盈度和延迟时间对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
本步骤在具体实施过程中，首先基于各所述充盈度确定目标充盈度：所述目标充盈度包括：基于各所述充盈度计算获取的平均充盈度，或满足预设条件的充盈度。即可以根据各充盈度计算获得以平均值，以此来作为目标充盈度。也可以通过计算方差的方式来从各充盈度中挑选出至少一个充盈度作为目标充盈度，即将充盈度中误差较大的充盈度剔除掉，由此能够使得确定目标充盈度更加合理准确、为后续基于目标充盈度准确的进行侧枝循环状态的评估奠定了基础。在确定出目标充盈度之后，就可以采用预定的充盈度与分值的映射关系确定该目标充盈度所对应的分值，获得第一分值，同理利用预定的延迟时间与分值的映射关系确定该延迟时间所对应的分值，获得第二分值，最后将第一分值和第二分值进行加权来确定目标分值，根据目标分值所对应的分值区间确定侧枝循环状态，即不同的分值区间对应不同的侧枝循环状态。
本实施例中，通过确定各血管掩膜图像中各第二目标供血区(位于患侧的供血区)的充盈度，进而能够更加全面、准确的确定出各第二目标供血区的目标充盈度，使得后续基于目标充盈度对侧支循环状态的评估更加准确。并且通过根据各血管掩膜图像中患侧与正常侧供血区的灰度值，来确定供血区的延迟时间，同时基于延迟时间和充盈度来评估侧枝循环状态，能使得最终获得的评估结果更加精准。
相关技术中提出了两种基于非创伤影响技术的侧枝状态评估方法，其中，一种侧枝状态评估方法是通过比较单期的CTA图像中患侧和正常侧侧支血管的充盈范围给出侧支得分，另一种评估方法是利用患者的动脉期、静脉期和静脉晚期数据，不仅比较了患侧和正常侧血管的充盈范围，而且给出了患侧侧支充盈延迟的时相。该评分方案有效避免了单期CTA数据低估侧支循环的缺点。但三期CTA数据采集对于设备和技师的要求都比较高且扫描时间相对较长，不适用于急性期患者。在一些实施方式中，也有根据时间密度曲线选取动脉血流达峰值时对应的一期CTP数据作为动脉期，此后每间隔8S分别选取静脉期、静脉晚期数据。然而，由于CTP各期数据扫描间隔是固定值，根据灌注曲线选择动脉期数据后直接按照固定时间间隔选取静脉期和静脉晚期数据这一做法存在数据选择不准确的问题。本申请上述实施例提出的侧支状态评估方法利用完整的CTP数据，省去了从CTP数据中手动选取动脉期、静脉期和静脉晚期数据这一复杂耗时的过程，提升侧枝状态评估的准确性和效率。此外，完整的CTP数据包含完整的血流动力学信息，且相比三期CTA数据在临床上容易获得，时间信息也更加丰富，相比单期CTA增加了时间信息。
本申请另一实施例提供一种侧枝循环状态评估装置，包括：
获取模块，获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征
提取模块，提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
获得模块，用于分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
评估模块，用于至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
本实施例中，所述侧枝循环状态评估装置还包括用于获取各大脑血管图像的获取模块，所述获取模块具体用于：基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描图像；对各所述扫描图像进行预处理，获得预处理后的第一图像；使用血管增强算法对所述预处理后的第一图像进行血管增强处理，以获取特征增强的目标图像。
所述提取模块具体用于：利用最大密度投影算法对各目标图像分别进行提取；采用非局部均值滤波的方式对各所述目标图像进行处理，获得各所述血管掩膜图像。
所述获得模块具体用于：分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度，具有包括：将各所述血管掩膜图像与标准供血区分区模板图像进行匹配，以确定各所述血管掩膜图像中正常侧的第一目标供血区以及第二目标供血区；所述第二目标供血区为患侧的供血区；基于各所述血管掩膜图像中正常侧的第一目标供血区的灰度值以及与所述第一目标供血区对应的第二目标供血区的灰度值，确定所述第二目标供血区的充盈度，以获得所述若干充盈度。
所述评估模块具体用于：至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果，具体包括：基于各所述充盈度确定目标充盈度，所述目标充盈度包括：基于各所述充盈度计算获取的平均充盈度，或满足预设条件的充盈度；至少基于所述目标充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
具体的，所述评估模块具体用于：基于各所述血管掩膜图像中第一目标供血区的灰度值以及第二目标供血区的灰度值，确定第二目标供血区的延迟时间；所述至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果，具体包括：基于各所述充盈度以及所述延迟时间对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
具体的实施过程中，所述评估模块用于，所述基于各所述血管掩膜图像中第一目标供血区的灰度值以及第二目标供血区的灰度值，确定第二目标供血区的延迟时间，具体包括：基于各所述血管掩膜图像中正常侧侧枝血管的灰度值，确定充盈度符合第一预设条件的至少一个第一血管掩膜图像；基于各所述血管掩膜图像中患侧侧枝血管的灰度值，确定充盈度符合第二预设条件的至少一个第二血管掩膜图像；基于各所述第一血管掩膜图像对应的扫描时间以及各所述第二血管掩膜图像对应扫描时间确定延迟时间。
本实施例中的方法，通过基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描数据(各期CTP数据)，然后通过确定侧枝循环的充盈度和延迟时间，基于充盈度和延迟时间来确定侧枝循环状态为良好状态或者非良好状态，能够使得评估结果更加准确，评估时间较短、效率较高。
本申请又一实施例提供一种存储介质，所述存储介质存储有计算机程序，所述计算机程序被处理器执行时实现如下方法步骤：
步骤一、获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征
步骤二、提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
步骤三、分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
步骤四、至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
上述方法步骤的具体实施例过程可参见上述任意侧枝循环状态评估方法的实施例，本实施例在此不再重复赘述。
本实施例中的方法，通过基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描数据(各期CTP数据)，然后通过确定侧枝循环的充盈度，基于充盈度来确定侧枝循环状态为良好状态或者非良好状态，能够使得评估结果更加准确，评估时间较短、效率较高。
本申请又一实施例提供一种电子设备，至少包括存储器、处理器，所述存储器上存储有计算机程序，所述处理器在执行所述存储器上的计算机程序时实现如下方法步骤：
步骤一、获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；
步骤二、提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；；
步骤三、用于分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
步骤四、用于至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
上述方法步骤的具体实施过程可参见上述侧枝循环状态评估方法的实施例，本实施例在此不再重复赘述。
本实施例中的方法，通过基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描数据(各期CTP数据)，然后通过确定侧枝循环的充盈度，基于充盈度来确定侧枝循环状态为良好状态或者非良好状态，能够使得评估结果更加准确，评估时间较短、效率较高。
以上实施例仅为本申请的示例性实施例，不用于限制本申请，本申请的保护范围由权利要求书限定。本领域技术人员可以在本申请的实质和保护范围内，对本申请做出各种修改或等同替换，这种修改或等同替换也应视为落在本申请的保护范围内。
Claims (10)
1.一种侧枝循环状态评估方法，其特征在于，包括：
获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；
提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
2.如权利要求1所述的方法，其特征在于，所述获取包含感兴趣区域的多个目标图像，具体包括：
基于预定的扫描次数进行CT灌注成像扫描，获得与各扫描次数对应的扫描图像；
对各所述扫描图像进行预处理，获得与各所述扫描图像对应的第一图像；
对各所述第一图像进行血管增强处理，以获得各侧枝血管特征增强的各所述目标图像。
3.如权利要求1所述的方法，其特征在于，所述提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，具体包括：
利用最大密度投影算法对各所述目标图像分别进行处理，获得处理后的各目标图像；
采用非局部均值滤波的方式对处理后的各所述目标图像进行处理，获得各所述血管掩膜图像。
4.如权利要求1所述的方法，其特征在于，所述分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度，具有包括：
将各所述血管掩膜图像与标准供血区分区模板图像进行匹配，以确定各所述血管掩膜图像中正常侧的第一目标供血区以及与所述第一目标供血区对应的第二目标供血区，所述第二目标供血区为患侧的供血区；
基于各所述血管掩膜图像中正常侧的第一目标供血区的灰度值以及所述第二目标供血区的灰度值，确定所述第二目标供血区的充盈度，以获得与所述各所述血管掩膜图像对应的所述若干充盈度。
5.如权利要求1所述的方法，所述至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果，具体包括：
基于各所述充盈度确定目标充盈度；所述目标充盈度包括：基于各所述充盈度计算获取的平均充盈度，或满足预设条件的充盈度；
至少基于所述目标充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
6.如权利要求1所述的方法，其特征在于，所述方法还包括：
基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值以及患侧侧枝血管的灰度值，确定患侧的侧枝血管充盈的延迟时间；
基于所述延迟时间对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
7.如权利要求6所述的方法，其特征在于，所述基于各所述血管掩膜图像中正常侧的侧枝血管的灰度值以及患侧侧枝血管的灰度值，确定患侧的侧枝血管充盈的延迟时间，具体包括：
基于各所述血管掩膜图像中正常侧侧枝血管的灰度值，确定充盈度符合第一预设条件的至少一个第一血管掩膜图像；
基于各所述血管掩膜图像中患侧侧枝血管的灰度值，确定充盈度符合第二预设条件的至少一个第二血管掩膜图像；
基于各所述第一血管掩膜图像对应的扫描时间以及各所述第二血管掩膜图像对应扫描时间确定延迟时间。
8.一种侧枝循环状态评估装置，其特征在于，包括：
获取模块，获取包含感兴趣区域的多个目标图像，各所述目标图像中包括正常侧和患侧的侧枝血管特征；
提取模块，提取各所述目标图像的侧枝血管特征，获得与各所述目标图像对应的血管掩膜图像，各所述血管掩膜图像包括正常侧和患侧的侧枝血管的灰度特征；
获得模块，用于分别基于同一血管掩膜图像中正常侧的灰度值以及患侧的灰度值，确定患侧的充盈度，以获得若干充盈度；
评估模块，用于至少基于各所述充盈度对侧枝循环状态进行评估，得到侧枝循环状态的评估结果。
9.一种存储介质，其特征在于，所述存储介质存储有计算机程序，所述计算机程序被处理器执行时实现上述权利要求1-7任一项所述侧枝循环状态评估方法的步骤。
10.一种电子设备，其特征在于，至少包括存储器、处理器，所述存储器上存储有计算机程序，所述处理器在执行所述存储器上的计算机程序时实现上述权利要求1-7任一项目标对象检测方法的步骤。
Priority Applications (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111308504.6A CN114209344A (zh) | 2021-11-05 | 2021-11-05 | 一种侧枝循环状态评估方法、装置、存储介质及电子设备 |
Applications Claiming Priority (1)
| Application Number | Priority Date | Filing Date | Title |
|---|---|---|---|
| CN202111308504.6A CN114209344A (zh) | 2021-11-05 | 2021-11-05 | 一种侧枝循环状态评估方法、装置、存储介质及电子设备 |
Publications (1)
| Publication Number | Publication Date |
|---|---|
| CN114209344A true CN114209344A (zh) | 2022-03-22 |
Family
ID=80696570
Family Applications (1)
| Application Number | Title | Priority Date | Filing Date |
|---|---|---|---|
| CN202111308504.6A Pending CN114209344A (zh) | 2021-11-05 | 2021-11-05 | 一种侧枝循环状态评估方法、装置、存储介质及电子设备 |
Country Status (1)
| Country | Link |
|---|---|
| CN (1) | CN114209344A (zh) |
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117291281A (zh) * | 2023-08-29 | 2023-12-26 | 强联智创（北京）科技有限公司 | 评估脑部侧支循环的评估模型进行训练的方法及相关产品 |
-
2021
- 2021-11-05 CN CN202111308504.6A patent/CN114209344A/zh active Pending
Cited By (1)
| Publication number | Priority date | Publication date | Assignee | Title |
|---|---|---|---|---|
| CN117291281A (zh) * | 2023-08-29 | 2023-12-26 | 强联智创（北京）科技有限公司 | 评估脑部侧支循环的评估模型进行训练的方法及相关产品 |
Similar Documents
| Publication | Publication Date | Title |
|---|---|---|
| US11615894B2 (en) | Diagnostically useful results in real time | |
| US11406337B2 (en) | Calculating a fractional flow reserve | |
| US10470730B2 (en) | Creating a vascular tree model | |
| CN106682636B (zh) | 血管提取方法及其系统 | |
| JP6542129B2 (ja) | リアルタイムの診断上有用な結果 | |
| EP3061015A2 (en) | Vascular characteristic determination with correspondence modeling of a vascular tree | |
| EP3537977B1 (en) | Method and system for modelling a human heart and atria | |
| CN117036530B (zh) | 基于跨模态数据的冠状动脉血流储备分数预测方法及装置 | |
| CN114209344A (zh) | 一种侧枝循环状态评估方法、装置、存储介质及电子设备 | |
| EP3658031B1 (en) | Motion compensated cardiac valve reconstruction | |
| EP4107697B1 (en) | Methods and systems for cardiac chamber imaging | |
| CN113538403A (zh) | 一种虚拟dsa影像获取装置、方法、介质及电子设备 |
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 |