[發(fā)明專利]一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法

有效

• 申請?zhí)枺?/td>: 201810507458.4
• 申請日: 2018-05-24
• 公開（公告）號: CN109523458B
• 公開（公告）日: 2022-12-02
• 發(fā)明（設(shè)計(jì)）人: 余維;王成祥;聶小英;曾德輝;湛慧苗;陳淑文
• 申請（專利權(quán)）人: 湖北科技學(xué)院
• 主分類號: G06T3/00
• 分類號: G06T3/00;G06T5/00
• 代理公司: 浙江千克知識產(chǎn)權(quán)代理有限公司 33246
• 代理人: 裴金華

1.一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于，包括如下步驟：

S1：通過利用濾波反投影方法對采集的投影數(shù)據(jù)進(jìn)行解析重建得到的圖像作為初始化圖像并進(jìn)行POCS（凸集投影）得到濾波輸入圖像；

S2：通過對所述濾波輸入圖像進(jìn)行圖像梯度稀疏正則化約束得到引導(dǎo)圖像；

S3：將所述引導(dǎo)圖像和所述濾波輸入圖像融合得到引導(dǎo)濾波圖像；

S4：將所述引導(dǎo)濾波圖像作為初始化圖像進(jìn)行迭代重建。

2.如權(quán)利要求1所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于，步驟S1之前包括：將X射線進(jìn)行掃描得到投影數(shù)據(jù)；采集所述投影數(shù)據(jù)。

3.如權(quán)利要求1或2所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于，步驟S1中具體包括：

S11：通過利用濾波反投影方法對采集的投影數(shù)據(jù)進(jìn)行解析重建得到的圖像作為初始化圖像并進(jìn)行數(shù)據(jù)一致性得到第一處理圖像；

S12：通過對第一處理圖像進(jìn)行非負(fù)約束得到濾波輸入圖像。

4.如權(quán)利要求3所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S11中數(shù)據(jù)一致性利用公式實(shí)現(xiàn)，其中是正松弛參數(shù)，是迭代指數(shù)。

5.如權(quán)利要求3所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S12中非負(fù)約束利用公式得到濾波輸入圖像。

6.如權(quán)利要求1或2所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于，步驟S2中具體包括：

S21：通過對濾波輸入圖像進(jìn)行圖像梯度稀疏約束得到正則化圖像；

S22：將得到的正則化圖像作為引導(dǎo)圖像。

7.如權(quán)利要求6所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S22中根據(jù)公式和，表示軟閾值濾波的參數(shù)，通過濾波輸入圖像計(jì)算引導(dǎo)圖像，其中，。

8.如權(quán)利要求3所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S12篩選出第一處理圖像中像素灰度值小于零的值置為零，確保濾波輸入圖像像素灰度值均大于或等于零。

9.如權(quán)利要求1或2所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S3中引導(dǎo)濾波圖像通過公式求出，其中，，，，是窗口中的像素總數(shù)，是濾波輸入圖像的平均值，是引導(dǎo)圖像的平均值，是引導(dǎo)圖像的方差。

10.如權(quán)利要求1或2所述的一種結(jié)合稀疏誘導(dǎo)動態(tài)引導(dǎo)濾波的高精度稀疏角度CT重建方法，其特征在于：步驟S4中設(shè)有預(yù)設(shè)重復(fù)步驟次數(shù)，當(dāng)重復(fù)的步驟次數(shù)超過預(yù)設(shè)重復(fù)步驟次數(shù)時(shí)則不論引導(dǎo)濾波圖像是否滿足預(yù)設(shè)終止值均輸出引導(dǎo)濾波圖像。