[發(fā)明專利]一種數(shù)字PET的時間校正方法及系統(tǒng)

在審

• 申請?zhí)枺?/td>: 201510065331.8
• 申請日: 2015-02-09
• 公開（公告）號: CN104644204A
• 公開（公告）日: 2015-05-27
• 發(fā)明（設計）人: 陳源寶;陳忻;呂云飛;佘亞軍
• 申請（專利權(quán)）人: 中國船舶重工集團公司第七一九研究所
• 主分類號: A61B6/03
• 分類號: A61B6/03
• 代理公司: 北京理工大學專利中心 11120
• 代理人: 仇蕾安;李愛英

1.一種數(shù)字PET的時間校正方法，其特征在于，包括：

步驟一、獲取各探測器模組生成的列表式數(shù)據(jù)，并對所述列表式數(shù)據(jù)進行能量篩選，剔除散射符合事件，形成用于時間校正的原始數(shù)據(jù)集；

步驟二、利用所述原始數(shù)據(jù)集進行時間符合判斷，計算得到每條響應線兩端探測單元的通道時延差Δt_i,j，同時獲得每條響應線上的符合事件數(shù)B_i,j；其中，i和j分別為響應線兩端探測單元的編號，i和j的取值范圍為1～n，n為探測單元總數(shù)，且i≠j；所述探測單元為模組級探測單元或晶體條級探測單元；

接著，利用計算得到的響應線兩端探測單元的通道時延差Δt_i,j，構(gòu)建探測單元的通道時延方程組：

其中，t_i為第i個探測單元的通道時延量；

將式(1)表示為式(2)：

即，A·x＝Δt_LOR?????????????(2)?

步驟三、利用步驟二獲得的每條響應線上的符合事件數(shù)B_i,j，計算每條響應線上的符合事件數(shù)占所述原始數(shù)據(jù)集總符合事件數(shù)的權(quán)重權(quán)重W_i,j越大即表征相應的響應線對時間分辨率的貢獻越大；

步驟四、從矩陣A中優(yōu)選出n個行向量，優(yōu)選準則為：1、該n個行向量可以組成矩陣A的一個極大線性無關(guān)組，2、該n個行向量對應的權(quán)重值加和最大；?從而獲得對時間分辨率貢獻最大的極大線性無關(guān)組A^*，實現(xiàn)所述探測單元的通道時延方程組的降維運算,獲得降維后的探測單元的通道時延方程組A^*·x＝Δt_LOR^*；利用降維后的探測單元的通道時延方程組進行計算，獲得探測單元的通道時延量；

步驟五、利用步驟四獲得的探測單元的通道時延量，對所述原始數(shù)據(jù)集進行時間校正。

2.如權(quán)利要求1所述的方法，其特征在于，所述列表式數(shù)據(jù)由探測器模組生成，該列表式數(shù)據(jù)包含所探測到的γ光子的時間信息、能量信息、位置信息、晶體單元編號、探測器模塊編號和探測器模組編號。

3.如權(quán)利要求1所述的方法，其特征在于，步驟四中，所述利用降維后的探測單元的通道時延方程組進行計算，獲得探測單元的通道時延量的具體步驟為：

求解降維后的探測單元的通道時延方程組，得到一組通解，從通解中選擇某一探測單元的通道時延量為參考基準探測單元，求解通道時延方程組的特解，即獲得探測器中其它探測單元相對于參考基準探測單元的通道時延量；其中，所述參考基準探測單元選擇為對時間分辨率貢獻最大的探測單元。