一種ct金屬偽影處理方法和裝置制造方法
【專利摘要】本發(fā)明提供了一種CT金屬偽影處理方法及裝置，所述處理方法包括：對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像COR1；提取所述第一處理圖像COR1的高頻部分，得到所述第一處理圖像的高頻圖像COR1HF；對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像COR2；提取所述第二處理圖像COR2的高頻部分得到所述第二處理圖像的高頻圖像COR2HF；提取金屬物圖像Metal，根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W；將所述第一處理圖像COR1、第一處理圖像COR1的高頻圖像COR1HF以及第二處理圖像COR2的高頻圖像COR2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。該處理方法提高了圖像質(zhì)量。
【專利說(shuō)明】一種CT金屬偽影處理方法和裝置

【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及醫(yī)學(xué)影像的圖像處理【技術(shù)領(lǐng)域】，尤其涉及一種CT金屬偽影處理方法和裝置。

【背景技術(shù)】
[0002]CT (Computed Tomography)，即電子計(jì)算機(jī)斷層掃描，它是利用精確準(zhǔn)直的X線束、γ射線、超聲波等，與靈敏度極高的探測(cè)器一同圍繞人體的某一部位作一個(gè)接一個(gè)的斷面掃描，具有掃描時(shí)間快，圖像清晰等特點(diǎn)，可用于多種疾病的檢查。
[0003]通過(guò)CT系統(tǒng)對(duì)被掃描對(duì)象重建斷層圖像時(shí)，若被掃描對(duì)象內(nèi)有金屬物體，則會(huì)產(chǎn)生嚴(yán)重的射線硬化、增大噪聲、增加散射和部分容積效應(yīng)，從而會(huì)在CT圖像中產(chǎn)生強(qiáng)烈的圖像偽影，也就是金屬偽影。該金屬偽影的存在會(huì)降低CT圖像的質(zhì)量。
[0004]現(xiàn)有技術(shù)中已經(jīng)有一些金屬偽影的處理方法可以減少或去除這種金屬偽影，大體過(guò)程如下:
[0005]先重建出未校正的原始圖像，再對(duì)該原始圖像進(jìn)行圖像分割，分割出只含有金屬物的圖像，對(duì)金屬物圖像進(jìn)行正投影獲得金屬在投影數(shù)據(jù)中的范圍，再對(duì)原始投影數(shù)據(jù)中金屬物所在范圍內(nèi)的數(shù)據(jù)進(jìn)行修正，金屬物以外的數(shù)據(jù)不變。修正后再重建出圖像，再將金屬物圖像加入重建的圖像當(dāng)中，就得到了校正后的圖像。
[0006]這種方法得到的校正后的圖像不含有金屬偽影，但是，在校正后的圖像中會(huì)損失金屬物附近區(qū)域的被掃描對(duì)象的信息，如骨骼信息，從而嚴(yán)重降低了圖像質(zhì)量。


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

[0007]有鑒于此，本發(fā)明提供了一種CT金屬偽影處理方法和裝置，以實(shí)現(xiàn)在去除金屬偽影的同時(shí)保留金屬物附近的被掃描對(duì)象的信息，進(jìn)而提高圖像質(zhì)量。
[0008]為了實(shí)現(xiàn)上述發(fā)明目的，本發(fā)明采用了如下技術(shù)方案:
[0009]一種CT金屬偽影處理方法，包括:
[0010]對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像CORl ;所述原始圖像中包括金屬偽影，所述第一金屬偽影處理用于去除所述金屬偽影；提取所述第一處理圖像CORl的高頻部分，得到所述第一處理圖像的高頻圖像CORlHF ;
[0011]對(duì)所述原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2 ;所述第二金屬偽影處理用于去除所述金屬偽影，在所述第二金屬偽影處理過(guò)程中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息；提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像C0R2HF ；
[0012]提取金屬物圖像Metal，根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W ;
[0013]將所述第一處理圖像C0R1、第一處理圖像CORl的高頻圖像CORlHF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。
[0014]可選地，所述加權(quán)組合的具體公式如下:
[0015]CORImp = CORl+ (C0R2HF-C0R1HF)轉(zhuǎn)。
[0016]可選地，所述金屬偽影包括硬化偽影，所述對(duì)所述原始圖像數(shù)據(jù)進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2，具體包括:
[0017]對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以到所述第二處理圖像C0R2。
[0018]可選地，所述金屬偽影還包括條狀偽影，所述對(duì)所述原始圖像進(jìn)行硬化校正之前，還包括:
[0019]對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中的條狀偽影，得到Or1-Air圖像；
[0020]所述對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以得到所述第二處理圖像C0R2，具體為:對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2。
[0021]可選地，所述對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2，具體包括:
[0022]對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air-Smo圖像；
[0023]計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像CORl進(jìn)行平滑濾波處理得到的；
[0024]計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
[0025]可選地，所述獲取所述金屬物圖像Metal的權(quán)重函數(shù)W，具體包括:
[0026]對(duì)所述金屬物圖像Metal進(jìn)行低通濾波，得到低頻金屬物圖像；
[0027]對(duì)所述低頻金屬物圖像進(jìn)行歸一化處理，得到所述權(quán)重函數(shù)W。
[0028]一種CT金屬偽影處理裝置，包括:
[0029]第一處理單元，用于對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像C0R1，所述原始圖像中包括金屬偽影；所述第一金屬偽影處理用于去除所述金屬偽影；
[0030]第一提取單元，用于提取所述第一處理圖像CORl的高頻部分，得到所述第一處理圖像的高頻圖像CORlHF ；
[0031]第二處理單元，用于對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2 ;所述第二金屬偽影處理用于去除所述金屬偽影，在所述第二金屬偽影處理過(guò)程中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息；
[0032]第二提取單元，用于提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像C0R2HF ；
[0033]第三提取單元，用于提取金屬物圖像Metal ;
[0034]獲取單元，用于根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W ；
[0035]加權(quán)組合單元，用于將所述第一處理圖像C0R1、第一處理圖像CORl的高頻圖像CORlHF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。
[0036]可選地，所述第二處理單元具體包括:
[0037]第一處理子單元，用于對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以得到所述第二處理圖像C0R2。
[0038]可選地，所述第二處理單元還包括:
[0039]第二處理子單元，用于對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中條狀偽影，得到Or1-Air圖像；
[0040]所述第一處理子單元包括對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2的子單元。
[0041]可選地，所述第一處理子單元具體包括:
[0042]平滑濾波處理子單元，用于對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air-Smo 圖像；
[0043]第一計(jì)算子單元，用于計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像CORl進(jìn)行平滑濾波處理得到的；
[0044]第二計(jì)算子單元，用于計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
[0045]相較于現(xiàn)有技術(shù)，本發(fā)明具有以下有益效果:
[0046]本發(fā)明實(shí)施例提供的CT金屬偽影處理方法中，將第一處理圖像C0R1、第一處理圖像CORl的高頻圖像CORlHF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。由于第二處理圖像C0R2中去除了金屬偽影但保留有金屬物附近的被掃描對(duì)象的信息，所以，第二處理圖像C0R2的高頻圖像C0R2HF中不含有金屬偽影，但包含有金屬物附近的被掃描對(duì)象的信息，而第一處理圖像以及其高頻圖像CORlHF中不包括金屬偽影，所以，C0R2HF與CORlHF的差值為金屬物附近的被掃描對(duì)象的信息，該金屬物附近的被掃描對(duì)象的信息進(jìn)行加權(quán)計(jì)算后，并與第一處理圖像CORl進(jìn)行組合后得到的處理后的圖像CORImp中就不含有金屬偽影，但含有金屬物附近的被掃描對(duì)象的信息。
[0047]所以，通過(guò)本發(fā)明提供的CT金屬偽影處理方法，既去除了原始圖像中的金屬偽影，又不會(huì)損失圖像中的金屬物附近的被掃描對(duì)象的信息，因此，該CT金屬偽影處理方法提高了圖像質(zhì)量。

【專利附圖】

【附圖說(shuō)明】
[0048]為了清楚地理解本發(fā)明的技術(shù)方案，下面將描述本發(fā)明的【具體實(shí)施方式】時(shí)用到的附圖進(jìn)行簡(jiǎn)要說(shuō)明。顯而易見(jiàn)地，這些附圖僅是本發(fā)明一部分實(shí)施例的附圖，本領(lǐng)域普通技術(shù)人員在不付出創(chuàng)造性勞動(dòng)的前提下，還可以獲得其它附圖。
[0049]圖1是本發(fā)明實(shí)施例一提供的CT金屬偽影處理方法的流程示意圖；
[0050]圖2是原始的CT圖像示意圖；
[0051]圖3是第一處理圖像CORl的示意圖；
[0052]圖4是第二處理圖像C0R2的示意圖；
[0053]圖5金屬物圖像Metal的示意圖；
[0054]圖6金屬物圖像Metal的權(quán)重函數(shù)示意圖；
[0055]圖7是本發(fā)明實(shí)施例二提供的CT金屬偽影處理裝置的結(jié)構(gòu)示意圖。

【具體實(shí)施方式】
[0056]為了使本發(fā)明的發(fā)明目的、技術(shù)效果更加清楚，技術(shù)方案更加完整，下面結(jié)合附圖對(duì)本發(fā)明的【具體實(shí)施方式】進(jìn)行描述。
[0057]實(shí)施例一
[0058]圖1是本發(fā)明實(shí)施例一提供的CT金屬偽影處理方法的流程示意圖。如圖1所示，該CT金屬偽影處理方法包括以下步驟:
[0059]S101、對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像CORl:
[0060]需要說(shuō)明的是，所述原始圖像中包括金屬偽影。其中，原始圖像是根據(jù)原始數(shù)據(jù)重建出的CT圖像，原始數(shù)據(jù)是在CT掃描過(guò)程中得到的掃描數(shù)據(jù)。其中，原始圖像的示意圖如圖2所示，其包括放射狀金屬偽影，該金屬偽影包括強(qiáng)烈的條狀偽影和嚴(yán)重的硬化偽影，這是由于金屬物體吸收射線后產(chǎn)生的射線硬化、噪聲和散射增大造成的。
[0061]所述第一金屬偽影處理用于去除金屬偽影和金屬偽影數(shù)據(jù)。本步驟所述的第一金屬偽影處理方法可以為本領(lǐng)域公知的一些金屬偽影校正方法。
[0062]作為一個(gè)示例，所述第一金屬偽影處理的方法可以如下:首先根據(jù)CT掃描得到的原始數(shù)據(jù)重建出未校正的原始圖像，再對(duì)該原始圖像進(jìn)行圖像分割，分割出只含有金屬物的圖像，對(duì)金屬物圖像進(jìn)行正投影獲得金屬在投影數(shù)據(jù)中的范圍，再對(duì)原始投影數(shù)據(jù)中金屬物所在范圍內(nèi)的數(shù)據(jù)進(jìn)行修正，金屬物以外的數(shù)據(jù)不變。修正后再重建出圖像，再將金屬物圖像加入重建的圖像當(dāng)中，就得到了校正后的圖像。
[0063]該第一金屬偽影處理方法得到的第一處理圖像CORl非常接近最終需要的校正圖像。但是，其不足之處是，在去除金屬偽影和金屬偽影數(shù)據(jù)的過(guò)程中，金屬物附近區(qū)域的被掃描對(duì)象的信息也損失掉了。因而，通過(guò)第一金屬偽影處理得到的第一處理圖像CORl相較于原始圖像，去除了原始圖像中的金屬偽影，并且金屬物附近區(qū)域的被掃描對(duì)象的信息通常情況下也被損失掉了。這些信息主要是骨骼信息。但是，第一處理圖像CORl中含有除金屬物附近的其它區(qū)域的被掃描對(duì)象的信息，該第一處理圖像CORl的示意圖如圖3所示。
[0064]S102、提取所述第一處理圖像CORl的高頻部分，得到所述第一處理圖像的高頻圖像 CORlHF:
[0065]具體地，可以利用高通濾波器對(duì)第一處理圖像CORl進(jìn)行濾波，提取第一處理圖像CORl的高頻部分，得到所述第一處理圖像的高頻圖像C0R1HF。
[0066]需要說(shuō)明的是，金屬偽影圖像為低頻圖像，被掃描對(duì)象的信息圖像為高頻圖像，所以，提取第一處理圖像得到的第一處理圖像的高頻圖像CORlHF為被掃描對(duì)象的信息圖像。
[0067]S103、對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2:
[0068]第二金屬偽影處理用于去除原始圖像中的金屬偽影。相較于第一金屬偽影處理，第二金屬偽影處理方法在對(duì)原始圖像進(jìn)行處理的過(guò)程中，只對(duì)原始圖像中的金屬偽影造成影響，而不會(huì)影響被掃描對(duì)象的信息圖像，因此，第二金屬偽影處理過(guò)程中不會(huì)損失金屬物附近區(qū)域的被掃描對(duì)象的信息，也就是說(shuō)，得到的第二處理圖像C0R2中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息。此外，不在金屬物附近的被掃描對(duì)象的信息在第二處理圖像中也保留下來(lái)。該第二處理圖像C0R2的示意圖如圖4所示。
[0069]所以，第一金屬偽影處理過(guò)程和第二金屬偽影處理過(guò)程的異同點(diǎn)具體如下:
[0070]兩者均能去除原始圖像中的金屬偽影，得到的第一處理圖像和第二處理圖像中均不再含有金屬偽影，且均含有除金屬物附近區(qū)域以外的其它區(qū)域的被掃描對(duì)象的信息。其區(qū)別點(diǎn)是:第一金屬偽影處理有可能損失了金屬物附近區(qū)域的被掃描對(duì)象的信息，而第二金屬偽影處理過(guò)程卻保留了金屬物附近區(qū)域的被掃描對(duì)象的信息。另外，第二金屬偽影處理的圖像整體效果沒(méi)有第一金屬偽影處理的圖像整體效果好，因此，可以利用圖像整體效果較好的第一處理圖像和第二處理圖像中的金屬物附近區(qū)域的信息來(lái)組合成圖像校正效果更好的處理圖像。本發(fā)明就是利用上述構(gòu)思來(lái)處理帶有金屬偽影的CT圖像的。
[0071]簡(jiǎn)單地說(shuō)，本發(fā)明進(jìn)行第二金屬偽影處理的主要目的就是為了得到金屬物附近區(qū)域的被掃描對(duì)象的信息，從而對(duì)后續(xù)的加權(quán)組合作準(zhǔn)備。
[0072]如上所述，金屬偽影可以包括硬化偽影，所述對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像COR2，具體包括:
[0073]對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以得到所述第二處理圖像COR2。
[0074]另外，所述金屬偽影除了包括硬化偽影外，還有可能包括條狀偽影，所述對(duì)初始圖像數(shù)據(jù)進(jìn)行第二金屬偽影處理，得到第二處理圖像COR2，具體包括消減條狀偽影和硬化偽影的步驟，其具體步驟如下:
[0075]A、對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中條狀偽影，得到Or1-Air圖像:
[0076]所述預(yù)處理方法可以為生數(shù)據(jù)迭代方法或自適應(yīng)濾波方法。由于在預(yù)處理的過(guò)程中，消減了原始圖像中的條狀偽影，所以，得到的Or1-Air圖像中不含有條狀偽影。
[0077]B、對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2:
[0078]由于Or1-Air圖像中不包含有條狀偽影，通過(guò)硬化校正又消減了 Or1-Air圖像中的硬化偽影，所以，處理得到的第二處理圖像C0R2中既不包括條狀偽影，也不包括硬化偽影，又由于金屬偽影由條狀偽影和硬化偽影組成，所以，得到的第二處理圖像C0R2中去除了金屬偽影。
[0079]作為本發(fā)明的一個(gè)優(yōu)選實(shí)施例，所述步驟B可以通過(guò)以下方法實(shí)現(xiàn):
[0080]B1、對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air-Smo圖像；
[0081]B2、計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像CORl進(jìn)行平滑濾波處理得到的；
[0082]B3、計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
[0083]S104、提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像 C0R2HF:
[0084]具體地，利用高通濾波器對(duì)第二處理圖像C0R2進(jìn)行濾波，提取第二處理圖像C0R2的高頻部分，得到所述第二處理圖像的高頻圖像C0R2HF。
[0085]需要說(shuō)明的是，與第一處理圖像CORl相同，金屬偽影圖像為低頻圖像，被掃描對(duì)象的信息圖像為高頻圖像，所以，提取第二處理圖像C0R2得到的第二處理圖像的高頻圖像C0R2HF為被掃描對(duì)象的信息圖像。
[0086]S105、提取金屬物圖像Metal，根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W:
[0087]具體地，將含有金屬物圖像Metal的CT圖像進(jìn)行圖像閾值分割，將CT值在預(yù)設(shè)閾值以上的像素點(diǎn)分割提取出來(lái)，以得到所述金屬物圖像Metal。將CT值在預(yù)設(shè)閾值以上的像素點(diǎn)定義為金屬點(diǎn)，CT值在預(yù)設(shè)閾值以下的像素點(diǎn)定義為無(wú)金屬點(diǎn)。這樣由分割出來(lái)的金屬點(diǎn)組成金屬物圖像Metal。得到的金屬圖像物圖像Metal的示意圖如圖5所示。在本發(fā)明實(shí)施例中，所述含有金屬物圖像Metal的CT圖像可以為原始圖像，也可以為第一處理圖像CORl，還可以為第二處理圖像C0R2。
[0088]然后，根據(jù)提取到的金屬物圖像Metal得到金屬物圖像Metal的權(quán)重函數(shù)W。其中，一個(gè)具體的方法可以為:
[0089]對(duì)所述金屬物圖像Metal進(jìn)行低通濾波，得到低頻金屬物圖像；
[0090]對(duì)所述低頻金屬物圖像進(jìn)行歸一化處理，得到所述權(quán)重函數(shù)W。
[0091]權(quán)重函數(shù)W的示意圖如圖6所述。需要說(shuō)明的是，權(quán)重函數(shù)W的特點(diǎn)是在金屬物及附近區(qū)域權(quán)重較高而在遠(yuǎn)離金屬物的區(qū)域權(quán)重較低。也就是說(shuō)，隨著相距金屬物距離的增大，權(quán)重逐漸變小。
[0092]S106、將所述第一處理圖像COR1、第一處理圖像CORl的高頻圖像CORlHF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp:
[0093]具體的加權(quán)組合的公式可以如下:
[0094]CORImp = CORl+(C0R2HF-C0R1HF)轉(zhuǎn)。
[0095]由于第二處理圖像的高頻圖像C0R2HF中不含有金屬偽影，但是含有金屬物附近的信息，而第一處理圖像的高頻圖像CORlHF中不含有金屬偽影，并且也不含有金屬物附近的信息，所以，兩者的差值即為金屬物附近的信息，該金屬物附近的信息與權(quán)重函數(shù)W的乘積即為金屬物附近的信息在CT圖像中的具體位置分布。將該金屬物附近的信息與權(quán)重函數(shù)的乘積與第一處理圖像CORl進(jìn)行加和得到的圖像即為本發(fā)明提供的CT金屬偽影處理方法最終得到的處理圖像CORImp。
[0096]由于第一處理圖像CORl中去除了金屬偽影，且整體處理效果較好，除了不含有金屬物附近的信息外，其它信息均非常接近理想的圖像。所以，利用該第一處理圖像CORl與金屬物附近的信息與權(quán)重函數(shù)的乘積進(jìn)行加和得到的處理圖像CORImp中既去除了原始圖像中的金屬偽影，也包括了金屬物附近的信息，彌補(bǔ)了現(xiàn)有技術(shù)中的處理方法損失金屬物附近的信息的不足。因此，通過(guò)本發(fā)明提供的CT金屬偽影處理方法得到的處理圖像既去除了金屬偽影，又保留了金屬物附近的信息，因此，相較于現(xiàn)有技術(shù)中的CT金屬偽影處理方法，本發(fā)明提供的CT金屬偽影處理方法提高了處理圖像的質(zhì)量。
[0097]此外，需要說(shuō)明的是，作為實(shí)施例的擴(kuò)展，本發(fā)明提供的CT金屬偽影處理方法不限于上述實(shí)施例一所述的各個(gè)步驟的執(zhí)行順序，其還可以先對(duì)原始圖像進(jìn)行第二金屬偽影處理，再對(duì)原始圖像進(jìn)行第一金屬偽影處理。但是，如果第二金屬偽影處理利用了第一處理圖像時(shí)，此時(shí)不能顛倒兩者的處理順序，此時(shí)，必須先執(zhí)行第一金屬偽影處理，再執(zhí)行第二金屬偽影處理。
[0098]以上為本發(fā)明實(shí)施例一提供的CT金屬偽影處理方法?；谏鲜鯟T金屬偽影處理方法，本發(fā)明實(shí)施例還提供了一種CT金屬偽影處理裝置。具體參見(jiàn)實(shí)施例二。
[0099]實(shí)施例二
[0100]圖7是本發(fā)明實(shí)施例二提供的CT金屬偽影處理裝置的結(jié)構(gòu)示意圖。如圖7所示，實(shí)施例二所述的CT金屬偽影處理裝置包括以下單元:
[0101]第一處理單元71，用于對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像COR1，所述原始圖像中包括金屬偽影；所述第一金屬偽影處理用于去除所述金屬偽影；
[0102]第一提取單元72，用于提取所述第一處理圖像CORl的高頻部分，得到所述第一處理圖像的高頻圖像CORlHF ；
[0103]第二處理單元73，用于對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2 ;所述第二金屬偽影處理用于去除所述金屬偽影，在所述第二金屬偽影處理過(guò)程中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息；
[0104]第二提取單元74，用于提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像C0R2HF ；
[0105]第三提取單元75，用于提取金屬物圖像Metal ；
[0106]獲取單元76，用于根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W ;
[0107]加權(quán)組合單元77，用于將所述第一處理圖像C0R1、第一處理圖像CORl的高頻圖像CORlHF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。
[0108]進(jìn)一步地，所述加權(quán)組合的具體公式如下:
[0109]CORImp = CORl+(C0R2HF-C0R1HF)轉(zhuǎn)。
[0110]通過(guò)上述實(shí)施例二提供的CT金屬偽影處理裝置得到的處理圖像，既去除了 CT圖像中的金屬偽影，又保留了金屬物附近的被掃描對(duì)象的信息，相較于現(xiàn)有技術(shù)中得到的處理圖像，圖像質(zhì)量得到了提高。
[0111]進(jìn)一步地，所述金屬偽影包括條狀偽影和硬化偽影，所以，用于去除金屬偽影的第二處理單元73既可以去除條狀偽影，又可以去除硬化偽影，因此，所述第二處理單元73可以具體包括:
[0112]第一處理子單元731，用于對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像或數(shù)據(jù)中的硬化偽影，以得到所述第二處理圖像C0R2。
[0113]進(jìn)一步地，所述第二處理單元73還可以包括:
[0114]第二處理子單元732，用于對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中條狀偽影，得到Or1-Air圖像；
[0115]在這種情形下，第二處理子單元731包括對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2的子單元。
[0116]為了較為簡(jiǎn)單地實(shí)現(xiàn)對(duì)Or1-Air圖像進(jìn)行硬化校正，上述第一處理子單元731可以具體包括:
[0117]平滑濾波處理子單元7311，用于對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air-Smo 圖像；
[0118]第一計(jì)算子單元7312，用于計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像CORl進(jìn)行平滑濾波處理得到的；
[0119]第二計(jì)算子單元7313，用于計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
[0120]另外，作為本發(fā)明的一個(gè)具體實(shí)施例，所述獲取單元76可以具體包括以下單元:
[0121]低通濾波子單元761，用于對(duì)所述金屬物圖像Metal進(jìn)行低通濾波，得到低頻金屬物圖像；
[0122]歸一化處理子單元762，用于對(duì)所述低頻金屬物圖像進(jìn)行歸一化處理，得到所述權(quán)重函數(shù)W。
[0123]對(duì)于裝置實(shí)施例而言，由于其基本對(duì)應(yīng)于方法實(shí)施例，所以相關(guān)之處參見(jiàn)方法實(shí)施例的部分說(shuō)明即可。以上所描述的裝置實(shí)施例僅僅是示意性的，其中所述作為分離部件說(shuō)明的單元可以是或者也可以不是物理上分開(kāi)的?？梢愿鶕?jù)實(shí)際的需要選擇其中的部分或者全部模塊來(lái)實(shí)現(xiàn)本實(shí)施例方案的目的。本領(lǐng)域普通技術(shù)人員在不付出創(chuàng)造性勞動(dòng)的情況下，即可以理解并實(shí)施。
[0124]以上所述僅是本發(fā)明的優(yōu)選實(shí)施方式，應(yīng)當(dāng)指出，對(duì)于本【技術(shù)領(lǐng)域】的普通技術(shù)人員來(lái)說(shuō)，在不脫離本發(fā)明原理的前提下，還可以做出若干改進(jìn)和潤(rùn)飾，這些改進(jìn)和潤(rùn)飾也應(yīng)視為本發(fā)明的保護(hù)范圍。
【權(quán)利要求】
1.一種CT金屬偽影處理方法，其特征在于，包括: 對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像C0R1 ;所述原始圖像中包括金屬偽影，所述第一金屬偽影處理用于去除所述金屬偽影；提取所述第一處理圖像C0R1的高頻部分，得到所述第一處理圖像的高頻圖像C0R1HF ； 對(duì)所述原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2 ;所述第二金屬偽影處理用于去除所述金屬偽影，在所述第二金屬偽影處理過(guò)程中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息；提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像 C0R2HF ； 提取金屬物圖像Metal，根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W ; 將所述第一處理圖像C0R1、第一處理圖像C0R1的高頻圖像C0R1HF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。
2.根據(jù)權(quán)利要求1所述的處理方法，其特征在于，所述加權(quán)組合的具體公式如下:
CORImp = C0R1+(C0R2HF-C0R1HF)轉(zhuǎn)。
3.根據(jù)權(quán)利要求1或2所述的處理方法，其特征在于，所述金屬偽影包括硬化偽影，所述對(duì)所述原始圖像數(shù)據(jù)進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2，具體包括: 對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以到所述第二處理圖像C0R2。
4.根據(jù)權(quán)利要求3所述的處理方法，其特征在于，所述金屬偽影還包括條狀偽影，所述對(duì)所述原始圖像進(jìn)行硬化校正之前，還包括: 對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中的條狀偽影，得到Or1-Air圖像； 所述對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以得到所述第二處理圖像C0R2，具體為:對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Ori_Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2。
5.根據(jù)權(quán)利要求4所述的處理方法，其特征在于，所述對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2，具體包括: 對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air_Smo圖像； 計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像C0R1進(jìn)行平滑濾波處理得到的； 計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
6.根據(jù)權(quán)利要求1或2所述的處理方法，其特征在于，所述獲取所述金屬物圖像Metal的權(quán)重函數(shù)W，具體包括: 對(duì)所述金屬物圖像Metal進(jìn)行低通濾波，得到低頻金屬物圖像； 對(duì)所述低頻金屬物圖像進(jìn)行歸一化處理，得到所述權(quán)重函數(shù)W。
7.—種CT金屬偽影處理裝置，其特征在于，包括: 第一處理單元，用于對(duì)原始圖像進(jìn)行第一金屬偽影處理，得到第一處理圖像C0R1，所述原始圖像中包括金屬偽影；所述第一金屬偽影處理用于去除所述金屬偽影； 第一提取單元，用于提取所述第一處理圖像C0R1的高頻部分，得到所述第一處理圖像的高頻圖像C0R1HF ； 第二處理單元，用于對(duì)原始圖像進(jìn)行第二金屬偽影處理，得到第二處理圖像C0R2 ;所述第二金屬偽影處理用于去除所述金屬偽影，在所述第二金屬偽影處理過(guò)程中保留有金屬物附近區(qū)域的被掃描對(duì)象的信息； 第二提取單元，用于提取所述第二處理圖像C0R2的高頻部分得到所述第二處理圖像的高頻圖像C0R2HF ； 第三提取單元，用于提取金屬物圖像Metal ； 獲取單元，用于根據(jù)所述金屬物圖像Metal獲取權(quán)重函數(shù)W ； 加權(quán)組合單元，用于將所述第一處理圖像C0R1、第一處理圖像C0R1的高頻圖像C0R1HF以及第二處理圖像C0R2的高頻圖像C0R2HF用所述權(quán)重函數(shù)W進(jìn)行加權(quán)組合，得到處理后的圖像CORImp。
8.根據(jù)權(quán)利要求7所述的處理裝置，其特征在于，所述第二處理單元具體包括: 第一處理子單元，用于對(duì)所述原始圖像進(jìn)行硬化校正，以消減所述原始圖像中的硬化偽影，以得到所述第二處理圖像C0R2。
9.根據(jù)權(quán)利要求8所述的處理裝置，其特征在于，所述第二處理單元還包括: 第二處理子單元，用于對(duì)原始數(shù)據(jù)進(jìn)行預(yù)處理，以消減所述原始圖像中條狀偽影，得到Or1-Air 圖像； 所述第一處理子單元包括對(duì)所述Or1-Air圖像進(jìn)行硬化校正，以消減所述Or1-Air圖像中的硬化偽影，以得到所述第二處理圖像C0R2的子單元。
10.根據(jù)權(quán)利要求9所述的處理裝置，其特征在于，所述第一處理子單元具體包括: 平滑濾波處理子單元，用于對(duì)所述Or1-Air圖像進(jìn)行平滑濾波處理，得到Or1-Air_Smo圖像； 第一計(jì)算子單元，用于計(jì)算所述Or1-Air-Smo圖像與CORl-Smo圖像的差值，所述差值為只含有硬化偽影的圖像HardenArti ;所述CORl-Smo圖像是對(duì)所述第一處理圖像C0R1進(jìn)行平滑濾波處理得到的； 第二計(jì)算子單元，用于計(jì)算所述Or1-Air圖像與所述HardenArti的差值，所述差值為第二處理圖像C0R2。
【文檔編號(hào)】G06T7/00GK104504655SQ201410741819
【公開(kāi)日】2015年4月8日 申請(qǐng)日期:2014年12月4日 優(yōu)先權(quán)日:2014年12月4日
【發(fā)明者】韓俊龍, 逄嶺, 樓珊珊 申請(qǐng)人:沈陽(yáng)東軟醫(yī)療系統(tǒng)有限公司