TW201346835A - 影像模糊程度之估測方法及影像品質之評估方法 - Google Patents
影像模糊程度之估測方法及影像品質之評估方法 Download PDFInfo
- Publication number
- TW201346835A TW201346835A TW101117289A TW101117289A TW201346835A TW 201346835 A TW201346835 A TW 201346835A TW 101117289 A TW101117289 A TW 101117289A TW 101117289 A TW101117289 A TW 101117289A TW 201346835 A TW201346835 A TW 201346835A
- Authority
- TW
- Taiwan
- Prior art keywords
- image
- input
- input image
- blur
- parameter
- Prior art date
Links
- 238000000034 method Methods 0.000 title claims abstract description 48
- 238000013441 quality evaluation Methods 0.000 title claims abstract description 16
- 238000012545 processing Methods 0.000 claims abstract description 37
- 230000008569 process Effects 0.000 claims abstract description 4
- 239000002131 composite material Substances 0.000 claims description 8
- 238000011156 evaluation Methods 0.000 claims description 7
- 238000013459 approach Methods 0.000 claims description 2
- 230000001419 dependent effect Effects 0.000 claims 1
- 230000015572 biosynthetic process Effects 0.000 abstract 3
- 238000003786 synthesis reaction Methods 0.000 abstract 3
- 238000003708 edge detection Methods 0.000 description 12
- 230000003287 optical effect Effects 0.000 description 10
- 238000003384 imaging method Methods 0.000 description 4
- 238000012805 post-processing Methods 0.000 description 4
- 238000006243 chemical reaction Methods 0.000 description 2
- 238000010586 diagram Methods 0.000 description 2
- 238000005315 distribution function Methods 0.000 description 2
- 230000000694 effects Effects 0.000 description 2
- 238000001914 filtration Methods 0.000 description 2
- 238000003702 image correction Methods 0.000 description 2
- 238000012986 modification Methods 0.000 description 2
- 230000004048 modification Effects 0.000 description 2
- 238000012634 optical imaging Methods 0.000 description 2
- 238000007781 pre-processing Methods 0.000 description 2
- 238000004364 calculation method Methods 0.000 description 1
- 230000015556 catabolic process Effects 0.000 description 1
- 230000008859 change Effects 0.000 description 1
- 239000003086 colorant Substances 0.000 description 1
- 238000012937 correction Methods 0.000 description 1
- 230000007547 defect Effects 0.000 description 1
- 230000003412 degenerative effect Effects 0.000 description 1
- 238000006731 degradation reaction Methods 0.000 description 1
- 238000013461 design Methods 0.000 description 1
- 238000001514 detection method Methods 0.000 description 1
- 229920013655 poly(bisphenol-A sulfone) Polymers 0.000 description 1
- 238000001303 quality assessment method Methods 0.000 description 1
- 230000009467 reduction Effects 0.000 description 1
- 230000007480 spreading Effects 0.000 description 1
- 230000001629 suppression Effects 0.000 description 1
- 230000002194 synthesizing effect Effects 0.000 description 1
- 230000000007 visual effect Effects 0.000 description 1
Classifications
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T7/00—Image analysis
- G06T7/0002—Inspection of images, e.g. flaw detection
-
- G—PHYSICS
- G06—COMPUTING; CALCULATING OR COUNTING
- G06T—IMAGE DATA PROCESSING OR GENERATION, IN GENERAL
- G06T2207/00—Indexing scheme for image analysis or image enhancement
- G06T2207/30—Subject of image; Context of image processing
- G06T2207/30168—Image quality inspection
Landscapes
- Engineering & Computer Science (AREA)
- Quality & Reliability (AREA)
- Computer Vision & Pattern Recognition (AREA)
- Physics & Mathematics (AREA)
- General Physics & Mathematics (AREA)
- Theoretical Computer Science (AREA)
- Image Processing (AREA)
- Studio Devices (AREA)
Abstract
本發明為一種影像模糊程度之估測方法及影像品質之評估方法,其先由一影像處理裝置傳送一輸入影像,依據該輸入影像之一輸入像素強度分布參數而產生一合成影像,且該合成影像包含一非線性影像感測函數;依據該輸入影像與該合成影像進行像素強度分布匹配處理,以產生一模糊程度參數;該影像處理裝置依據該模糊程度參數計算出該輸入影像之一模糊估測結果,且本發明之影像模糊程度之估測方法更可應用於針對複數輸入影像之複數有效區域進行模糊分布之估測,而獲得該些輸入影像之模糊分布之估測,進一步針對複數輸入影像進行影像品質之評估。
Description
本發明係有關於一種影像運算方法,特是指一種影像模糊程度之估測方法及影像品質之評估方法。
影像處理可分為影像前處理以及影像後處理。我們可以使用影像擷取裝置拍攝景物以獲得景物的數位影像資料,此時所獲得的影像資料一般稱為原始資料(Raw data)。原始資料會再被處理以便產生特定的影像效果。其中,使用影像擷取裝置拍攝取得原始資料影像的過程被稱為影像前處理。在此之後所進行的影像處理程序則稱為影像後處理。前處理程序例如為自動對焦(Auto focus)、自動曝光(Auto exposure)等在影像擷取時的控制。一般常見的影像後處理程序包括:對原始資料進行減少雜訊(Noise reduction)、白平衡(White balancing)、色彩內插法(Interpolation)、色彩校正(Color calibration)、γ補償(Gamma correction)、RGB轉換為YCbCr(Color space conversion)、邊緣加強(Edge enhancement)、飽和度加強(Saturation enhancement)以及錯色壓制(False color suppression)等程序,則可獲得良好的影像。
再者,隨著數位相機、數位攝影機等數位光學裝置的普及,人們對於影像品質的要求也隨著日益提升。然而,數位光學裝置在於光學設計殘存的成像缺陷、鏡片加工的誤差、影像感測器的非線性特性與雜訊等因素都將嚴重影響成像品質,特別是影像感測器的失焦容易導致模糊影像的產生。但,上述之影像後處理方式未能針對模糊影像的模糊程度進行估測,因此影像處理裝置在無法得知模糊影像之模糊程度的情況下,無法改善失焦之模糊影像的模糊情況。為了解決此一問題,影像處理技術領域中的研發人員因而發展出一種點擴散函數,以用於解決失焦問題。
一般來說,數位光學裝置可利用點擴散函數(Point Spreading Function,PSF)來表示一光路(Optical Path)或成像系統(Imaging System),其中每一物距就會有一組對應的PSF用以表示該光路在該物距下的光學特性,且影像感測器所擷取到的影像則是目標影像與PSF之卷積(Convolution)結果。但現今PSF係依循高斯函數所建立之高斯PSF,而高斯PSF係用以呈現線性感測元件下的光學成像系統,也就是高斯PSF為一理想狀態下的點擴散函數,然,數位光學成像系統實際上為建立於非線性感測元件之影像擷取系統,因此高斯PSF所得到之成像無法貼切表達實際的目標影像。
因此,即便在聚焦(Focus)的情形下,對焦的景物也無法完美地成像於影像感測器,而對於其它物距下的物體,更會因為失焦使得成像品質嚴重劣化。未聚焦影像通常在彼此靠近之像素間僅具有略微差異,因而讓目標影像與PSF之卷積結果造成複數像素之間的影像強度均勻分配,所以影像感測器所擷取到的影像為模糊且強度值均勻的影像。由於現今針對失焦的模糊影像仍以高斯PSF為修飾模糊的主要函數,且因高斯PSF並未能建構使用非線性感測元件光學系統之成像模型,而高斯PSF僅能針對線性光學系統進行模糊修飾,以及現今影像修正之技術並未能確實估測出模糊影像之模糊程度,所以現今影像修正之技術並未能有效地消除影像感測器因失焦所產生之模糊影像。
因此,本發明即針對上述問題而提出一種影像模糊程度之估測方法及影像品質之評估方法,其係利用非均勻影像相關聯之合成影像而與輸入影像作像素強度分布匹配處理,藉此,獲得模糊估測結果,並可進一步應用於影像品質評估。
再者,隨著數位相機、數位攝影機等數位光學裝置的普及,人們對於影像品質的要求也隨著日益提升。然而,數位光學裝置在於光學設計殘存的成像缺陷、鏡片加工的誤差、影像感測器的非線性特性與雜訊等因素都將嚴重影響成像品質,特別是影像感測器的失焦容易導致模糊影像的產生。但,上述之影像後處理方式未能針對模糊影像的模糊程度進行估測,因此影像處理裝置在無法得知模糊影像之模糊程度的情況下,無法改善失焦之模糊影像的模糊情況。為了解決此一問題,影像處理技術領域中的研發人員因而發展出一種點擴散函數,以用於解決失焦問題。
一般來說,數位光學裝置可利用點擴散函數(Point Spreading Function,PSF)來表示一光路(Optical Path)或成像系統(Imaging System),其中每一物距就會有一組對應的PSF用以表示該光路在該物距下的光學特性,且影像感測器所擷取到的影像則是目標影像與PSF之卷積(Convolution)結果。但現今PSF係依循高斯函數所建立之高斯PSF,而高斯PSF係用以呈現線性感測元件下的光學成像系統,也就是高斯PSF為一理想狀態下的點擴散函數,然,數位光學成像系統實際上為建立於非線性感測元件之影像擷取系統,因此高斯PSF所得到之成像無法貼切表達實際的目標影像。
因此,即便在聚焦(Focus)的情形下,對焦的景物也無法完美地成像於影像感測器,而對於其它物距下的物體,更會因為失焦使得成像品質嚴重劣化。未聚焦影像通常在彼此靠近之像素間僅具有略微差異,因而讓目標影像與PSF之卷積結果造成複數像素之間的影像強度均勻分配,所以影像感測器所擷取到的影像為模糊且強度值均勻的影像。由於現今針對失焦的模糊影像仍以高斯PSF為修飾模糊的主要函數,且因高斯PSF並未能建構使用非線性感測元件光學系統之成像模型,而高斯PSF僅能針對線性光學系統進行模糊修飾,以及現今影像修正之技術並未能確實估測出模糊影像之模糊程度,所以現今影像修正之技術並未能有效地消除影像感測器因失焦所產生之模糊影像。
因此,本發明即針對上述問題而提出一種影像模糊程度之估測方法及影像品質之評估方法,其係利用非均勻影像相關聯之合成影像而與輸入影像作像素強度分布匹配處理,藉此,獲得模糊估測結果,並可進一步應用於影像品質評估。
本發明之一主要目的,在於提供一種影像模糊程度之估測方法及影像品質之評估方法,其提供較佳之模糊估測結果,以供後續應用。
本發明之一次要目的,在於提供一種影像模糊程度之估測方法及影像品質之評估方法,其利用像素強度分布匹配處理,以提高模糊估測精確度。
本發明係提供一種影像模糊程度之估測方法,其先傳送一輸入影像至一影像處理裝置,接續由影像處理裝置依據一非線性影像感測函數與該輸入影像之一輸入像素強度分布參數而產生一合成模糊影像,再由影像處理裝置依據該輸入影像與該合成模糊影像進行像素強度分布匹配處理,以產生一模糊程度參數,最後由該影像處理裝置依據該輸入影像與該合成影像之該模糊程度參數計算出該輸入影像之一模糊估測結果,因而得知該輸入影像之模糊程度。
本發明係提供一種影像品質之評估方法,其係先傳送複數輸入影像至一影像處理裝置,接續由該影像處理裝置依據該些輸入影像之複數輸入像素強度分布參數而獲得複數邊緣影像,再由該影像處理裝置依據該些邊緣影像選取複數有效區域,然後由該影像處理裝置依據本發明所述之影像模糊程度之估測方法而針對該些有效區域進行模糊分布估測,以獲得該些有效區域之複數模糊估測結果,最後由該影像處理裝置依據該些模糊估測結果評估該些輸入影像之影像品質。
本發明之一次要目的,在於提供一種影像模糊程度之估測方法及影像品質之評估方法,其利用像素強度分布匹配處理,以提高模糊估測精確度。
本發明係提供一種影像模糊程度之估測方法,其先傳送一輸入影像至一影像處理裝置,接續由影像處理裝置依據一非線性影像感測函數與該輸入影像之一輸入像素強度分布參數而產生一合成模糊影像,再由影像處理裝置依據該輸入影像與該合成模糊影像進行像素強度分布匹配處理,以產生一模糊程度參數,最後由該影像處理裝置依據該輸入影像與該合成影像之該模糊程度參數計算出該輸入影像之一模糊估測結果,因而得知該輸入影像之模糊程度。
本發明係提供一種影像品質之評估方法,其係先傳送複數輸入影像至一影像處理裝置,接續由該影像處理裝置依據該些輸入影像之複數輸入像素強度分布參數而獲得複數邊緣影像,再由該影像處理裝置依據該些邊緣影像選取複數有效區域,然後由該影像處理裝置依據本發明所述之影像模糊程度之估測方法而針對該些有效區域進行模糊分布估測,以獲得該些有效區域之複數模糊估測結果,最後由該影像處理裝置依據該些模糊估測結果評估該些輸入影像之影像品質。
茲為使 貴審查委員對本發明之結構特徵及所達成之功效有更進一步之瞭解與認識,謹佐以較佳之實施例及配合詳細之說明,說明如後:
請參閱第一圖,其為本發明之一實施例之流程圖。如圖所示,本發明影像模糊程度之估測方法係應用於一影像處理裝置,以針對一失焦影像進行估測,以進一步了解失焦影像的模糊分布程度。本發明之估測方法係包含:
步驟S10:傳輸一輸入影像;
步驟S20:依據該輸入影像產生一合成影像;
步驟S30:選擇至少一邊緣區域;
步驟S40:過濾影像;
步驟S50:像素強度分布匹配;以及
步驟S60:獲得模糊程度估測結果。
於步驟S10中,利用一影像擷取裝置擷取一輸入影像並傳送至一影像處理裝置,本實施例之影像擷取裝置為一相機系統之一影像感測器,感測該輸入影像並傳送至該影像處理裝置,例如:中央處理單元(Central Processing Unit,CPU)、圖形處理單元(Graphic Processing Unit,GPU),其中該輸入影像為一失焦模糊影像,其影像函數為,而一對焦影像之影像函數為,其亦即表示該對焦影像之像素強度分布,因此,依據該失焦模糊影像與對焦影像之間的關聯性,藉由影像函數之卷積(convolution)可求得方程式(1),如下列方程式所示:
─────────────(1)
其中,為一影像雜訊之影像函數,一般利用影像感測器擷取影像時,皆會獲得該影像雜訊;為影像感測器之一點擴散函數(Point Spread Function,PSF),基於空間不變性退化的特性,點擴散函數與影像雜訊之影像函數實際上的分布特性相當於對焦影像之影像函數的分佈特性,一般影像感測器之點擴散函數為非線性函數。
於步驟S20中,影像處理裝置依據該輸入影像之像素強度分布而對應產生一合成模糊影像,其中本發明之影像處理裝置為藉由該輸入影像之一模糊像素分布參數的二相鄰像素分布參數而合成該合成模糊影像,例如:一失焦模糊之黑白輸入影像由不同灰階亮度之像素所構成,而失焦模糊部分為一影像之輪廓模糊,亦即較亮之像素與較暗之像素混合排列,未亮暗分明,本發明之影像處理裝置係針對失焦模糊部分二側的較亮像素分布與較暗像素分布進行疊合而合成對應於該黑白輸入影像,再者,本發明之合成模糊影像的影像函數可被表示如上方程式(1),但其中點擴散函數則如下方程式(2)所示:
──────────(2)
其中為一光學函數,該光學函數係關聯於光圈大小與等效焦長,藉由影像強度值的非線性調整與方程式(2)提供卷積,以讓合成模糊影像趨近實際上的模糊影像,例如:該輸入影像。如第二A圖至第二C圖所示,其分別為高斯合成影像、輸入影像與本發明之合成模糊影像的直方圖,模糊影像之像素強度分布為大多數像素之影像強度坐落於影像強度值25-225之間,其中第二A圖為習知合成影像之像素強度分布的直方圖,其為一般影像還原技術所採用之高斯點擴散函數,以形成一高斯合成影像,因此第二A圖之高斯合成影像的像素強度分布會較為均勻分布,例如:以影像強度值125為中心點,而呈左右對稱;第二B圖為實際擷取之輸入影像的像素強度分布的直方圖,由第二B圖可知,輸入影像的像數強度分布並未對稱,由於高斯點擴散函數為理想化參數,使高斯合成影像的像數強度分布呈現均勻分布,因此影像處理裝置依據高斯點擴散函數所產生之高斯合成影像無法適切地趨近於輸入影像之像素強度分布曲線;第二C圖為影像處理裝置依據本實施例所述之合成模糊影像,由於本實施例之合成影像受到非線性影像感測函數與點擴散函數的影響亦呈現非對稱分布,所以第二C圖之像素強度分布曲線趨近於輸入影像之像素強度分布曲線,因此本發明之合成模糊影像較適合搭配輸入影像進行像素強度分布匹配處理。
為了凸顯像素強度分布之差異,請一併參閱第三A圖至第三C圖,其分別為高斯合成影像、輸入影像與本發明之合成模糊影像之邊緣的影像對亮度之曲線圖。如圖所示,第三A圖之高斯合成影像會在影像之邊緣部分呈現線性結果,亦即趨近一直線,但第三B圖之輸入影像相較於第三A圖可知,實際上的輸入影像於影像之邊緣部分並非呈現線性結果,因而呈現一曲線,而本發明之之合成模糊影像卻是如第三C圖所示,在影像之邊緣呈現一曲線,因此本發明之合成模糊影像可較趨近於輸入影像之邊緣部分的曲線。
此外,本發明為了讓合成模糊影像更逼近輸入影像之像素強度分布,更可將一比例因子加入上述之方程式(2)中,以讓本發明之合成模糊影像更加以逼近輸入影像,該比例因子之表示式如下列方程式(3)所示:
────────(3)
,且 ────────(4)
因此方程式(2)可改寫為如下方程式(5)所示,以清楚定義出點擴散函數之邊緣。
───────────(5)
於步驟S30中,依據輸入影像與合成模糊影像選取出有效之邊緣,一般影像處理裝置以一邊緣偵測法偵測有效之邊緣,例如:Canny邊緣偵測法、高斯邊緣偵測法等之偵測法;
於步驟S40中,基於像素強度分布在高頻、低頻訊號中表現出的變異性較不明顯,而中頻訊號較能呈現出細緻的變異性,因此針對中頻訊號加入雜訊變異度可更為突顯出像素強度分布之變異,因此方程式(1)將改寫成方程式(6),其表示式如下所示:
───────────(6)
其中即為將方程式(1)之雜訊函數改寫為零均值(zero-mean)的高斯分布區線函數,而為影像之標準差,為了接續影像匹配,而讓合成模糊影像所呈現之變異性可與輸入影像相關,因此將標準差定義為,如此影像在亮部與暗部所呈現之雜訊強度會低於在亮、暗之間的雜訊強度。
於步驟S50中,依據步驟S40過濾完成之輸入影像與合成模糊影像進行像素強度分布匹配處理,因而獲得一影像模糊程度參數,其為平均值、與變異數、,其關係如下方程式(7)所示:
────────────(7)
於步驟S60中,依據方程式(7)獲得模糊程度之估測結果,其表示式如下列方程式(8)所示:
────────────(8)
藉由方程式(8)可求得最低模糊程度與最大模糊程度,其分別為下列方程式(9)與方程式(10)。
──────(9)
───(10)
此外,方程式(9)與方程式(10)更可在去除變數、之後,求得方程式(11)與方程式(12),其分別為失焦模糊影像之模糊程度的初始值與最大估計值。
─────────────────(11)
────────────────(12)
由以上所述可知,本發明之模糊分布之估測方法係針對失焦模糊影像可有效地提供客觀且逼近於實施需求的模糊估測,而非利用輸入影像之理想條件進行模糊分布之估測,因此本發明之估測方法方便使用者可確實得知影像之實際模糊程度,且針對後續應用可獲得更加精確的模糊估測結果,例如:影像品質評估、模糊影像還原等應用可提高效能。
請參閱第四圖,其為本發明之另一實施例之流程圖。如圖所示,本發明之影像品質之評估方法為採用前一實施例之模糊程度估測方法,而應用影像品質評估。該評估方法之步驟包含:
步驟S100:傳輸複數輸入影像;
步驟S110:邊緣偵測;
步驟S120:選取有效區域;
步驟S130:進行模糊估測,以獲得模糊估測結果;以及
步驟S140:依據輸入影像之模糊估測結果而獲得影像品質評估結果。
於步驟S100中,利用影像擷取裝置(如數位相機、網路監視器等)或儲存裝置將複數輸入影像傳送至影像處理裝置;於步驟S110中,其為基於仿人類視覺系統而針對該些輸入影像濾除非有效邊緣,亦即篩選出該些輸入影像之有效邊緣,例如:利用Canny邊緣偵測法偵測該些輸入影像之邊緣而過濾該些輸入影像,使該些輸入影像剩下的邊緣相對較少,或者採用Sobel邊緣偵測法針對仿人類視覺系統之RGB三原色中的G色系進行邊緣偵測,如第五A圖與第五B圖所示,本實施例之邊緣偵測係針對輸入影像之邊緣進行偵測,而獲得第五A圖之邊緣影像,再經由設置邊緣偵測之偵測條件而濾除難以運算之非有效邊緣,而獲得如第五B圖所示之邊緣偵測結果;於步驟S120,針對步驟S110之邊緣偵測結果選取出有效區域(Region Of Interest,ROI),本實施例之選取方式係以標準差作為一選取門檻值,以保持被選取之輸入影像的標準差大於或等於該選取門檻值,而選取有效區域,以接續進行模糊估測。
於步驟S130所示,採用前一實施例之模糊估測方法而將步驟S120之選取結果進行模糊估測,而獲得該些輸入影像之模糊估測結果,由於模糊估測方式已於前一實施例說明,因此本實施例不再詳述;於步驟S140中,依據步驟S130之估測結果而進行影像品質評估,其中影像品質評估的判斷依據可為前一實施例之最低模糊程度與最大模糊程度,或者為模糊程度之初始估計值與最大估計值,因而獲得評估結果。
由以上所述可知,使用者無需在每一次獲得擷取影像後,隨即人工判斷是否影像品質過差,可藉由本發明之評估方法獲得影像品質之評估,在做決定是否重新擷取影像,或者針對影像品質進行排序等等。
綜上所述,本發明之影像模糊程度之估測方法,其為針對失焦模糊影像提供較為逼近之合成模糊影像,以進行像素強度分布匹配處理,而獲得較逼近輸入影像之模糊估測結果,以用於後續之應用。此外,本發明之影像品質之評估方法即為上述影像模糊程度之估測方法的應用,以用於評估影像品質。藉此提供較為方便之影像處理技術,讓使用者減少針對影像品質或模糊程度的人工判斷,並方便使用者快速得知失焦模糊影像的模糊程度以及影像品質,更可提供配套軟體可獲得較精確之模糊程度的估測結果。
雖然本發明已以較佳實施例揭露如上,然其並非用以限定本發明,任何熟習此技藝者,在不脫離本發明之精神和範圍內,當可作些許之更動與潤飾,因此本發明之保護範圍當視後附之申請專利範圍所界定者為準。
請參閱第一圖,其為本發明之一實施例之流程圖。如圖所示,本發明影像模糊程度之估測方法係應用於一影像處理裝置,以針對一失焦影像進行估測,以進一步了解失焦影像的模糊分布程度。本發明之估測方法係包含:
步驟S10:傳輸一輸入影像;
步驟S20:依據該輸入影像產生一合成影像;
步驟S30:選擇至少一邊緣區域;
步驟S40:過濾影像;
步驟S50:像素強度分布匹配;以及
步驟S60:獲得模糊程度估測結果。
於步驟S10中,利用一影像擷取裝置擷取一輸入影像並傳送至一影像處理裝置,本實施例之影像擷取裝置為一相機系統之一影像感測器,感測該輸入影像並傳送至該影像處理裝置,例如:中央處理單元(Central Processing Unit,CPU)、圖形處理單元(Graphic Processing Unit,GPU),其中該輸入影像為一失焦模糊影像,其影像函數為,而一對焦影像之影像函數為,其亦即表示該對焦影像之像素強度分布,因此,依據該失焦模糊影像與對焦影像之間的關聯性,藉由影像函數之卷積(convolution)可求得方程式(1),如下列方程式所示:
其中,為一影像雜訊之影像函數,一般利用影像感測器擷取影像時,皆會獲得該影像雜訊;為影像感測器之一點擴散函數(Point Spread Function,PSF),基於空間不變性退化的特性,點擴散函數與影像雜訊之影像函數實際上的分布特性相當於對焦影像之影像函數的分佈特性,一般影像感測器之點擴散函數為非線性函數。
於步驟S20中,影像處理裝置依據該輸入影像之像素強度分布而對應產生一合成模糊影像,其中本發明之影像處理裝置為藉由該輸入影像之一模糊像素分布參數的二相鄰像素分布參數而合成該合成模糊影像,例如:一失焦模糊之黑白輸入影像由不同灰階亮度之像素所構成,而失焦模糊部分為一影像之輪廓模糊,亦即較亮之像素與較暗之像素混合排列,未亮暗分明,本發明之影像處理裝置係針對失焦模糊部分二側的較亮像素分布與較暗像素分布進行疊合而合成對應於該黑白輸入影像,再者,本發明之合成模糊影像的影像函數可被表示如上方程式(1),但其中點擴散函數則如下方程式(2)所示:
其中為一光學函數,該光學函數係關聯於光圈大小與等效焦長,藉由影像強度值的非線性調整與方程式(2)提供卷積,以讓合成模糊影像趨近實際上的模糊影像,例如:該輸入影像。如第二A圖至第二C圖所示,其分別為高斯合成影像、輸入影像與本發明之合成模糊影像的直方圖,模糊影像之像素強度分布為大多數像素之影像強度坐落於影像強度值25-225之間,其中第二A圖為習知合成影像之像素強度分布的直方圖,其為一般影像還原技術所採用之高斯點擴散函數,以形成一高斯合成影像,因此第二A圖之高斯合成影像的像素強度分布會較為均勻分布,例如:以影像強度值125為中心點,而呈左右對稱;第二B圖為實際擷取之輸入影像的像素強度分布的直方圖,由第二B圖可知,輸入影像的像數強度分布並未對稱,由於高斯點擴散函數為理想化參數,使高斯合成影像的像數強度分布呈現均勻分布,因此影像處理裝置依據高斯點擴散函數所產生之高斯合成影像無法適切地趨近於輸入影像之像素強度分布曲線;第二C圖為影像處理裝置依據本實施例所述之合成模糊影像,由於本實施例之合成影像受到非線性影像感測函數與點擴散函數的影響亦呈現非對稱分布,所以第二C圖之像素強度分布曲線趨近於輸入影像之像素強度分布曲線,因此本發明之合成模糊影像較適合搭配輸入影像進行像素強度分布匹配處理。
為了凸顯像素強度分布之差異,請一併參閱第三A圖至第三C圖,其分別為高斯合成影像、輸入影像與本發明之合成模糊影像之邊緣的影像對亮度之曲線圖。如圖所示,第三A圖之高斯合成影像會在影像之邊緣部分呈現線性結果,亦即趨近一直線,但第三B圖之輸入影像相較於第三A圖可知,實際上的輸入影像於影像之邊緣部分並非呈現線性結果,因而呈現一曲線,而本發明之之合成模糊影像卻是如第三C圖所示,在影像之邊緣呈現一曲線,因此本發明之合成模糊影像可較趨近於輸入影像之邊緣部分的曲線。
此外,本發明為了讓合成模糊影像更逼近輸入影像之像素強度分布,更可將一比例因子加入上述之方程式(2)中,以讓本發明之合成模糊影像更加以逼近輸入影像,該比例因子之表示式如下列方程式(3)所示:
,且
因此方程式(2)可改寫為如下方程式(5)所示,以清楚定義出點擴散函數之邊緣。
於步驟S30中,依據輸入影像與合成模糊影像選取出有效之邊緣,一般影像處理裝置以一邊緣偵測法偵測有效之邊緣,例如:Canny邊緣偵測法、高斯邊緣偵測法等之偵測法;
於步驟S40中,基於像素強度分布在高頻、低頻訊號中表現出的變異性較不明顯,而中頻訊號較能呈現出細緻的變異性,因此針對中頻訊號加入雜訊變異度可更為突顯出像素強度分布之變異,因此方程式(1)將改寫成方程式(6),其表示式如下所示:
其中即為將方程式(1)之雜訊函數改寫為零均值(zero-mean)的高斯分布區線函數,而為影像之標準差,為了接續影像匹配,而讓合成模糊影像所呈現之變異性可與輸入影像相關,因此將標準差定義為,如此影像在亮部與暗部所呈現之雜訊強度會低於在亮、暗之間的雜訊強度。
於步驟S50中,依據步驟S40過濾完成之輸入影像與合成模糊影像進行像素強度分布匹配處理,因而獲得一影像模糊程度參數,其為平均值、與變異數、,其關係如下方程式(7)所示:
於步驟S60中,依據方程式(7)獲得模糊程度之估測結果,其表示式如下列方程式(8)所示:
藉由方程式(8)可求得最低模糊程度與最大模糊程度,其分別為下列方程式(9)與方程式(10)。
此外,方程式(9)與方程式(10)更可在去除變數、之後,求得方程式(11)與方程式(12),其分別為失焦模糊影像之模糊程度的初始值與最大估計值。
由以上所述可知,本發明之模糊分布之估測方法係針對失焦模糊影像可有效地提供客觀且逼近於實施需求的模糊估測,而非利用輸入影像之理想條件進行模糊分布之估測,因此本發明之估測方法方便使用者可確實得知影像之實際模糊程度,且針對後續應用可獲得更加精確的模糊估測結果,例如:影像品質評估、模糊影像還原等應用可提高效能。
請參閱第四圖,其為本發明之另一實施例之流程圖。如圖所示,本發明之影像品質之評估方法為採用前一實施例之模糊程度估測方法,而應用影像品質評估。該評估方法之步驟包含:
步驟S100:傳輸複數輸入影像;
步驟S110:邊緣偵測;
步驟S120:選取有效區域;
步驟S130:進行模糊估測,以獲得模糊估測結果;以及
步驟S140:依據輸入影像之模糊估測結果而獲得影像品質評估結果。
於步驟S100中,利用影像擷取裝置(如數位相機、網路監視器等)或儲存裝置將複數輸入影像傳送至影像處理裝置;於步驟S110中,其為基於仿人類視覺系統而針對該些輸入影像濾除非有效邊緣,亦即篩選出該些輸入影像之有效邊緣,例如:利用Canny邊緣偵測法偵測該些輸入影像之邊緣而過濾該些輸入影像,使該些輸入影像剩下的邊緣相對較少,或者採用Sobel邊緣偵測法針對仿人類視覺系統之RGB三原色中的G色系進行邊緣偵測,如第五A圖與第五B圖所示,本實施例之邊緣偵測係針對輸入影像之邊緣進行偵測,而獲得第五A圖之邊緣影像,再經由設置邊緣偵測之偵測條件而濾除難以運算之非有效邊緣,而獲得如第五B圖所示之邊緣偵測結果;於步驟S120,針對步驟S110之邊緣偵測結果選取出有效區域(Region Of Interest,ROI),本實施例之選取方式係以標準差作為一選取門檻值,以保持被選取之輸入影像的標準差大於或等於該選取門檻值,而選取有效區域,以接續進行模糊估測。
於步驟S130所示,採用前一實施例之模糊估測方法而將步驟S120之選取結果進行模糊估測,而獲得該些輸入影像之模糊估測結果,由於模糊估測方式已於前一實施例說明,因此本實施例不再詳述;於步驟S140中,依據步驟S130之估測結果而進行影像品質評估,其中影像品質評估的判斷依據可為前一實施例之最低模糊程度與最大模糊程度,或者為模糊程度之初始估計值與最大估計值,因而獲得評估結果。
由以上所述可知,使用者無需在每一次獲得擷取影像後,隨即人工判斷是否影像品質過差,可藉由本發明之評估方法獲得影像品質之評估,在做決定是否重新擷取影像,或者針對影像品質進行排序等等。
綜上所述,本發明之影像模糊程度之估測方法,其為針對失焦模糊影像提供較為逼近之合成模糊影像,以進行像素強度分布匹配處理,而獲得較逼近輸入影像之模糊估測結果,以用於後續之應用。此外,本發明之影像品質之評估方法即為上述影像模糊程度之估測方法的應用,以用於評估影像品質。藉此提供較為方便之影像處理技術,讓使用者減少針對影像品質或模糊程度的人工判斷,並方便使用者快速得知失焦模糊影像的模糊程度以及影像品質,更可提供配套軟體可獲得較精確之模糊程度的估測結果。
雖然本發明已以較佳實施例揭露如上,然其並非用以限定本發明,任何熟習此技藝者,在不脫離本發明之精神和範圍內,當可作些許之更動與潤飾,因此本發明之保護範圍當視後附之申請專利範圍所界定者為準。
無。
第一圖為本發明之一實施例之流程圖;
第二A圖為習知高斯合成影像之像素強度分布的曲線圖;
第二B圖為輸入影像之像素強度分布的曲線圖;
第二C圖為本發明之合成模糊影像之像素強度分布的曲線圖;
第三A圖為習知高斯合成影像之像素強度分布的曲線圖;
第三B圖為輸入影像之像素強度分布的曲線圖;
第三C圖為本發明之合成模糊影像之像素強度分布的曲線圖;
第四圖為本發明之另一實施例之流程圖;
第五A圖為本發明之另一實施例之邊緣偵測的示意圖;以及
第五B圖為本發明之另一實施例之邊緣濾除的示意圖。
第二A圖為習知高斯合成影像之像素強度分布的曲線圖;
第二B圖為輸入影像之像素強度分布的曲線圖;
第二C圖為本發明之合成模糊影像之像素強度分布的曲線圖;
第三A圖為習知高斯合成影像之像素強度分布的曲線圖;
第三B圖為輸入影像之像素強度分布的曲線圖;
第三C圖為本發明之合成模糊影像之像素強度分布的曲線圖;
第四圖為本發明之另一實施例之流程圖;
第五A圖為本發明之另一實施例之邊緣偵測的示意圖;以及
第五B圖為本發明之另一實施例之邊緣濾除的示意圖。
Claims (9)
- 一種影像模糊程度之估測方法,其包含:
傳送一輸入影像至一影像處理裝置;
依據該輸入影像之一模糊像素分布的二相鄰像素強度分布參數而產生一合成模糊影像,該合成影像包含一非線性影像感測函數;以及
比對該輸入影像與該合成模糊影像,以進行像素強度分布匹配處理而產生該輸入影像之一模糊程度參數模糊程度參數。 - 如申請專利範圍第1項所述之估測方法,模糊程度參數更包含:
依據該模糊程度參數計算出該輸入影像之一初始模糊估測值;以及
依據該模糊程度參數計算出該輸入影像之一最大模糊估測值。 - 如申請專利範圍第1項所述之估測方法,其中該合成影像包含一影像卷積參數與一雜訊參數,該影像卷積參數為該輸入影像與一點擴散函數之卷積結果。
- 如申請專利範圍第3項所述之估測方法,其中該點擴散函數係與一像素強度相依的比例因子相關聯。
- 如申請專利範圍第1項所述之估測方法,其中於依據一非線性影像感測函數與該輸入影像之一輸入像素強度分布參數而產生一合成影像之步驟前,更包含:
依據一比例因子過濾該合成影像,以讓該合成影像之一合成像素強度分布參數趨近於該輸入影像之該輸入像素強度分布參數。 - 一種影像品質之評估方法,其包含:
傳送一第一輸入影像與一第二輸入影像至一影像處理裝置;
依據該第一輸入影像與該第二輸入影像之複數輸入像素強度分布參數而獲得一第一邊緣影像與一第二邊緣影像;
依據該第一邊緣影像與該第二邊緣影像選取複數有效區域;
該影像處理裝置依據該些有效區域產生對應之一第一合成影像與一第二合成影像;
該影像處理裝置依據該第一合成影像與該第二合成影像分別比對該第一輸入影像與該第二輸入影像之該些有效區域,以進行模糊分布估測,而獲得該第一輸入影像與該第二輸入影像之一第一模糊程度參數與一第二模糊成度參數;以及
比對該第一模糊程度參數與該第二模糊成度參數,以評估該第一輸入影像與該第二輸入影像之影像品質。 - 如申請專利範圍第6項所述之評估方法,其中該些有效區域之一亮部標準差與一暗部標準差皆小於一分布門檻值。
- 如申請專利範圍第6項所述之評估方法,其中該些有效區域之一影像標準差差異大於一差異門檻值。
- 如申請專利範圍第6項所述之評估方法,其中該些有效區域之一影像標準差大於等於一影像標準門檻值。
Priority Applications (2)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW101117289A TWI462054B (zh) | 2012-05-15 | 2012-05-15 | Estimation Method of Image Vagueness and Evaluation Method of Image Quality |
US13/684,324 US8917938B2 (en) | 2012-05-15 | 2012-11-23 | Method for estimating blur degree of image and method for evaluating image quality |
Applications Claiming Priority (1)
Application Number | Priority Date | Filing Date | Title |
---|---|---|---|
TW101117289A TWI462054B (zh) | 2012-05-15 | 2012-05-15 | Estimation Method of Image Vagueness and Evaluation Method of Image Quality |
Publications (2)
Publication Number | Publication Date |
---|---|
TW201346835A true TW201346835A (zh) | 2013-11-16 |
TWI462054B TWI462054B (zh) | 2014-11-21 |
Family
ID=49581350
Family Applications (1)
Application Number | Title | Priority Date | Filing Date |
---|---|---|---|
TW101117289A TWI462054B (zh) | 2012-05-15 | 2012-05-15 | Estimation Method of Image Vagueness and Evaluation Method of Image Quality |
Country Status (2)
Country | Link |
---|---|
US (1) | US8917938B2 (zh) |
TW (1) | TWI462054B (zh) |
Families Citing this family (12)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
TWI462054B (zh) * | 2012-05-15 | 2014-11-21 | Nat Univ Chung Cheng | Estimation Method of Image Vagueness and Evaluation Method of Image Quality |
US9235711B1 (en) * | 2014-06-24 | 2016-01-12 | Voxience S.A.R.L. | Systems, methods and devices for providing visual privacy to messages |
WO2016203282A1 (en) * | 2015-06-18 | 2016-12-22 | The Nielsen Company (Us), Llc | Methods and apparatus to capture photographs using mobile devices |
US10181208B2 (en) | 2016-02-10 | 2019-01-15 | Microsoft Technology Licensing, Llc | Custom heatmaps |
CN107784645A (zh) * | 2016-08-26 | 2018-03-09 | 广州康昕瑞基因健康科技有限公司 | 图像清晰度评价方法及系统、自动聚焦方法 |
EP3668077B1 (en) * | 2017-08-09 | 2023-08-02 | FUJIFILM Corporation | Image processing system, server device, image processing method, and image processing program |
DE102018103714A1 (de) * | 2018-02-20 | 2019-08-22 | Volume Graphics Gmbh | Verfahren zur Bestimmung von Fehlern von aus digitalen Objektdarstellungen abgeleiteten Parametern |
US11195257B2 (en) * | 2018-08-24 | 2021-12-07 | Canon Kabushiki Kaisha | Image processing method, image processing apparatus, imaging apparatus, lens apparatus, storage medium, and image processing system |
CN111798421A (zh) * | 2020-06-29 | 2020-10-20 | 浙江同善人工智能技术有限公司 | 一种图像质量判断方法、设备及存储介质 |
CN111861963B (zh) * | 2020-07-28 | 2023-08-01 | 平安科技(深圳)有限公司 | 图像合成方法及相关设备 |
CN111968052B (zh) * | 2020-08-11 | 2024-04-30 | 北京小米松果电子有限公司 | 图像处理方法、图像处理装置及存储介质 |
CN112396050B (zh) * | 2020-12-02 | 2023-09-15 | 度小满科技(北京)有限公司 | 图像的处理方法、设备以及存储介质 |
Family Cites Families (13)
Publication number | Priority date | Publication date | Assignee | Title |
---|---|---|---|---|
US7015954B1 (en) * | 1999-08-09 | 2006-03-21 | Fuji Xerox Co., Ltd. | Automatic video system using multiple cameras |
GB0017039D0 (en) * | 2000-07-11 | 2000-08-30 | King S College London | Image enhancement |
CN101297545B (zh) * | 2005-10-28 | 2012-05-02 | 株式会社尼康 | 摄影装置、图像处理装置 |
KR101035824B1 (ko) * | 2006-05-09 | 2011-05-20 | 샤프 가부시키가이샤 | 촬상 장치 및 촬상 방법 |
EP2053557B1 (en) * | 2006-08-17 | 2013-11-06 | Fujitsu Limited | Image processing device, image processing program, and image processing method |
WO2008144180A1 (en) * | 2007-05-18 | 2008-11-27 | Samsung Electronics Co., Ltd. | Image color balance adjustment for display panels with 2d subpixel layouts |
JP2009110137A (ja) * | 2007-10-29 | 2009-05-21 | Ricoh Co Ltd | 画像処理装置、画像処理方法及び画像処理プログラム |
KR101574733B1 (ko) * | 2008-11-19 | 2015-12-04 | 삼성전자 주식회사 | 고화질 컬러 영상을 획득하기 위한 영상 처리 장치 및 방법 |
US8395642B2 (en) * | 2009-03-17 | 2013-03-12 | Mitsubishi Electric Research Laboratories, Inc. | Method for virtual image synthesis |
US8823873B2 (en) * | 2009-10-10 | 2014-09-02 | Thomson Licensing | Method and device for calculating blur in video images |
TWI405147B (zh) * | 2010-03-16 | 2013-08-11 | Novatek Microelectronics Corp | 單張影像的階層式去動態模糊方法 |
TWI421618B (zh) * | 2010-04-09 | 2014-01-01 | Ind Tech Res Inst | 景深擴展之投影系統及影像處理方法 |
TWI462054B (zh) * | 2012-05-15 | 2014-11-21 | Nat Univ Chung Cheng | Estimation Method of Image Vagueness and Evaluation Method of Image Quality |
-
2012
- 2012-05-15 TW TW101117289A patent/TWI462054B/zh not_active IP Right Cessation
- 2012-11-23 US US13/684,324 patent/US8917938B2/en active Active
Also Published As
Publication number | Publication date |
---|---|
US8917938B2 (en) | 2014-12-23 |
US20130308866A1 (en) | 2013-11-21 |
TWI462054B (zh) | 2014-11-21 |
Similar Documents
Publication | Publication Date | Title |
---|---|---|
TWI462054B (zh) | Estimation Method of Image Vagueness and Evaluation Method of Image Quality | |
US10997696B2 (en) | Image processing method, apparatus and device | |
KR101662846B1 (ko) | 아웃 포커싱 촬영에서 빛망울 효과를 생성하기 위한 장치 및 방법 | |
JP5213670B2 (ja) | 撮像装置及びぶれ補正方法 | |
US20110128422A1 (en) | Image capturing apparatus and image processing method | |
CN105046655B (zh) | 一种视频图像自动锐化方法和装置 | |
JP6833415B2 (ja) | 画像処理装置、画像処理方法、及びプログラム | |
JP4454657B2 (ja) | ぶれ補正装置及び方法、並びに撮像装置 | |
US20120307009A1 (en) | Method and apparatus for generating image with shallow depth of field | |
JP2022036179A (ja) | 画像処理装置およびその制御方法及びプログラム | |
KR20140118031A (ko) | 영상 처리 장치 및 방법 | |
US11836900B2 (en) | Image processing apparatus | |
CN114584700B (zh) | 对焦标记方法、标记装置以及电子设备 | |
WO2012008116A1 (en) | Image processing apparatus, image processing method, and program | |
JP6624785B2 (ja) | 画像処理方法、画像処理装置、撮像装置、プログラム、および、記憶媒体 | |
JP2012085083A (ja) | 画像処理装置、撮像装置および画像処理プログラム | |
JP6075835B2 (ja) | 距離情報取得装置、撮像装置、距離情報取得方法、及び、プログラム | |
JP5115297B2 (ja) | 画像処理装置、撮像装置、画像処理方法およびプログラム | |
JP2018160024A (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP6348883B2 (ja) | 画像撮像装置、画像撮像方法及びコンピュータプログラム | |
van Zwanenberg et al. | Camera system performance derived from natural scenes | |
JP7297566B2 (ja) | 画像処理装置、撮像装置、画像処理方法およびプログラム | |
JP7408353B2 (ja) | 画像処理装置、画像処理方法及びプログラム | |
JP2017219737A (ja) | 撮像装置及びその制御方法、プログラム、記憶媒体 | |
JP2016170637A (ja) | 画像処理装置および画像処理方法 |
Legal Events
Date | Code | Title | Description |
---|---|---|---|
MM4A | Annulment or lapse of patent due to non-payment of fees |