[發明專利]多層容積CT灌注成像源影像的自動化減影處理方法

有效

• 申請號: 201210592536.8
• 申請日: 2012-12-29
• 公開（公告）號: CN103908276A
• 公開（公告）日: 2014-07-09
• 發明（設計）人: 楊朝輝;楊秀軍;汪雪林;佟慶彬;雷大偉
• 申請（專利權）人: 北京握奇數據系統有限公司
• 主分類號: A61B6/03
• 分類號: A61B6/03
• 代理公司: 北京天悅專利代理事務所(普通合伙) 11311
• 代理人: 田明;任曉航

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。