颅内血管血流延迟的影像标注方法及系统技术方案

本发明专利技术公开了一种颅内血管血流延迟的影像标注方法，包括以下步骤：步骤1，血管分割及延迟时间的计算；步骤2，计算得到血管显影的三维视频；步骤3，量化；利用动态灌注数据分割出低灌注区和梗死核心区，找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。本发明专利技术基于磁共振或计算机断层扫描的动态脑灌注影像数据，按照造影剂抵达的时间延迟，对动态灌注数据中的颅内血管进行定量标注，对延迟时间进行彩色方式编码，并与结构影像进行融合显示，所得到的图像可以作为中间结果，为医生对侧支循环的建立程度以及建立通路进行评估和判读提供参考。本发明专利技术还公开了一种颅内血管血流延迟的影像标注系统。

Image Labeling Method and System for Delayed Intracranial Blood Flow

The invention discloses an image labeling method for intracranial blood flow delay, which includes the following steps: step 1, calculation of blood vessel segmentation and delay time; step 2, calculation of three-dimensional video of blood vessel development; step 3, quantification; using dynamic perfusion data to segment the low perfusion area and infarction core area, finding the part located in the low perfusion area for statistics, and obtaining the low perfusion area. Quantitative value of blood flow delay. The invention is based on dynamic cerebral perfusion image data of magnetic resonance or computed tomography. According to the time delay of arrival of contrast agent, the intracranial blood vessels in dynamic perfusion data are quantitatively labeled, the time delay is coded in color mode, and fused with the structural image to display. The obtained image can be used as an intermediate result to establish the collateral circulation for doctors. It also provides reference for the establishment of pathways for evaluation and interpretation. The invention also discloses an image labeling system for intracranial blood flow delay.

【技术实现步骤摘要】
颅内血管血流延迟的影像标注方法及系统
本专利技术涉及一种灌注成像数据处理技术，具体涉及一种颅内血管血流延迟的影像标注方法。本专利技术还涉及一种颅内血管血流延迟的影像标注系统。
技术介绍
脑动脉的急性或慢性闭塞使得部分脑组织得不到充足的血液供应而坏死，从而形成缺血型卒中。闭塞发生后，脑内血管网侧支开放，经由willis环、眼动脉、软脑膜吻合支或新生血管代偿，为缺血部位供给血液。定量评估这些侧支循环的建立情况，对静脉溶栓、血管内治疗、动静脉桥接治疗的患者预后具有显著的预测价值。目前，急性卒中对侧支循环的影像学检查的金标准是有创血管成像DSA(Digitalsubtractionangiography，血管造影的影像数字化处理)，其需要侵入患者血管，注射造影剂；同时将DSA的C型臂调整至各观察角度以追踪观察血管内造影剂的动态充盈情况，并录制视频。而进行有创DSA检查前，一般会首先进行计算机断层成像(ComputedTomography，CT)扫描，包括无造影剂增强的平扫、无创断层扫描血管成像(ComputedTomographyAngiography，CTA)、以及动态脑灌注成像。现有的对于侧支循环建立情况的无创评估方法，如CTA-ASPECTS、ASTIN/SIR、Miteff评分等，均为基于CTA的定性评分，不仅需要耗费人力进行判读，而且结果还会受判读医生主观影响，评估结论不够客观准确。另外还可以利用动态脑灌注成像进行间接判断，即注射对比剂后，使用CT或磁共振对脑部进行动态连续扫描(CTP/MR-PWI)。这种成像方法可以用于计算脑血流量(CBF)、脑血容量(CBV)、平均通过时间(MTT)、达峰时间(TTP)、最大增强时间(Tmax)等脑血流动力学参数，对脑组织的血流供给情况进行量化从而间接推测侧支循环的建立。如美国专利技术专利文献US6650928B1所公开的基于动态灌注数据计算出“脑组织”血流供给的参数图，以rainbow、three-color、hot-iron或Puh-Thallium伪彩编码方式与结构影像进行融合叠加。然而该血供分析结果无法描述供血来源，无法确定血液经由血管网络的何种通路抵达缺血区域。而且对于超急性期或部分慢性闭塞的缺血性卒中病人，采用上述血供分析并不能观察到脑血供异常。另外，还有少量医疗单位选择了多时态CTA扫描并按照美国专利技术专利文献US20170265829A1所公开的方法进行处理。即分多个时间点扫描CTA，并按照不同时态的血管显影以彩色进行标注。这种方法由于也需要在注射造影剂后的各时间点进行成像，因而不能和动态脑灌注成像同时存在。其较低的时间分辨率无法准确计算动态脑灌注成像所得到的各参数图，从而难以反应脑组织的缺血状态。
技术实现思路
本专利技术所要解决的技术问题是通过重新利用无创检查中的动态脑灌注成像数据进行再次处理，从而得到与有创DSA检查相仿的效果，生成血管的三维的充盈过程视频，并且其可以在计算机上进行多角度的观察，在有创DSA检查前，就能得到颅内血管侧支循环的评估参考。为解决上述技术问题，本专利技术颅内血管血流延迟的影像标注方法的技术解决方案为，基于脑部的动态灌注数据，得到颅内血管的血流延迟时间；包括以下步骤：步骤1，血管分割及延迟时间的计算；步骤1.1，血管分割；找到血管体素；步骤1.2，计算所有血管体素的延迟时间；在另一实施例中，所述步骤1.2计算延迟时间的方法为：根据血管内位于处体素的信号时间曲线得到其最大峰值的峰位,即以作为血管内处体素的延迟时间；或者，计算与动脉输入曲线的相对时间，以大脑闭塞侧中动脉根部体素的时间信号曲线作为动脉输入曲线SA(t),计算其达到最大值时的时间LA＝argmaxt(SA(t))，再以作为血管内处体素的延迟时间；其中步骤2，计算得到血管显影的三维视频；根据步骤1得到的血管体素的不同延迟时间，对不同延迟时间点的血管进行三维建模，并将血管的所有延迟时间点的三维模型拼接为一个三维血管的视频，该视频含有表明造影剂充盈各血管显影的先后顺序的信息；步骤3，量化；利用动态灌注数据分割出低灌注区和梗死核心区，再根据步骤1中得到的已知血流延迟时间的血管，找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。在另一实施例中，对所述步骤3中所得到的低灌注区和梗死核心区进行三维建模，并标注在步骤2所得到的三维视频模型上。在另一实施例中，对所述步骤2所得到的三维视频中的三维血管模型进行颜色标注。进一步地，所述进行颜色标注的方法为：首先计算血管内位于处体素的时间信号曲线与动脉输入曲线SA(t)的反卷积结果的最大值，以作为血管内处的颜色标注根据。在另一实施例中，所述步骤3包括以下步骤：步骤3.1，根据动态灌注数据的参数图得到低灌注区和梗死核心区；步骤3.2，对步骤1中得到的已知延迟时间的血管找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。进一步地，所述步骤3.2中的低灌注区内血流延迟的量化值包括具有血流延迟时间的平均值、中位数、最大值、最小值以及血管体积。在另一实施例中，所述步骤1之前执行步骤0，读入磁共振或计算机断层扫描的动态脑灌注影像数据，对动态灌注数据进行预处理。在另一实施例中，所述步骤0中的预处理方法包括：逐时间点进行运动矫正、背景去除、图像去噪。本专利技术还提供一种颅内血管血流延迟的影像标注系统，其技术解决方案为：包括血管分割及延迟时间计算模块、三维建模模块和量化模块；所述血管分割及延迟时间计算模块包括血管分割单元、延迟时间计算单元：血管分割单元，被配置为利用动态灌注数据中的血管体素的时间信号曲线峰值高、峰宽窄两个特点，找到动态灌注数据中的血管体素；延迟时间计算单元，被配置为计算所述血管分割单元找到的所有血管体素的延迟时间；所述三维建模模块，被配置为根据延迟时间计算单元得到的血管体素的不同延迟时间，对不同延迟时间点的血管进行三维建模，并将所有延迟时间点组合成视频，得到一个血管的三维视频，该视频含有表明造影剂充盈各血管显影的先后顺序的信息；所述量化模块，被配置为利用动态灌注数据分割出低灌注区和梗死核心区，再根据延迟时间计算单元所得到的已知血流延迟时间的血管，找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。在另一实施例中，所述量化模块包括低灌注区获取单元、梗死核心区获取单元、统计单元；低灌注区获取单元，被配置为根据动态灌注数据的参数图得到低灌注区；梗死核心区获取单元，被配置为动态灌注数据的参数图得到梗死核心区；统计单元，被配置为根据所述延迟时间计算单元所得到的已知延迟时间的血管找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。在另一实施例中，还包括三维模型颜色标注单元，被配置为对所述三维建模模块所得到的三维视频中的三维血管模型进行颜色标注。本专利技术可以达到的技术效果是：本专利技术提供了磁共振(MRI)或计算机断层扫描(CT)的动态脑灌注数据的另一种处理方式，即利用动态脑灌注数据时造影剂抵达血管的动态信息对血管三维模型进行各时间点的重建，从而得到与有创检查DSA相仿的结果，在计算机上还可以在任意截面进行观察，能够直观地为判断颅内血管的闭塞程度、血管网中血液流向、侧支循环代偿通路、血管膨胀情况提供判断依据。而且，使用灌本文档来自技高网...

【技术保护点】
1.一种颅内血管血流延迟的影像标注方法，其特征在于，基于脑部的动态灌注数据，得到颅内血管的血流延迟时间；包括以下步骤：步骤1，血管分割及延迟时间的计算；步骤1.1，血管分割；找到血管体素；步骤1.2，计算所有血管体素的延迟时间；步骤2，计算得到血管显影的三维视频；根据步骤1得到的血管体素的不同延迟时间，对不同延迟时间点的血管进行三维建模，并将血管的所有延迟时间点的三维模型拼接为一个三维血管的视频，该视频含有表明造影剂充盈各血管显影的先后顺序的信息；步骤3，量化；利用动态灌注数据分割出低灌注区和梗死核心区，再根据步骤1中得到的已知血流延迟时间的血管，找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。

【技术特征摘要】
1.一种颅内血管血流延迟的影像标注方法，其特征在于，基于脑部的动态灌注数据，得到颅内血管的血流延迟时间；包括以下步骤：步骤1，血管分割及延迟时间的计算；步骤1.1，血管分割；找到血管体素；步骤1.2，计算所有血管体素的延迟时间；步骤2，计算得到血管显影的三维视频；根据步骤1得到的血管体素的不同延迟时间，对不同延迟时间点的血管进行三维建模，并将血管的所有延迟时间点的三维模型拼接为一个三维血管的视频，该视频含有表明造影剂充盈各血管显影的先后顺序的信息；步骤3，量化；利用动态灌注数据分割出低灌注区和梗死核心区，再根据步骤1中得到的已知血流延迟时间的血管，找到位于低灌注区的部分进行统计，得到低灌注区内血流延迟的量化值。2.根据权利要求1所述的颅内血管血流延迟的影像标注方法，其特征在于：对所述步骤3中所得到的低灌注区和梗死核心区进行三维建模，并标注在步骤2所得到的三维视频模型上。3.根据权利要求1所述的颅内血管血流延迟的影像标注方法，其特征在于：对所述步骤2所得到的三维视频中的三维血管模型进行颜色标注。4.根据权利要求3所述的颅内血管血流延迟的影像标注方法，其特征在于：所述进行颜色标注的方法为：首先计算血管内位于处体素的时间信号曲线与动脉输入曲线SA(t)的反卷积结果的最大值，以作为血管内处的颜色标注根据。5.根据权利要求1所述的颅内血管血流延迟的影像标注方法，其特征在于：所述步骤1.2计算延迟时间的方法为：根据血管内位于处体素的信号时间曲线得到其最大峰值的峰位,即以作为血管内处体素的延迟时间；或者，计算与动脉输入曲线的相对时间，以大脑闭塞侧中动脉根部体素的时间信号曲线作为动脉输入曲线SA(t),计算其达到最大值时的时间LA＝argmaxt(SA(t))，再以作为血管内处体素的延迟时间；其中6.根据权利要求1所述的颅内血管血流延迟的影像标注方法，其特征在于：所述步骤3包括以下步骤：步骤3.1，根据动态灌注数据的参数图得到低...

【专利技术属性】
技术研发人员：杨鸣，李鑫，丁晶，汪昕，方文星，
申请(专利权)人：脑玺上海智能科技有限公司，
类型：发明
国别省市：上海,31