[发明专利]多层容积CT灌注成像源影像的自动化减影处理方法

权利要求书

1.一种多层容积CT灌注成像源影像的自动化减影处理方法，包括以下步骤：

（1）进行容积层CT灌注成像扫描，获取容积成像数据，根据扫描时间和容积成像数据计算CT成像设备的球管每转一周所扫描的层数n；所述的容积成像数据包括扫描得到的CT灌注影像序列中的DICOM图像；

（2）确定扫描时间中每一扫描时刻容积成像数据的n幅图像,获取每一扫描时刻n幅图像的容积CT值集合；

（3）确定CT减影的被减时刻，根据被减时刻的容积CT值集合得到其它扫描时刻的减影容积CT值集合；

（4）根据减影容积CT值集合与其对应的文件元信息集合，获取容积成像数据的减影DICOM文件。

2.如权利要求1所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（1）中，计算CT成像设备的球管每转一周所扫描的层数n的具体方式为：

（1-1）计算CT灌注影像序列中所有图像的图像平面与原点的最小距离，并存储在一个一维数组中；

（1-2）删除所述一维数组中的重复的元素，数组中剩余元素的个数为CT成像设备的球管每转一周所扫描的层数n。

3.如权利要求2所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（2）中，获取每一时刻所有图像的容积CT值集合的具体方式为：

（2-1）根据所获取的CT灌注影像序列中DICOM图像的图像序列号Instance Nunber确定某一扫描时刻timei容积成像数据的n幅图像；

（2-2）根据n幅图的图像序列号Instance Nunber的顺序依次读入对应图像的像素数值，并转换成对应图像的CT值，得到n个一维CT值数组；

（2-3）将所得到的n个一维CT值数组按照其对应图像的Instance Nunber串接得到timei时刻的容积CT集合。

4.如权利要求3所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（3）中，所述的其它扫描时刻的减影容积CT值集合包括被减时刻之前的扫描时刻的减影容积CT值集合和被减时刻之后的扫描时刻的减影容积CT值集合。

5.如权利要求4所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（3）中，获取被减时刻之后的扫描时刻的减影容积CT值集合的具体方式为：

（3-1）确定被减时刻subtractedi，获取subtractedi时刻的容积CT集合HUset0；

（3-2）获取被减时刻subtractedi之后的某一扫描时刻t的容积容积CT集合HUset；

（3-3）计算某一扫描时刻t的减影容积CT值集合deltaHUset，计算公式为:deltaHUset=HUset-HUset0。

6.如权利要求5所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（3）中，被减时刻subtractedi之前的所有扫描时刻的减影容积CT值集合等于subtractedi+1时刻的减影容积CT集合。

7.如权利要求6所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（3）中，如果减影容积CT值集合中某一元素的数值小于0，将该元素的数值置为0。

8.如权利要求7所述的一种多层容积CT灌注成像源影像的自动化减影处理方法，其特征在于：步骤（4）中，获取容积成像数据的减影DICOM文件的具体方式为：

假设需要获取t时刻的容积成像数据中第j幅图像的DICOM文件，首先获取t时刻的减影容积CT值集合deltaHUset，对集合deltaHUset中的第j个元素deltaHUset[j]进行线性变换得到第j幅图像的像素值，根据转换得到的第j幅图像的像素值与t时刻容积成像数据中第j幅图像的DICOM文件元信息得到subtractedi时刻容积成像数据第j幅图像的DICOM文件；其中，1≤j≤n。