專利名稱：圖像處理設(shè)備和方法
技術(shù)領(lǐng)域：
本公開涉及圖像處理領(lǐng)域，具體而言，涉及用于對圖像(如醫(yī)學(xué)圖像)降噪的圖像處理設(shè)備和方法
背景技術(shù)：
圖像降噪是圖像處理的一個重要環(huán)節(jié)。例如，在醫(yī)學(xué)圖像處理中，由于設(shè)備的信號檢測性能、成像算法不適當(dāng)、掃描劑量低等諸如此類的原因，利用計算機(jī)斷層掃描(Computed Tomography, CT)設(shè)備、核磁共振成像(Magnetic Resonance Imaging, MRI)設(shè)備、超聲波(Ultrasound，UL)診斷裝置等醫(yī)療設(shè)備獲取的醫(yī)學(xué)圖像常常會包含噪聲。這些圖像中的噪聲會嚴(yán)重影響醫(yī)生的診斷，并且會影響后續(xù)圖像處理(如圖像分割、測量和體繪制(Volume Rending)等)的正確性。

發(fā)明內(nèi)容
在下文中給出關(guān)于本公開的一些方面的簡要概述，以便提供對于本公開的基本理解。該概述并非意圖確定本公開的關(guān)鍵或重要部分，也不是意圖限定本公開的范圍。其目的僅僅是以簡化的形式給出某些概念，以此作為稍后論述的更詳細(xì)描述的前序。根據(jù)本公開的一個方面，提供了一種圖像處理設(shè)備。該圖像處理設(shè)備可以包括 濾波裝置，用于對輸入圖像進(jìn)行迭代濾波；以及迭代濾波停止裝置，用于根據(jù)所述迭代濾波過程中每次濾波后得到的經(jīng)濾波圖像隨迭代次數(shù)的變化速率來確定是否停止所述迭代濾波。根據(jù)本公開的另一方面，提供了一種圖像處理方法。該圖像處理方法可以包括對輸入圖像進(jìn)行迭代濾波；以及根據(jù)所述迭代濾波過程中每次濾波后得到的經(jīng)濾波圖像隨迭代次數(shù)的變化速率來確定是否停止所述迭代濾波。另外，本公開的實施例還提供了用于實現(xiàn)上述方法的計算機(jī)程序。此外，本公開的實施例還提供了至少計算機(jī)可讀介質(zhì)形式的計算機(jī)程序產(chǎn)品，其上記錄有用于實現(xiàn)上述方法的計算機(jī)程序代碼。


參照下面結(jié)合附圖對本公開實施例的說明，會更加容易地理解本公開的以上和其它目的、特點和優(yōu)點。附圖中的部件只是為了示出本公開的原理。在附圖中，相同的或類似的技術(shù)特征或部件將采用相同或類似的附圖標(biāo)記來表示。圖I是示出根據(jù)本公開的一個實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖；圖2是示出根據(jù)該實施例的圖像處理方法的示意性流程圖；圖3是示出根據(jù)本公開的一個具體實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖；圖4是示出根據(jù)該具體實施例的圖像處理方法的示意性流程圖；圖5是示出用于確定停止迭代濾波的時間的一個具體示例的示意性流程圖；圖6(A)、6(B)和6(C)是分別示出對不同類型的圖像進(jìn)行迭代濾波的過程中經(jīng)濾波圖像與噪聲圖像的相關(guān)系數(shù)的變化曲線的示意圖；圖7是示出根據(jù)本公開的實施例的應(yīng)用場景的一個示例的示意圖；圖8是示出根據(jù)本公開的另一具體實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖；圖9是示出根據(jù)本公開的另一具體實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖；圖10是示出圖8所示的實施例的圖像處理方法的示意性流程圖；圖11是示出圖9所示的實施例的圖像處理方法的示意性流程圖；圖12(A)、12(B)是分別示出在迭代濾波過程中經(jīng)濾波圖像和噪聲圖像的方差的變化曲線的示意圖；圖13是示出根據(jù)本公開的另一具體實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框 圖；圖14是示出根據(jù)該具體實施例的圖像處理方法的示意性流程圖；圖15是示出該具體實施例的一個具體應(yīng)用示例的示意圖；圖16是示出根據(jù)本公開的另一實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖；圖17是示出根據(jù)該實施例的圖像處理方法的示意性流程圖；圖18是示出了圖像中的各個圖像單元的梯度值分布的示例的示意圖；及圖19是示出用于實現(xiàn)本公開的計算機(jī)的結(jié)構(gòu)的示例性框圖。
具體實施例方式下面參照附圖來說明本公開的實施例。在本公開的一個附圖或一種實施方式中描述的元素和特征可以與一個或更多個其它附圖或?qū)嵤┓绞街惺境龅脑睾吞卣飨嘟Y(jié)合。應(yīng)當(dāng)注意，為了清楚的目的，附圖和說明中省略了與本公開無關(guān)的、本領(lǐng)域普通技術(shù)人員已知的部件和處理的表示和描述。在進(jìn)行圖像處理時，通常可以對圖像進(jìn)行迭代濾波，以盡量濾除圖像中的噪聲。如果迭代的次數(shù)太少，則降噪效果可能不理想，而另一方面，如果迭代次數(shù)過多，則圖像中的有用信息也會被平滑，從而導(dǎo)致圖像模糊，影響后續(xù)的處理。本公開的一些實施例提供了能夠在圖像的迭代濾波過程中確定停止迭代濾波的適當(dāng)時機(jī)的圖像處理設(shè)備和方法。應(yīng)理解，根據(jù)本公開的實施例的設(shè)備和方法可用于對各種類型的圖像進(jìn)行降噪處理。例如，根據(jù)本公開的實施例的設(shè)備和方法可用于醫(yī)學(xué)圖像的降噪處理。這里所述的醫(yī)學(xué)圖像可以是根據(jù)利用醫(yī)療診斷成像裝置獲得的被檢測者的數(shù)據(jù)而形成的圖像。這里所述的醫(yī)療診斷裝置包括但不限于X射線成像診斷裝置、超聲波(UL)診斷成像裝置、計算機(jī)斷層掃描(CT)裝置、磁共振成像(MRI)診斷成像裝置、正電子發(fā)射斷層掃描(PositronEmission Tomography, PET)裝置等。圖7是示出了利用醫(yī)療診斷成像裝置獲取醫(yī)學(xué)圖像的示意圖。如圖7所示，對于每種成像設(shè)備(例如上述醫(yī)療診斷裝置)均有一些預(yù)定的或用戶自定義的掃描協(xié)議。以CT為例,可以有諸如心電同步掃描(cardiac synchronization)、腹部4相增強(qiáng)掃描(abdominal 4 phases enhance)、動態(tài)掃描等協(xié)議,如圖7中所示的協(xié)議
I、協(xié)議2........協(xié)議N(NSl);利用某個協(xié)議，可以獲得多個圖像序列(如圖7中所示的
序列I、序列2........序列N等)。因此，待處理的圖像可以是根據(jù)利用醫(yī)療診斷裝置、基
于相應(yīng)的協(xié)議獲取的患者數(shù)據(jù)所形成的醫(yī)學(xué)圖像。圖I是示出根據(jù)本公開的一個實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖，而圖2是示出根據(jù)該實施例的圖像處理方法的示意性流程圖。如圖I所示，圖像處理設(shè)備100可以包括濾波裝置101和迭代濾波停止裝置103。圖像處理設(shè)備100可以采用如圖2所示的方法進(jìn)行圖像濾波。下面參考圖2所示的方法流程來描述圖像處理設(shè)備100中的各個裝置的功能。濾波裝置101用于對輸入到該圖像處理設(shè)備的輸入圖像進(jìn)行迭代濾波(圖2的步驟202)。根據(jù)實際應(yīng)用場景，濾波裝置101 (即步驟202)可以采用任何適當(dāng)?shù)臑V波方法，如非線性擴(kuò)散(Non-Linear Diffusion, NLD)濾波方法、線性濾波方法、各向異性擴(kuò)散(Anisotropic Diffusion, ATD)濾波方法等等,這里不作限定。
濾波裝置101在每次濾波后將得到的經(jīng)濾波圖像輸出到迭代濾波停止裝置103，而迭代濾波停止裝置103用于根據(jù)每次濾波后得到的經(jīng)濾波圖像來判斷是否停止迭代濾波。具體地，迭代濾波停止裝置103可以根據(jù)迭代濾波過程中經(jīng)濾波圖像隨著迭代次數(shù)的增加而變化的速率來判斷是否停止迭代濾波(圖2的步驟206)。例如，如果迭代濾波停止裝置103判斷經(jīng)濾波圖像的變化速率下降為小于或等于預(yù)定的閾值，這就說明此時的濾波操作對圖像的影響變得足夠小，也就是說可被濾除的噪聲變得足夠少了，則可以指示濾波裝置101停止迭代濾波，否則，則指示濾波裝置101繼續(xù)進(jìn)行下一次迭代。這里所述的經(jīng)濾波圖像的“變化”是指當(dāng)前濾波得到的經(jīng)濾波圖像與之前濾波(“之前濾波”是指當(dāng)前迭代之前進(jìn)行的某一次迭代濾波(如上一次迭代等)，在此不作限定)得到的經(jīng)濾波圖像之間的差異。在對圖像的迭代濾波過程中，隨著迭代次數(shù)的增加，每次濾波得到的經(jīng)濾波圖像與之前濾波得到的經(jīng)濾波圖像相比的變化(即差異)會隨著迭代次數(shù)而變小。也就是說，隨著迭代次數(shù)的增加，經(jīng)濾波圖像的變化速率會變慢。圖I和圖2所示的設(shè)備和方法基于這一原理，根據(jù)經(jīng)濾波圖像的變化速率來確定停止迭代濾波的時間。利用這種方法或設(shè)備，可以有效地確定停止迭代濾波的適當(dāng)時間，從而在濾除圖像中的大部分噪聲的同時保證圖像中的有用信息不被模糊，另外，還可避免手工設(shè)置停止迭代濾波時間的麻煩，從而提高圖像處理的自動化程度。經(jīng)濾波圖像在迭代過程中的變化可以反映于其不同特征或特性，如每次迭代得到的經(jīng)濾波圖像與噪聲圖像的相關(guān)性在迭代過程中的變化、經(jīng)濾波圖像或噪聲圖像的某些圖像特征(如方差或標(biāo)準(zhǔn)差)的值在迭代過程中的變化等等，但不限于此，例如還可以是上述各種變化的數(shù)學(xué)組合等。下文中將參考圖3至圖15來描述在迭代濾波過程中根據(jù)經(jīng)濾波圖像的變化速率來判斷是否停止迭代濾波的圖像處理設(shè)備和方法的一些更為具體的實施例。圖3是示出根據(jù)一個具體實施例的圖像處理設(shè)備的結(jié)構(gòu)的示意性框圖，而圖4是示出根據(jù)該具體實施例的圖像處理方法的示意性流程圖。如上所述，在迭代濾波過程中，每次濾波得到的經(jīng)濾波圖像與之前濾波得到經(jīng)濾波圖像相比，其變化會變小，即變化速率會變慢。在迭代濾波的開始階段，每次濾波濾除的噪聲較多，而隨著迭代次數(shù)的增加，圖像中的噪聲越來越少，而圖像本身也逐漸平滑。濾波圖像的變化速率與濾波圖像的梯度成正比，隨著圖像本身逐漸平滑，圖像的梯度隨之變小，因此，經(jīng)濾波圖像的變化速率會變慢。相應(yīng)地，隨著經(jīng)濾波圖像中噪聲含量的變化，其與所濾除的噪聲圖像之間的相關(guān)性也會發(fā)生變化。也就是說，所述相關(guān)性的變化也反映了經(jīng)濾波圖像的變化，相關(guān)性的變化速率也反映了經(jīng)濾波圖像的變化速率。圖6 (A)、6⑶和6 (C)分別示出了這種相關(guān)性在迭代濾波過程中的示例性變化曲線。圖6(A)、6(B)和6(C)中的曲線分別表示在對三個不同圖像(在所圖示的示例中，三個圖像分別為采用UL、CT和MRI設(shè)備針對不同人體部位而獲取的醫(yī)學(xué)圖像)進(jìn)行迭代濾波過程中每次濾波得到的經(jīng)濾波圖像與相應(yīng)的噪聲圖像的相關(guān)系數(shù)隨時間(迭代次數(shù))的變化趨勢。如圖6(A)和6 (B)所示，在迭代濾波的開始階段，經(jīng)濾波圖像與噪聲圖像的相關(guān)系數(shù)的值呈下降趨勢，這表示越來越多的噪聲被濾除。經(jīng)過一段時間之后，相關(guān)系數(shù)的曲線到達(dá)谷底(圖中用“X”表示的谷點T1、T2)。在谷點之后，經(jīng)濾波圖像與噪聲的相關(guān)系數(shù)又呈現(xiàn)緩慢上升的趨勢，這表示在谷點之后，隨著迭代次數(shù)的增加，圖像中的有用信息也被逐漸平滑掉，也就是說一些有用信息也被濾除，或者說進(jìn)入了噪聲圖像，從而使得經(jīng)濾波圖像與濾除的噪聲圖像之間的相關(guān)性又有所增大。一種已有的方法即通過尋找上述谷點來確定停止迭代濾波的時間。在這種已有的方法中，將經(jīng)濾波圖像與噪聲圖像的相關(guān)系數(shù)最小的時間點確定為停止迭代濾波的時間點，相關(guān)描述可以參見 Pavel Mrdzek 的論文 “Selection of Optimal Stopping Timefor Nonlinear Diffusion Filtering，，(干丨J于“Scale-Space and Morphology in ComputerVision”，第三次國際會議，2001年)(下文中將該論文稱為相關(guān)文獻(xiàn)I)。在圖6(A)所示的情況下，存在一個谷點Tl，且該谷點的位置對應(yīng)于正確的停止時間，因此，在這種情況下，利用上述已有方法可以確定正確的停止時間。但是，上述已有方法并不適用于圖6 (B)和6(C)所示的情況例如，在圖6(B)所示的曲線中，出現(xiàn)谷點T2的位置明顯滯后于正確的停止時間，且在谷點T2之后，相關(guān)系數(shù)曲線保持平穩(wěn)(即在谷點T2之后，相關(guān)系數(shù)值幾乎沒有變化)，因此，在這種情況下，如果利用上述已有方法，則所確定的時間點為谷點T2，而不是正確的停止時間；又如，圖6(C)所示的曲線是單調(diào)下降的，沒有谷點，因此，在這種情況下，如果利用上述已有方法，則無法找到相關(guān)系數(shù)最小的時間點，從而無法確定正確的停止時間。本申請的發(fā)明人發(fā)現(xiàn)，如6(A)、6 (B)和6 (C)的三個曲線所示，不管是否有谷點，這三個曲線在迭代濾波的開始階段都是單調(diào)下降的，且下降的速率越來越慢。圖3和圖4所示的設(shè)備和方法即利用經(jīng)濾波圖像與噪聲圖像之間的相關(guān)性在迭代濾波過程中下降的速率(而不是相關(guān)性最小的時間點)來確定停止迭代濾波的時間。如圖3所示，圖像處理設(shè)備300可以包括濾波裝置301和迭代濾波停止裝置303，其中迭代濾波停止裝置303可以包括噪聲圖像獲取單元303-1、相關(guān)性獲取單元303-2和停止判斷單元303-3。圖像處理設(shè)備300可以采用如圖4所示的方法流程進(jìn)行圖像濾波。下面參考圖4所示的方法流程來描述圖像處理設(shè)備300中的各個裝置和/或單元的功能。濾波裝置301與圖I所示的濾波裝置101相似，用于采用任何適當(dāng)?shù)姆椒▽斎氲皆搱D像處理設(shè)備的輸入圖像進(jìn)行迭代濾波(圖4的步驟402)，這里不再重復(fù)。濾波裝置301在每次濾波后將得到的經(jīng)濾波圖像輸出到噪聲圖像獲取單元303-1。噪聲圖像獲取單元303-1用于在每次濾波后根據(jù)所得到的經(jīng)濾波圖像和輸入圖像來估算噪聲圖像(圖4中的步驟406-1)。作為一個示例，可以計算每次濾波后得到的經(jīng)濾波圖像與輸入圖像之間的差值圖像，作為對應(yīng)的噪聲圖像。設(shè)I。表示輸入圖像，Ist表示當(dāng)前濾波后得到的經(jīng)濾波圖像，則對應(yīng)的噪聲圖像Int可以通過下式來估算
Int = I0-Ist(I)應(yīng)理解，每次濾波后，可以采用任何適當(dāng)?shù)钠渌椒?如取根據(jù)式(I)計算得到的差值圖像的絕對值，即Int = IIt5-IstI，這里不一一列舉)來估算噪聲圖像，包括對上式(I)進(jìn)行任何形式的數(shù)學(xué)變化，因此本公開并不局限于上述示例。噪聲圖像獲取單元303-1將得到的噪聲圖像輸出到相關(guān)性獲取單元303-2，而相關(guān)性獲取單元303-2可以估算該噪聲圖像與當(dāng)前迭代得到的經(jīng)濾波圖像的相關(guān)性(圖4的步驟406-2)，并將估算結(jié)果輸出到停止判斷單元303-3。作為一個示例，可以計算噪聲圖像與經(jīng)濾波圖像之間的相關(guān)系數(shù)，以表征兩個圖像之間的相關(guān)性。例如，可以采用下面的公式來計算兩個圖像之間的相關(guān)系數(shù)
權(quán)利要求
1.一種圖像處理設(shè)備，包括 濾波裝置，用于對輸入圖像進(jìn)行迭代濾波；以及 迭代濾波停止裝置，用于根據(jù)所述迭代濾波過程中每次濾波后得到的經(jīng)濾波圖像隨迭代次數(shù)的變化速率來確定是否停止所述迭代濾波。
2.如權(quán)利要求I所述的圖像處理設(shè)備，其中，所述迭代濾波停止裝置包括 噪聲圖像獲取單元，用于在所述濾波裝置每次進(jìn)行濾波后根據(jù)所得到的經(jīng)濾波圖像和所述輸入圖像來估算噪聲圖像； 相關(guān)性獲取單元，用于估算所述噪聲圖像與所述經(jīng)濾波圖像的相關(guān)性；及停止判斷單元，用于根據(jù)所述相關(guān)性在所述迭代濾波過程中下降的速率來確定是否停止所述迭代濾波。
3.如權(quán)利要求2所述的圖像處理設(shè)備，其中，所述停止判斷單元被配置為 計算表示所述相關(guān)性在所述迭代濾波過程中下降的速率的特征值，并判斷所述特征值與預(yù)定閾值之間是否滿足預(yù)定的關(guān)系，若是，則停止所述迭代濾波。
4.如權(quán)利要求3所述的圖像處理設(shè)備，其中，所述相關(guān)性獲取單元被配置為計算所述噪聲圖像與所述經(jīng)濾波圖像之間的相關(guān)系數(shù)，以表征所述相關(guān)性。
5.如權(quán)利要求4所述的圖像處理設(shè)備，其中，所述特征值為表示所述相關(guān)系數(shù)在所述迭代濾波過程中隨迭代次數(shù)的增加而下降的曲線的傾斜度的值。
6.如權(quán)利要求3到5中任一項所述的圖像處理設(shè)備，其中，所述預(yù)定閾值為根據(jù)輸入圖像的特性而預(yù)先確定的值。
7.如權(quán)利要求I所述的圖像處理設(shè)備，其中，所述迭代濾波停止裝置包括 圖像特征計算單元，用于在每次濾波后計算經(jīng)濾波圖像的預(yù)定圖像特征值；及 停止判斷單元，用于根據(jù)所述預(yù)定圖像特征值隨著迭代次數(shù)的增加而變化的速率來確定是否停止所述迭代濾波。
8.如權(quán)利要求7所述的圖像處理設(shè)備，其中，所述迭代濾波停止裝置還包括 噪聲圖像獲取單元，用于在所述濾波裝置每次濾波后根據(jù)所得到的經(jīng)濾波圖像和所述輸入圖像來估算噪聲圖像，并且 其中，所述圖像特征計算單元被配置為計算所述噪聲圖像的所述預(yù)定圖像特征值，用于反映所述經(jīng)濾波圖像的預(yù)定圖像特征。
9.如權(quán)利要求7或8所述的圖像處理設(shè)備，其中，所述預(yù)定圖像特征為圖像的方差或標(biāo)準(zhǔn)差。
10.如權(quán)利要求7或8所述的圖像處理設(shè)備，其中，所述停止判斷單元被配置為 計算表示所述預(yù)定圖像特征在所述迭代濾波過程中變化的速率的特征值，并判斷所述特征值與預(yù)定閾值之間是否滿足預(yù)定的關(guān)系，若是，則停止所述迭代濾波。
11.如權(quán)利要求10所述的圖像處理設(shè)備，其中，所述特征值為表示所述預(yù)定圖像特征在所述迭代濾波過程中隨迭代次數(shù)的增加而變化的曲線的傾斜度的值。
12.如權(quán)利要求2所述的圖像處理設(shè)備，其中，所述迭代濾波停止裝置還包括 圖像特征計算單元，用于計算所述噪聲圖像的預(yù)定圖像特征值； 組合單元，用于組合所述相關(guān)性與所述噪聲圖像的預(yù)定圖像特征值，得到組合特征值，并且其中，所述停止判斷單元被配置為當(dāng)所述組合特征值在所述迭代濾波過程中達(dá)到極小值時停止所述迭代濾波。
13.如權(quán)利要求12所述的圖像處理設(shè)備，其中，所述預(yù)定圖像特征為圖像的方差或標(biāo)準(zhǔn)差。
14.如權(quán)利要求12或13所述的圖像處理設(shè)備，其中，所述圖像特征計算單元還被配置用于對計算得到的所述噪聲圖像的預(yù)定圖像特征值進(jìn)行變換，得到經(jīng)變換的預(yù)定圖像特征值，并且 其中，所述組合單元被配置為計算所述相關(guān)性與經(jīng)變換的預(yù)定圖像特征值的和。
15.如權(quán)利要求12或13所述的圖像處理設(shè)備，其中，所述相關(guān)性獲取單元被配置為計算所述噪聲圖像與所述經(jīng)濾波圖像之間的相關(guān)系數(shù)，以表征所述相關(guān)性。
16.如權(quán)利要求I到5中任一項所述的圖像處理設(shè)備，其中，所述濾波裝置被配置為基于非線性擴(kuò)散方法進(jìn)行濾波。
17.如權(quán)利要求16所述的圖像處理設(shè)備，還包括 梯度閾值確定裝置，用于在濾波裝置進(jìn)行每次濾波之前計算所述輸入圖像或上一次濾波得到的經(jīng)濾波圖像中的每一圖像單元的梯度值，并根據(jù)所述圖像中的所有圖像單元的梯度分布來確定梯度閾值， 其中，所述濾波裝置被配置為利用所述梯度閾值確定裝置確定的梯度閾值來基于非線性擴(kuò)散方法進(jìn)行濾波。
18.如權(quán)利要求17所述的圖像處理設(shè)備，其中，所述梯度閾值確定裝置被配置為 獲取與所述輸入圖像對應(yīng)的比例值，并按照所述比例值、根據(jù)所述梯度分布來計算所述梯度閾值，其中，圖像中梯度值小于所述梯度閾值的圖像單元的個數(shù)與所有圖像單元的個數(shù)的比接近所述比例值。
19.如權(quán)利要求18所述的圖像處理設(shè)備，其中，所述比例值為根據(jù)所述輸入圖像的特性而預(yù)先確定的值。
20.如權(quán)利要求I到5中任一項所述的圖像處理設(shè)備，其中，所述輸入圖像為根據(jù)通過醫(yī)療診斷裝置獲得的數(shù)據(jù)而形成的醫(yī)學(xué)圖像。
21.—種圖像處理方法，包括 對輸入圖像進(jìn)行迭代濾波；以及 根據(jù)所述迭代濾波過程中每次濾波后得到的經(jīng)濾波圖像隨迭代次數(shù)的變化速率來確定是否停止所述迭代濾波。
全文摘要
本公開提供了圖像處理設(shè)備和方法。該圖像處理設(shè)備可以包括濾波裝置，用于對輸入圖像進(jìn)行迭代濾波；以及迭代濾波停止裝置，用于根據(jù)所述迭代濾波過程中每次濾波后得到的經(jīng)濾波圖像隨迭代次數(shù)的變化速率來確定是否停止所述迭代濾波。
文檔編號G06T5/00GK102646265SQ20111004524
公開日2012年8月22日 申請日期2011年2月22日 優(yōu)先權(quán)日2011年2月22日
發(fā)明者山形仁, 張玭, 王艷麗, 馬紅楓 申請人:東芝醫(yī)療系統(tǒng)株式會社, 株式會社東芝