回歸擬合NR函數(shù)及GPDR先驗的圖像霧濃度檢測

中圖分類號：TP391.9 文獻標志碼：A 文章編號：1000-582X（2025）04-115-12

Inspection of image fog Concentration using regression-fitting NR function and GPDR prior

WEN Limin*b， WANG Huifenga， JU Yongfenga （a.School of Electronic amp; Engineering； b. Experimental Teaching Center of Electronics and Electronics， Chang'an University， Xi'an 710064， P.R.China）

Abstract：Addressing the limitations of fog concentration inspection in image defogging，an algorithm based on the scatterplot prior of the generalized pixel diference-ratio（GPDR）and the Naka-Rushton（NR）fitting function was proposed.First，the GPDR prior for gray scaterplots in standard foggy image sets across various scenes was extracted.Next， the NR function，constrained by the prior， was introduced， and a lookup table of parameters （n，k） corresponding to fog concentration levels was established by calculating the parameters （n，k） of NR function for standard image sets. Regression analysis was then used to calculate the parameters for real foggy images， and the comprehensive correlation coefficient between （n，k） and was calculated. Parameters （n，k） with correlation coeffcients exceeding a set threshold were considered indicative of the fog concentration level. Simulations show that the algorithm accurately reflect changes in fog concentration across images with varying densities.Additionally，correlation coefficients between the algorithm's results and PM2.5 measurements reached up to 0.95，both within the same andacrossdiffrent scenes.This shows thatthealgorithm can be effectively used for fog concentration rating in visual field.Horizontal comparison tests show that the inspection accuracy of the proposed algorithm can reach up to 4.8% ， making it suitable for field fog concentration detection. Keywords： fog concentration； LIVE library； GPDR； Naka-Rushton function； inspection

霧是生活中常見的自然現(xiàn)象，其存在會影響人們的生產(chǎn)和生活，通常這種影響是負面的，需要減緩或消除。中國是霧的多發(fā)區(qū)，特別是在南方很多山區(qū)，霧給人們生活帶來的隱患到了不容忽視的地步\"。因此，研究如何預(yù)警和預(yù)報霧是當前亟待解決的問題。然而預(yù)警和預(yù)報的前提是正確評定霧濃度的等級，目前有多種用于霧濃度等級評定方法，如基于激光、紅外技術(shù)檢測等[但該類方法是通過研究傳播介質(zhì)布朗運動對激光或紅外光路的遮斷或衰減特性來實現(xiàn)，通常在較長物距且需要主動發(fā)射或接收裝置才能實現(xiàn)，適合特定場合，因此，研究基于機器視覺或圖像處理的霧濃度檢測具有現(xiàn)實意義。

目前國內(nèi)外很多學者從事基于機器視覺或圖像處理技術(shù)霧濃度檢測研究，并取得了豐碩成果，Miital等是較早開展相關(guān)研究的學者，其主要是采用歸一化亮度系數(shù)法（mean subtracted contrast normalized，MSCN）實現(xiàn)圖像對比度的評定，雖然該方法只是針對圖像對比度評價，并未與視場霧濃度檢測聯(lián)系起來，但為霧濃度檢測評價提供了參考依據(jù)。溫立民等采用對比度、信息建立多維向量，通過求解空間馬氏距離的方法映射視場霧濃度，取得較好效果，但檢測精度略顯不足。Cui等利用暗盒方法，通過控制盒內(nèi)水蒸氣（代替霧）濃度，采集視場圖像并研究影響霧濃度的參數(shù)因素，但其測試是在限定背景參數(shù)的條件下獲得，不具有一般性。李春明等通過統(tǒng)計霧濃度與圖像色度、亮度和飽和度之間的關(guān)系，建立回歸方程并求解視場霧濃度指數(shù)，雖然該方法在特定環(huán)境下的霧濃度檢測具有較高精度，但因?qū)舛取⑸取⒘炼燃帮柡投戎g的統(tǒng)計具有片面性，使其在應(yīng)用上不具有通用性。Guo等基于暗原色先驗方法建立邊緣信息索引并利用索引值解算霧濃度參數(shù)，雖在中低濃度檢測中取得較好效果，但因其是采用暗原色先驗方法實現(xiàn)霧濃度檢測，難以避免該方法的弊端，即在高濃度或高亮背景下失效。Ling等將暗原色先驗、統(tǒng)計離差、隨機深林等作為霧圖像相關(guān)統(tǒng)計參量，構(gòu)建線性組合方程求解視場霧濃度。Jiang等研究了視場景深與圖像色彩、色度等參數(shù)之間的關(guān)系，通過回歸分析實現(xiàn)霧濃度計算，雖取得一定效果，但需要大量能反映現(xiàn)場景深的參考霧圖。因此，研究提出基于標準圖像庫的廣義灰度差-比（generalized pixel diffrence-ratio，GPDR）先驗，建立基于 Naka-Rushton（NR）函數(shù)的3段擬合曲線，計算標準圖像散點GPDR值的擬合NR函數(shù)參數(shù) （n，k） 與真實圖像擬合 的相關(guān)系數(shù)，建立參數(shù) （n，k） 與霧濃度對應(yīng)查找表和聞值門限，通過遍歷查找表并判斷綜合相關(guān)系數(shù)是否高于門限來確定視場霧濃度等級。經(jīng)相同場景不同濃度霧圖像測試，證明圖像集擬合Naka-Rushton函數(shù) （n，k） 與真實霧圖像 高度相關(guān)，可以作為視場霧濃度等級評定。

1基于標準圖像庫灰度差-比散點先驗

1.1 影響圖像對比度的因素

通常霧的濃度是反映在圖像的對比度和信息上，對比度和信息隨著霧濃度增高而降低。霧濃度大小不僅與視場中傳播的介質(zhì)有關(guān)，還與光照條件有關(guān)，傳播介質(zhì)是勻質(zhì)的還是非勻質(zhì)的，光照條件是各向同性還是各向異性都會影響霧濃度的大小[-13]。要通過圖像實現(xiàn)視場霧濃度檢測，首先要研究對比度和信息隨霧的變化規(guī)律，這需要以獲得具有不同屬性的有霧和無霧圖像樣本為前提，而現(xiàn)實中獲得各種不同性質(zhì)、不同光照條件的圖像較為困難，本文在標準圖像數(shù)據(jù)庫（LIVE）中選取參考圖像作為測試樣本，建立參照標準實現(xiàn)霧濃度檢測。當前通用標準圖像庫有LIVE、FRIDA和 。其中，LIVE庫擁有豐富的資源，收集了各種圖片資源，僅有霧圖像就分為各向同性、各向異性介質(zhì)和光照等不同條件下圖像，適合各種圖像處理和機器視覺應(yīng)用，具有權(quán)威性。

1.2基于標準有霧和無霧圖像的GPDR先驗

研究選取4組標準圖像集（記為I）如圖1（a）～（d）所示，每組圖像集場景和光照條件相同或相近，每個場景包括有霧（下排）和無霧圖像（上排），圖1（a）為霧和光照都是各向異性， PM2.5=35 μg/m³"；圖1（b）為霧是各向異性，光照是各向同性， PM2.5=67μg/m³ ；圖1（c）為霧是各向同性，光照為各向異性， PM2.5=80 μg/m³ 圖1（d）為霧和光照均為各向同性， PM2.5=120μg/m³ 。每個場景有霧和無霧圖像各100張，其霧濃度（a） lt; （以每組標準圖像的 值為準）。圖 1（e）～（h） 分別為（a）～（d）對應(yīng)的散點圖，散點圖以圖像灰度級[0，255]作為 x 軸，以同場景有霧圖像與無霧圖像歸一化廣義灰度差-比（GPDR）作為y軸（廣義灰度差-比為有霧圖像灰度值與對應(yīng)的無霧圖像灰度值的差與廣義有霧圖像灰度值之比，其中廣義有霧圖像灰度值為有霧圖像灰度與天空光的差值），藍色曲線為100幅圖像對應(yīng)的散點分布，紅色直線表示過渡階段對應(yīng)的線性區(qū)域。由圖 1 （e）～（h）知，無論是各向同性還是各向異性的霧圖，其GPDR都呈S形曲線，即在低灰度區(qū)其值較低且緩慢上升，呈現(xiàn)下凹型態(tài)；在高灰度區(qū)卻是呈現(xiàn)緩慢上凸形態(tài)（如圖紅色箭頭所示），且隨著霧濃度的增大，下凹和上凸越明顯。中間段為近似的線性，隨著霧濃度的增加，線性度越明顯，通常下凹與上凸之間存在過渡點，且過渡點的斜率隨著霧濃度的增加呈上升趨勢（如圖 。

從以上分析可以得到一個關(guān)于散點圖先驗：霧的濃度與S型的散點圖一一對應(yīng)，即一定形狀的S型散點圖對應(yīng)一定的視場霧濃度。因此，根據(jù)此散點圖先驗，可通過改變凹凸及線性度等參數(shù)確定一個S型函數(shù)，表征視場霧濃度大小。為此研究引入Naka-Rushton函數(shù)來表征這種S型散點形狀，確定對應(yīng)標準圖像的霧濃度。

對于真實的圖像1（j），同樣可繪制有霧圖像與無霧圖像間GPDR值S型曲線，如圖1（i）所示。根據(jù)標準圖像集獲得的S型函數(shù)可確定一個濃度的參照標準，比如某參數(shù)下確定的S型函數(shù)對應(yīng) PM2.5=50 μg/m³，如果真實霧圖對應(yīng)的S曲線與之匹配程度達到一定級別時，則認為此時的真實霧圖的濃度即為 PM2.5=50 μg/m³。

1.3廣義灰度差-比先驗證明

廣義灰度差-比的先驗是通過統(tǒng)計200幅標準圖像得出的先驗規(guī)律，對該先驗證明其具有一般性。在圖像去霧領(lǐng)域中，應(yīng)用較為廣泛的為Narasimhan、Nayar等提出的大氣散射模型，如下所示

其中： I（x，λ） 為采集器采集到的有霧圖像；為人射光波長； 為同場景下的無霧圖像； 為天空光，對于攝像頭采集到的圖像，天空光可用最大像素值或最高 10% 的像素均值代替，對于式（1）可認為 為已知常數(shù)； β（x，λ） 為傳播介質(zhì)的衰減系數(shù)； d（x） 為視場景深函數(shù)； k 為常數(shù)。對式（1）整理得

進一步有

若定義式（3）等號左邊為灰度差-比函數(shù)，則由式（3）等號右邊可知，該函數(shù)等價于過原點的超越函數(shù)（或稱類指數(shù)函數(shù)），該函數(shù)與指數(shù)函數(shù)類似，這是散點圖出現(xiàn)類線性的原因。而式（3）與指數(shù)函數(shù)不同的是隨著x 的增大而不是無限增大，是有界函數(shù)。證明如下：由文獻[14]知

衰減系數(shù)是波長與景深的函數(shù)，衰減系數(shù)與景深成反比，則有

因此式（3）等號右邊 為有界函數(shù)，廣義灰度差-比函數(shù)為S型函數(shù)。若令

其中 I（d（x）） 即為廣義像素灰度差-比函數(shù)，可表示為有霧圖像灰度 I（x，λ） 與無霧圖像灰度 的差比與廣義有霧圖像灰度差 （等價于有霧圖像灰度與常數(shù)之差，稱為廣義灰度）。因為式（3）是個超越方程，所以筆者采用擬合分析方法計算霧濃度，即首先建立廣義灰度散點圖，選取合適的參數(shù)點建立擬合函數(shù)。因為圖像灰度值與景深存在正相關(guān)，所以圖1以圖像灰度值為橫軸，廣義像素灰度-比為縱軸建立散點圖。由于衰減系數(shù)的冪指數(shù)是景深的函數(shù)，致使在景深較小（灰度較低）時曲線出現(xiàn)不同程度的下凹，在景深較大（灰度較高）時出現(xiàn)不同程度的上凸。

2 改進的Naka-Rushton函數(shù)擬合

2.1 基本Naka-Rushton擬合函數(shù)

從分析可知，S型散點圖為100幅圖像在灰度級上的GPDR分布，由于這種散點圖分布呈現(xiàn)如圖1（e）～（h）粗線條軌跡，而研究是采用回歸分析法擬合曲線，在每個灰度級上選擇的回代樣本將影響最終擬合結(jié)果，不同的回代樣本點將擬合出不同曲線，本文以每個灰度級GPDR均值作為回代樣本，擬合NR函數(shù)表征視場霧濃度的S型曲線，計算霧濃度等級。而擬合曲線的選擇，需要滿足以下約束條件：

1）擬合曲線需整體呈現(xiàn)S形狀，并存在大范圍的近似線性區(qū)域；

2）中間線性區(qū)過渡點可調(diào)整，且過渡點的切線斜率可調(diào)整；

3）在低值區(qū)域函數(shù)值很小或近似為0且呈現(xiàn)下凹特性，在高值區(qū)域近似滿幅且呈現(xiàn)上凸特性。

引人Naka-Rushton（NR）擬合函數(shù)如圖2所示，一是因為NR函數(shù)曲線圖形為S型；二是因為該函數(shù)滿足上述3條約束。如圖2（a）所示，該函數(shù)可分為3段區(qū)間，即低值區(qū)，線性區(qū)和滿幅區(qū)。低值區(qū)幅值可近似為0，在這段區(qū)域內(nèi)，是靠近采集設(shè)備區(qū)域；滿幅值區(qū)間是灰度值最高的區(qū)間，是無窮遠端或天空光區(qū)域（滿足約束條件3））；函數(shù)存在大范圍線性區(qū)，且線性區(qū)范圍和斜率可根據(jù)需要調(diào)整（滿足約束條件1），2））。基于上述2點原因，采用NR函數(shù)擬合標準圖集GPDR散點圖估計視場霧濃度。

Naka-Rushton（NR）函數(shù)形如

式中： 為NR函數(shù)參數(shù)， 為NR函數(shù)上界； k 可實現(xiàn)對凹凸過渡點的調(diào)整； n 為速降因子，用于調(diào)整凹凸的幅度。當用NR函數(shù)擬合散點圖時，不同霧濃度可計算出不同的 n，k 值。如果通過回歸分析法計算出 ，也就確定了當前霧圖像對應(yīng)的視場霧濃度大小，即一個 （n，k） 對應(yīng)一個視場濃度，不同 （n，k） 對應(yīng)不同的霧濃度（如圖2（b）所示），圖 2（b）x"軸為圖像灰度級， y 軸為歸一化廣義灰度差-比（GPDR），其中的藍、黑、綠、紅色曲線表示在保持 k 為定值時，改變 n 值為2.5，3，5，7時的NR曲線，由圖可知隨著 n 值的增大，曲線的過渡區(qū)域增加，上凸和下凹的幅度增大。

若將式（7）等號兩邊同除 ，有

變形后

將等號兩邊取對數(shù)

式（10）中， 分別為散點圖中灰度樣本值，該樣本值可從散點圖中直接獲得，若將樣本值代入式（10）即可估計出參數(shù) 。

2.2 Naka-Rushton擬合函數(shù)的改進

根據(jù)式（1）估計出的 （n，k） 值對應(yīng)的擬合NR函數(shù)與散點區(qū)在局部上有時會存在較大偏差，主要表現(xiàn)在2個方面： ① 當灰度 xk 時呈現(xiàn)上凸特性。因為， x 是在灰度[0，255]范圍內(nèi)變化，導(dǎo)致曲線下凹和上凸的變化趨勢過快，這不符合視場霧濃度的變化規(guī)律； ② 如果 取圖像像素的最大差值，將導(dǎo)致GPDR值整體偏高，需要對Naka-Rushton函數(shù)作以改進。對于問題 ① ，改進如下

改進方法如式（11）所示，將整個灰度區(qū)域[0，255]分為3段區(qū)間，分別為：低值區(qū) 近似線性區(qū) 和高值區(qū) ，如式（11）所示在低值區(qū)將分母中的 x 變?yōu)?x/2 ；在近似線性區(qū)保持不變；在高值區(qū)將 x 變?yōu)?倍。經(jīng)此改進后在低值區(qū)的下凹和高值區(qū)的上凸變緩，更符合灰度散點的變化規(guī)律。

對于問題 ② ，需要對 根據(jù)下式進行補償。

式（12）中： 為圖像的最大GPDR值點； ε 為補償量，該值需根據(jù)當前散點圖所處的灰度級范圍確定，補償經(jīng)驗公式如下

式（13）中： x 為灰度； ε 為補償量，則改進后NR函數(shù)為

當確定修正量ε后，可根據(jù)式（12）計算最終 ，并將 代人公式（14），可實現(xiàn) 參數(shù)估計。

圖3為改進前后的對比測試結(jié)果，紅色線為濃度 PM2.5=135 μg/m³的散點圖，綠色線為改進前擬合生成的NR函數(shù)，由圖3可知擬合函數(shù)GPDR值在低值區(qū)下凹處幅值過低，而在高值區(qū)上凸處幅值過高，導(dǎo)致這2段擬合誤差較大。改進后的擬合曲線為藍色線，由圖3可知，紅色線和藍色線基本重合，表明擬合曲線與散點圖具有較高的重合度，消除了下回處幅值過低和上凸處過高的不足。因此，經(jīng)過改進的NR函數(shù)，符合S型散點圖擬合約束要求。

3NR函數(shù)霧濃度的參數(shù)估計

由Naka-Rushton函數(shù)擬合GPDR均值散點圖曲線原理可知，擬合的關(guān)鍵是如何確定參數(shù) （n，k） 。由式（14）和參數(shù) （n，k） 通過回歸分析法估計出來，即在標準圖像散點中選取關(guān)鍵點坐標代入式（14）中，估計參數(shù) （n，k） 。具體步驟如下，流程如圖4所示。

1）計算每個灰度級對應(yīng)的標準散點的GPDR均值 ， x 為灰度級， 為灰度級對應(yīng)的GPDR均值。遍歷整個灰度級范圍，找到中最大均值 ，將此最大均值按公式（12）（13）進行補償，得到最大估計

2）在步驟1）均值坐標 中確定關(guān)鍵坐標樣本點。因為研究采用回歸分析法估計參數(shù) （n，k） ，則選擇代入公式（14）的散點坐標十分重要，如果坐標樣本點選擇不合適，估計的參數(shù) （n，k） 偏差大，則擬合出的Naka-Rushton函數(shù)就不能正確反映視場散點規(guī)律，因此，散點樣本坐標應(yīng)該在下凹處、上凸處和過渡點處選取。

3）將步驟1）2）確定的 和散點均值坐標 分別代人公式（14）計算 （n，k） ，進一步計算標準圖像集散點擬合NR函數(shù)。不同濃度的圖集散點對應(yīng)不同的擬合函數(shù)和 （n，k） ，即霧濃度等級與 （n，k） 一一對應(yīng)，如果將每個濃度所對應(yīng)的 （n，k） 通過上述方法全部估計出來，可建立一個查找表（如表1所示，表1只給出部分（n，k） 值），以這個查找表作為尺度衡量真實有霧圖像的霧濃度。表1是在標準圖像集中選取一定數(shù)量的霧圖，并將霧圖按PM2.5值等間隔分組，間隔值約為0.5，分別計算出每組的 （n，k） 值，建立PM2.5與 （n，k） 查找表。

4）查找表建立后，可以將查找表作為標尺衡量真實圖像霧濃度，步驟如下：

a.參照步驟1）建立真實圖像（記為 有霧和無霧圖像GPDR散點圖如圖1（i）（j）所示b.采用回歸分析法計算 的Naka-Rushton擬合曲線，并計算相應(yīng)的參數(shù)估計，記為 ：

c.分別計算標準圖像集 （n，k） 和真實圖像參數(shù) 的相關(guān)系數(shù)

式（15）（16）中， 為 n 和 k 各自的相關(guān)系數(shù)，cov 為協(xié)方差， 和 為標準差；d.計算綜合相關(guān)系數(shù) ρ

式（17）中 為權(quán)重系數(shù)，因為在Naka-Rushton函數(shù)中，決定曲線形狀的參數(shù) n 起的作用要大些，所以經(jīng)實驗測試 ；

e.建立遍歷參照表1搜索進程，搜索綜合相關(guān)系數(shù) ρ 大于門限的 （n，k）

f.查找表1確定當前霧濃度指數(shù)。

4仿真實驗

4.1 實驗條件

研究實驗仿真上位機平臺為CPUIntel酷睿i7內(nèi)存16G，獨立顯卡，64位MicrosoftWindows7；測軟件：Matlab2016。選取的測試圖像為標準LIVE圖像和實際隨機采集的霧圖，圖像大小從 400×300～1024×768 不等。

4.2霧濃度與NR函數(shù)相關(guān)性測試

為驗證基于NR函數(shù)擬合能否用于視場霧濃度計算，需要對標準圖像NR函數(shù)擬合的 （n，k） 與真實圖像擬合 進行相關(guān)性測試。研究選擇某地連續(xù)14d采集的相同場景不同霧濃度圖像作為測試樣本，將這14幅采集樣本按PM2.5值由低到高的順序排列（如圖5所示）[7-18。以氣象預(yù)報的PM2.5濃度歸一化值為橫坐標，將圖像按PM2.5由小到大的順序排列，從左至右從上到下編號1～14。因為視場霧濃度大小在圖像上體現(xiàn)在對比度上，以對比度歸一化值作為縱坐標。對比度計算公式°為

式（18）中： ，即相鄰元素幅度差值； 為幅度差值為 δ 的概率；

和 為行列最大。

測試數(shù)據(jù)如表2所示， PM2.5 和測試結(jié)果都是歸一化值，若以PM2.5為 x 軸，對比度值為y軸做散點圖，如圖6所示。14d的 和測試結(jié)果整體上呈現(xiàn)線性，經(jīng)計算二者的綜合相關(guān)系數(shù)可達0.95，說明散點斜率與視場濃度高度相關(guān)，因此，算法可以作為視場霧濃度的度量。

4.3不同場景不同濃度霧測試

為測試算法在不同場景下的檢測效果，隨機選取了8組不同濃度的場景圖像（如圖7所示），包括無霧（上排）和有霧（下排）圖像，采用研究算法對有霧和無霧圖像分別進行計算。從左至右按有霧濃度由小到大順序排列，分別編號為場景序號1～8，測試濃度（PM2.5）如表3所示。圖8中為不同濃度不同場景對比曲線，圖 8：x 軸為場景序號，y軸為PM2.5。無霧圖像的測試結(jié)果通常較低，濃度值大都在70以下，按PM2.5空氣質(zhì)量標準都在優(yōu)和良，這一點從無霧圖像也能反映出來，圖像質(zhì)量清晰；而有霧圖像測試值則隨著對比度減小（霧濃度增大）呈遞增趨勢，因此，實驗證明算法對不同濃度、不同場景同樣適用。

4.4 橫向比較

當前，采用圖像處理的方法實現(xiàn)霧濃度檢測的途徑是基于圖像對比度的評價，為更好地評估算法有效性，將對比度評價較好的成果與研究作對比。目前用于對比度或霧濃度評價效果較好的有MSCN、文獻[9]和[15]，文獻[9]是基于自然統(tǒng)計特性實現(xiàn)霧濃度檢測，文獻[15]是基于空間距離計算視場霧濃度。為比較時具有共同的衡量尺度，采用圖像的MOS（mean opinion scores）值3作為參照標準。MOS法是對比度評價中普遍采用的圖像主觀評價法，其值能較為準確反映圖像對比度大小，因此，以MOS值作為標準尺度。由于MSCN和文獻[15]針對的是圖像對比度，在比較時統(tǒng)一化為對比度。

眾所周知，不同背景、不同紋理信息對于霧濃度檢測結(jié)果會產(chǎn)生影響，為消除干擾，保證測試一致性，選用同場景、同光照、同大小不同濃度的霧圖作為測試樣本，將測試樣本從左到右從上到下按濃度由小到大的順序排列1～14#（如圖9所示），分別采用文獻[9，15]、MSCN及研究方法進行計算，測試結(jié)果如圖10（a）所示。圖10（a）中： x 軸為測試樣本序號；y軸為測試對比度值；黑色十字星型曲線為標準MOS測試結(jié)果；紅色星型曲線為文獻[15]測試結(jié)果；紫色曲虛線為MSCN的測試結(jié)果；藍色點劃畫線為文獻[9]的測試結(jié)果；綠色叉實線為測試結(jié)果。從圖10（a）可知MSCN和文獻[15]相比于MOS測試結(jié)果相差較大，整個測試時間段內(nèi)對比度值都小于MOS方法，表明MSCN與文獻[15]對比度評價結(jié)果相比文獻[9]和研究結(jié)果要差些；比較研究與文獻[9]，從整體上看研究波動小些，文獻[9]波動范圍較大，如1～10幅圖，偏離MOS值較大，特別是第10幅達到最大，因此，研究測試結(jié)果要優(yōu)于文獻[9]。

圖10（b）為文獻[9]與研究測量誤差對比， x 軸為樣本序號，y軸為對比誤差。為保證一致性，測試樣本與文獻[9]相同，測試結(jié)果為連續(xù)14d某地PM2.5值（以當天天氣預(yù)報為準）。其中：紅色星線為文獻[9]的誤差；綠色的叉實線為研究的測試誤差。由圖10（b）可知，文獻[9]與研究誤差都隨著濃度增而有所增加，但各點的誤差文獻[9]要高于研究結(jié)果。在濃度最大時誤差也達到最大，研究結(jié)果只有 4.8% ，而文獻[9]誤差達 10% ，因此，無論是從檢測效果上還是誤差上對比，研究都要優(yōu)于其他方法。

5結(jié)論

研究提出無霧圖像集與有霧圖像集灰度差-比散點先驗原理，引人NR函數(shù)擬合散點圖，建立霧濃度查找表，實現(xiàn)視場霧濃度等級評定。仿真實驗表明：研究算法測試結(jié)果與視場濃度高度相關(guān)，相關(guān)系數(shù)達0.95，通過同場景、不同濃度，不同場景、不同濃度樣本測試，及與現(xiàn)有的測試效果較好的MSCN、文獻[9]等對比，誤差小于 5% ，結(jié)果表明所提算法能夠用于視場霧濃度檢測。

研究采用的是基于圖像散點斜率先驗作為對比度測量，采用有參圖像統(tǒng)計估計，使檢測誤差較大，特別是在高濃度檢測中測試精度不能滿足要求，在惡劣氣象環(huán)境中（如雨、雪），檢測精度也會受到影響。另外本算法存在的另一個不足是需要同場景的有霧和無霧圖像才能解算，上述不足需要進一步研究。

參考文獻

[1］高濤，劉夢尼，陳婷，等.結(jié)合暗亮通道先驗的遠近景融合去霧算法[J].西安交通大學學報，2021，55（10）：78-86.Gao T，LiuMNChenTetal.Afarandnearscenefusiondefogging algorithmbasedonthe priorofdark-lightchanel[].

［2］俞任姣，李振華，來建成，等.雨霧共存天氣 1.55μm 波長激光動態(tài)傳輸特性[J].中國激光，2021，48（13）：153-160. Yu R J， Li Z H， Lai J C， et al. Dynamic propagation characteristics of 1.55-μm wavelength laser in rain and fog coexisting weather[J]. Chinese Journalof Lasers， 2021， 48（13）： 153-160. （in Chinese）

[3］王亞平，周裕豐，張寶華.基于去霧增強和張量恢復(fù)的紅外小目標檢測[J].紅外與激光工程，2022，51（4）：429-437. Wang YP， Zhou YF， Zhang BH.Infrared smaltarget detectionbased on dehazing enhancement and tensor recovery[]. Infrared and Laser Engineering， 2022， 51（4）： 429-437.（in Chinese）

[4]MitalA，Moorthy AK，Bovik ACNo-reference image qualityasessment inthe spatialdomain[]EEETransactionson Image Processing， 2012， 21（12）： 4695-4708.

［5］溫立民，巨永鋒，閆茂德.基于自然統(tǒng)計特征分布的交通圖像霧濃度檢測[J].電子學報，2017，45（8）：1888-1895. WenLM，JuFanMDnspectnoffgenityforaffcmagebasedonstbutnhaactersticsofatuatatstics[]. Acta Electronica Sinica， 2017， 45（8）： 1888-1895. （in Chinese）

[6]CuiGG，TangLY，LiuMQetal.Quantativeresponseofsubjectivevisualrecognitiontofogconcentrationattuation based on image standard deviation[J]. Optik， 2021，232： 166446.

［7］李春明，姜雨彤，宋海平，等.基于光學厚度代理模型的霧濃度估計及圖像去霧[J].兵工學報，2019，40（7）：1425-1433. LiCM，JiangYSong HP，etal.Researchonpticaldepthsurrogatemodel-basedmethodforestimating fogdensityand removing fog effect from images[J]. Acta Armamentarii， 2019， 40（7）： 1425-1433. （in Chinese）

[8]GuoH，Wang XC，LiHJ.Densityestimationoffog inimagebasedondarkchannelprior[J].Atmosphere，22，3（5）：710.

[9]Ling ZGGongJWFanGL，etal.Optimaltransmissonestimationviafogdensityperceptionforefficient single image defogging[J]. IEEE Transactions on Multimedia， 2018， 20（7）： 1699-1711.

[10]JiangYTSunCZhaoYetal.Foenityestiatiandimagefoggingasedonogateelingfortaet[]. IEEE Transactions on Image Processing， 2017， 26（7）： 3397-3409.

[11]DuminilAarelJBremndRSingleimageatosphericeilemalusingewprorsfotergenericity[]Atmh， 2021，12（6）： 772.

[12]LiuD，SunJF，GaoS，etal.Single-paameterestimationconstructionalgorithmforGm-APDladarimaging through fog[] Optics Communications， 2021， 482： 126558.

[13]Wang TZhaoLHuang PCetal.Hazeconcentratioadaptivenetwork frimagedehazing[].Neurocomputing，439： 75-85.

[14］賈之娟.激光在霧環(huán)境下的傳輸衰減特性研究[D].西安：西安理工大學，2021. Jia ZJ.tudyontransmissionateuationharacteristicsoflaserinfogenvironmentDXi'an：Xi'anUniversityofTechnology， 2021.（in Chinese）

[15]陳勇，帥鋒，樊強.基于自然統(tǒng)計特征分布的無參考圖像質(zhì)量評價[J].電子與信息學報，2016，38（7）：1645-1653. ChenY，ShuaiFFanQ.ANreferenceimagequalityasessmentbasedondistributioncharacteristicsofnaturalstatistics[]. Journal of Electronics amp; Information Technology， 2016， 38（7）： 1645-1653.（in Chinese）

[16]Narasimhan SG，Nayar SK.Contrastrestorationofweatherdegradedimages[]IEEETransactions on PaternAnalysisand Machine Intelligence， 2003， 25（6）： 713-724.

[17]HautiereNarelJPHalmaouiHetalnhancedfogdetectionandfree-space segmentationforcaravigation[achine Vision and Applications，2014，25（3）：667-679.

[18]付曉薇，代蕓，陳黎，等.基于局部的量子衍生醫(yī)學超聲圖像去斑[J].電子與信息學報，2015，37（3）：560-566. Fu XW，DaiY，Chen L，etal.Quantum-inspireddespeckling ofmedicalultrasound images basedonlocalentropy[]Joual of Electronics amp; Information Technology， 2015， 37（3）： 560-566.（in Chinese）

[19]WangCJZhangQY，SunYetal.Anewumericalcrectinmethodforgammaspectabasedonthesystemtrafoatin theory of random signals[J].Applied Radiation and Isotopes， 2021，172： 109671.

[20］楊愛萍，楊蘇輝，田鑫，等.基于統(tǒng)計特性和亮度估計的夜晚圖像去霧[J].天津大學學報，2018，51（3）：299-307. Yang AP，Yang SH，TianXetal.Nightimeimagefogremovalbasedonstatistical propertiesandintensityestimatin[]. Journal of Tianjin University （Science and Technology）， 2018， 51（3）： 299-307. （in Chinese）

（編輯侯湘）