1.基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,包括以下步驟:
步驟1:獲得被檢測(cè)對(duì)象的PET/CT掃描的體模原始數(shù)據(jù)和CT圖像數(shù)據(jù);
步驟2:對(duì)原始數(shù)據(jù)進(jìn)行門控校正處理,得到新的PET掃描數(shù)據(jù),即門控?cái)?shù)據(jù);
步驟3:分析處理軟件讀入CT圖像數(shù)據(jù),根據(jù)CT圖像得到CT衰減系數(shù),將CT衰減系數(shù)轉(zhuǎn)化為PET能量下的衰減校正系數(shù),求出各響應(yīng)線上的衰減校正因子;
步驟4:將門控?cái)?shù)據(jù)導(dǎo)入到分析處理軟件中,并將門控?cái)?shù)據(jù)與得到的衰減校正因子相乘,從而得到新的PET掃描數(shù)據(jù),即衰減校正后的數(shù)據(jù);
步驟5:對(duì)進(jìn)行門控校正和衰減校正后的數(shù)據(jù)進(jìn)行重建,得到清晰的PET圖像。
2.根據(jù)權(quán)利要求1所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟2中的門控校正的具體方法如下:
步驟a1:將原始數(shù)據(jù)以固定間隔時(shí)間劃分為N幀;
步驟a2:在每一幀掃描數(shù)據(jù)中,統(tǒng)計(jì)每一幀各個(gè)探測(cè)環(huán)中的真光子數(shù);
步驟a3:選取其中一幀為參考幀,分別與其他所有幀進(jìn)行比較,通過絕對(duì)誤差公式計(jì)算其他幀探測(cè)環(huán)真光子數(shù)與參考幀探測(cè)環(huán)真光子數(shù)的絕對(duì)誤差;
步驟a4:設(shè)定誤差閾值,通過所有的絕對(duì)誤差值與設(shè)定閾值的比較,將絕對(duì)誤差值小于設(shè)定閾值的所有幀的數(shù)據(jù)與參考幀的數(shù)據(jù)進(jìn)行整合,得到新的PET掃描數(shù)據(jù),即門控?cái)?shù)據(jù)。
3.根據(jù)權(quán)利要求2所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟a1中的N為:
N=fix(Y/T)
式中:fix為朝零方向取整;
Y為PET掃描的總時(shí)間;
T為固定間隔時(shí)間。
4.根據(jù)權(quán)利要求2所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟a3中的絕對(duì)誤差公式為:
式中DRTCE(M,I)為第M幀第I探測(cè)環(huán)真光子數(shù),DRTCE(N,I)為第N幀第I探測(cè)環(huán)真光子數(shù),I為PET探測(cè)器探測(cè)環(huán)編號(hào),M、N為幀的編號(hào)。
5.根據(jù)權(quán)利要求1所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟3的衰減校正因子的算法如下:
ACF表示衰減校正因子;
a表示探測(cè)環(huán)上探測(cè)器A、探測(cè)器B之間的距離;
P點(diǎn)表示A點(diǎn)到B點(diǎn)這條響應(yīng)線上的湮滅點(diǎn);
μ(x)表示這條響應(yīng)線上離P點(diǎn)x距離處的衰減校正系數(shù)。
6.根據(jù)權(quán)利要求1所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟4的將門控?cái)?shù)據(jù)導(dǎo)入到分析處理軟件中的具體方法為:
通過ROOT軟件輸出門控校正后的數(shù)據(jù),輸出的形式為列表模型數(shù)據(jù),將列表模型的文本數(shù)據(jù)進(jìn)行分割,將分割后的數(shù)據(jù)導(dǎo)入分析處理軟件中,再將導(dǎo)入的數(shù)據(jù)進(jìn)行整合得到原數(shù)據(jù)。
7.根據(jù)權(quán)利要求1所述的基于PET/CT成像的呼吸運(yùn)動(dòng)門控校正和衰減校正方法,其特征在于,所述步驟5中的重建采用MLEM迭代算法。