白丝美女被狂躁免费视频网站,500av导航大全精品,yw.193.cnc爆乳尤物未满,97se亚洲综合色区,аⅴ天堂中文在线网官网

首頁 / 專利庫 / 數(shù)學(xué)與統(tǒng)計(jì) / 濾波反投影 / 一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法

一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法

閱讀:762發(fā)布:2020-05-13

專利匯可以提供一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法專利檢索,專利查詢,專利分析的服務(wù)。并且一種基于 卷積神經(jīng)網(wǎng)絡(luò) 的CT圖像重建方法,包括對原始弦圖數(shù)據(jù)IK進(jìn)行反投影操作,得到反投影圖像數(shù)據(jù)I’K,對反投影圖像數(shù)據(jù)I’K進(jìn)行歸一化處理,得到歸一化反投影圖像數(shù)據(jù)PK,將歸一化反投影圖像PK通過卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行卷積神經(jīng)網(wǎng)絡(luò)濾波,生成待處理圖像P’K,對待處理圖像P’K進(jìn)行反歸一化處理,得到最終的重建圖像Pfinal。不需要對圖像濾波進(jìn)行特別的設(shè)計(jì),通過對卷積神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練自動(dòng)完成圖像濾波的學(xué)習(xí),本 發(fā)明 重建方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和偽影,同時(shí)較好地保持原有圖像的 分辨率 ,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。,下面是一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法專利的具體信息內(nèi)容。

1.一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:包括如下步驟:
A1、對原始弦圖數(shù)據(jù)IK進(jìn)行反投影操作,得到反投影圖像數(shù)據(jù)I’K;
A2、對反投影圖像數(shù)據(jù)I’K進(jìn)行歸一化處理,得到歸一化反投影圖像數(shù)據(jù)PK;
A3、將歸一化反投影圖像PK通過卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行卷積神經(jīng)網(wǎng)絡(luò)濾波,生成待處理圖像P’K;
A4、對待處理圖像P’K進(jìn)行反歸一化處理,得到最終的重建圖像Pfinal。
2.根據(jù)權(quán)利要求1所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A1中反投影操作是通過CT掃描機(jī)對原始弦圖數(shù)據(jù)IK進(jìn)行幾何成像處理。
3.根據(jù)權(quán)利要求2所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A2中歸一化處理的方法步驟如下:
T1、計(jì)算反投影圖像數(shù)據(jù)I’K的均值XI’和方差SI’;
T2、根據(jù)式(1)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK:
PK=(I’K-XI’)/SI’······式(1)。
4.根據(jù)權(quán)利要求3所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A3中濾波的具體操作如下:
設(shè)第s層卷積層有ns卷積核,第s層的輸出特征圖記為Fs,第s層的輸入特征圖記為Fs-1,輸出特征圖Fs的第i個(gè)通道記為Fis-1(i=1,2,...,ns),輸入特征圖Fs-1的第j個(gè)通道記為Fjs-1(j=1,2,...,ns),n為自然數(shù);
卷積神經(jīng)網(wǎng)絡(luò)的輸入記為F0,卷積神經(jīng)網(wǎng)絡(luò)中第s層卷積層輸出Fs與輸入Fs-1之間的關(guān)系如下:
式(2)中,δ(·)表示輸入矩陣中元素的非線性激活函數(shù),*表示二維卷積操作, 為第s層卷積層中輸出特征圖第i個(gè)通道與輸入特征圖中第j個(gè)通道進(jìn)行卷積的二維卷積核, 為第s層卷積層中用于計(jì)算輸出特征圖第i個(gè)通道的偏置項(xiàng),S表示卷積網(wǎng)絡(luò)中的卷積層的數(shù)目,S為自然數(shù)。
5.根據(jù)權(quán)利要求4所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A3將步驟A2中得到的歸一化反投影圖像數(shù)據(jù)PK作為卷積神經(jīng)網(wǎng)絡(luò)的輸入F0。
6.根據(jù)權(quán)利要求5所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A4中的反歸一化處理為步驟A2的逆操作。
7.根據(jù)權(quán)利要求6所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:步驟A4中的反歸一化處理的具體方法如下:
M1、計(jì)算反待處理圖像P’K的均值XP’和方差SP’;
M2、根據(jù)式(3)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK:
Pfinal=P’K*SP’+XP’······式(3)。
8.根據(jù)權(quán)利要求7所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:所述非線性激活函數(shù)δ(·)為sigmoid函數(shù),根據(jù)式(4)計(jì)算:
δ(x)=1/(1+e-x),e為自然底數(shù),······式(4)。
9.根據(jù)權(quán)利要求7所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:所述非線性激活函數(shù)δ(·)為雙曲線正切函數(shù),根據(jù)式(5)計(jì)算:
δ(x)=(ex-e-x)/(ex+e-x),其中e為自然底數(shù),······式(5)。
10.根據(jù)權(quán)利要求7所述的一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其特征在于:所述非線性激活函數(shù)δ(·)為ReLU函數(shù),根據(jù)式(6)計(jì)算:

說明書全文

一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法

技術(shù)領(lǐng)域

[0001] 本發(fā)明涉及醫(yī)學(xué)影像的圖像重建方法技術(shù)領(lǐng)域,特別是涉及一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法。

背景技術(shù)

[0002] X射線CT已經(jīng)廣泛應(yīng)用于臨床醫(yī)學(xué)影像診斷,但CT掃描中過高的X射線輻射劑量存在致癌風(fēng)險(xiǎn)。如何最大限度地降低X射線使用劑量已經(jīng)成為醫(yī)學(xué)CT成像領(lǐng)域研究的關(guān)鍵技術(shù)。
[0003] 當(dāng)前,降低CT掃描過程中的管電流和掃描時(shí)間是實(shí)現(xiàn)低劑量CT成像的最簡便且常用的途徑。然而,由于降低管電流和掃描時(shí)間,使得投影數(shù)據(jù)中含有大量的噪聲,基于傳統(tǒng)的濾波反投影方法重建的圖像質(zhì)量存在嚴(yán)重的退化現(xiàn)象,難以滿足臨床診斷需要。為了在保證圖像質(zhì)量的前提下大幅降低X射線輻射劑量,諸多基于降低管電流和掃描時(shí)間的低劑量CT圖像重建方法相繼提出,例如基于統(tǒng)計(jì)模型的迭代重建方法和基于投影數(shù)據(jù)濾波的解析重建方法。其中,基于統(tǒng)計(jì)模型的迭代重建方法,通過對采集的投影數(shù)據(jù)的噪聲以及成像系統(tǒng)進(jìn)行圖像重建模型構(gòu)建,可以實(shí)現(xiàn)低劑量CT圖像優(yōu)質(zhì)重建;基于投影數(shù)據(jù)濾波的解析重建方法,同樣通過對采集的投影數(shù)據(jù)的噪聲以及成像系統(tǒng)進(jìn)行數(shù)據(jù)濾波建模,再通過解析重建方法實(shí)現(xiàn)快速且優(yōu)質(zhì)的低劑量CT圖像重建。
[0004] 基于統(tǒng)計(jì)模型的迭代重建方法的缺點(diǎn)及導(dǎo)致原因:在重建同樣大小的CT圖像時(shí),基于統(tǒng)計(jì)模型的迭代重建方法所花費(fèi)的時(shí)間遠(yuǎn)遠(yuǎn)大于傳統(tǒng)的解析重建方法,不能滿足臨床的所需要的CT實(shí)時(shí)顯像要求,其原因在于基于統(tǒng)計(jì)模型的迭代重建方法,需要對目標(biāo)函數(shù)進(jìn)行幾十甚至上百次的反復(fù)迭代求解,導(dǎo)致圖像重建時(shí)間大幅增加。
[0005] 基于投影數(shù)據(jù)濾波的解析重建方法的缺點(diǎn)及導(dǎo)致原因:傳統(tǒng)的基于投影數(shù)據(jù)濾波的解析重建方法,在投影數(shù)據(jù)降噪過程中不可避免地會(huì)導(dǎo)致圖像原有細(xì)節(jié)信息的丟失,從而導(dǎo)致相應(yīng)CT圖像分辨率的下降。。
[0006] 因此,針對現(xiàn)有技術(shù)不足,提供一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法以解決現(xiàn)有技術(shù)不足甚為必要。

發(fā)明內(nèi)容

[0007] 本發(fā)明的目的在于避免現(xiàn)有技術(shù)的不足之處而提供一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,該CT圖像重建方法基于級聯(lián)方式X射線CT圖像重建方法,處理方法操作簡單,處理方便,圖像精度高。
[0008] 本發(fā)明的上述目的通過如下技術(shù)手段實(shí)現(xiàn)。
[0009] 提供一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,包括如下步驟:
[0010] A1、對原始弦圖數(shù)據(jù)IK進(jìn)行反投影操作,得到反投影圖像數(shù)據(jù)I’K;
[0011] A2、對反投影圖像數(shù)據(jù)I’K進(jìn)行歸一化處理,得到歸一化反投影圖像數(shù)據(jù)PK;
[0012] A3、將歸一化反投影圖像PK通過卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行卷積神經(jīng)網(wǎng)絡(luò)濾波,生成待處理圖像P’K;
[0013] A4、對待處理圖像P’K進(jìn)行反歸一化處理,得到最終的重建圖像Pfinal。
[0014] 具體而言的,步驟A1中反投影操作是通過CT掃描機(jī)對原始弦圖數(shù)據(jù)IK進(jìn)行幾何成像處理。
[0015] 優(yōu)選的,步驟A2中歸一化處理的方法步驟如下:
[0016] T1、計(jì)算反投影圖像數(shù)據(jù)I’K的均值XI’和方差SI’;
[0017] T2、根據(jù)式(1)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK:
[0018] PK=(I’K-XI’)/SI’······式(1)。
[0019] 具體而言的,步驟A3中濾波的具體操作如下:
[0020] 設(shè)第s層卷積層有ns卷積核,第s層的輸出特征圖記為Fs,第s層的輸入特征圖記為Fs-1,輸出特征圖Fs的第i個(gè)通道記為 輸入特征圖Fs-1的第j個(gè)通道記為 n為自然數(shù);
[0021] 卷積神經(jīng)網(wǎng)絡(luò)的輸入記為F0,卷積神經(jīng)網(wǎng)絡(luò)中第s層卷積層輸出Fs與輸入Fs-1之間的關(guān)系如下:
[0022]
[0023] 式(2)中,δ(·)表示輸入矩陣中元素的非線性激活函數(shù),*表示二維卷積操作,為第s層卷積層中輸出特征圖第i個(gè)通道與輸入特征圖中第j個(gè)通道進(jìn)行卷積的二維卷積核, 為第s層卷積層中用于計(jì)算輸出特征圖第i個(gè)通道的偏置項(xiàng),S表示卷積網(wǎng)絡(luò)中的卷積層的數(shù)目,S為自然數(shù)。
[0024] 進(jìn)一步的,步驟A3將步驟A2中得到的歸一化反投影圖像數(shù)據(jù)PK作為卷積神經(jīng)網(wǎng)絡(luò)的輸入F0。
[0025] 進(jìn)一步的,步驟A4中的反歸一化處理為步驟A2的逆操作。
[0026] 優(yōu)選的,步驟A4中的反歸一化處理的具體方法如下:
[0027] M1、計(jì)算反待處理圖像P’K的均值XP’和方差SP’;
[0028] M2、根據(jù)式(3)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK:
[0029] Pfinal=P’K*SP’+XP’······式(3)。
[0030] 進(jìn)一步的,所述非線性激活函數(shù)δ(·)為sigmoid函數(shù),根據(jù)式(4)計(jì)算:
[0031] δ(x)=1/(1+e-x),e為自然底數(shù),······式(4)。
[0032] 具體而言的,所述非線性激活函數(shù)δ(·)為雙曲線正切函數(shù),根據(jù)式(5)計(jì)算:
[0033] δ(x)=(ex-e-x)/(ex+e-x),其中e為自然底數(shù),······式(5)。
[0034] 優(yōu)選的,所述非線性激活函數(shù)δ(·)為ReLU函數(shù),根據(jù)式(6)計(jì)算:
[0035]
[0036] 本發(fā)明不需要對圖像濾波進(jìn)行特別的設(shè)計(jì),通過對卷積神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練自動(dòng)完成圖像濾波的學(xué)習(xí),本發(fā)明重建方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和偽影,同時(shí)較好地保持原有圖像的分辨率,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。附圖說明
[0037] 利用附圖對本發(fā)明作進(jìn)一步的說明,但附圖中的內(nèi)容不構(gòu)成對本發(fā)明的任何限制。
[0038] 圖1是本發(fā)明一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法的操作流程示意圖。
[0039] 圖2是圖1卷積神經(jīng)網(wǎng)絡(luò)濾波中的卷積神經(jīng)網(wǎng)絡(luò)構(gòu)架示意圖。
[0040] 圖3為圖1中部分操作的效果圖。

具體實(shí)施方式

[0041] 結(jié)合以下實(shí)施例對本發(fā)明作進(jìn)一步描述。
[0042] 實(shí)施例1。
[0043] 如圖1-3所示,一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,包括如下步驟:
[0044] A1、對原始弦圖數(shù)據(jù)IK進(jìn)行反投影操作,得到反投影圖像數(shù)據(jù)I’K。
[0045] 步驟A1中反投影操作具體是通過CT掃描機(jī)對原始弦圖數(shù)據(jù)IK進(jìn)行幾何成像處理。
[0046] A2、對反投影圖像數(shù)據(jù)I’K進(jìn)行歸一化處理,得到歸一化反投影圖像數(shù)據(jù)PK。
[0047] 步驟A2中歸一化處理的方法步驟如下:
[0048] T1、計(jì)算反投影圖像數(shù)據(jù)I’K的均值XI’和方差SI’。
[0049] T2、根據(jù)式(1)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK。
[0050] PK=(I’K-XI’)/SI’······式(1)。
[0051] A3、將歸一化反投影圖像PK通過卷積神經(jīng)網(wǎng)絡(luò)進(jìn)行卷積神經(jīng)網(wǎng)絡(luò)濾波,生成待處理圖像P’K。
[0052] 卷積神經(jīng)網(wǎng)絡(luò)具有多個(gè)卷積層,卷積神經(jīng)網(wǎng)絡(luò)的卷積層的輸入和輸出均稱為特征圖,卷積層的輸入對應(yīng)于輸入特征圖,卷積層的輸出對應(yīng)于輸出特征圖;每個(gè)特征圖具有多個(gè)通道,通道數(shù)目取決于對應(yīng)卷積層中的卷積核的數(shù)目。
[0053] 步驟A3中濾波的具體操作如下:
[0054] 設(shè)第s層卷積層有ns卷積核,第s層的輸出特征圖記為Fs,第s層的輸入特征圖記為Fs-1,輸出特征圖Fs的第i個(gè)通道記為 輸入特征圖Fs-1的第j個(gè)通道記為 n為自然數(shù)。
[0055] 卷積神經(jīng)網(wǎng)絡(luò)的輸入記為F0,卷積神經(jīng)網(wǎng)絡(luò)中第s層卷積層輸出Fs與輸入Fs-1之間的關(guān)系如下:
[0056]
[0057] 式(2)中,δ(·)表示輸入矩陣中元素的非線性激活函數(shù),*表示二維卷積操作,為第s層卷積層中輸出特征圖第i個(gè)通道與輸入特征圖中第j個(gè)通道進(jìn)行卷積的二維卷積核, 為第s層卷積層中用于計(jì)算輸出特征圖第i個(gè)通道的偏置項(xiàng),S表示卷積網(wǎng)絡(luò)中的卷積層的數(shù)目,S為自然數(shù)。
[0058] 步驟A3將步驟A2中得到的歸一化反投影圖像數(shù)據(jù)PK作為卷積神經(jīng)網(wǎng)絡(luò)的輸入F0。
[0059] A4、對待處理圖像P’K進(jìn)行反歸一化處理,得到最終的重建圖像Pfinal。
[0060] 步驟A4中的反歸一化處理為步驟A2的逆操作。
[0061] 步驟A4中的反歸一化處理的具體方法如下:
[0062] M1、計(jì)算反待處理圖像P’K的均值XP’和方差SP’。
[0063] M2、根據(jù)式(3)計(jì)算得出歸一化反投影圖像數(shù)據(jù)PK。
[0064] Pfinal=P’K*SP’+XP’······式(3)。
[0065] 本發(fā)明不需要對圖像濾波進(jìn)行特別的設(shè)計(jì),通過對卷積神經(jīng)網(wǎng)絡(luò)模型的訓(xùn)練自動(dòng)完成圖像濾波的學(xué)習(xí),本發(fā)明重建方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和偽影,同時(shí)較好地保持原有圖像的分辨率,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。
[0066] 實(shí)施例2。
[0067] 一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其它特征與實(shí)施例1相同,不同之處在于:非線性激活函數(shù)δ(·)為sigmoid函數(shù),根據(jù)式(4)計(jì)算:
[0068] δ(x)=1/(1+e-x),e為自然底數(shù),······式(4)。
[0069] 需要說明的是,非線性激活函數(shù)δ(·)可以根據(jù)卷積神經(jīng)網(wǎng)絡(luò)訓(xùn)練的目標(biāo)函數(shù)和優(yōu)化算法選取非線性激活函數(shù)的類型。
[0070] 處理方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和和偽影的同時(shí)較好地保持原有圖像的分辨率,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。
[0071] 實(shí)施例3。
[0072] 一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其它特征與實(shí)施例1相同,不同之處在于:非線性激活函數(shù)δ(·)為雙曲線正切函數(shù),根據(jù)式(5)計(jì)算:
[0073] δ(x)=(ex-e-x)/(ex+e-x),其中e為自然底數(shù),······式(5)。
[0074] 需要說明的是,非線性激活函數(shù)δ(·)可以根據(jù)卷積神經(jīng)網(wǎng)絡(luò)訓(xùn)練的目標(biāo)函數(shù)和優(yōu)化算法選取非線性激活函數(shù)的類型。
[0075] 處理方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和和偽影的同時(shí)較好地保持原有圖像的分辨率,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。
[0076] 實(shí)施例4。
[0077] 一種基于卷積神經(jīng)網(wǎng)絡(luò)的CT圖像重建方法,其它特征與實(shí)施例1相同,不同之處在于:非線性激活函數(shù)δ(·)為ReLU函數(shù),根據(jù)式(6)計(jì)算:
[0078]
[0079] 需要說明的是,非線性激活函數(shù)δ(·)可以根據(jù)卷積神經(jīng)網(wǎng)絡(luò)訓(xùn)練的目標(biāo)函數(shù)和優(yōu)化算法選取非線性激活函數(shù)的類型。
[0080] 處理方法操作簡單,處理方便,能夠在大幅度減少圖像噪聲和和偽影的同時(shí)較好地保持原有圖像的分辨率,最終實(shí)現(xiàn)CT圖像的優(yōu)質(zhì)重建。
[0081] 最后應(yīng)當(dāng)說明的是,以上實(shí)施例僅用以說明本發(fā)明的技術(shù)方案而非對本發(fā)明保護(hù)范圍的限制,盡管參照較佳實(shí)施例對本發(fā)明作了詳細(xì)說明,本領(lǐng)域的普通技術(shù)人員應(yīng)當(dāng)理解,可以對本發(fā)明的技術(shù)方案進(jìn)行修改或者等同替換,而不脫離本發(fā)明技術(shù)方案的實(shí)質(zhì)和范圍。
高效檢索全球?qū)@?/div>

專利匯是專利免費(fèi)檢索,專利查詢,專利分析-國家發(fā)明專利查詢檢索分析平臺(tái),是提供專利分析,專利查詢,專利檢索等數(shù)據(jù)服務(wù)功能的知識產(chǎn)權(quán)數(shù)據(jù)服務(wù)商。

我們的產(chǎn)品包含105個(gè)國家的1.26億組數(shù)據(jù),免費(fèi)查、免費(fèi)專利分析。

申請?jiān)囉?/a>

QQ群二維碼
意見反饋