基于自適應(yīng)霧濃度系數(shù)的暗通道先驗法能見度測量

陳鐘榮 張炎 張瑤

摘 ?要： 針對傳統(tǒng)暗通道先驗法處理圖像計算能見度時，透射率優(yōu)化計算復(fù)雜以及透射率估測準(zhǔn)確性差的問題，提出一種引導(dǎo)濾波的改進(jìn)方法，以灰度圖代替RGB圖作為引導(dǎo)圖，極大地提高了濾波效率。同時，重點分析霧濃度系數(shù)[ω]和圖像灰度均值差異對透射率估測準(zhǔn)確性的影響，總結(jié)出灰度均值的合理范圍，提出一種自適應(yīng)霧濃度系數(shù)的新思路。通過構(gòu)建BP神經(jīng)網(wǎng)絡(luò)，以灰度均值及目標(biāo)區(qū)域計算值作為輸入，霧濃度系數(shù)作為輸出，實現(xiàn)了能見度測量系統(tǒng)對于霧濃度系數(shù)的自適應(yīng)修正。將大量圖像處理的能見度結(jié)果與CJY散射式能見度儀測量結(jié)果作比較，二者基本一致，97.1%的相對誤差在20%以內(nèi)，滿足能見度測量要求。

關(guān)鍵詞： 霧濃度系數(shù); 能見度測量; 暗通道先驗法; 引導(dǎo)濾波; BP神經(jīng)網(wǎng)絡(luò); 灰度圖

中圖分類號： TN911.73?34; P427.2 ? ? ? ? ? ? ? ? ? ? 文獻(xiàn)標(biāo)識碼： A ? ? ? ? ? ? ? ?文章編號： 1004?373X（2019）09?0039?07

Visibility measurement using dark channel prior method based on adaptive fog concentration coefficient

CHEN Zhongrong1， 2， ZHANG Yan1， 2， ZHANG Yao1， 2

（1. Key Laboratory for Aerosol?Cloud?Precipitation of China Meteorological Administration， Nanjing 210044， China;

2. School of Atmospheric Physics， Nanjing University of Information Science and Technology， Nanjing 210044， China）

Abstract： The conventional dark channel prior method used to process the image and calculate the visibility has the problems of complicated transmittance optimization calculation and poor transmittance estimation accuracy. Therefore， an improved guided filtering method taking grey?scale map as the guiding picture to replace RGB image is proposed to greatly improve the filtering efficiency. The influence of the difference of fog concentration coefficient [ω] and image gray mean value on the accuracy of transmittance estimation is analyzed emphatically， and the appropriate range of gray mean value is summarized. A new idea of adaptive fog concentration coefficient is proposed. By constructing the BP neural network， the gray mean value and calculated value of the target area are taken as inputs， and the fog concentration coefficient is taken as the output to realize the adaptive correction of the fog concentration coefficient of the visibility measurement system. The visibility results processed by a large number of images are compared with the measured results of CJY scattering visibility meter. The comparison results show that the values of two results are basically the same， and the 97.1% relative error is within 20%， and the method can satisfy the requirements of visibility measurement.

Keywords： fog concentration coefficient; visibility measurement; dark channel prior method; guided filtering; BP neural network; grey?scale map

0 ?引 ?言

大氣能見度是氣象觀測的重要指標(biāo)，對輪渡、民航、高速公路等交通運(yùn)輸及市民日常生活都有重大影響[1]。目前，其測量途徑主要分為目測、儀器測量和基于圖像的能見度測量。隨著生活中的圖像采集設(shè)備越發(fā)常見，以數(shù)字圖像測量能見度的優(yōu)勢不斷突顯。相比能見度測量儀器的高昂價格和高架設(shè)要求，生活中的數(shù)字圖像不僅獲取途徑豐富，且成本也更加低廉。因此，以圖像測量能見度的技術(shù)有較好的經(jīng)濟(jì)優(yōu)勢和發(fā)展前景。

國際上，文獻(xiàn)[2]提出使用圖像的視覺特征測量能見度，但受制于技術(shù)水平并未推廣應(yīng)用;文獻(xiàn)[3]采用小波變換對圖像邊緣進(jìn)行檢測，進(jìn)而估算大氣能見度，但適用場景受限;文獻(xiàn)[4]提出基于目標(biāo)物的亮度特征點測量能見度，但需預(yù)設(shè)目標(biāo)物，成本高且繁瑣。

在國內(nèi)，文獻(xiàn)[5]利用雙亮度差法估測日間能見度，但該方法依賴人工標(biāo)記且對觀測條件要求高;文獻(xiàn)[6]使用自標(biāo)定、小波變換等方法，以道路監(jiān)控視頻檢測能見度，效果尚佳但測量距離較短，且不適用道路結(jié)冰的情況。文獻(xiàn)[7]提出暗通道先驗理論，并廣泛應(yīng)用于圖像去霧領(lǐng)域，而利用該方法可以得到圖像每個像素點的透射率，進(jìn)一步求得每點的消光系數(shù)后，即可獲得大氣能見度信息。

1 ?測量原理

1.1 ?能見度測量原理

氣象能見度定義為：標(biāo)準(zhǔn)視力眼睛觀察水平方向以天空為背景的黑體目標(biāo)物（視角[0.5°～5°]），能從背景上分辨出目標(biāo)物輪廓的最大水平距離。設(shè)目標(biāo)物亮度為[IT]，背景視亮度為[IB]，則視亮度對比（[C]）可以表示為：[C=IB-ITIB]。世界氣象組織采用氣象光學(xué)視程（MOR）來衡量大氣光學(xué)狀態(tài)，定義為白熾燈發(fā)出色溫為2 700 K的平行光柱，通過大氣，亮度減少到初始的0.05時的路徑長度。根據(jù)朗伯?布格定律：

式中：[F]是經(jīng)過路徑[L]后的光通量;[F0]為[L=0]時的光通量;[σ]為平均消光系數(shù)。根據(jù)透射率的概念：[T=FF0]，則[T=e-σL]，令[T=0.05]，則光學(xué)視程[1]為：

所以，能見度可以通過測量消光系數(shù)來計算。

1.2 ?基于暗通道先驗理論的能見度測量

在計算機(jī)圖形領(lǐng)域中，霧霾圖像的物理模型為[8]：

式中：[I（x）]是待去霧圖像;[J（x）]是無霧圖像;參數(shù)[A]是全球大氣光成分;[t（x）]為透射率。假定大氣光均勻分布，則透射率[t（x）]和大氣散射系數(shù)[α]滿足如下關(guān)系：

式中[d（x）]表示目標(biāo)物到攝像頭的距離。

在可見光和紅外波段，大氣對光的吸收作用可忽略不計，因此可將散射系數(shù)近似為消光系數(shù)，即[σ≈α]。將式（4）代入式（2）得到能見度與透射率的關(guān)系式：

故測量能見度的關(guān)鍵是通過圖像處理獲得準(zhǔn)確的透射率。

透射率計算中，對圖像分窗口處理，設(shè)窗口小鄰域中的透射率[t（x）]為常數(shù)[t（x）]，[A]在圖像中可以解出，對式（3）變換后兩邊求最小值得：

式中：[Jc]表示圖像[J]中RGB三色通道其中的一個;[Ω（x）]是以[x]為中心的局部區(qū)域;[y]表示濾波窗口中的像素點。

而暗通道先驗原理指在絕大多數(shù)非天空的局部區(qū)域內(nèi)，某些像素總會有至少一個顏色通道的灰度值很低。即[Jdark（x）=miny∈Ω（x）[mincJc（y）]=0]，代入式（6）就可估測透射率：

在現(xiàn)實生活中，盡管是晴天，空氣中依然會含有雜質(zhì)分子。所以引入一個霧濃度系數(shù)[9][ω]（[0<ω≤1]）。最終，由式（8）求出透射率后代入式（5）即得能見度。

1.3 ?透射率優(yōu)化

由于圖像一定區(qū)域塊內(nèi)的透射率[t（x）]不總是常量，所以上述方式獲得的透射率較粗糙，導(dǎo)致透射率圖存在塊狀效應(yīng)。文獻(xiàn)[9?10]先后采用軟摳圖法和引導(dǎo)濾波法精細(xì)化透射率圖，獲得了不錯的效果。本文同樣采用引導(dǎo)濾波，但在原有基礎(chǔ)上作出改進(jìn)，在保證濾波效果的前提下，進(jìn)一步提升了濾波效率。

1.3.1 ?引導(dǎo)濾波

引導(dǎo)濾波是邊緣保留的平滑算子，濾波后透射率[t（x）]與霧霾圖像[I（x）]滿足線性關(guān)系：

1.3.2 ?濾波改進(jìn)

引導(dǎo)濾波時，通常取原霧霾圖像為引導(dǎo)圖。但以RGB圖（三維矩陣）濾波時，運(yùn)算量遠(yuǎn)大于基于灰度圖（一維矩陣）對透射率處理時的情況。文獻(xiàn)[12]預(yù)處理壓縮圖片大小，通過減小數(shù)據(jù)量提高引導(dǎo)濾波的運(yùn)算速度。不過，壓縮圖片減少了有效的像素數(shù)量，必然會造成信息量的損失。所以，本文提出一種以灰度圖代替RGB圖進(jìn)行濾波的新思路，既保持像素數(shù)量，又提高了濾波效率。表1是分別采用灰度圖像和彩色RGB圖像作為引導(dǎo)圖，對不同天氣條件下觀測圖像進(jìn)行濾波實驗的對比結(jié)果。

表1 ?RGB圖濾波及灰度圖濾波的誤差和時間比較

如表1所示，在三種不同能見度下，當(dāng)基于灰度圖濾波時，時間大大縮短，且與基于RGB圖濾波的能見度估測結(jié)果相對誤差極?。ㄐ∮?%）。后對2018年4月28日—5月3日，每天7：00～18：00的整點觀測數(shù)據(jù)進(jìn)行統(tǒng)計分析，基于灰度圖濾波的結(jié)果較RGB濾波結(jié)果的平均相對偏差為0.434%，而濾波效率提升95.93%。圖1為采用兩種濾波方法測量能見度結(jié)果的最小二乘法擬合，擬合函數(shù)為[y=1.006 ?2x-12.761 ?9]，二者相關(guān)系數(shù)[R]接近于1。因此，以灰度圖進(jìn)行引導(dǎo)濾波切實可行，且極大提高了圖像處理效率。

圖1 ?兩種濾波結(jié)果的擬合曲線

2 ?暗通道法透射率估測影響因素分析

2.1 ?目標(biāo)區(qū)域選取

由于暗通道理論存在假設(shè)前提，適用范圍是非天空的絕大多數(shù)局部區(qū)域，因此，符合理論適用條件目標(biāo)區(qū)域的選取，對透射率的估測準(zhǔn)確性影響極大[9]。實驗中，目標(biāo)區(qū)域采取人工篩選的方式進(jìn)行確定。

本實驗地點位于南京信息工程大學(xué)氣象樓10樓，向東側(cè)進(jìn)行圖像采集。前期對觀測方向上的各目標(biāo)物進(jìn)行詳細(xì)分析，發(fā)現(xiàn)如圖2所示的4個區(qū)域暗通道的灰度值較低，滿足暗通道理論條件。依據(jù)文獻(xiàn)[9]提出的目標(biāo)物選取標(biāo)準(zhǔn)對4個區(qū)域進(jìn)行篩選。其中，1區(qū)域距離遠(yuǎn)，在低能見度天氣條件下受霧霾影響，與天空的區(qū)分度小，測量效果較差。4區(qū)域經(jīng)測發(fā)現(xiàn)其觀測角度不滿足視角0.5°～5°的范圍標(biāo)準(zhǔn)[9]，對水平能見度的代表性不足。3區(qū)域為305 m處的綠色植被，大風(fēng)天氣時，樹冠晃動會帶來目標(biāo)區(qū)域不穩(wěn)定，導(dǎo)致能見度測量誤差偏大。綜上考慮，實驗中選取545 m外固定建筑的2區(qū)域作為研究對象進(jìn)行能見度的測量。

圖2 ?觀測場景

2.2 ?霧濃度系數(shù)

霧濃度系數(shù)[ω]值固定為0.95，但文獻(xiàn)[13]研究發(fā)現(xiàn)，霧濃度系數(shù)的變化對透射率計算影響很大，不同霧化程度的圖像計算透射率時[ω]取值不可絕對固定。同樣，在測量能見度中，若是不考慮天氣條件變化，而僅選取固定的[ω]會出現(xiàn)較大誤差。

式（13）為式（12）的導(dǎo)函數(shù)。由式（13）可發(fā)現(xiàn)，當(dāng)能見度一定時，處理圖像獲得的計算值[miny∈Ω（x）mincIc（y）Ac]也將固定不變，式（13）可看作關(guān)于自變量[ω]和因變量[MOR]的函數(shù)關(guān)系，則能見度與霧濃度系數(shù)滿足對數(shù)倒數(shù)的關(guān)系，因此[ω]的小幅變化將引起能見度的大幅計算誤差。

以真實能見度9 035 m和3 834 m的圖像實驗結(jié)果為例，其誤差分布如圖3所示。結(jié)果表明，[ω]變化確實對能見度計算結(jié)果存在很大影響。而霧濃度系數(shù)主要受天氣條件變化的影響，在晴朗天氣下能見度較好時，[ω]取值偏小;而當(dāng)空氣中顆粒物濃度大，能見度較差時，[ω]的取值應(yīng)相應(yīng)擴(kuò)大。

2.3 ?圖像質(zhì)量對[ω]的選取影響

2.3.1 ?圖像灰度均值對[ω]的影響

實際圖像采樣中發(fā)現(xiàn)，即使是同一天氣條件下，由于不同觀測時間下環(huán)境光變化以及硬件設(shè)備響應(yīng)速率等多方面噪聲因素的綜合影響，會造成圖像質(zhì)量差異，這種差異主要表現(xiàn)為圖像視覺上的亮度差異（以灰度均值來衡量）。且若霧濃度系數(shù)[ω]僅按天氣條件選取確定，該差異將引起較大的能見度估計誤差。

圖3 ?霧濃度系數(shù)變化的能見度誤差分布

如圖4所示，圖4a）～圖4c）是在1 min內(nèi)連續(xù)拍攝的真實能見度為7 133 m的圖像，通過手動調(diào)節(jié)相機(jī)的曝光率模擬實際由環(huán)境光等噪聲因素引起的圖像亮度差異，得到三者的灰度均值分別為161，128和82。圖4d）～圖4f）為圖4a）～圖4c）對應(yīng)的灰度分布直方圖。

觀察發(fā)現(xiàn)，當(dāng)圖像的灰度均值增大時，整張圖像的灰度分布會向高一側(cè)拉伸，即各像素的灰度值都將升高，導(dǎo)致[miny∈Ω（x）mincIc（y）Ac]增大，但此時能見度（透射率）是一定的，由式（9）可知，霧濃度系數(shù)[ω]需要調(diào)整才能準(zhǔn)確修正透射率。若控制[ω]為0.3～0.95，以0.05間隔變化，按三種情況分別計算能見度的結(jié)果繪制曲線如圖5所示?？梢园l(fā)現(xiàn)，為完成透射率修正，得到準(zhǔn)確的能見度，隨著圖像灰度均值增大，選取的霧濃度系數(shù)需減小。

圖5 ?不同灰度均值圖像能見度估測結(jié)果

由上述結(jié)果可知，計算透射率時，同一天氣條件下，不同灰度均值的圖像對應(yīng)的[ω]不同，而灰度均值范圍一定時，不同天氣條件的霧濃度系數(shù)[ω]對應(yīng)的規(guī)律保持穩(wěn)定。且天氣條件一定時，隨灰度均值遞減，[ω]取值將遞增，二者保持穩(wěn)定的負(fù)相關(guān)性。

2.3.2 ?灰度均值合理范圍

實驗發(fā)現(xiàn)，當(dāng)灰度均值過小時，若能見度小于3 km，[ω]在定義域內(nèi)無法完成對透射率的修正;當(dāng)灰度均值過大時，若能見度大于11 km，能見度的估測誤差波動會增大。為保證能見度測量的動態(tài)范圍較為理想，圖像灰度均值存在一定合理范圍。

引起該現(xiàn)象的原因是灰度均值變化會導(dǎo)致目標(biāo)區(qū)域計算值[miny∈Ω（x）mincIc（y）Ac]變化。當(dāng)能見度小時，對應(yīng)實際透射率應(yīng)小，由式（8）可知，[ωminy∈Ω（x）mincIc（y）Ac]應(yīng)該較大，而此時若圖像整體灰度均值過小，會導(dǎo)致[miny∈Ω（x）mincIc（y）Ac]自身過小。[ω]取值范圍為0～1，無法彌補(bǔ)當(dāng)前修正需要，使得最終結(jié)果表現(xiàn)為[ω]對透射率無法完成修正。

同理，當(dāng)能見度大時，對應(yīng)實際透射率大，則[ωminy∈Ω（x）mincIc（y）Ac]應(yīng)較小，而此時若圖像整體灰度均值過大，會導(dǎo)致[miny∈Ω（x）mincIc（y）Ac]自身過大，使得[ω]單位變化引起的透射率變化量增大。而由式（8）的導(dǎo)數(shù)是非線性增函數(shù)（值域大于0）可知，能見度關(guān)于透射率的函數(shù)遞增速度在自變量取大值時變快。所以，當(dāng)[miny∈Ω（x）mincIc（y）Ac]過大時，[ω]單位變化引起的能見度變化量更大，最終結(jié)果表現(xiàn)為[ω]修正的能見度誤差波動劇烈。

綜上，不同灰度均值的圖像對于能見度測量的敏感響應(yīng)范圍是有限的，理想灰度均值范圍應(yīng)使[ω]既可完成全范圍的透射率修正，又不至于讓能見度對[ω]的變化過于敏感。因此，經(jīng)實驗統(tǒng)計發(fā)現(xiàn)，當(dāng)圖像灰度均值處于80～139時，其對于能見度的估測響應(yīng)效果較好，是較合理的灰度范圍。

3 ?自適應(yīng)[ω]的能見度測量系統(tǒng)

鑒于能見度測量中圖像灰度均值直接影響[ω]的選取，而不同天氣條件下[ω]的取值又決定了所得透射率的準(zhǔn)確性。因此，本文提出一種基于BP神經(jīng)網(wǎng)絡(luò)的自適應(yīng)[ω]的能見度測量系統(tǒng)。

3.1 ?能見度測量系統(tǒng)

能見度測量系統(tǒng)基本流程如圖6所示。首先計算圖像的灰度均值及暗通道圖像，然后取暗通道圖中最亮的前0.1%像素平均值為[A]，接著計算目標(biāo)區(qū)域的暗通道與全球大氣光成分[A]的比值，以灰度均值和計算值作為輸入，經(jīng)BP神經(jīng)網(wǎng)絡(luò)預(yù)測輸出霧濃度系數(shù)，完成透射率求解。最后通過改進(jìn)的引導(dǎo)濾波優(yōu)化透射率，即可獲得能見度信息。

3.2 ?BP神經(jīng)網(wǎng)絡(luò)

BP神經(jīng)網(wǎng)絡(luò)算法是一種有監(jiān)督的學(xué)習(xí)算法，即在輸入訓(xùn)練樣本后，運(yùn)用反向傳播算法對網(wǎng)絡(luò)的權(quán)值和偏差不斷調(diào)整訓(xùn)練，使得輸出結(jié)果逼近期望輸出。當(dāng)網(wǎng)絡(luò)輸出層的誤差平方達(dá)到或小于指定誤差時完成訓(xùn)練，保存網(wǎng)絡(luò)的權(quán)值和偏差[14]。

圖6 ?能見度測量流程圖

BP神經(jīng)網(wǎng)絡(luò)設(shè)計方案和訓(xùn)練樣本的選擇都會影響最終結(jié)果。訓(xùn)練樣本過少不足以反映輸入量總體信息，過多則會導(dǎo)致網(wǎng)絡(luò)規(guī)模過大及映射關(guān)系復(fù)雜化。隱層層數(shù)由映射函數(shù)的連續(xù)性決定，連續(xù)函數(shù)采用單隱層前饋網(wǎng)絡(luò)，非連續(xù)函數(shù)用多隱層。隱層節(jié)點數(shù)要求滿足最優(yōu)化前提，具體根據(jù)經(jīng)驗公式推導(dǎo)和嘗試得出[m=n+l+a]，其中，[m]為隱層節(jié)點個數(shù)，[n]為輸入節(jié)點個數(shù)，[l]為輸出節(jié)點個數(shù)，[a]為可調(diào)用常數(shù)（1～10）。

3.3 ?基于BP神經(jīng)網(wǎng)絡(luò)算法的[ω]自適應(yīng)選取

基于以上特性，本文選用經(jīng)典三層BP神經(jīng)網(wǎng)絡(luò)結(jié)構(gòu)，輸入層設(shè)置兩個神經(jīng)元，分別為圖像灰度均值和目標(biāo)區(qū)域計算值[miny∈Ω（x）mincIc（y）Ac]，輸出層設(shè)置一個神經(jīng)元，為對應(yīng)灰度均值及天氣條件下的霧濃度系數(shù)[ω]。隱含層使用Tansig函數(shù)，輸出層使用Purelin函數(shù)。根據(jù)經(jīng)驗公式，確定隱含層神經(jīng)元節(jié)點的個數(shù)為10。

使用2018年4月18日—5月6日能見度整點觀測實驗所得的715組數(shù)據(jù)作為訓(xùn)練樣本，該訓(xùn)練集包含不同天氣條件下能見度范圍為500～16 000 m的情況，具備普遍代表性。圖7a）為訓(xùn)練完成網(wǎng)絡(luò)的誤差收斂情況，訓(xùn)練誤差為3.24%。圖7b）是網(wǎng)絡(luò)輸出結(jié)果與期望值的擬合結(jié)果，擬合曲線為線性函數(shù)，相關(guān)系數(shù)[R]為0.941 4，接近1，結(jié)果表明該神經(jīng)網(wǎng)絡(luò)性能較好，滿足實驗要求。

圖7 ?神經(jīng)網(wǎng)絡(luò)訓(xùn)練結(jié)果

4 ?實驗結(jié)果與分析

為進(jìn)一步檢驗自適應(yīng)霧濃度系數(shù)能見度測量方法的效果。本文自2018年5月7日—5月16日的10天內(nèi)，對圖2場景在不同天氣條件下，不同時間段進(jìn)行能見度觀測實驗，每日11：00，12：00，14：00，15：00，16：00，17：00和18：00分7個整點觀測，共計210組數(shù)據(jù)。實驗設(shè)備為索尼WW328261型數(shù)碼相機(jī)，搭載FE28?70 mm變焦鏡頭，采集圖像大小為4 000×6 000 pixel。經(jīng)前期實驗得知，不同時間及天氣條件下，當(dāng)相機(jī)自動感光值一定時，所得圖像的灰度均值范圍穩(wěn)定不變。因此，拍攝時采用相機(jī)的手動曝光M模式，依據(jù)其自動感光指示值，在每個整點分別采集感光值為-0.3，0，+0.3的3張不同亮度圖像（分別代表欠曝光、正常曝光、過度曝光），模擬實際中由于各種環(huán)境噪聲因素導(dǎo)致的灰度均值變化的情況。這三種圖像灰度均值均滿足80～139的合理范圍。對三者分別計算能見度，將計算值與安裝在南京信息工程大學(xué)綜合觀測培訓(xùn)實習(xí)基地的CJY散射式能見度儀測量值進(jìn)行對比。

圖8a）～8c）分別表示高灰度均值圖像（過度曝光）、中灰度均值圖像（正常曝光）和低灰度均值圖像（欠曝光）的能見度測量結(jié)果與能見度儀測量結(jié)果的對比情況。可以看出，4～16 km能見度范圍內(nèi)，不同天氣條件下各灰度均值圖像的能見度計算結(jié)果和實測數(shù)據(jù)基本一致。對三種圖像的估測結(jié)果和實測數(shù)據(jù)進(jìn)行最小二乘擬合后，擬合函數(shù)為[y=0.81x+1 ?664]，相關(guān)系數(shù)[R]為0.88，接近于1。其中，僅2.9%的估測相對誤差超出20%，但最大誤差不超過22.6%，而97.1%的能見度估測相對誤差落在20%以內(nèi)。因此，各灰度均值圖像的計算結(jié)果與真實能見度基本符合，BP神經(jīng)網(wǎng)絡(luò)能夠完成不同天氣條件和圖像灰度均值下對霧濃度系數(shù)[ω]的自動修正，獲得準(zhǔn)確的能見度。據(jù)世界氣象組織規(guī)定，在探測范圍內(nèi)，能見度儀的相對誤差小于等于±20%即達(dá)標(biāo)，允許投入使用。而上述利用神經(jīng)網(wǎng)絡(luò)自適應(yīng)修正霧濃度系數(shù)的暗通道先驗?zāi)芤姸葴y量方法達(dá)到誤差標(biāo)準(zhǔn)，滿足能見度的測量要求[9]。

圖8 ?不同曝光率圖像能見度測量結(jié)果

5 ?結(jié) ?論

基于暗通道先驗原理測量能見度時，本文提出一種引導(dǎo)濾波改進(jìn)思路和自適應(yīng)[ω]的能見度測量方法。引導(dǎo)濾波時，以灰度圖像作為引導(dǎo)圖像，在保證圖像透射率優(yōu)化效果良好的前提下，相較RGB圖像引導(dǎo)方式效率更高。測量能見度時，不同天氣條件下霧濃度系數(shù)[ω]的選擇對透射率準(zhǔn)確性影響很大，不可取定值。且在同等天氣條件下，由于環(huán)境噪聲等因素引起的圖像灰度均值變化也對霧濃度系數(shù)選取存在影響，霧濃度系數(shù)的確定需同時考慮天氣條件和噪聲影響。因此，本文采用三層經(jīng)典結(jié)構(gòu)BP神經(jīng)網(wǎng)絡(luò)算法，實現(xiàn)了霧濃度系數(shù)的自適應(yīng)修正，獲得了準(zhǔn)確的能見計算結(jié)果。計算結(jié)果較CJY散射能見度儀測量值，97.1%的計算結(jié)果相對誤差在20%以內(nèi)，滿足能見度的測量要求，僅有2.9%的數(shù)據(jù)誤差超出20%，多出現(xiàn)在極高能見度時，原因可能是前期訓(xùn)練樣本中高能見度的樣本比重不足，影響網(wǎng)絡(luò)性能。綜上，自適應(yīng)霧濃度系數(shù)的暗通道能見度測量方法操作簡單，實施成本小，運(yùn)用范圍廣泛，有較好的使用前景。當(dāng)然，本文研究的目標(biāo)區(qū)域是手動選擇，后期的研究重點應(yīng)關(guān)注目標(biāo)區(qū)域的自動選擇，并擴(kuò)大觀測周期以增強(qiáng)訓(xùn)練樣本的代表性來進(jìn)一步優(yōu)化神經(jīng)網(wǎng)絡(luò)性能。

參考文獻(xiàn)

[1] 苗苗.視頻能見度檢測算法綜述[J].現(xiàn)代電子技術(shù)，2012，35（15）：72?75.

MIAO Miao. Survey on video?based visibility detection method [J]. Modern electronics technique， 2012， 35（15）： 72?75.

[2] STEFFENS C. Measurement of visibility by photographic photometry [J]. Industrial engineering chemistry， 1949， 41（11）： 2396?2399.

[3] BUSH C， DEDES C. Wavelet transform for analyzing fog visibility [J]. IEEE intelligent systems， 1998， 13（6）： 66?71.

[4] KWON T M. An automatic visibility measurement system based on video cameras [R]. Minnesota： Minnesota Department of Transportation Duluth Office and Office of Research Administration， 1998.

[5] 呂偉濤，陶善昌，劉亦風(fēng)，等.基于數(shù)字?jǐn)z像技術(shù)測量氣象能見度：雙亮度差方法和試驗研究[J].大氣科學(xué)，2004，28（4）：559?570.

L? Weitao， TAO Shanchang， LIU Yifeng， et al. Measuring meteorological visibility based on digital photography： dual differential luminance method and experimental study [J]. Chinese journal of atmospheric science， 2004， 28（4）： 559?570.

[6] 陳釗正，陳啟美.基于攝像機(jī)自標(biāo)定的視頻對比度能見度檢測算法與實現(xiàn)[J].電子與信息學(xué)報，2010，32（12）：2907?2912.

CHEN Zhaozheng， CHEN Qimei. Video contrast visibility detection algorithm and its implementation based on camera self calibration [J]. Journal of electronics & information technology， 2010， 32（12）： 2907?2912.

[7] HE K， SUN J， TANG X. Single image haze removal using dark channel prior [J]. IEEE transaction on pattern analysis & machine intelligence， 2011， 34（12）： 2341?2353.

[8] 郭艷穎，楊國慶.基于暗原色和加權(quán)形態(tài)學(xué)濾波的圖像去霧算法[J].現(xiàn)代電子技術(shù)，2012，35（20）：71?73.

GUO Yanying， YANG Guoqing. Auto image haze removal algorithm based on dark channel prior and weighted morphological filtering [J]. Modern electronics technique， 2012， 35（20）： 71?73.

[9] 趙瑞，胡遼林，趙江濤.基于暗原色先驗理論的大氣能見度測量[J].激光與光電子學(xué)進(jìn)展，2016，53（6）：78?85.

ZHAO Rui， HU Liaolin， ZHAO Jiangtao. Atmospheric visibility measurement based on dark channel prior [J]. Laser & optoelectronics progress， 2016， 53（6）： 78?85.