[发明专利]获取脑部CT灌注成像参数图的方法、装置、系统和计算机存储介质

说明书

本申请公开了获取脑部CT灌注成像参数图的方法、装置、系统和计算机存储介质，所述方法包括：获取CT灌注影像、对所述CT灌注影像进行前处理，获得离散的组织密度变化曲线C(n)；读取CT灌注影像的采集时间信息，获得采集时刻数组T(n)；截取采集时刻数组T(n)，得到相对采集时刻数组t(n)；将组织密度变化曲线C(n)与相对应的相对采集时刻数组t(n)结合，得到各像素点的组织密度时间曲线C(t_n)；对离散的组织密度时间曲线C(t_n)进行拟合或插值后，按照相同的时间间隔重新离散，得到离散后的组织密度时间曲线C(n)’。对动脉输入变化曲线AIF(n)和静脉输出变化曲线VOF(n)进行同样的处理。本申请提高了组织密度时间曲线的可用性，有助于降低求解难度、提高求解的准确度。

技术领域

本申请涉及图像处理和医疗工程交叉技术领域，特别是涉及一种获取脑部 CT灌注成像参数图的方法、装置、系统和计算机存储介质。

背景技术

脑卒中，即中风，是一种急性脑血管疾病，是由于脑部血管突然破裂或因 血管阻塞导致血液不能流入大脑而引起脑组织损伤的一组疾病，包括缺血性和 出血性卒中，具有极高的致残率和较高的致死率。人类生活结构的改变导致急 性缺血性脑卒中发病率越来越高。急性缺血性脑卒中(脑梗死)是最常见的脑卒 中类型,占全部脑卒中的60％～80％。

灌注成像在急性脑卒中和肿瘤学中应用最为广泛。卒中诊断时，灌注成像 的目的是确定病变组织的范围，并描绘出可以再灌注的缺血组织。灌注测量越 来越多地用于评估、分期和治疗后监测。

脑CT灌注成像(CTP)是一种评价脑实质血流灌注状态的成像技术，也 是当前检查急性缺血性脑卒中的重要影像学方法。这项技术包括在静脉注射造 影剂(如碘化造影剂)后对头部进行快速的连续CT扫描。基于示踪动力学模 型及示踪剂稀释理论，通过基于去卷积或非去卷积的算法对获取的对比增强和 时间分辨CT图像进行处理，以生成灌注参数(如脑血容量(cerebral blood volume, CBV)、脑血流量(cerebral blood flow,CBF)、平均传输时间(mean transient time, MTT)等)图。

现有技术在CT灌注影像采集的过程中，由于影像采集设备或采集方式的 不同，同一影像序列的不同帧的影像可能具有不同的采集时间。去卷积方法是 目前最常用的CT灌注图像处理方法，参阅公式(1)：c(t)＝AIF(t)*k(t)； 其中，C(t)为组织时间密度曲线；AIF(t)为动脉输入函数；k(t)为残余函数，或 称流量标度残余函数。由于不同帧的影像具有不同的采集时间，所以直接采用 采集的原始影像计算残余函数k(t)，将导致后续灌注参数图的准确性下降、求 解难度困难、残余函数k(t)时间分辨率低等问题。

发明内容

为解决现有技术中灌注参数图的准确性低、求解灌注参数图的难度大、残 余函数k(t)时间分辨率低的问题，本申请提供一种获取脑部CT灌注成像参数图 的方法。

本申请获取脑部CT灌注成像参数图的方法，包括：

获取CT灌注影像、对所述CT灌注影像进行前处理，获得离散的组织密度 变化曲线C(n)，其中n为第n个采集时刻，C(n)为第n个时刻对应的组织密度；

读取CT灌注影像的采集时间信息，获得采集时刻数组T(n)，其中n为第n 个时刻，T(n)为采集第n个影像时，影像采集设备的物理时刻，T(n)包括与离散 的组织密度变化曲线C(n)相对应的Tc(n)；

截取采集时刻数组T(n)，得到相对采集时刻数组t(n)，t(n)包括与离散的组 织密度变化曲线C(n)相对应的tc(n)；

将组织密度变化曲线C(n)与相对应的相对采集时刻数组t(n)结合，得到各 像素点的组织密度时间曲线C(t_n)；