專利名稱::圖像處理裝置及其方法
技術(shù)領(lǐng)域:
:本發(fā)明涉及在收集被攝體的圖像的情況下,提高被收集的圖像數(shù)據(jù)的信號對噪聲之比(S/N比),并且有助于抑制上述收集的圖像的分辨率的降低的技術(shù)。
背景技術(shù):
:醫(yī)用圖像為,對被攝體(被檢查者)進(jìn)行包括X射線、γ射線的放射線、磁、超聲波等的用于測量的能量照射,對通過上述能量照射得到的上述被撿查者的透射信號或反射信號等由計(jì)算機(jī)進(jìn)行數(shù)據(jù)處理,在該被數(shù)據(jù)處理的結(jié)果中,上述被檢查者的內(nèi)部被圖像化。近年來,在上述醫(yī)用圖像的測量時,一方面由于X射線輻射等的測量而從保護(hù)被檢查者的角度要求降低對檢查者照射的能量。另一方面,即使在能量照射被降低也要求高像質(zhì)的圖像。雙方的要求有矛盾的關(guān)系,用于滿足上述要求的一個方法為通過抑制用于上述測量的能量而除去增加的噪聲信號,改善S/N比的技術(shù)?,F(xiàn)有的圖像處理方法在“專利文獻(xiàn)1”中被公開?!皩@墨I(xiàn)1”中具備下述階段,關(guān)于第一圖像數(shù)據(jù)的各像素,取得強(qiáng)調(diào)清晰度的清晰化圖像數(shù)據(jù)的階段;對上述第一圖像數(shù)據(jù)的各像素取得已平滑化的平滑化圖像數(shù)據(jù)的階段;對上述第一圖像數(shù)據(jù)的各像素算出邊緣度的階段;對所算出的上述邊緣度的最頻值按照提高上述平滑化圖像數(shù)據(jù)的融合比率的方式設(shè)定邊緣度-融合比率相關(guān)的階段;基于該邊緣度-融合比率相關(guān)對每個像素融合上述清晰化圖像數(shù)據(jù)和上述平滑化圖像數(shù)據(jù)并得到第二圖像數(shù)據(jù)的圖像數(shù)據(jù)融合階段。由此,可具有圖像數(shù)據(jù)的輪廓部強(qiáng)調(diào)和粒子狀噪聲降低的兩方面效果。專利文獻(xiàn)1特開2003-1312352號公報(bào)但是,由于“專利文獻(xiàn)1”所公開的技術(shù),考慮了算出上述邊緣度的階段在該算出時將尖峰噪聲(spikenoise)以外的本來的圖像的像素信息也作為噪音進(jìn)行識別,由此也會除去識別出的本來的像素信息的情況,因此從有可能導(dǎo)致圖像的分辨率降低的角度來說考慮不充分。
發(fā)明內(nèi)容本發(fā)明的圖像處理裝置,具備圖像生成機(jī)構(gòu)和混合圖像生成機(jī)構(gòu),其中所述圖像生成機(jī)構(gòu)對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像,所述混合圖像生成機(jī)構(gòu)生成將由該圖像生成機(jī)構(gòu)生成的平滑化圖像和清晰化圖像混合的混合圖像,所述圖像處理裝置具備分析量計(jì)算機(jī)構(gòu),其對分別包圍所述圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;控制機(jī)構(gòu),其按照通過該分析量計(jì)算機(jī)構(gòu)算出的分析量控制上述混合圖像生成機(jī)構(gòu)的平滑化圖像和清晰化圖像的混合比率;和顯示機(jī)構(gòu),其顯示通過由該控制機(jī)構(gòu)控制的混合圖像生成機(jī)構(gòu)輸出的混合圖像。由此,可提供一種抑制圖像的分辨率降低并可除去噪聲的圖像處理裝置。此外,本發(fā)明的圖像處理方法包括圖像生成步驟和混合圖像生成步驟,其中所述圖像生成步驟對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像,所述混合圖像生成步驟生成將由該圖像生成機(jī)構(gòu)生成的平滑化圖像和清晰化圖像混合的混合圖像,所述圖像處理方法包括分析量計(jì)算步驟,其對分別包圍所述圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;和控制步驟,其按照通過該分析量計(jì)算機(jī)構(gòu)算出的分析量控制上述混合圖像生成機(jī)構(gòu)的平滑化圖像和清晰化圖像的混合比率;和顯示步驟,其顯示通過由該控制機(jī)構(gòu)控制的混合圖像生成機(jī)構(gòu)輸出的混合圖像。由此,提供一種抑制圖像的分辨率的降低并可除去噪聲的圖像處理方法。發(fā)明效果通過本發(fā)明能夠抑制圖像的分辨率的降低并能除去噪聲。此外,由于能夠適當(dāng)?shù)貨Q定混合平滑化圖像和清晰化圖像的比率,因此能夠有效地實(shí)現(xiàn)圖像數(shù)據(jù)的輪廓的維持和噪聲降低兩方面,能夠生成具有期望的分析量的圖像。圖1為本發(fā)明的圖像處理裝置的示意結(jié)構(gòu)圖。圖2為實(shí)施方式1的從圖像讀入到圖像保存為止的流程圖。圖3為方差或標(biāo)準(zhǔn)偏差與混合權(quán)重之間的關(guān)系的說明圖。圖4為各部位的威納頻譜的說明圖。圖5為方差或標(biāo)準(zhǔn)偏差與混合權(quán)重之間的關(guān)系的說明圖。圖6為平滑化濾波器的流程圖。圖7為平滑化濾波器的示意圖。圖8為表示標(biāo)準(zhǔn)偏差和矩陣大小的關(guān)系的圖。圖9為實(shí)施方式2的從圖像讀入到圖像保存為止的流程圖。圖10為設(shè)定某ROI的源圖像的威納頻譜和期望的威納頻譜的說明圖。圖11為表示分析量為方差或標(biāo)準(zhǔn)偏差時的第二混合權(quán)重的初始值決定用查表的一例。圖12為對某ROI的源圖像的威納頻譜和期望的威納頻譜的差值進(jìn)行規(guī)一化的一例的圖。圖13為表示分析量為威納頻譜時的第二混合權(quán)重的初始值決定用查表的一例的圖。圖14為表示圖像噪聲顯眼的圖像例的圖。圖15為表示在圖14的圖像中應(yīng)用現(xiàn)有技術(shù)的圖像的一例的圖。圖16為表示在圖14的圖像中應(yīng)用本發(fā)明的圖像的一例的圖。圖17為表示具有某標(biāo)準(zhǔn)偏差值的現(xiàn)有圖像例的圖。圖18為表示在圖17的圖像中應(yīng)用本發(fā)明的圖像的一例的圖。具體實(shí)施例方式以下參照附圖對本發(fā)明相關(guān)的圖像診斷支援裝置的優(yōu)選實(shí)施方式進(jìn)行詳細(xì)說明。(實(shí)施方式1)圖1為表示本發(fā)明相關(guān)的圖像處理裝置的結(jié)構(gòu)例的框圖。本發(fā)明相關(guān)的圖像處理裝置具有圖像取得機(jī)構(gòu)1、與圖像取得機(jī)構(gòu)1可傳遞信號地連接的運(yùn)算處理機(jī)構(gòu)2;與運(yùn)算處理機(jī)構(gòu)2可傳遞信號地連接的圖像顯示機(jī)構(gòu)3。圖像取得機(jī)構(gòu)1根據(jù)作為處理對象的圖像種類而不同,例如如果成為處理對象的圖像為X射線透視圖像、CT圖像、MR圖像或超聲波圖像等的醫(yī)療用數(shù)字圖像數(shù)據(jù),則這些圖像取得裝置為X射線裝置、CT裝置、MR裝置或超聲波裝置等的醫(yī)療用圖像診斷裝置?;蛘?,在成為處理對象的圖像為由數(shù)碼相機(jī)所拍攝的圖像或由掃描儀讀取的圖像等的數(shù)字圖像數(shù)據(jù)的情況下,這些圖像取得裝置為數(shù)碼相機(jī)或掃描儀等的圖像取得機(jī)器。由于在這些圖像取得機(jī)器中以較少的入射光量取得的圖像的噪聲成分較多,因此本發(fā)明公開的方法在除去該噪聲的情況下也能采用。運(yùn)算處理機(jī)構(gòu)2為與圖像取得機(jī)構(gòu)1一體化地計(jì)算機(jī)或另外獨(dú)立地計(jì)算機(jī),具有平滑化圖像生成機(jī)構(gòu)4、清晰度圖像生成機(jī)構(gòu)5、分析量計(jì)算機(jī)構(gòu)6、圖像混合機(jī)構(gòu)7、保存機(jī)構(gòu)8、存儲機(jī)構(gòu)9以及輸入機(jī)構(gòu)10。平滑化圖像生成機(jī)構(gòu)4進(jìn)行處理對象源圖像上的噪聲去除。清晰度圖像生成機(jī)構(gòu)5對處理對象源圖像的邊緣進(jìn)行強(qiáng)化。分析量計(jì)算機(jī)構(gòu)6算出圖像的方差值、標(biāo)準(zhǔn)偏差值、威納頻譜值等。圖像混合機(jī)構(gòu)7基于分析量計(jì)算機(jī)構(gòu)6算出的分析量混合平滑化圖像和清晰化圖像。保存機(jī)構(gòu)8具有保存平滑化圖像、清晰度化圖像和混合圖像的功能和保存各種程序的功能,為硬盤等。存儲機(jī)構(gòu)9具有暫時保存平滑化圖像、清晰化圖像和混合圖像的功能,為存儲器等。輸入機(jī)構(gòu)10具有操作者輸入用于使圖像處理裝置動作的各種參數(shù)的功能,為鼠標(biāo)或鍵盤等。此外,圖像取得機(jī)構(gòu)1具有作為計(jì)算機(jī)的功能,因此至少包括數(shù)字信號處理器(DSP)、微處理器(MPU)以及中央處理單元(CPU)中的一種。圖像顯示機(jī)構(gòu)3為圖像取得機(jī)構(gòu)1、運(yùn)算處理機(jī)構(gòu)2或該任一個一體化或獨(dú)立的顯示器等的顯示裝置。圖2為表示從本發(fā)明相關(guān)的圖像處理裝置的實(shí)施方式1的數(shù)據(jù)讀入到輸出圖像的顯示的流程圖。步驟201中,圖像取得機(jī)構(gòu)讀入成為處理對象的圖像,運(yùn)算處理機(jī)構(gòu)2將上述所讀入的圖像存儲在存儲機(jī)構(gòu)9中?;蛘撸\(yùn)算處理機(jī)構(gòu)2將已保存在保存機(jī)構(gòu)8中的圖像讀入,運(yùn)算處理機(jī)構(gòu)2也可將上述讀入的圖像保存在保存機(jī)構(gòu)9中。此時讀入的圖像也可為通過X射線裝置、MR裝置、CT裝置、超聲波診斷裝置、數(shù)碼相機(jī)等新生成的數(shù)字圖像數(shù)據(jù),并不限于上述種類。在步驟202中,操作者(未圖示)選擇用于算出混合平滑化圖像和清晰化圖像的比率的源的分析量的各種參數(shù),采用輸入機(jī)構(gòu)10將該所選擇的各種參數(shù)輸入到運(yùn)算處理機(jī)構(gòu)2中。該參數(shù)為分析量的種類和矩陣大小。所謂分析量指統(tǒng)計(jì)學(xué)的物理量,用將矩陣內(nèi)的像素值作為基礎(chǔ)計(jì)算得到的值反映圖像的粒狀性。分析量的種類例如為矩陣內(nèi)的像素值的方差、標(biāo)準(zhǔn)偏差、維也納光譜以及差值中的至少一個。在此所述的差值表示將某關(guān)注的像素值和位于其周邊的像素的像素值相減后的值。矩陣尺寸指算出該分析量的像素的范圍,例如指定某像素周邊的3×3矩陣或5×5矩陣。在此,這些分析量的種類和矩陣大小的參數(shù)也可由操作者任意指定,也可由程序預(yù)先設(shè)定。另外,步驟202能夠省略預(yù)先設(shè)定上述參數(shù)的情況。在步驟203中,分析量計(jì)算機(jī)構(gòu)6基于由步驟202決定的參數(shù)算出該分析量。在此矩陣內(nèi)的像素值的方差、標(biāo)準(zhǔn)偏差以及威納頻譜的計(jì)算方法是公知的方法。步驟204中,圖像混合機(jī)構(gòu)7按照由步驟203算出的分析量計(jì)算混合平滑化圖像和傾斜化圖像的權(quán)重。混合權(quán)重的計(jì)算方法在后面敘述。對每一個像素決定混合權(quán)重W(x,y),對所有的像素重復(fù)步驟203~204直到該計(jì)算終止為止。在步驟205中,如果運(yùn)算處理機(jī)構(gòu)2對所有的像素結(jié)束運(yùn)算,則轉(zhuǎn)移到下一個步驟。在步驟206中,平滑化圖像生成機(jī)構(gòu)4對處理對象源圖像進(jìn)行噪聲除去處理。在噪聲除去處理中可使用移動平均濾波器或中間值濾波器等的公知的平滑化濾波器,但也可基于后述的分析量使用平滑化濾波器。在步驟207中,清晰化生成機(jī)構(gòu)5對處理對象的原圖像進(jìn)行邊緣強(qiáng)化處理。在邊緣強(qiáng)化處理中,也可使用采用拉普拉斯算子(二次微分)的清晰化濾波器等的公知的邊緣強(qiáng)化濾波器。在步驟208中,圖像混合機(jī)構(gòu)8基于由步驟204算出的混合權(quán)重W(x,y)根據(jù)式1對平滑化圖像AIMG(x,y)和清晰化圖像BIMG(x,y)進(jìn)行混合,生成混合圖像OIMG(x,y)。OIMG(x,y)=W(x,y)·AIMG(x,y)+(1-W(x,y))·BIMG(x,y)(式1)在步驟209中,運(yùn)算處理機(jī)構(gòu)2將步驟208中生成的混合圖像OIMG(x,y)顯示在圖像顯示機(jī)構(gòu)3中。在步驟210中,運(yùn)算處理機(jī)構(gòu)2讓操作者判斷是否需要保存混合圖像。由于該判斷,運(yùn)算處理機(jī)構(gòu)2在圖像顯示機(jī)構(gòu)3中顯示保存確認(rèn)信息,操作者按照該顯示判斷要否保存。在需要保存混合圖像的情況下,運(yùn)算處理機(jī)構(gòu)2在步驟211中保存于保存機(jī)構(gòu)8中。在上述動作說明中,以混合權(quán)重計(jì)算、平滑化圖像生成、清晰化圖像生成的順序記述了處理內(nèi)容,但這些處理的順序可不同。即步驟202~205之前的處理和步驟206、207也可任意地替換處理順序。該處理順序的替換,操作者通過輸入機(jī)構(gòu)10將該處理順序顯示在圖像顯示機(jī)構(gòu)3中并任意地進(jìn)行設(shè)定。此外,如果運(yùn)算處理機(jī)構(gòu)2的處理能力高,則對上述混合權(quán)重計(jì)算、上述平滑化圖像生成以及上述清晰化圖像生成進(jìn)行并行處理或流水線處理,也可進(jìn)一步實(shí)現(xiàn)運(yùn)算的高速化。接下來,對分析量為方差或標(biāo)準(zhǔn)偏差的情況的混合權(quán)重的計(jì)算方法的一例進(jìn)行說明。式2、式3對以某像素P(x,y)為中心的M×N矩陣的方差D(x,y)和標(biāo)準(zhǔn)偏差σ(x,y)進(jìn)行定義。D(x,y)=Σi=x-M/2x+M/2Σj=y-N/2y+N/2(P(i,j)-p-)2/(M·N)]]>(式2)σ(x,y)=D(x,y)=Σi=x-M/2x+M/2Σj=y-N/2y+N/2(P(i,j)-P‾)2/(M·N)]]>(式3)在此由下式定義以像素P(x,y)為中心的M×N矩陣內(nèi)的像素值的平均值。(式4)p-=Σi=x-M/2x+M/2Σj=y-N/2y+N/2P(i,j)/(M·N)]]>(式4)圖3為表示方差或標(biāo)準(zhǔn)偏差與混合權(quán)重的關(guān)系的曲線圖。該曲線圖對方差或標(biāo)準(zhǔn)偏差從0到THσ為止的范圍決定閾值THY~1的范圍。其決定方法也可為圖3中的黑線11那樣的線性或如圖3中的虛線12那樣的三角函數(shù)型的任一個。一般來說,矩陣內(nèi)的像素值的方差D(x,y)和標(biāo)準(zhǔn)偏差σ(x,y)有如果圖像噪聲大則變大的傾向。在此,如果方差D(x,y)和標(biāo)準(zhǔn)偏差σ(x,y)比規(guī)定值越大,越提高平滑化圖像的混合比率,則可發(fā)揮噪聲降低效果。另一方面,由于有急劇的濃度變化,因此在邊緣附近的矩陣的方差和標(biāo)準(zhǔn)偏差極端地變高。在此,邊緣附近的矩陣能夠解釋存儲某些構(gòu)造物、例如臟器的邊界。在這種情況下,如果降低平滑化圖像的混合比率,則使分辨率優(yōu)先。閾值THσ也可為根據(jù)經(jīng)驗(yàn)得到的邊緣附近的標(biāo)準(zhǔn)偏差或方差。另外,圖3中的THY的起點(diǎn)也可為0,但在THY為0的情況下,在存在急劇的濃度變換的區(qū)域中會有在清晰化圖像的部分中產(chǎn)生不自然的邊界線(二重線)等的情況。因此,為了防止該不自然的邊界線的產(chǎn)生,優(yōu)選THY不為0,為漸近于0的值。此外,矩陣的大小(M×N)也可為任意值,但優(yōu)選為維持邊緣或結(jié)構(gòu)物的清晰度并在除去噪聲時基本上為3×3或5×5的大小。如果矩陣尺寸的大小為3×3或5×5左右,則從運(yùn)算時間縮短的角度來看優(yōu)選。接下來,對分析量為威納頻譜時的混合權(quán)重的計(jì)算方法進(jìn)行說明。威納頻譜通過自相關(guān)函數(shù)的傅立葉變換求出,為在圖像的粒狀性的評價(jià)中一般使用的分析量,威納頻譜越低則粒狀度越好。作為例子在圖4中表示腹部CT圖像的各部位的威納頻譜。圖4的威納頻譜用32×32矩陣算出。部位21的矩陣內(nèi)存在急劇的濃度變化,因此解釋為矩陣內(nèi)包括骨等的輪廓。由于部位22和部位23在矩陣內(nèi)不存在急劇的濃度變化,因此解釋為在矩陣內(nèi)不包括輪廓或結(jié)構(gòu)物。此外,部位22與部位23相比,由于威納頻譜低,粒狀度較好,因此解釋為矩陣內(nèi)的噪聲較小。威納頻譜的形狀,由圖4可知,根據(jù)急劇的濃度變化的有無和噪聲的強(qiáng)弱而變化。由于在輪廓附近有急劇的濃度變化,因此低頻區(qū)域的威納頻譜變地極端高。在此,在低頻區(qū)域的威納頻譜比閾值WSTH高的矩陣的情況下,優(yōu)先分辨率,如圖5所示給平滑化圖像的混合權(quán)重賦予最低值THY。在低頻區(qū)域的威納頻譜比閾值WSTH低的情況下,優(yōu)先除去噪聲。具體地來說,威納頻譜如圖5所示,按照越大的像素越提高平滑化圖像的混合比率的方式?jīng)Q定混合權(quán)重。閾值WSTH為根據(jù)經(jīng)驗(yàn)得到的輪廓附近的威納頻譜。圖5的橫軸也可是例如圖4的F1所示那樣特定的頻帶的威納頻譜,或者也可是F2~F3的區(qū)間的平均值?;旌蠙?quán)重相對威納頻譜在THY~1的范圍中決定。該決定的方法也可是圖5中的黑線13所示那樣為線性或圖5中的虛線14的三角函數(shù)狀的任一個形狀。圖5中的THY也可為零,但在THY為0的情況下,在存在急劇的濃度變化的區(qū)域中會有清晰化圖像具有不自然的邊界線(二重線)的情況。因此,優(yōu)選THY不為0且為漸進(jìn)于0的值。由圖4可知,如果噪聲小,則威納頻譜變低,輪廓等的急劇的濃度剃度所引起的威納頻譜的變化具有在低頻區(qū)域中出現(xiàn)的傾向。因此,在將特定的頻帶的威納頻譜作為圖5的橫軸的情況下,確切地判定輪廓的有無,因此優(yōu)選將該特定的頻帶設(shè)定為高頻區(qū)域。此外,在特定的區(qū)間的威納頻譜的平均值為圖5的橫軸的情況下,能夠確切地判定有無輪廓,因此優(yōu)選特定區(qū)間設(shè)定為低~中頻區(qū)域。矩陣的大小(M×N)越小,威納頻譜的頻率分辨率越低,矩陣的大小越大則越增強(qiáng)使細(xì)小的輪廓不清楚的效果。因此,在分析量為威納頻譜的情況下,優(yōu)選矩陣的大小為16×16、32×32程度的矩陣大小。接下來,基于圖6對圖2的步驟206的平滑化圖像的生成方法進(jìn)行說明。圖6為基于本實(shí)施方式的分析量的平滑化濾波器的示意圖和處理流程圖。在步驟601中,分析量計(jì)算機(jī)構(gòu)6求出標(biāo)準(zhǔn)偏差。求出標(biāo)準(zhǔn)偏差的范圍優(yōu)選以像素P(x,y)為中心的3×3至5×5矩陣,但并不限于該例示的矩陣大小。在此所述的范圍與后面的矩陣大小不同,由于求得標(biāo)準(zhǔn)偏差,因此為暫定的范圍。步驟602中,運(yùn)算處理機(jī)構(gòu)2基于由步驟601算出的標(biāo)準(zhǔn)偏差,根據(jù)圖7和圖8決定作為用于平滑化的單位的矩陣大小。運(yùn)算處理機(jī)構(gòu)2如圖8所示,按照標(biāo)準(zhǔn)偏差越大越增大矩陣大小的方式?jīng)Q定。其中,在輪廓部等中標(biāo)準(zhǔn)偏差極端地高,因此在標(biāo)準(zhǔn)偏差比閾值THSD大的情況下,優(yōu)先分辨率,也可為最小矩陣大小(3×3)。在此,THSD為根據(jù)經(jīng)驗(yàn)得到的輪廓附近的標(biāo)準(zhǔn)偏差值。運(yùn)算處理機(jī)構(gòu)2也可將任意值賦予TH1~TH5,也可根據(jù)下式?jīng)Q定。(式5)TH1=THSD/NTH2=2·THSD/NTH3=3·THSD/NTH4=4·THSD/NTH5=5·THSD/N]]>(式5)式5中的N為0到THSD的領(lǐng)域的分割數(shù),圖7的情況下,該分割數(shù)為6。另外,在圖7中最小矩陣大小為3×3、最大矩陣大小為13×13,設(shè)從0到THSD的分割數(shù)為6,但最小以及最大矩陣大小、分割數(shù)并不限于這些,也可任意地設(shè)定。由此,包括噪聲較多的標(biāo)準(zhǔn)偏差大的部位以大的矩陣大小進(jìn)行分析量的運(yùn)算處理,噪聲少的標(biāo)準(zhǔn)偏差較小的像素以小的矩陣大小進(jìn)行分析量的運(yùn)算處理,此外存在輪廓等的極端的濃度變化的部位以最小的矩陣大小進(jìn)行分析量的運(yùn)算處理。由此,可實(shí)現(xiàn)分辨率的降低抑制的效果和噪聲除去效果兩方面。另外,在像素P(x,y)接近圖像的端部的情況下,在按照圖7決定的矩陣大小中會有矩陣超出圖像范圍的情況,但也可將該情況置換為不超出范圍的大小而進(jìn)行處理。在步驟603中,運(yùn)算處理機(jī)構(gòu)2在以由步驟602決定的大小規(guī)定的像素P(x,y)為中心的矩陣內(nèi)乘以基于分析量的權(quán)重并對像素值進(jìn)行平均化,決定濾波后的像素值。接下來,對基于分析量的濾波后的像素值的決定方法進(jìn)行說明。設(shè)某坐標(biāo)(x,y)的濾波前的像素值為P(x,y),濾波后的像素值為P’(x,y),則根據(jù)下式算出P’(x,y)。(式6)p'(x,y)=Σi=x-(N-1)/2x+(N-1)/2Σj=y-(N-1)/2y+(N-1)/2p(i,j)·W(i,j)Σi=x-(N-1)/2x+(N-1)/2Σj=y-(N-1)/2y+(N-1)/2W(i,j)]]>(式6)在此,N為矩陣大小,在上述的方法中對每一個像素進(jìn)行決定。此外,坐標(biāo)(i,j)為矩陣內(nèi)的各坐標(biāo),P(i,j)為坐標(biāo)(i,j)的像素值。在此i為到x-(N-1)/2~x+(N-1)/2的整數(shù),j為到y(tǒng)-(N-1)/2~y+(N-1)/2的整數(shù)。從而例如如果矩陣大小為3×3,則i、j為-1、0、1,如果矩陣大小為5×5,則i、j為-2、-1、0、1、2。還有,W(i,j)為對像素值進(jìn)行平均化時的權(quán)重,按照下式?jīng)Q定。(式7)W(i,j)={1-(P(i,j)-P(x,y)α·σ)}2]]>(式7)在式7中,σ為標(biāo)準(zhǔn)偏差,該值也可采用在圖6的步驟601中算出的值。此外,α為決定平滑化的程度的任意常數(shù),α越大平滑化的程度越強(qiáng),噪聲除去效果增強(qiáng)。在式7中,平均化的權(quán)重W(i,j)按照像素值的差或偏差(標(biāo)準(zhǔn)偏差)、即像素值的類似性決定,類似性越強(qiáng)權(quán)重越大。因此,根據(jù)式6像素值越為較相似的值,則稱以越大的權(quán)重并進(jìn)行平均化,從而與以往的平滑化濾波器相比能夠期待平均化所引起的分辨率降低的抑制效果。因此,運(yùn)算處理機(jī)構(gòu)2按照能夠起到充足的噪聲除去效果的方式?jīng)Q定α并由式7決定權(quán)重,如果按照式6進(jìn)行處理,則可達(dá)到分辨率的降低抑制的效果和充足的噪聲除去效果兩方面。在式7中,W(i,j)在像素值中按照類似性決定。此外,運(yùn)算處理機(jī)構(gòu)2也可根據(jù)式8那樣按照像素值的類似性和中心像素之間的距離來決定。在這種情況下也能期待與由式7決定權(quán)重的情況相同的效果。(式8)W(i,j)=1(x-i)2+(y-i)2{1-P(i,j)-P(x,y)α·σ}2]]>(式8)以上,如圖2所示的流程圖,例如通過圖1的運(yùn)算處理機(jī)構(gòu)2可自動地執(zhí)行。接下來,在適用本實(shí)施方式1的情況下,邊提高信號對噪聲比(SN比),邊將抑制分辨率的降低的混合圖像與現(xiàn)有技術(shù)對比來表示。圖14為信號噪聲比不良好且圖像噪聲顯眼的CT圖像的采樣圖像。圖15為對圖14的圖像適用現(xiàn)有技術(shù)(移動平均濾波器)并進(jìn)行噪聲減低的情況的采樣圖像。與此相對,圖16為對圖14的圖像適用本發(fā)明的實(shí)施方式1并進(jìn)行噪聲降低時的采樣圖像。對圖15和16進(jìn)行比較,在圖15中除去了圖14中所有的噪聲,但分辨率降低且圖像的清晰度也降低。與此相對,可知在圖16中,圖14中的噪聲明顯減少,分辨率也維持與原圖像相同的分辨率。從而,通過本實(shí)施方式,可達(dá)到除去噪聲和維持分辨率兩方面,能夠?qū)υ搱D像的再現(xiàn)提供最佳狀態(tài)。(實(shí)施方式2)圖9為表示從本發(fā)明相關(guān)的圖像處理裝置的實(shí)施方式2的數(shù)據(jù)讀入到輸出圖像的表示為止的流程圖。與實(shí)施方式1的不同點(diǎn)為根據(jù)對圖像內(nèi)的某關(guān)心區(qū)域內(nèi)的圖像的混合權(quán)重和對圖像全體的第二混合權(quán)重那樣適用多個混合權(quán)重。在步驟901中,運(yùn)算處理機(jī)構(gòu)2,讀入作為處理對象的圖像,保存在保存機(jī)構(gòu)9中。此時,運(yùn)算處理機(jī)構(gòu)2也可讀入己保存在保存機(jī)構(gòu)8中的圖像并保存,也可讀入通過CT裝置或數(shù)碼相機(jī)新生成的數(shù)字圖像數(shù)據(jù)這一情況與實(shí)施方式1中所說明的相同。在步驟902中,運(yùn)算處理機(jī)構(gòu)2將用于算出作為混合平滑化圖像和清晰化圖像的比率的基礎(chǔ)的分析量的各種參數(shù)顯示在圖像顯示機(jī)構(gòu)3中,操作者經(jīng)由輸入機(jī)構(gòu)10從該顯示的各種參數(shù)中進(jìn)行選擇。與實(shí)施方式1中所說明的相同,任一個都為反映圖像的粒狀性的分析量。矩陣大小指算出上述的分析量的范圍,例如運(yùn)算處理機(jī)構(gòu)2在用于算出分析量的參數(shù)中指定某像素周邊的3×3矩陣或5×5矩陣。此外,不需要矩陣的形狀為正方形這一點(diǎn)也與實(shí)施方式1相同。在此操作者也可任意指定用于算出分析量的參數(shù),也可預(yù)先設(shè)定。此外,在預(yù)先設(shè)定的情況下,可省略步驟902。在步驟903中,操作者經(jīng)由輸入機(jī)構(gòu)10在濾波前的圖像上設(shè)定關(guān)心區(qū)域(以下稱作ROI)。此時運(yùn)算處理機(jī)構(gòu)2也計(jì)算ROI內(nèi)的分析量,優(yōu)選在圖像顯示機(jī)構(gòu)3上與圖像一起顯示。在步驟904中,操作者經(jīng)由輸入機(jī)構(gòu)10指定應(yīng)取得濾波后的圖像的ROI內(nèi)的分析量的期望值。例如在分析量為方差或標(biāo)準(zhǔn)偏差的情況下,操作者介由輸入機(jī)構(gòu)10經(jīng)由用戶接口輸入任意的方差或標(biāo)準(zhǔn)偏差的值。該輸入例,作為處理對象的圖像為CT圖像,如果通過濾波生成相當(dāng)于以k倍的線量攝影的圖像的圖像,則也可指定當(dāng)前的1/k倍的方差值或當(dāng)前倍的標(biāo)準(zhǔn)偏差值。此外,在分析量為威納頻譜時,將如圖10所示的濾波前的威納頻譜25顯示在圖像顯示機(jī)構(gòu)3上,操作者也可經(jīng)由輸入機(jī)構(gòu)10在圖形上指定期望形狀的威納頻譜26。此時,操作者經(jīng)由輸入機(jī)構(gòu)10通過圖形用戶界面(GUI)對ROI內(nèi)的威納頻譜的形狀進(jìn)行變更或設(shè)定。在此所指定的值按照圖像中最應(yīng)重視的ROI內(nèi)具有期望的分析量的方式指定。在步驟905中,分析量計(jì)算機(jī)構(gòu)6基于以步驟902決定的參數(shù)算出以各像素為中心的矩陣內(nèi)的分析量。在此作為分析量的方差、標(biāo)準(zhǔn)偏差、威納頻譜的計(jì)算方法也可為公知的方法。在步驟906中,圖像混合機(jī)構(gòu)7基于由步驟905算出的分析量計(jì)算第一混合權(quán)重W1(x,y)。第一混合權(quán)重的計(jì)算方法也可與上述實(shí)施方式1的混合權(quán)重的計(jì)算方法相同。對每一個像素決定第一混合權(quán)重,運(yùn)算處理機(jī)構(gòu)2在步驟907中重復(fù)步驟905~906直到所有的像素計(jì)算結(jié)束為止。運(yùn)算處理機(jī)構(gòu)2,如果所有的像素計(jì)算結(jié)束,則進(jìn)入到下一個步驟。在步驟908中,運(yùn)算處理機(jī)構(gòu)2采用清晰化圖像生成機(jī)構(gòu)5,對處理對象源圖像進(jìn)行邊緣強(qiáng)化處理。邊緣強(qiáng)化處理與實(shí)施方式1中所述的方法相同。在步驟909中,運(yùn)算處理機(jī)構(gòu)2采用平滑化圖像生成機(jī)構(gòu)4,對作為處理對象源圖像進(jìn)行噪聲除去處理。噪聲除去處理也可與實(shí)施方式1所示的方法相同。在步驟910中,操作者基于由步驟903指定的源圖像的ROI內(nèi)的分析量,運(yùn)算處理機(jī)構(gòu)2基于由步驟904決定的混合圖像的ROI內(nèi)的期望的分析量值來決定第二混合權(quán)重W2的初始值。對第二混合權(quán)重W2的初始值的決定方法后述。在步驟911中,圖像混合機(jī)構(gòu)7基于由步驟906得到的第一混合權(quán)重和由步驟910得到的第二混合權(quán)重根據(jù)式9混合平滑化圖像AIMG(x,y)和清晰化圖像BIMG(x,y),生成混合圖像OIMG(x,y)。(式9)OIMG(x,y)=W2·W1(x,y)·AIMG(x,y)+(1-W2·W1(x,y))·BIMG(x,y)(式9)在步驟912中,運(yùn)算處理機(jī)構(gòu)2算出混合圖像OIMG(x,y)中的ROI內(nèi)的分析量。在步驟913中,運(yùn)算處理機(jī)構(gòu)2判定由步驟912算出的ROI內(nèi)的分析量與步驟904所指定的值是否一致,如果一致,則進(jìn)入到下一個步驟。在步驟914中,不一致的情況下,如果混合權(quán)重沒有達(dá)到極限(0或1),則運(yùn)算處理機(jī)構(gòu)2邊適當(dāng)變更第二混合權(quán)重的值,邊重復(fù)步驟910~912的處理直到與步驟904所指定的期望的值一致為止,在混合權(quán)重達(dá)到極限的情況下,重復(fù)步驟909~912的處理直到與步驟904所指定的期望的值一致為止。在此所述的平滑化水平(level)是指在平滑化濾波器中采用移動平均濾波器或中間值濾波器的情況下矩陣大小的事,在采用基于上述的分析量的平滑化濾波器的情況下指式7中的α。運(yùn)算處理機(jī)構(gòu)2重復(fù)步驟910~912或909~912的處理,也有與步驟904所指定的規(guī)定值不完全一致的情況。此時,運(yùn)算處理機(jī)構(gòu)2也可重復(fù)步驟910~912或909~912的處理直到由步驟904指定的期望值之間的誤差收斂在任意的范圍內(nèi)(例如5%)。此外,運(yùn)算處理機(jī)構(gòu)2也可預(yù)先固定重復(fù)步驟910~912或909~912的次數(shù),按照其中與期望的值之間的誤差最小的方式生成混合圖像。在步驟915中,運(yùn)算處理機(jī)構(gòu)2將由步驟912生成的混合圖像OIMG(x,y)顯示在圖像顯示機(jī)構(gòu)3中。在步驟916中,運(yùn)算處理機(jī)構(gòu)2盤點(diǎn)是否需要保存混合圖像。在需要保存的情況下,運(yùn)算處理機(jī)構(gòu)2采用步驟917保存在圖像保存機(jī)構(gòu)8中。接下來,對第二混合權(quán)重的初始值的決定方法進(jìn)行說明。首先對分析量為分散或標(biāo)準(zhǔn)偏差的情況進(jìn)行說明。在對具有某方差或標(biāo)準(zhǔn)偏差的源圖像以某比率混合由特定的平滑化濾波器生成的平滑化圖像和由特定的清晰化濾波器生成的清晰化圖像的情況下,混合圖像中的方差或標(biāo)準(zhǔn)偏差可唯一決定。將上述情況總結(jié)在表格中,能夠生成如圖11所示的查表(lookuptable)。對采用該查表的第二混合權(quán)重的初始值的決定方法進(jìn)行說明。例如,在圖9的步驟903中得到源圖像中的ROI內(nèi)的方差或標(biāo)準(zhǔn)偏差的值,在圖9的步驟904中指定混合圖像的ROI內(nèi)的期望的方差或標(biāo)準(zhǔn)偏差的值。將在圖11的查表的縱軸上ROI內(nèi)的方差或標(biāo)準(zhǔn)偏差的值一致的位置和橫軸上混合圖像的ROI內(nèi)的期望的方差或標(biāo)準(zhǔn)偏差的值一致的位置的交點(diǎn)的某值設(shè)為第二混合權(quán)重的初始值。另外,圖11中采用相對方差或標(biāo)準(zhǔn)偏差標(biāo)本間隔0.5的查表,但標(biāo)本間隔并不限于此,也可為任意的大小。接下來,對分析量為威納頻譜時的第二混合權(quán)重的初始值的決定方法進(jìn)行說明。將源圖像的威納頻譜和期望的威納頻譜相減,如果用1進(jìn)行規(guī)一化,則變?yōu)閳D12那樣。在對具有某威納頻譜的源圖像,以某比率混合由特定的平滑化濾波器生成的平滑化圖像和由特定的清晰化濾波器生成的清晰化圖像的情況下,源圖像和混合圖像的威納頻譜的差值的各頻率成分值可唯一地決定。因此,能夠生成圖13那樣的查表。在圖9的步驟803中,可得到源圖像的ROI內(nèi)的威納頻譜,在圖9的步驟804中指定混合圖像的期望的威納頻譜,因此將對它們進(jìn)行減法運(yùn)算并規(guī)一化的兩個頻帶27、28適用于圖13的查表中來決定第二混合權(quán)重的初始值。在圖13中采用相對威納頻譜的差標(biāo)本間隔0.1的查表,但標(biāo)本間隔并不限于此,可為任意的大小。以上圖9中所示的流程圖例如可通過圖1的運(yùn)算處理機(jī)構(gòu)2自動地執(zhí)行。接下來,與現(xiàn)有例的情況對比表示適用本實(shí)施方式2,一邊提高信號噪聲比(S/N比)一邊抑制分辨率的降低的結(jié)果。圖17為ROI31、32內(nèi)的標(biāo)準(zhǔn)偏差為約8的CT圖像的樣品圖像。圖18為在圖17的圖像中適用本發(fā)明的實(shí)施方式2的樣品圖像且將標(biāo)準(zhǔn)偏差設(shè)定為4。根據(jù)圖18可知,能夠得到一邊抑制分辨率的降低一邊降低噪聲的期望的分析量的圖像。通過本發(fā)明,在提高數(shù)字圖像數(shù)據(jù)的信號對噪聲比的情況下,具有一邊抑制分辨率的降低一邊可降低圖像噪聲的效果。此外,也具有在期望的關(guān)心區(qū)域內(nèi)能夠生成具有期望的分析量的圖像的效果。以上所說明的實(shí)施方式只是用于說明本發(fā)明的一例,分別單獨(dú)實(shí)施它們、組合實(shí)施都屬于本發(fā)明的權(quán)利范圍。產(chǎn)業(yè)上的利用可能性本發(fā)明的圖像處理裝置及其方法,在收集被攝體的圖像的情況下,提高被收集的圖像數(shù)據(jù)的信號對噪聲之比(S/N比),并且抑制上述收集的圖像的分辨率的降低。權(quán)利要求1.一種圖像處理裝置,具備圖像生成機(jī)構(gòu)和混合圖像生成機(jī)構(gòu),其中所述圖像生成機(jī)構(gòu)對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像,所述混合圖像生成機(jī)構(gòu)生成將該圖像生成機(jī)構(gòu)所生成的平滑化圖像和清晰化圖像混合的混合圖像,所述圖像處理裝置具備分析量計(jì)算機(jī)構(gòu),其對分別包圍所述圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;和控制機(jī)構(gòu),其按照通過該分析量計(jì)算機(jī)構(gòu)算出的分析量控制上述混合圖像生成機(jī)構(gòu)的平滑化圖像和清晰化圖像的混合比率;和顯示機(jī)構(gòu),其顯示通過由該控制機(jī)構(gòu)控制的混合圖像生成機(jī)構(gòu)輸出的混合圖像。2.根據(jù)權(quán)利要求1所述的圖像處理裝置,其特征在于,所述圖像數(shù)據(jù),為包括X線透視圖像、CT圖像、MR圖像、超聲波圖像的醫(yī)療用數(shù)字圖像數(shù)據(jù)或由包括數(shù)碼相機(jī)、掃描儀的圖像取得設(shè)備取得的數(shù)字圖像數(shù)據(jù)。3.根據(jù)權(quán)利要求1所述的圖像處理裝置,其特征在于,通過所述分析量計(jì)算機(jī)構(gòu)算出的分析量為圖像的方差值、標(biāo)準(zhǔn)偏差值、威納頻譜、關(guān)注像素值和位于其周邊的像素的像素值的差值中的至少一種。4.根據(jù)權(quán)利要求1所述的圖像處理裝置,其特征在于,還具備任意地設(shè)定所述圖像混合處理的平滑化圖像和清晰化圖像之間的混合比率、所述平滑化圖像生成以及所述清晰化圖像生成的處理順序的機(jī)構(gòu)。5.根據(jù)權(quán)利要求4所述的圖像處理裝置,其特征在于,所述處理順序生成機(jī)構(gòu)按照并行進(jìn)行所述混合比率計(jì)算、所述平滑化圖像生成以及所述清晰化圖像生成的方式進(jìn)行設(shè)定。6.根據(jù)權(quán)利要求1所述的圖像處理裝置,其特征在于,所述圖像混合處理的平滑化圖像和清晰化圖像之間的混合比率為線性或三角函數(shù)型的任一種。7.根據(jù)權(quán)利要求1所述的圖像處理裝置,其特征在于,所述圖像平滑化機(jī)構(gòu),按照圖像數(shù)據(jù)內(nèi)的標(biāo)準(zhǔn)偏差對所述圖像內(nèi)的各像素的每一個變更作為平滑化處理的單位的平滑化矩陣大小,按照該平滑化矩陣內(nèi)的像素值的類似性以及距離中心像素的距離的至少一個變更平均化權(quán)重,并對平滑化矩陣內(nèi)的多個像素值進(jìn)行平均化。8.根據(jù)權(quán)利要求7所述的圖像處理裝置,其特征在于,還具備自動調(diào)整的機(jī)構(gòu),其按照所述混合圖像的某關(guān)心區(qū)域內(nèi)的分析量自動調(diào)整由所述圖像平滑化機(jī)構(gòu)的平滑化矩陣大小、平均化權(quán)重的系數(shù)構(gòu)成的平滑化水平以及圖像混合機(jī)構(gòu)混合平滑圖像和清晰化圖像的比率中的至少一方,以使混合圖像的期望的關(guān)心區(qū)域內(nèi)的分析量達(dá)到期望的值。9.根據(jù)權(quán)利要求8所述的圖像處理裝置,其特征在于,在所述顯示機(jī)構(gòu)中顯示所述關(guān)心區(qū)域的威納頻譜,還具備可任意地變更并設(shè)定威納頻譜的形狀的圖形用戶界面。10.一種圖像處理方法,包括圖像生成步驟和混合圖像生成步驟,其中所述圖像生成步驟對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像,所述混合圖像生成步驟生成將由該圖像生成機(jī)構(gòu)生成的平滑化圖像和清晰化圖像混合的混合圖像,所述圖像處理方法包括分析量計(jì)算步驟,其對分別包圍所述圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;和控制步驟,其按照通過該分析量計(jì)算機(jī)構(gòu)算出的分析量控制上述混合圖像生成機(jī)構(gòu)的平滑化圖像和清晰化圖像的混合比率;和顯示步驟,其顯示通過由該控制機(jī)構(gòu)控制的混合圖像生成機(jī)構(gòu)輸出的混合圖像。11.根據(jù)權(quán)利要求10所述的圖像處理方法,其特征在于,所述圖像數(shù)據(jù),為包括X線透視圖像、CT圖像、MR圖像、超聲波圖像的醫(yī)療用數(shù)字圖像數(shù)據(jù)或由包括數(shù)碼相機(jī)、掃描儀的圖像取得設(shè)備取得的數(shù)字圖像數(shù)據(jù)。12.根據(jù)權(quán)利要求10所述的圖像處理方法,其特征在于,通過所述分析量計(jì)算機(jī)構(gòu)算出的分析量為圖像的方差值、標(biāo)準(zhǔn)偏差值、威納頻譜、關(guān)注像素值和位于其周邊的像素的像素值的差值中的至少一種。13.根據(jù)權(quán)利要求10所述的圖像處理方法,其特征在于,還具備任意地設(shè)定所述圖像混合處理的平滑化圖像和清晰化圖像之間的混合比率、所述平滑化圖像生成以及所述清晰化圖像生成的處理順序的步驟。14.根據(jù)權(quán)利要求13所述的圖像處理方法,其特征在于,所述處理順序生成步驟按照并行進(jìn)行所述混合比率計(jì)算、所述平滑化圖像生成以及所述清晰化圖像生成的方式進(jìn)行設(shè)定。15.根據(jù)權(quán)利要求10所述的圖像處理方法,其特征在于,所述圖像混合處理的平滑化圖像和清晰化圖像之間的混合比率為線性或三角函數(shù)型的任一種。16.根據(jù)權(quán)利要求10所述的圖像處理方法,其特征在于,所述圖像平滑化步驟,按照圖像數(shù)據(jù)內(nèi)的標(biāo)準(zhǔn)偏差對所述圖像內(nèi)的各像素的每一個變更作為平滑化處理的單位的平滑化矩陣大小,按照該平滑化矩陣內(nèi)的像素值的類似性以及距離中心像素的距離的至少一個變更平均化權(quán)重,并對平滑化矩陣內(nèi)的多個像素值進(jìn)行平均化。17.根據(jù)權(quán)利要求16所述的圖像處理方法,其特征在于,還具備自動調(diào)整的步驟,其按照所述混合圖像的某關(guān)心區(qū)域內(nèi)的分析量自動調(diào)整由所述圖像平滑化機(jī)構(gòu)的平滑化矩陣大小、平均化權(quán)重的系數(shù)構(gòu)成的平滑化水平以及圖像混合機(jī)構(gòu)混合平滑圖像和清晰化圖像的比率中的至少一方,以使混合圖像的期望的關(guān)心區(qū)域內(nèi)的分析量達(dá)到期望的值。18.根據(jù)權(quán)利要求17所述的圖像處理方法,其特征在于,在所述顯示機(jī)構(gòu)中顯示所述關(guān)心區(qū)域的威納頻譜,還具備可任意地變更并設(shè)定威納頻譜的形狀的圖形用戶界面。19.一種圖像處理方法,其特征在于,包括分析量計(jì)算步驟,其對分別包圍圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;按照由該分析量計(jì)算機(jī)構(gòu)算出的分析量對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像的步驟;混合所述平滑化圖像和清晰化圖像的步驟;和顯示通過該步驟混合的混合圖像的顯示步驟。全文摘要本發(fā)明的圖像處理裝置具備對圖像數(shù)據(jù)內(nèi)的至少一部分分別生成平滑化圖像和清晰化圖像的圖像生成機(jī)構(gòu)和生成將由該圖像生成機(jī)構(gòu)生成的平滑化圖像和清晰化圖像混合的混合圖像的混合圖像生成機(jī)構(gòu),所述圖像處理裝置具備分析量計(jì)算機(jī)構(gòu),其對分別包圍所述圖像數(shù)據(jù)內(nèi)的各像素的任意大小的每一個矩陣算出該矩陣內(nèi)的多個像素值相關(guān)的分析量;控制機(jī)構(gòu),其按照通過該分析量計(jì)算機(jī)構(gòu)算出的分析量控制上述混合圖像生成機(jī)構(gòu)的平滑化圖像和清晰化圖像的混合比率;和顯示機(jī)構(gòu),其顯示通過由該控制機(jī)構(gòu)控制的混合圖像生成機(jī)構(gòu)輸出的混合圖像。由此,能夠抑制圖像的分辨率的降低并能除去噪聲。文檔編號G06T5/00GK1953709SQ200580015769公開日2007年4月25日申請日期2005年5月17日優(yōu)先權(quán)日2004年5月19日發(fā)明者尾見康夫,宮崎靖,三竹毅申請人:株式會社日立醫(yī)藥