肖堯，易凱

（四川警察學(xué)院 四川瀘州 646000）

現(xiàn)代科學(xué)技術(shù)的發(fā)展，帶動(dòng)了公安部門(mén)查案辦案技術(shù)和方式的改變。隨著公安機(jī)關(guān)視頻監(jiān)控系統(tǒng)以及民用、商用監(jiān)控設(shè)備的完善和普及，圖像為公安偵破各類(lèi)案件的過(guò)程中帶來(lái)了極大便利。由圖像記錄的重要信息常被收集作為案件證據(jù)，還原案件事實(shí)的重要突破口。圖像中復(fù)原的運(yùn)動(dòng)模糊車(chē)牌圖像能夠使公安機(jī)關(guān)快速鎖定嫌疑車(chē)輛，快速分析嫌疑人行動(dòng)軌跡，快速偵破案件。

當(dāng)前，復(fù)原運(yùn)動(dòng)模糊車(chē)牌圖像的技術(shù)有很多，例如：朱一帆運(yùn)用逆濾波算法和維納濾波算法復(fù)原運(yùn)動(dòng)模糊圖像，逆濾波算法需要在不存在噪聲干擾情況下才有較好的復(fù)原效果，而維納濾波算法能夠有效地抑制噪聲，在有先驗(yàn)信息的情況下也有很好的復(fù)原效果[1]；Richardson 和Lucy 二人基于貝葉斯理論提出的Richardson-Lucy（簡(jiǎn)稱(chēng)RL）算法[2][3]，該算法利用最大似然估計(jì)對(duì)模糊圖像進(jìn)行迭代求解；也有最小二乘方算法[4]和劉湘武改進(jìn)最小二乘方算法復(fù)原運(yùn)動(dòng)模糊圖像[5]；周志強(qiáng)等人提出的使用Butterworth 帶通濾波器對(duì)模糊圖像頻譜圖進(jìn)行預(yù)處理的方法，也能夠較好地求出重建模糊核參數(shù)[6]；陳娟提出的基于MAP 的多尺度復(fù)原運(yùn)動(dòng)模糊圖像算法[7]能夠提高運(yùn)動(dòng)去模糊的效果，但振鈴效應(yīng)仍然無(wú)法完全消除；李桐提出的基于圖像中奇異值的復(fù)原模糊圖像算法[8]能夠避免的以往算法中所帶來(lái)的振鈴效應(yīng)，去模糊的效果也相對(duì)較好，但是算法任過(guò)于復(fù)雜。

目前公安機(jī)關(guān)復(fù)原運(yùn)動(dòng)模糊車(chē)牌主要使用一些專(zhuān)業(yè)圖像處理軟件，但這些軟件大多數(shù)都存在使用費(fèi)用高、專(zhuān)業(yè)性強(qiáng)，普及率等問(wèn)題，不能滿(mǎn)足公安一線的實(shí)際需求[9]。所以本文提出將人工估計(jì)運(yùn)動(dòng)模糊核及點(diǎn)擴(kuò)散函數(shù)PSF[10]與Matlab 軟件中基于維納濾波、正則化濾波、Richardson-Lucy、盲反褶積的四種算法相結(jié)合的方法來(lái)解決公安領(lǐng)域運(yùn)動(dòng)模糊車(chē)牌的復(fù)原問(wèn)題。

一、運(yùn)動(dòng)模糊車(chē)牌復(fù)原的原理

（一）運(yùn)動(dòng)模糊車(chē)牌的物理形成原因

運(yùn)動(dòng)模糊車(chē)牌形成的原因有二：一種是成像設(shè)備每秒幀數(shù)的原因；另一種是監(jiān)控系統(tǒng)在拍攝時(shí)車(chē)輛運(yùn)動(dòng)速度的原因。車(chē)輛在道路上短時(shí)間內(nèi)的運(yùn)動(dòng)可以視為直線運(yùn)動(dòng)，在每幀的曝光時(shí)間內(nèi)，車(chē)輛因快速運(yùn)動(dòng)在成像設(shè)備的傳感器上留下一段位移圖像，從而形成了運(yùn)動(dòng)模糊的車(chē)牌[11]。

圖1 是成像設(shè)備形成運(yùn)動(dòng)模糊的簡(jiǎn)單原理圖。由圖1 我們知道，當(dāng)拍被攝物體在成像設(shè)備曝光時(shí)間內(nèi)有多次位移的情況下，成像設(shè)備的傳感器上也會(huì)相應(yīng)的留下多重影像，最終導(dǎo)致成像設(shè)備輸出圖像是產(chǎn)生運(yùn)動(dòng)模糊。

（二）運(yùn)動(dòng)模糊車(chē)牌的數(shù)字退化過(guò)程

將運(yùn)動(dòng)模糊車(chē)牌的物理形成過(guò)程轉(zhuǎn)化為數(shù)字退化過(guò)程（如圖2所示）。

圖1 成像設(shè)備形成運(yùn)動(dòng)模糊的簡(jiǎn)單原理

圖2 物理形成過(guò)程轉(zhuǎn)化為數(shù)字退化過(guò)程

在數(shù)字退化過(guò)程中監(jiān)控設(shè)備受到天氣、溫度、磁場(chǎng)等因素干擾，使數(shù)字圖像存在噪聲[11]。復(fù)原運(yùn)動(dòng)模糊車(chē)牌是指根據(jù)運(yùn)動(dòng)模糊核以及加入車(chē)牌圖像中的噪聲，盡可能的還原原始車(chē)牌中的有用信息。

其中：f（x，y）為理想的無(wú)運(yùn)動(dòng)模糊原始車(chē)牌，b（x，y）為運(yùn)動(dòng)模糊形成的模糊核，h（x，y）為無(wú)噪聲的運(yùn)動(dòng)模糊車(chē)牌，N（x，y）為加入各種因素所引起的噪聲，g（x，y）為加入噪聲待復(fù)原的車(chē)牌Matlab中運(yùn)動(dòng)模糊圖像復(fù)原算法。

本節(jié)主要介紹本課題用于研究運(yùn)動(dòng)模糊車(chē)牌復(fù)原的四種算法，具體如下：

基于維納濾波的算法：J=deconvwnr（I，psf，nsr），其中運(yùn)動(dòng)模糊圖像I 利用維納濾波算法復(fù)原為去運(yùn)動(dòng)模糊圖像J，psf是點(diǎn)擴(kuò)散函數(shù)及模糊核，nsr是加性噪聲與信號(hào)功率的比值。

基于正則化濾波的算法：J=deconvreg（I，psf，np），其中模糊圖像I利用正則化濾波算法復(fù)原為去模糊圖像J；psf是指點(diǎn)擴(kuò)散函數(shù)及模糊核；np是指定噪聲的功率。

基于Richardson-Lucy 的算法：J=deconvlucy（I，psf，iter），其中模糊圖像I利用R-L 算法復(fù)原為去模糊圖像J；psf是指點(diǎn)擴(kuò)散函數(shù)及模糊核；iter是指迭代次數(shù)。

基于盲反褶積的算法：[J，psfr]=deconvblind（I，psf，iter），其中模糊圖像I 利用盲反褶積的算法復(fù)原為去模糊圖像J；psf是指點(diǎn)擴(kuò)散函數(shù)及模糊核；iter是指迭代次數(shù)。

二、運(yùn)動(dòng)模糊車(chē)牌復(fù)原算法仿真實(shí)驗(yàn)

（一）不加噪復(fù)原運(yùn)動(dòng)模糊車(chē)牌算法比較

本組實(shí)驗(yàn)中，為原始車(chē)牌（圖3）添加長(zhǎng)度60 像素，方向角45 度的點(diǎn)擴(kuò)散函數(shù)在未添加噪聲干擾的情況下模擬無(wú)加噪運(yùn)動(dòng)模糊車(chē)牌（圖4），將已知的點(diǎn)擴(kuò)散函數(shù)（L=60，α=45°）與Matlab 中基于維納濾波（圖5）、正則化濾波（圖6）、Richardson-Lucy（圖7）、盲反褶積（圖8）的四種算法相結(jié)合復(fù)原運(yùn)動(dòng)模糊車(chē)牌。

圖3 原始車(chē)牌

圖4 無(wú)加噪運(yùn)動(dòng)模糊車(chē)牌

圖6 正則化濾波復(fù)原

圖7 R-L算法復(fù)原

圖8 盲反褶積算法復(fù)原

從實(shí)驗(yàn)結(jié)果可以看出基于維納濾波的算法（圖5），基于正則化濾波的算法（圖6），基于Rich?ardson-Lucy 的算法（圖7）能夠在已知點(diǎn)擴(kuò)散函數(shù)的情況下較好的復(fù)原運(yùn)動(dòng)模糊車(chē)牌，基于盲反褶積的算法（圖8）雖然也能復(fù)原圖像但是其效果不如前三種算法。

加噪復(fù)原運(yùn)動(dòng)模糊車(chē)牌算法比較：本組實(shí)驗(yàn)研究噪聲對(duì)上述四類(lèi)算法的影響，噪聲采用高斯噪聲，噪聲添加值為：0.01。

圖9 維納濾波復(fù)原

圖10 正則化濾波復(fù)原

圖11 R-L算法復(fù)原

圖12 盲反褶積算法復(fù)原

由實(shí)驗(yàn)結(jié)果可以看出，在加入噪聲后，基于維納濾波的算法（圖9）和基于正則化濾波的算法（圖10）受噪聲的影響較強(qiáng)，而后兩種算法（圖11、12）受噪聲的影響較弱。下面對(duì)基于維納濾波的算法和基于正則化濾波的算法添加預(yù)估噪聲值0.01。如圖13、14 所示，兩種算法都能較好的通過(guò)抑制噪聲復(fù)原運(yùn)動(dòng)模糊車(chē)牌，并且可以通過(guò)人員無(wú)阻力的識(shí)別車(chē)牌信息，較圖9、10 對(duì)應(yīng)的兩種算法的復(fù)原效果圖有很大提升。

圖13 維納濾波復(fù)原

圖14 正則化濾波復(fù)原

三、運(yùn)功模糊車(chē)牌復(fù)原算法實(shí)驗(yàn)

（一）人工估計(jì)運(yùn)動(dòng)模糊PSF的數(shù)學(xué)模型

數(shù)字圖像是由一個(gè)一個(gè)像素點(diǎn)緊密排列形成，運(yùn)動(dòng)模糊是由于像素點(diǎn)發(fā)生位移而形成。因?yàn)楫?dāng)車(chē)輛快速?gòu)谋O(jiān)控設(shè)備下通過(guò)時(shí)，可以將車(chē)輛看作在做直線運(yùn)動(dòng)，所以運(yùn)動(dòng)模糊車(chē)牌的運(yùn)動(dòng)軌跡也可以視為直線。因此估計(jì)運(yùn)動(dòng)模糊長(zhǎng)度的數(shù)學(xué)模型就可以理解為以1像素為長(zhǎng)度單位的直角三角形求解斜邊。

創(chuàng)建一個(gè)大小為50 像素*50 像素的白色背景，輸入字母A，通過(guò)Matlab 模擬運(yùn)動(dòng)模糊，其運(yùn)動(dòng)模糊長(zhǎng)度為8像素，運(yùn)動(dòng)模糊方向?yàn)?6度。如圖15所示：

圖15 運(yùn)動(dòng)模糊像素位移

圖11最右側(cè)的圖像是運(yùn)動(dòng)模糊后圖像的局部放大圖，放大軟件為Windows電腦自帶的繪圖軟件，使用此軟件的原因：一是因?yàn)槊颗_(tái)Windows 電腦都自帶，便于以后公安一線直接使用；二是因?yàn)橛么塑浖糯髨D像可以放到的單位像素級(jí)別，便于運(yùn)動(dòng)模糊車(chē)牌的運(yùn)動(dòng)模糊長(zhǎng)度（l）以及運(yùn)動(dòng)模糊角度（α）的計(jì)算。

通過(guò)使用Windows自帶繪圖軟件進(jìn)行估計(jì)運(yùn)動(dòng)模糊長(zhǎng)度（l）及角度（α），其具體的計(jì)算公式為：

將人工估計(jì)的運(yùn)動(dòng)模糊長(zhǎng)度l 和角度α 代入基于維納濾波、正則化濾波、Richardson-Lucy、盲反褶積的四種算法中，根據(jù)結(jié)果（如圖16 所示）可以看出，人工估計(jì)的運(yùn)動(dòng)模糊長(zhǎng)度及角度對(duì)于復(fù)原運(yùn)動(dòng)模糊圖像有一定的效果，此種方法效率高，簡(jiǎn)單易學(xué)能夠廣泛適用于公安一線，從而避免了繁雜的自動(dòng)識(shí)別運(yùn)動(dòng)模糊長(zhǎng)度及角度的算法，降低了算法使用的難度。

圖16 像素位移復(fù)原

復(fù)原運(yùn)動(dòng)模糊車(chē)牌的方法步驟：

第一步：使用Windows繪圖軟件打開(kāi)待復(fù)原圖圖像，勾選標(biāo)尺和網(wǎng)格線。

第二步：放大圖像至最大及每方格大小為1 像素。使用形狀工具選擇直角三角形并調(diào)整粗細(xì)為1像素。

第三部：根據(jù)待復(fù)原圖像的模糊長(zhǎng)度和方向調(diào)整直角三角形使三角形斜邊方向與模糊方向平行，與模糊長(zhǎng)度大致一致。

第五步：將求出的模糊長(zhǎng)度l 和模糊角度α 代入Matlab 軟件中的四種經(jīng)典的運(yùn)動(dòng)模糊復(fù)原算法，復(fù)原運(yùn)動(dòng)模糊圖像。

估計(jì)待原車(chē)牌的運(yùn)動(dòng)模糊核待復(fù)原的三幅運(yùn)動(dòng)模糊車(chē)牌（如圖17）。

圖17 待復(fù)原車(chē)牌（從左至右：a、b、c）

將三幅待運(yùn)動(dòng)模糊復(fù)原的車(chē)牌使用Windows 自帶的繪圖軟件打開(kāi)，并估算出其運(yùn)動(dòng)模糊長(zhǎng)度以及運(yùn)動(dòng)模糊角度（如圖18、19、20）。

圖18 待復(fù)原車(chē)牌a運(yùn)動(dòng)模糊長(zhǎng)度和角度

圖19 待復(fù)原車(chē)牌b運(yùn)動(dòng)模糊長(zhǎng)度和角度

圖20 待復(fù)原車(chē)牌c運(yùn)動(dòng)模糊長(zhǎng)度和角度

圖18 中估計(jì)的角度為運(yùn)動(dòng)模糊的相對(duì)角度，其相對(duì)的是垂直參考線，而在Matlab 的運(yùn)算中角度參考線為水平參考線，所以最左側(cè)圖中的真實(shí)的運(yùn)動(dòng)模糊估計(jì)角度應(yīng)為106°。人工估計(jì)運(yùn)動(dòng)模糊噪聲在視頻監(jiān)控系統(tǒng)采集圖像，采集圖像往往會(huì)受到外界諸多因素形成噪聲點(diǎn)，形成噪聲的原因有很多，每幅圖像噪聲點(diǎn)的數(shù)量也各不相同，這導(dǎo)致了人工無(wú)法估計(jì)噪聲值，在復(fù)原運(yùn)動(dòng)模糊圖像時(shí)，經(jīng)驗(yàn)性的輸入預(yù)估噪聲值如0.001，0.005，0.01等等，具體值根據(jù)算法復(fù)原運(yùn)動(dòng)模糊車(chē)牌具體效果而定。

基于人工估計(jì)PSF復(fù)原運(yùn)動(dòng)模糊車(chē)牌將人工估計(jì)的運(yùn)動(dòng)模糊長(zhǎng)度l1、l2、l3和角度α1、α2、α3代入基于維納濾波、正則化濾波、Richardson-Lucy、盲反褶積的四種算法中。如圖21、22、23 所示，實(shí)驗(yàn)證明通過(guò)人工估計(jì)點(diǎn)擴(kuò)散函數(shù)（PSF）是行之有效的，這樣以來(lái)，避免了大量且復(fù)雜的自動(dòng)估計(jì)運(yùn)動(dòng)模糊模糊核的算法，通過(guò)簡(jiǎn)單易學(xué)的數(shù)學(xué)知識(shí)即可完成點(diǎn)擴(kuò)散函數(shù)的估計(jì)。

圖21 真實(shí)運(yùn)動(dòng)模糊車(chē)牌復(fù)原結(jié)果之一

圖22 真實(shí)運(yùn)動(dòng)模糊車(chē)牌復(fù)原結(jié)果之二

圖23 真實(shí)運(yùn)動(dòng)模糊車(chē)牌復(fù)原結(jié)果之三

基于公安專(zhuān)業(yè)圖像處理軟件復(fù)原運(yùn)動(dòng)模糊車(chē)牌將圖17 待復(fù)原的三幅車(chē)牌圖像使用專(zhuān)業(yè)圖像處理軟件Amped FIVE 打開(kāi)，使用去運(yùn)動(dòng)模糊功能將三幅車(chē)牌圖像復(fù)原，復(fù)原效果和復(fù)原參數(shù)如圖24所示。

圖24 Amped FIVE復(fù)原運(yùn)動(dòng)模糊車(chē)牌

復(fù)原結(jié)果比較運(yùn)動(dòng)模糊參數(shù)比較如表1 所示，人工估計(jì)的運(yùn)動(dòng)模糊參數(shù)與專(zhuān)業(yè)圖像處理軟件處理時(shí)提供的參數(shù)相差不大，可見(jiàn)人工估計(jì)運(yùn)動(dòng)模糊參數(shù)方法可行。

表1 運(yùn)動(dòng)模糊參數(shù)比較

復(fù)原結(jié)果比較如圖25所示，表中選取本文四種復(fù)原算法中效果較好的兩種算法的復(fù)原圖與專(zhuān)業(yè)圖像處理軟件的復(fù)原圖相比較。

圖25 復(fù)原結(jié)果比較（a組）

圖25 復(fù)原結(jié)果比較（b組）

圖25 復(fù)原結(jié)果比較（c組）

四、結(jié)語(yǔ)

本文研究的是公安領(lǐng)域運(yùn)動(dòng)模糊車(chē)牌的復(fù)原問(wèn)題，在計(jì)算運(yùn)動(dòng)模糊核時(shí)采用人工估計(jì)的方法，人工估計(jì)運(yùn)動(dòng)模糊核存在一定的誤差，但從復(fù)原的效果來(lái)看，該方法行之有效，值得推薦。從公安實(shí)戰(zhàn)角度來(lái)看，監(jiān)控系統(tǒng)形成的運(yùn)動(dòng)模糊圖像還存在亮度、對(duì)比度不足，混有光學(xué)模糊等因素影響，若直接使用本文所提出的方法復(fù)原圖像可能達(dá)不到預(yù)期效果。接下來(lái)將結(jié)合圖像處理軟件調(diào)整亮度、對(duì)比度，處理光學(xué)模糊等影響因素做進(jìn)一步研究。