圖像處理設(shè)備、圖像處理方法和圖像處理程序的制作方法
【專利摘要】[目標(biāo)]為了恰當(dāng)確定用于抑制包括周期性圖樣(例如，柵格圖像)在內(nèi)的放射線照相圖像中的周期性圖樣的通帶特性。[解決方案]通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與周期性圖樣相對應(yīng)的頻譜的峰位置(pf)和峰寬(pw)。基于檢測到的峰位置(pf)和檢測到的峰寬(pw)，確定一維濾波器的通帶特性(F)，該一維濾波器對圖像信號中與周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。
【專利說明】圖像處理設(shè)備、圖像處理方法和圖像處理程序
【技術(shù)領(lǐng)域】
[0001]本發(fā)明涉及對圖像信號中與周期性圖樣相對應(yīng)的空間頻率分量加以抑制的圖像處理設(shè)備、方法和程序。具體的，本發(fā)明涉及對放射線照相中使用的柵格所引起的放射線照相圖像中的周期性圖樣加以抑制的圖像處理設(shè)備、方法和程序。
【背景技術(shù)】
[0002]傳統(tǒng)上，使用了蓄積性熒光體(輝盡性熒光體)。當(dāng)使用放射線(X射線、α射線、β射線、Y射線、電子束、紫外線等)來輻照蓄積性熒光體時，一部分放射線能量被存儲在蓄積性熒光體中。之后，當(dāng)向蓄積性熒光體輸出激發(fā)光(例如可見光和激光束)時，蓄積性熒光體發(fā)射與其中存儲的放射線能量相對應(yīng)的光致發(fā)光。例如，在CR(計算放射線照相)中廣泛使用采用蓄積性熒光體的放射線照相圖像讀出設(shè)備。使用已穿透對象(例如，人體)的放射線來輻照蓄積性熒光片，在該蓄積性熒光片中，在襯底上沉積蓄積性熒光體，且在蓄積性熒光片上臨時存儲和記錄放射線照相信息。向蓄積性熒光片輸出激發(fā)光(例如，激光束)，以引起光致發(fā)光。此外，針對光致發(fā)光執(zhí)行光電轉(zhuǎn)換，以獲得圖像信號。
[0003]此處，當(dāng)在蓄積性熒光片等上對對象的放射線照相圖像進(jìn)行成像和記錄時，在一些情況下通過在對象和片之間放置柵格來執(zhí)行成像，使得由對象所散射的放射線不輻照該片。例如，在柵格中以約4行/_的狹窄間距來交替布置鉛等(其不允許放射線從中穿過)和鋁(其傾向于允許放射線從中穿過)等。當(dāng)通過使用柵格來執(zhí)行放射線照相時，由對象散射的放射線不傾向于輻照該片。因此，有可能改進(jìn)對象的放射線照相圖像的對比度。然而，由于在圖像中包括柵格圖像，因此存在圖像質(zhì)量劣化的問題。此外，當(dāng)對包括柵格圖像在內(nèi)的圖像的大小進(jìn)行放大或縮小時，取決于放大或縮小比率，會出現(xiàn)由折疊所造成的混疊(aliasing)。此外，當(dāng)混疊與柵格圖像等的空間頻率重疊時，生成了狹窄的條紋圖樣(莫爾條紋)，且對重新生成的圖像的觀察變得困難。
[0004]專利文獻(xiàn)I公開了用于移除柵格圖像的方法，作為用于從圖像中移除這種周期性圖樣的方法。在該方法中，在布置柵格圖像的條紋圖樣的方向上，針對圖像執(zhí)行一維高通濾波處理。此外，在與柵格圖像的條紋圖樣相平行的方向上，針對圖像執(zhí)行一維低通濾波處理。因此，從原始圖像中提取出與柵格圖像相對應(yīng)的空間頻率分量，且從原始圖像中減去所提取的空間頻率分量。
[0005]此外，專利文獻(xiàn)2 公開了用于獲得圖像的方法，其中，抑制了由柵格引起的偽像。在該方法中，通過二維或一維傅立葉變換，對圖像數(shù)據(jù)執(zhí)行頻率分析，且獲得莫爾條紋的峰位置(頻率)、峰高、半值全寬(HWFM)、峰值處的總能量、柵格的方向等等。此外，以基于圖像的像素大小、峰位置和峰值處的能量確定的內(nèi)核大小，執(zhí)行平滑濾波。
[0006]此外，專利文獻(xiàn)3公開了圖像處理方法。在該方法中，針對圖像數(shù)據(jù)執(zhí)行二維傅立葉變換，以檢測與周期性圖樣相對應(yīng)的空間頻率分量。此外，通過使用僅移除檢測到的空間頻率分量的濾波器來執(zhí)行濾波。
[0007]相關(guān)技術(shù)文獻(xiàn)[0008]專利文獻(xiàn)
[0009]專利文獻(xiàn)1:日本未審專利公開N0.2003-150954
[0010]專利文獻(xiàn)2:美國專利N0.6269176
[0011]專利文獻(xiàn)3:日本未審專利公開N0.2009-195512

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

[0012]此處，需要在盡可能不移除表示對象的分量的情況下移除周期性圖樣(例如，莫爾條紋)，以獲得適用于基于圖像的診斷的高質(zhì)量圖像。然而，在專利文獻(xiàn)I所公開的方法中，與周期性圖樣相對應(yīng)的空間頻率分量(例如，柵格圖像的條紋圖樣)是通過執(zhí)行一維高通濾波處理來檢測的。因此，如果周期性圖樣未出現(xiàn)在高頻帶中，而是出現(xiàn)在低頻帶中，則不可能檢測到期望的周期性圖樣。因此，不可能采用專利文獻(xiàn)I中公開的方法。
[0013]此外，在專利文獻(xiàn)2中公開的方法中，可以通過執(zhí)行平滑濾波處理來抑制局部化的噪聲分量。然而，該方法不僅從圖像中移除檢測到的與周期性圖樣相對應(yīng)的空間頻率分量。因此，不可能集中移除包括周期性圖樣的空間頻率分量。
[0014]在專利文獻(xiàn)3所公開的方法中，由于通過二維傅立葉變換來檢測與周期性圖樣相對應(yīng)的空間頻率分量，因此計算的負(fù)荷很高。因此，需要使用較低計算負(fù)荷來確定與周期性圖樣相對應(yīng)的空間頻率分量的通帶特性的方法。 [0015]考慮到前述事實(shí)，本發(fā)明的目的是提供圖像處理設(shè)備、圖像處理方法以及圖像處理程序，即使在中頻帶或低頻帶中出現(xiàn)與周期性圖樣相對應(yīng)的頻率分量，它們也可以使用低計算負(fù)荷來恰當(dāng)?shù)卮_定通帶特性。
[0016]本發(fā)明的圖像處理設(shè)備是一種圖像處理設(shè)備，包括:
[0017]周期性圖樣檢測裝置，通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及
[0018]通帶特性確定裝置，基于檢測到的峰位置和檢測到的峰寬，來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)is抑制。
[0019]根據(jù)本發(fā)明的圖像處理方法是一種圖像處理方法，包括以下步驟:
[0020]通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及
[0021]基于檢測到的峰位置和檢測到的峰寬，來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。
[0022]本發(fā)明的圖像處理程序是一種圖像處理程序，其使得計算機(jī)擔(dān)當(dāng):
[0023]周期性圖樣檢測裝置，通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及
[0024]通帶特性確定裝置，基于檢測到的峰位置和檢測到的峰寬，來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)is抑制。
[0025]術(shù)語“周期性圖樣”指 在原始圖像中包括的具有周期性圖樣的噪聲。例如，周期性圖樣指當(dāng)通過使用柵格在蓄積性熒光片上對放射線照相圖像進(jìn)行成像時在原始圖像中包括的柵格圖像、莫爾條紋等。
[0026]此外，術(shù)語“通帶特性”指允許預(yù)定頻帶通過或?qū)ζ溥M(jìn)行移除的濾波器的特性。通帶特性包括帶通特性(帶通特性)和頻帶移除特性(帶阻特性)。
[0027]此外，只要峰位置表示頻譜的峰頻的代表值，則可以定義任意峰位置。類似地，只要峰寬表示頻譜的峰寬的代表值，則可以定義任意峰寬。
[0028]只要通帶特性包括與周期圖案相對應(yīng)的頻率分量的峰位置，且峰位置在被包括在要通過的頻帶中，通帶特性確定裝置就可以確定任何類型的通帶特性，且基于峰位置和峰寬可確定通帶特性。
[0029]期望帶通寬度寬于或等于峰寬，以移除與周期性圖樣相對應(yīng)的頻率分量。例如，期望通帶特性確定裝置以下述方式來確定通帶特性:包括峰位置且通帶寬度隨著峰寬更寬而更寬。
[0030]在通帶特性中，低頻帶包括很多對象分量。因此，在低頻帶中，將對象分量和與周期性圖樣相對應(yīng)的頻率分量一起移除的風(fēng)險隨著濾波器的帶通更寬而更高。因此，在低頻帶中，期望帶通特性包括周期性圖樣的峰且通帶的寬度較窄。例如，通帶特性確定裝置可以用下述方式來確定通帶特性:通帶寬度隨著峰位置位于更低頻帶而更窄。
[0031]此外，通帶特性確定裝置可以通過由峰位置和峰寬來限定的正態(tài)分布函數(shù)來確定通帶特性。例如，可以通過使用峰位置作為平均值，且使用半值全寬作為標(biāo)準(zhǔn)差，來定義正態(tài)分布函數(shù)?？梢詫⒁韵率龇绞椒糯蠡蚩s小的正態(tài)分布曲線確定為通帶特性:正態(tài)分布函數(shù)的峰對應(yīng)于頻率響應(yīng)中的I。
[0032]通帶特性確定裝置可以通過矩形函數(shù)來確定通帶特性，所述矩形函數(shù)具有比峰寬更寬或相等的寬度，且該矩形函數(shù)的中心是峰位置。
[0033]此外，通帶特性確 定裝置可以基于峰高來確定通帶特性。這是因?yàn)楫?dāng)頻譜的量值小時，即使不對包括與周期性圖樣相對應(yīng)的頻率分量在內(nèi)的整個頻帶進(jìn)行抑制(與周期性圖樣相對應(yīng)的頻率分量的一部分保留在圖像信號中)，依然有可能以下述方式來抑制周期性圖樣:質(zhì)量足以用作診斷用圖像。例如，可以用下述方式來確定通帶特性:帶通濾波器的帶通頻率響應(yīng)隨著峰高越高而越高,且?guī)V波器的帶通頻率響應(yīng)隨著峰高越低而越低。
[0034]期望根據(jù)本發(fā)明的圖像處理設(shè)備還包括:抑制裝置，通過基于所確定的通帶特性執(zhí)行一維濾波，來對所述圖像信號中的周期性圖樣進(jìn)行抑制。
[0035]此外，根據(jù)本發(fā)明的圖像處理設(shè)備中的通帶特性確定裝置可以在所述通帶特性中設(shè)置通帶寬度的上限或下限。這是因?yàn)橛蓛H通過(或阻止)極窄頻帶的帶通(或帶阻)濾波器進(jìn)行的濾波導(dǎo)致在圖像中生成偽像等，且這是不想要的。此外，通過過寬頻帶的帶通(帶阻)濾波器進(jìn)行的濾波是不想要的，因?yàn)榇嬖趯⒃诨趫D像的診斷中重要的信息(例如，對象分量)和與周期性圖樣相對應(yīng)的分量一起移除的風(fēng)險。
[0036]通帶特性確定裝置可以用下述方式來確定通帶特性:該通帶特性根據(jù)圖像的成像條件而不同。此處，可以通過使用成像菜單(例如，與成像條件和成像區(qū)域相關(guān)的信息)來獲得圖像的成像區(qū)域和圖像的成像條件，該成像菜單是在原始圖像的成像設(shè)備處指定的。備選地，可以使用由用 戶通過他/她的手動操作來輸入的成像區(qū)域信息，例如對象的器官和病變。此外，本發(fā)明的圖像處理設(shè)備還可以包括區(qū)域提取裝置，其從圖像中提取圖像的成像區(qū)域，且使用所提取的成像區(qū)域信息。例如，可以使用日本未審專利公開N0.2002-109548和日本未審專利公開N0.2003-6661中公開的方法，它們是由本申請的 申請人:所提出的。在該方法中，通過使用模板執(zhí)行模板匹配來自動檢測胸腔，該模板實(shí)質(zhì)上類似于作為參考的平均心胸腔的輪廓。
[0037]根據(jù)本發(fā)明的圖像處理設(shè)備、圖像處理方法和圖像處理程序，通過關(guān)于至少一個方向針對表示包括周期性圖樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬。此外，基于檢測到的峰位置和檢測到的峰寬來確定對圖像信號中與周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制的一維濾波器的通帶特性。因此，有可能基于峰位置和峰寬來恰當(dāng)確定一維濾波器的通帶特性，該一維濾波器可以在不過度抑制圖像中的對象分量的情況下恰當(dāng)?shù)匾种祁l率分量。因此，有可能以相對低的計算負(fù)荷來獲得高質(zhì)量圖像。
【專利附圖】

【附圖說明】
[0038]圖1是示出了放射線照相設(shè)備的示意圖；
[0039]圖2是示出了通過使用柵格成像來獲得的放射線照相圖像的圖；
[0040]圖3是示出了放射線照相圖像讀出設(shè)備的示例的立體圖；
[0041]圖4是示出了掃描方向和要讀出的圖像之間的關(guān)系的圖；
[0042]圖5是示出了根據(jù)本發(fā)明的實(shí)施例的圖像處理設(shè)備的配置的示意圖；
[0043]圖6是示出了根據(jù)本發(fā)明的實(shí)施例的由圖像處理設(shè)備進(jìn)行的處理的流程圖；
[0044]圖7是用于解釋本發(fā)明的實(shí)施例中的區(qū)域提取處理的示意圖；
[0045]圖8是用于解釋本發(fā)明的實(shí)施例中的歸一化處理的圖；
[0046]圖9A是用于解釋本發(fā)明的實(shí)施例中的通帶特性確定處理的圖；
[0047]圖9B是用于解釋本發(fā)明的實(shí)施例中的通帶特性確定處理的修改示例的圖；
[0048]圖1OA是用于解釋本發(fā)明的實(shí)施例中的抑制處理的圖；
[0049]圖1OB是用于解釋本發(fā)明的實(shí)施例中的抑制處理的修改示例的圖；
[0050]圖11是示出了本發(fā)明的實(shí)施例中的通帶特性的修改示例的圖(I號)；
[0051]圖12是示出了本發(fā)明的實(shí)施例中的通帶特性的修改示例的圖(2號)；以及
[0052]圖13是示出了本發(fā)明的實(shí)施例中的通帶特性的修改示例的圖(3號)。
【具體實(shí)施方式】
[0053]下文中，將參考附圖來描述本發(fā)明的實(shí)施例。在以下實(shí)施例中，將描述在放射線照相圖像讀出設(shè)備中使用根據(jù)本發(fā)明的周期性圖樣抑制處理設(shè)備的情況。備選地，可以在用于抑制放射線照相圖像中包括的周期性圖樣的圖像處理設(shè)備等中使用該周期性圖樣抑制處理設(shè)備，該放射線照相圖像是當(dāng)通過窗口屏幕、遮光物等來執(zhí)行放射線照相時在使用數(shù)字相機(jī)等的普通放射線照相中獲得的。此處，放射線照相圖像讀出設(shè)備通過使用激光束掃描放射線照相圖像，讀出在蓄積性熒光片上記錄的人體的放射線照相圖像作為數(shù)字圖像信號。
[0054]圖1是示出了放射線照相設(shè)備的示意圖。從放射線源I輸出放射線2，且放射線2穿過對象3，并到達(dá)靜態(tài)柵格(下文中，簡稱為“柵格”)4。在柵格4中，以例如約4行/mm的間距來交替布置鉛4a (其吸收放射線2)和鋁4b (其讓放射線2通過)。此外，以鉛4a的傾斜根據(jù)其位置而略微不同的方式來設(shè)置鉛4a，使得放射線2穿過鋁4b，并進(jìn)入蓄積性熒光片11。
[0055]因此，已穿過對象3的放射線2由鉛4a來吸收，且不到達(dá)蓄積性熒光片11。然而，放射線2穿過鋁4b，并到達(dá)蓄積性熒光片11。將具有4行/mm的條紋圖樣的柵格圖像與對象圖像一起記錄在蓄積性熒光片11上。同時，在對象3中被散射的散射放射線2a由鉛4a所吸收，或由柵格4的表面來反射，鉛4a是以根據(jù)其位置而傾斜的方式設(shè)置的。因此，散射的放射線2a不到達(dá)蓄積性熒光片11。因此，蓄積性熒光片11可以記錄具有少量散射放射線2a的銳利的放射線照相圖像，該少量散射放射線2a輻照蓄積性熒光片11。圖2是示出了對象圖像5和由柵格圖像引起的莫爾條紋圖樣6的放射線照相圖像的示例的圖。當(dāng)通過使用圖1所示的放射線照相設(shè)備來執(zhí)行放射線照相時，在蓄積性熒光片11上存儲和記錄放射線照相圖像。
[0056]圖3是示出了放射線照相圖像讀出設(shè)備的立體圖和功能框圖的組合的圖。在讀出單元10的預(yù)定位置上設(shè)置蓄積性熒光片11，并通過由驅(qū)動裝置(未示出)來驅(qū)動的片傳送裝置15 (例如環(huán)帶)來傳送。在箭頭方向Y上以例如掃描間距10行/mm來傳送(輔掃描)蓄積性熒光片11。同時，從激光束源16輸出光束17，且由旋轉(zhuǎn)多面鏡18向聚光透鏡19 (例如，f0透鏡)反射。旋轉(zhuǎn)多面鏡18由電機(jī)24驅(qū)動，并在箭頭方向上以高速旋轉(zhuǎn)。在光束17穿過聚光透鏡19之后，光束17由鏡子20向蓄積性熒光片11反射。由光束17在箭頭方向X上主掃描蓄積性熒光片11，該箭頭方向X與副掃描方向(箭頭方向Y)實(shí)質(zhì)上成直角。
[0057]當(dāng)使用光束17來照射蓄積性熒光片11時，從蓄積性熒光片11的被照射位置發(fā)射在量上與在該位置上存儲和記錄的放射線照相圖像信息相對應(yīng)的受激發(fā)射光21。受激發(fā)射光21從光導(dǎo)22的入射端面22a進(jìn)入光導(dǎo)22,并在光導(dǎo)22中重復(fù)全反射,并從光導(dǎo)22的輸出端面22b輸出。由光電倍增器23來接收輸出的受激發(fā)射光22,并通過光電轉(zhuǎn)換轉(zhuǎn)換為模擬圖像信號Sa。
.[0058]在由對數(shù)放大器26將模擬圖像信號Sa進(jìn)行對數(shù)放大之后，使用與空間頻率fs =10周期/mm相對應(yīng)的采樣間隔對模擬圖像信號Sa采樣，并由A/D轉(zhuǎn)換器28來數(shù)字化，且輸出數(shù)字圖像信號Sd(下文中，簡稱為“圖像信號Sd”)。圖像信號Sd表示通過對蓄積性熒光片11進(jìn)行二維掃描所獲得的放射線照相圖像信息。如圖4所示，通過在主掃描方向(水平方向)上使用光束17對蓄積性突光片11進(jìn)行掃描的同時在副掃描方向(垂直方向)上移動蓄積性熒光片11來獲得放射線照相圖像信息。除了與對象3相對應(yīng)的放射線照相圖像5相關(guān)的信息之外，以這種方式獲得的圖像信號Sd還包括由與柵格4相對應(yīng)的柵格圖像(未示出)所引起的莫爾條紋圖樣6相關(guān)的信息。
[0059]在存儲單元29中臨時存儲圖像信號Sd之后，向圖像信號處理單元30輸入圖像信號Sd。圖像信號處理單元30包括用于執(zhí)行本發(fā)明中的圖像處理方法的圖像處理設(shè)備40。
[0060]圖5是示出了根據(jù)本發(fā)明的實(shí)施例的圖像處理設(shè)備40的配置的示意圖。將參考圖5來描述根據(jù)本發(fā)明的實(shí)施例的圖像處理設(shè)備40。
[0061]根據(jù)本發(fā)明的實(shí)施例的圖像處理設(shè)備40包括周期性圖樣檢測裝置41、通帶特性確定裝置42、以及抑制裝置43。周期性圖樣檢測裝置41通過關(guān)于至少一個方向針對表示包括周期性圖樣在內(nèi)的圖像的圖像信號Sd執(zhí)行頻率分析，來檢測與周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬。通帶特性確定裝置基于檢測到的峰位置和檢測到的峰寬來確定一維濾波器的通帶特性，該一維濾波器對圖像信號中與周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。抑制裝置43通過基于所確定的通帶特性執(zhí)行一維濾波，來抑制圖像信號中的周期性圖樣。
[0062]當(dāng)安裝圖像處理程序時，在存儲單元29中存儲本發(fā)明的實(shí)施例中的圖像處理程序和該圖像處理程序所涉及的數(shù)據(jù)，并當(dāng)啟動圖像處理程序時，將其加載到存儲單元29中包括的存儲器。圖像處理程序定義了周期性圖樣檢測處理、通帶特性確定處理以及抑制處理，作為由圖像信號處理單元30的中央處理單元來執(zhí)行的處理，該圖像信號處理單元30構(gòu)成了圖像處理設(shè)備40。中央處理單元基于程序來執(zhí)行上述各個類型處理。因此，圖像信號處理單元30的中央處理單元擔(dān)當(dāng)周期性圖樣檢測裝置41、通帶特性確定裝置42、以及抑制裝置43。
[0063]圖6是示出了根據(jù)本發(fā)明的實(shí)施例中的圖像處理設(shè)備40中的圖像處理的流程的流程圖。圖7是用于解釋周期性圖樣檢測裝置41的線狀(line-shaped)采樣區(qū)域提取處理的圖。將參考圖6和圖7來詳細(xì)描述根據(jù)本發(fā)明的實(shí)施例的圖像處理設(shè)備40中的圖像處理。
[0064]首先，·如圖6所示，向周期性圖樣檢測裝置41輸入圖像數(shù)據(jù)Sd(SOl)。然后，周期性圖樣檢測裝置41在由圖像數(shù)據(jù)Sd表示的圖像中關(guān)于X方向和y方向中的每一個方向提取nxm個線狀區(qū)域，如圖7所示。由于用于關(guān)于X方向提取區(qū)域的方法和用于關(guān)于y方向提取區(qū)域的方法相同，因此此處將僅描述用于關(guān)于X方向提取區(qū)域的方法。
[0065]首先，周期性圖樣檢測裝置41提取η個(在本發(fā)明的實(shí)施例中，η = 9)矩形區(qū)域Al至An，它們在由圖像數(shù)據(jù)Sd表示的圖像Gl上沿著X方向延伸(S02)。然后，周期性圖樣檢測裝置41獲得η個區(qū)域圖像數(shù)據(jù)集合SRxi中每一個的m個(在本發(fā)明的實(shí)施例中，m =9)線狀區(qū)域圖像數(shù)據(jù)SRxij的集合(i = I至n，且j = I至m)。此處，區(qū)域Al在y方向上的中心位于遠(yuǎn)離圖像上邊緣1/30圖像高度的位置處。區(qū)域A9在Y方向上的中心位于遠(yuǎn)離圖像上邊緣1/30圖像高度的位置處。此外，區(qū)域A2到區(qū)域A7在區(qū)域Al和A9之間以相等距離來布置。在y方向上以三個像素間隔來布置m(在本發(fā)明的實(shí)施例中，m = 9)個線狀區(qū)域圖像數(shù)據(jù)的集合SRxij(i = I至n，且j = I至m)。在X方向上的像素數(shù)目是1024。
[0066]周期性圖樣檢測裝置41通過對每個SRxij (i = I至n，且j = I至m)執(zhí)行快速傅立葉變換來執(zhí)行頻率分析(S03)。具體地，針對每個區(qū)域圖像數(shù)據(jù)SRxij的像素串SRxii (k)(k = O至1023)來執(zhí)行快速傅立葉變換。換言之，當(dāng)頻譜是Pii (k) (k = O至1023)且通過快速傅立葉變換獲得的運(yùn)算結(jié)果的實(shí)部和虛部是Re (k)和Im(k)時，通過以下公式來計算頻譜 Pij (k):
[0067][式I]
[0068]
【權(quán)利要求】
1.一種圖像處理設(shè)備，包括: 周期性圖樣檢測裝置，通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及 通帶特性確定裝置，基于檢測到的峰位置和檢測到的峰寬來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。
2.根據(jù)權(quán)利要求1所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置以下述方式來確定所述通帶特性:包括所述峰位置且通帶寬度隨著所述峰寬更寬而更寬。
3.根據(jù)權(quán)利要求2所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置通過由所述峰位置和所述峰寬限定的正態(tài)分布函數(shù)來確定所述通帶特性。
4.根據(jù)權(quán)利要求1至3中任一項(xiàng)所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置以下述方式來確定所述通帶特性:通帶寬度隨著所述峰位置位于更低頻帶而更窄。
5.根據(jù)權(quán)利要求1至4中任一項(xiàng)所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置基于峰高來確定所述通帶特性。
6.根據(jù)權(quán)利要求1至5中任一項(xiàng)所述的圖像處理設(shè)備，所述圖像處理設(shè)備還包括: 抑制裝置，通過基于所確定的通帶特性執(zhí)行一維濾波，對所述圖像信號中的周期性圖樣進(jìn)行抑制。.
7.根據(jù)權(quán)利要求1至6中任一項(xiàng)所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置在所述通帶特性中設(shè)置通帶寬度的上限或下限。
8.根據(jù)權(quán)利要求1至7中任一項(xiàng)所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置以下述方式來確定所述通帶特性:所述通帶特性根據(jù)所述圖像的成像條件而不同。
9.根據(jù)權(quán)利要求1至8中任一項(xiàng)所述的圖像處理設(shè)備，其中，所述通帶特性確定裝置通過矩形函數(shù)來確定所述通帶特性，所述矩形函數(shù)具有比所述峰寬更寬或相等的寬度，且所述矩形函數(shù)的中心是所述峰位置。
10.一種圖像處理方法，包括以下步驟: 通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及 基于檢測到的峰位置和檢測到的峰寬來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。
11.一種圖像處理程序，使得計算機(jī)擔(dān)當(dāng): 周期性圖樣檢測裝置，通過關(guān)于至少一個方向?qū)Ρ硎景ㄖ芷谛詧D樣在內(nèi)的圖像的圖像信號執(zhí)行頻率分析，來檢測與所述周期性圖樣相對應(yīng)的頻譜的峰位置和峰寬；以及 通帶特性確定裝置，基于檢測到的峰位置和檢測到的峰寬來確定一維濾波器的通帶特性，所述一維濾波器對所述圖像信號中與所述周期性圖樣相對應(yīng)的空間頻率分量進(jìn)行抑制。
【文檔編號】G06T5/10GK103442641SQ201280014904
【公開日】2013年12月11日 申請日期:2012年3月23日 優(yōu)先權(quán)日:2011年3月24日
【發(fā)明者】今井睦朗 申請人:富士膠片株式會社