針對(duì)單一模態(tài)的醫(yī)學(xué)圖像無(wú)法為臨床診斷提供全面、互補(bǔ)信息的問(wèn)題,提出一種基于非下采樣剪切波變換(NSST)與改進(jìn)稀疏表示(ISR)的多模態(tài)醫(yī)學(xué)圖像融合方法.首先用NSST分解工具將待融合圖像分解為一個(gè)低頻子帶和若干個(gè)高頻子帶;其次,用ISR方法融合低頻子帶,通過(guò)Sobel算子和引導(dǎo)濾波器去除低頻子帶的細(xì)節(jié)特征,從而提高低頻子帶的融合效率,同時(shí)對(duì)高頻子帶采用絕對(duì)值最大的融合規(guī)則進(jìn)行融合;最后,將融合后的低頻子帶和高頻子帶進(jìn)行逆NSST變換得到最終的融合圖像.實(shí)驗(yàn)結(jié)果表明,該方法在主觀視覺(jué)性能和客觀評(píng)價(jià)上均優(yōu)于其他對(duì)比融合方法. 

【文章來(lái)源】：吉林大學(xué)學(xué)報(bào)(理學(xué)版). 2020年04期 北大核心

【文章頁(yè)數(shù)】：6 頁(yè)

【部分圖文】：

過(guò)完備字典訓(xùn)練集

算法流程

本文實(shí)驗(yàn)采用MATLAB R2016a, CPU為Intel(R) Core (TM) i7-7700, 內(nèi)存為8 GB. 為了驗(yàn)證本文算法的有效性, 使用兩種不同模態(tài)的醫(yī)學(xué)圖像, 其中包括20對(duì)MR-T1與MR-T2圖像, 20對(duì)MR與CT圖像. 所有圖像均來(lái)自美國(guó)哈佛大學(xué)創(chuàng)建的The Whole Brain Atlas醫(yī)學(xué)圖像數(shù)據(jù)庫(kù)(http://www.med.harvard.edu/AANLIB/), 所有源圖像的大小均為256×256, 且實(shí)驗(yàn)所用圖像都經(jīng)過(guò)了配準(zhǔn). 為了驗(yàn)證算法的有效性, 將本文算法與基于邊界尋找方法(BF)[13]、 基于密集尺度不變特征變換(DSITF)[14]、 基于卷積稀疏表示方法(CSR)[15]、 基于卷積稀疏性的形態(tài)成分分析方法(CSMCA)[16]和基于四叉樹的方法(QB)[17]進(jìn)行比較. 由于融合質(zhì)量評(píng)估尚未形成統(tǒng)一的評(píng)價(jià)指標(biāo), 在實(shí)際應(yīng)用中只能根據(jù)經(jīng)驗(yàn)進(jìn)行選擇, 這將導(dǎo)致評(píng)估結(jié)果缺乏說(shuō)服力. 因此, 本文采用4個(gè)不同方面的評(píng)價(jià)指標(biāo)進(jìn)行驗(yàn)證, 并對(duì)融合效果進(jìn)行綜合評(píng)估, 其中包括局部重要性質(zhì)量指數(shù)(Q0)[18]、 Piella的結(jié)構(gòu)相似性度量(Qw)、 熵(EN)和標(biāo)準(zhǔn)差(SD). 圖3為MR-T1與MR-T2圖像的融合結(jié)果.由圖3可見(jiàn), 利用BF算法得到的融合圖像只存在MR-T1圖像的信息, 缺少M(fèi)R-T2圖像的信息; 利用DSIFT算法得到的融合圖像出現(xiàn)多處塊效應(yīng); 利用CSR與QB算法得到的融合圖像也存在塊效應(yīng); 利用CSMCA算法得到的融合圖像雖然沒(méi)有塊效應(yīng), 但亮度較低; 利用本文算法得到的融合圖像既不存在塊效應(yīng), 也未損失源圖像的重要信息.

【參考文獻(xiàn)】：
期刊論文
[1]多級(jí)方向引導(dǎo)濾波器及其在多傳感器圖像融合中的應(yīng)用[J]. 陳廣秋,梁小偉,段錦,才華.  吉林大學(xué)學(xué)報(bào)(理學(xué)版). 2019(01)



本文編號(hào)：2941225