點擴(kuò)散函數(shù)的顏色影響性研究

宮本興

摘 要： 點擴(kuò)散函數(shù)是圖像處理的核心元素，是研究圖像復(fù)原，增強的基礎(chǔ)參數(shù)，也是構(gòu)建退化函數(shù)的主要部分。在此應(yīng)用Matlab程序，對擁有不同影響因素的實驗圖像進(jìn)行點擴(kuò)散函數(shù)的計算、歸納、比對和總結(jié)。探究圖像顏色對點擴(kuò)散函數(shù)的影響因素，圖像自身參數(shù)影響因素，從而為建立應(yīng)用更為廣泛的退化模型奠定基礎(chǔ)，也為圖像復(fù)原工作的準(zhǔn)確性提供更準(zhǔn)確的依據(jù)。

關(guān)鍵詞： 點擴(kuò)散函數(shù)； Matlab； 圖像復(fù)原； 退化函數(shù)

中圖分類號： TN919?34 文獻(xiàn)標(biāo)識碼： A 文章編號： 1004?373X（2015）10?0052?03

0 引 言

隨著人類對陸地資源的不斷開發(fā)，資源匱乏是人類在21世紀(jì)共同面對的嚴(yán)峻問題。豐富的海洋資源對面臨能源危機的人類有不可替代的重要作用，開發(fā)及利用海洋已然刻不容緩。了解掌握足夠有效的海洋信息是開發(fā)的前提，水下成像技術(shù)便是獲取海洋信息的有效途徑。但水底環(huán)境復(fù)雜多變，主要困難源自散射介質(zhì)對光的吸收和散射。目標(biāo)對照明光的反射在向傳感器傳輸?shù)倪^程中被散射而導(dǎo)致成像細(xì)節(jié)模糊或丟失，稱為前向散射，而照明光在到達(dá)目標(biāo)之前也受到水體的散射而被傳感器接收形成一種“霧化”背景，降低了圖像對比度，稱為后向散射。散射介質(zhì)的吸收以及前向散射和后向散射的同時作用導(dǎo)致水下圖像嚴(yán)重降質(zhì)，從而限制了水下成像系統(tǒng)的觀測距離。雖然提高照明功率可以增加目標(biāo)光強度，但是后向散射光的強度也同時增加，照明功率達(dá)到一定水平后，進(jìn)一步增大照明功率并不能提高圖像對比度和信噪比，所以后向散射在很大程度上成為限制觀測距離的重要因素。光的散射特性對成像質(zhì)量有著嚴(yán)重影響，突出的問題是水下成像清晰度差，分辨率低，傳統(tǒng)的圖像增強和復(fù)原技術(shù)都需要建立有效的點擴(kuò)散函數(shù)模型。有關(guān)水下圖像復(fù)原的研究基本上都集中在點擴(kuò)散函數(shù)（PSF）的描述上。點擴(kuò)散函數(shù)描述的是點光源在水下成像過程中因前向散射而形成的擴(kuò)散現(xiàn)象，取決于水體的固有光學(xué)特性以及成像距離。一般采用退化函數(shù)將退化過程模型化，它和加性噪聲項一起，作用于輸入圖像[fx，y]，產(chǎn)生一幅退化的圖像[gx，y]：

[gx，y=Hfx，y+ηx，y]

若H是線性的，空間不變的過程，則退化圖像在空間域形式為：

[gx，y=hx，y*fx，y+ηx，y]

空間域的卷積和頻域的乘法組成了一個傅里葉變換對，所以可以用等價的頻域表示：

[Gu，v=Hu，vFu，v+Nu，v]

式中大寫字母表示的項是卷積方程式中相應(yīng)的傅里葉變換。退化函數(shù)[Hu，v]有時稱為光學(xué)傳遞函數(shù)（OTF），在空間域，[hx，y]稱為點擴(kuò)散函數(shù)（PSF）。因此，復(fù)原圖像f可以通過線性系統(tǒng)反卷積得到。點擴(kuò)散函數(shù)在圖像復(fù)原及增強領(lǐng)域都有重要的作用，也是圖像退化的主要成因。目前常用的圖像復(fù)原方法，一般可分為參數(shù)法與迭代法，然而這兩種方式都要對點擴(kuò)散函數(shù)進(jìn)行估計和計算，要獲得水下清晰化的圖像，點擴(kuò)散函數(shù)有著極其重要的作用。在水下圖像的復(fù)原研究方面，研究的重點一直集中在 PSF 的描述上。點擴(kuò)散函數(shù)PSF的準(zhǔn)確性極大程度地影響著圖像復(fù)原質(zhì)量。探討點擴(kuò)散函數(shù)的影響因素對進(jìn)一步進(jìn)行水下圖像研究有十分重要的理論指導(dǎo)意義。與以往水下圖像復(fù)原技術(shù)的研究角度不同的是，本實驗從點擴(kuò)散函數(shù)本身的特性著手，主要以點擴(kuò)散函數(shù)的成像影響，及水下環(huán)境目標(biāo)對點擴(kuò)散函數(shù)的影響性為主要內(nèi)容，探討水下成像過程中顏色對于點擴(kuò)散函數(shù)的影響性。

1 實驗環(huán)境

散射介質(zhì)光傳輸特性的研究以實驗為主，現(xiàn)有的關(guān)于散射機理的理論，包括經(jīng)驗或半經(jīng)驗性結(jié)論作為分析的參考依據(jù)。實驗借助標(biāo)準(zhǔn)光學(xué)水池和自制實驗箱進(jìn)行。在光學(xué)水池內(nèi)添加泥沙或其他懸浮雜質(zhì)，待沉淀不同時間后，可以形成不同能見度的成像環(huán)境。使用已有的水下攝像系統(tǒng)（束光源照明，CCD與光源分離）采集視頻圖像。此項研究的主要目的是分析光在散射介質(zhì)中傳輸時所發(fā)生的前向和后向散射的一般性統(tǒng)計規(guī)律，以此作為散射介質(zhì)圖像復(fù)原算法中的普適性先驗知識。實驗將在不同能見度、不同照明強度、不同傳輸距離的可變條件下，針對圖像退化復(fù)原問題進(jìn)行探討，應(yīng)用水池基本數(shù)據(jù)為長4.5 m，寬2 m，高2 m，每次實驗用水約15 t，水池三面帶有4個觀測窗口，室內(nèi)有水源和下水道，并配有三相交流電源（見圖1）。通過自制帶孔圓桶創(chuàng)造散射條件，圓孔半徑5 cm，并應(yīng)用不同顏色的玻璃紙置于小孔前制造顏色差異（見圖2）。應(yīng)用水下攝像機拍攝特定條件下的動態(tài)視頻，根據(jù)需求截取相應(yīng)的實驗圖片。

[圖1 三相交流電源 圖2 帶孔圓桶]

2 計算依據(jù)及方法

2.1 點擴(kuò)散函數(shù)測量

點擴(kuò)散函數(shù)估計是圖像復(fù)原的敏感問題，與前向和后向散射的成像描述直接相關(guān)。已有的各種理論或經(jīng)驗?zāi)Ｐ停ㄔ谇捌谘芯恐刑岢龅拿枋鰡挝缓穸人w散射的高斯模型。其適用范圍并沒有經(jīng)過系統(tǒng)的實驗驗證，尤其在渾濁度較高的強散射環(huán)境中，PSF會呈現(xiàn)出更為復(fù)雜的形式。比如發(fā)現(xiàn)介質(zhì)渾濁程度越高，PSF分布中所包含的隨機成分就越明顯，PSF中隨機成分的存在，導(dǎo)致成像過程中產(chǎn)生除了后向散射以外新的圖像噪聲，也就是說，目標(biāo)反射信號中其實也存在噪聲源，其機理源于散射過程的隨機性。這一問題在以往的研究中都被忽略了。這里所關(guān)心的是：此種噪聲對成像有多大的影響以及在復(fù)原過程中如何進(jìn)行抑制。最終在實驗的基礎(chǔ)上，提出一個適用性強、計算復(fù)雜度低的更有效的點擴(kuò)散函數(shù)的描述模型。

本實驗的主要目的是對點擴(kuò)散函數(shù)數(shù)據(jù)進(jìn)行分析，得到點擴(kuò)散函數(shù)分布情況是實驗進(jìn)行的基礎(chǔ)。本實驗計算點擴(kuò)散函數(shù)的基本依據(jù)為：一般將輸入圖像[ft]和退化圖像[gt]應(yīng)用公式表示為：

[gt=Hft+ηt]

在空間域可以變化為：[gt=ht*ft+ηt]。其中：[ηt]表示加性噪聲項；[ht]稱為點擴(kuò)散函數(shù)（PSF）。在本次實驗中，假設(shè)加性噪聲項不存在，即忽略加性噪聲的影響，則有：[gt=ht*ft]，由卷積的基本定理可知，對[gt=ht*ft]兩邊同時求導(dǎo)可得：[g′t=f′t*ht=ft*h′t]；其中，[ft]為階躍函數(shù)，根據(jù)階躍函數(shù)的性質(zhì)可以知道，[ft]的導(dǎo)數(shù)為δ函數(shù)，則[g′t]可重新表示為[g′t=δt*ht]；由δ函數(shù)的基本性質(zhì)，δ函數(shù)與其他函數(shù)的卷積為其他函數(shù)本身，因此有結(jié)論：

[g′t=δt*ht=ht]

即對退化圖像求導(dǎo)可得到點擴(kuò)散函數(shù)。

2.2 差分法

在Matlab當(dāng)中，讀取的圖像是以矩陣的形式儲存的，矩陣中的元素是獨立的像素坐標(biāo)，在此次實驗中，要求取沿著一幅圖像某一確定方向上，多個點處的點擴(kuò)散函數(shù)，因此不能應(yīng)用一般的連續(xù)求導(dǎo)方式，在此引入了差分法近似計算導(dǎo)數(shù)。差分法將連續(xù)的定解區(qū)塊用有限個離散點構(gòu)成的方陣代替，當(dāng)中的所有離散點稱作為節(jié)點；將連續(xù)區(qū)間中的連續(xù)變量函數(shù)用方陣中的離散變量函數(shù)近似表示，將求解方程條件中的微分商用差商來近似求取，積分應(yīng)用積分和近似，所以原微分方程和定性條件可以近似的用代數(shù)方程組表示，解此方程組便可以得到離散點處的近似解。由于圖像像素間的間距很小，用差商近似表示導(dǎo)數(shù)的誤差在可控的范圍內(nèi)，對實驗數(shù)據(jù)的準(zhǔn)確性不會產(chǎn)生明顯影響。

3 實驗?zāi)康?/p>

目前國內(nèi)外關(guān)于散射介質(zhì)圖像處理的研究多集中在處理方法的研究上，而有關(guān)散射介質(zhì)成像的機理或規(guī)律等基礎(chǔ)性研究很不成熟。

本實驗進(jìn)行的相關(guān)實驗研究，不僅作為水下成像過程分析和處理的依據(jù)，也為散射介質(zhì)單幅圖像復(fù)原的普適性方法研究提供了有價值的參考。本次實驗將圖像的顏色作為假定影響因素，應(yīng)用Matlab處理相關(guān)數(shù)據(jù)結(jié)果。實驗選取圖像圓中心點坐標(biāo)作為計算起始點，沿一個固定方向求取圖像矩陣中不同點處的點擴(kuò)散函數(shù)分布情況。在實驗結(jié)果中，將數(shù)據(jù)做曲線圖形式表示，包含3種顏色圖像的點擴(kuò)散函數(shù)分布情況，通過比較，得到圖像自身顏色的影響性結(jié)論。

4 實驗結(jié)論

實驗應(yīng)用紅，綠，藍(lán)，3種顏色玻璃紙制造顏色差異，在清水環(huán)境下中拍攝水下圖像如圖3所示。

[（a） （b）]

（c）

圖3 清水環(huán)境中拍攝的水下圖像

實驗主要研究在圓邊緣處點擴(kuò)散函數(shù)的越變情況，并以此分析在同樣水域環(huán)境下的點擴(kuò)散函數(shù)的差異性，分析顏色對點擴(kuò)散函數(shù)的影響性。在計算過程中，選取了一個固定方向作為參考，本實驗選取了x軸正向45°作為計算方向，通過統(tǒng)計學(xué)方式根據(jù)差分法求取這一方向上多個像素點處的點擴(kuò)散函數(shù)，總結(jié)計算公式為：

[hm=bx-m+1，y+m-1-bx-m-1，y+m+1]

式中：x，y為選取的起始點坐標(biāo)；m為從1開始的正整數(shù)，計算的點擴(kuò)散函數(shù)通過plot函數(shù)將結(jié)果表示成曲線形式如圖4～圖6所示。

圖4 圖像（圖3（a））曲線分布

圖5 圖像（圖3（b））曲線分布

圖6 圖像（圖3（c））曲線分布

通過觀察3幅曲線圖的峰值，可以發(fā)現(xiàn)，在顏色不同的圖像中，點擴(kuò)散函數(shù)的分布形式大致相同，并且在圓邊緣處的越變數(shù)值基本接近，由此，可以得到結(jié)論：在外部環(huán)境條件相同的情況下，顏色對點擴(kuò)散函數(shù)沒有明顯影響。

5 結(jié) 語

點擴(kuò)散函數(shù)是圖像處理的一個核心問題，應(yīng)用于圖像處理的多個方面。對提高圖像清晰度，散射環(huán)境成像研究等都有重要的應(yīng)用價值。廣泛應(yīng)用于與光學(xué)成像系統(tǒng)，工業(yè)成像研究等圖像領(lǐng)域，是圖像復(fù)原研究的重要理論模型。散射是自然界普遍存在的光學(xué)現(xiàn)象，解決散射對成像系統(tǒng)的影響對現(xiàn)代圖像處理有深遠(yuǎn)的意義。國內(nèi)外的研究，由于各種光學(xué)系統(tǒng)采用的期間并不相同，點擴(kuò)散函數(shù)的形式必然擁有差異，同時PSF的構(gòu)造原理，方法和推導(dǎo)方式也不盡相同，目前，還沒有一套完整的確定PSF的基礎(chǔ)理論。

本次實驗從PSF的成因角度，探討顏色對PSF的影響，是初步對PSF自身理論的研究，根據(jù)實驗數(shù)據(jù)，得到顏色并沒有明顯影響PSF的結(jié)論。但是實驗有不足及改進(jìn)空間的，比如對環(huán)境，光照的要求，數(shù)據(jù)的精準(zhǔn)性，算法的改進(jìn)等方面，這是接下來工作的目標(biāo)和方向。

參考文獻(xiàn)

[1] CASTLEMAN K R.數(shù)字圖像處理中文版[M].北京：電子工業(yè)出版社，2002.

[2] 王風(fēng)鵬.用CCD測定光學(xué)系統(tǒng)的點擴(kuò)散函數(shù)[J].贛南師范學(xué)院學(xué)報，2005（6）：17?18.

[3] 李奇，馮華君.用于全數(shù)字對焦的點擴(kuò)散函數(shù)性能分析與評價[J].浙江大學(xué)學(xué)報：工學(xué)版，2006，40（6）：1093?1096.

[4] 孫傳東，陳良益，高立民，等.水的光學(xué)特性及其對水下成像的影響[J].應(yīng)用光學(xué)，2000，21（4）：39?46.

[5] 阮秋琦.數(shù)字圖像處理[M].北京：電子工業(yè)出版社，2001.

[6] 李建.水下圖像后向散射噪聲的去噪問題研究[D].青島：中國海洋大學(xué)，2009.