瀝青路面離析的數字圖像評價方法

張爭奇 徐耀輝++胡紅松 石偉 張苛

摘要：為快速有效評價瀝青路面的離析程度，文章采用數字圖像技術研究瀝青路面的非均勻性.文中分析了常用的瀝青路面離析的圖像處理評價方法的不足，提出了基于分析路面紋理中下凹區域的宏觀構造寬度評價瀝青路面離析的新方法.此法通過對宏觀構造中的下凹區域照相并獲得二值圖像，并將二值圖像中對應的下凹區域的平均寬度定義為路表宏觀構造寬度K，以構造寬度K和構造深度TD建立模型，來評價瀝青混合料的集料離析程度.研究結果表明：構造寬度K值與傳統的手工鋪砂法測出的構造深度TD具有良好的相關性，能較好地表征瀝青路面的集料離析程度.

關鍵詞：數字圖像；瀝青路面；宏觀構造寬度K值；構造深度；離析；道路工程

中圖分類號：U416 文獻標識碼：A

瀝青路面集料離析[1-4]是影響路面施工質量的重要因素，常用構造深度值TD來評價瀝青路面的集料離析程度.目前測量路面構造深度TD的方法有鋪砂法和激光法[5-7].傳統的鋪砂法測量路面構造深度效率低、操作繁瑣且主觀隨意性較大；激光法雖然測量構造深度快速、準確，但其所測構造深度值為路面一定長度范圍內的平均值，雖能評價路面的橫向離析，但無法評價路面的縱向和局部離析，且設備價格昂貴，不利于推廣應用.

隨著二維數字圖像技術[8-9]的日趨成熟，這一技術也被引用到瀝青路面質量控制[10]中.許多專家學者也對其進行了研究[11-15]，王端宜教授采用普通數碼相機記錄道路表面構造特征，由圖像的像素值分析得出路表面構造深度.張肖寧教授運用數字圖像處理技術對瀝青混合料內部結構組成的三維圖像重構及虛擬力學試驗等方面進行了研究.彭勇等運用基于數字圖像處理技術，直接對混合料中集料的分布狀態，即瀝青混合料均勻性進行了研究.以上研究采用圖像處理技術研究了瀝青路面表面及內部的材料分布情況，為圖像技術在瀝青路面質量控制研究中起到了促進作用.但鑒于離析問題的復雜性，如何快速、客觀、有效地判斷瀝青路面的離析程度還是一個研究熱點.

因此，本文嘗試以路面紋理中下凹區域的宏觀構造寬度作為研究對象，并在圖像處理技術基礎上，探討宏觀構造寬度與豎向深度之間的相關性，最終通過易于得到的構造寬度K來評價瀝青混合料的離析情況，以達到高效的目的.

1常用的瀝青路面離析的數字圖像評價方

法分析

數字圖像處理技術是指通過數碼相機或掃描儀將一幅連續的圖像離散化為計算機可以處理的信息，并對此信息進行各種處理，以達到某一效果的方法.其過程主要包括圖像采集、圖像處理及圖像分析等.由于數字圖像處理技術具有使用方便、經濟性好、可利用信息量大和形象化儲存等優點，其在道路領域的工程應用已成為國際上的研究熱點.數字圖像處理由兩部分組成：物理圖像轉化為對應的數字圖像和通過各種運算法則[16]從數字圖像中提取有用信息.

自數字圖像技術引用到對路面離析進行評價之后，有許多專家學者對其如何較準確地評價路面離析提出了不同的方法并對其進行了研究.比較常見的方法主要有灰度矩陣擬合構造曲面法、灰度均值法以及分形維數法.

在實驗過程中為了減少或避免圖像采集時原始圖像的信息誤差，需要規范瀝青路面數字圖像采集的操作過程.本文采用一種自制的圖像采集設備來獲取穩定性、效果均較好的原始圖像.此設備如圖1所示.它包括一個遮光暗箱，3個LED燈，一個點光控制器和一個相機.LED燈布設在暗箱里面，通過暗箱外的點光控制器可以調節LED燈的功率，從而達到調節箱內光線強度的作用.采集圖像時將試樣放于遮光暗箱內，調節點光控制器使暗箱內達到理想的光線強度，然后將照相機置于遮光暗箱上的開口處，調節相機焦距來采集圖像.

1.1灰度矩陣擬合構造曲面法

有不少研究利用MATLAB數學軟件將瀝青路面RGB數字圖像轉化為灰度圖像，并得到圖像對應的二維灰度矩陣，然后用灰度矩陣信息擬合構造曲面，最后應用曲面積分等公式計算得到表征路面構造深度的HP.

在操作過程中發現，這種圖像處理模型在實際應用中存在一定的局限性，原因是路面的構造紋理極不規則且復雜，當相機像素越高時，由灰度圖像得到的灰度矩陣數據量越大，且分布離散.很多情況下，由這些三維數據點擬合的空間曲面方程的相關系數偏低，難以保證精度要求，且在下一步的積分求體積中，偏差進一步擴大；同時由灰度矩陣擬合曲面方程涉及到一個龐大的數據量處理過程，要通過反復的調節擬合方式和擬合參數才能獲得相關性較好的曲面方程，這對計算機運行功能要求較高.因此，本文認為基于灰度矩陣擬合構造曲面模型在用于瀝青混合料離析評價的應用前還需要進一步的研究.

1.2灰度均值法

采用數字圖像方法評價瀝青路面離析比較常用的方法為灰度均值法.灰度均值法首先通過數碼相機拍攝路面不同離析狀態下的表面紋理構造圖像，然后利用MATLAB軟件繪制出其相應的路表灰度直方圖，接下來跳過擬合曲面的過程，直接考慮路面圖像灰度均值與構造深度的關系，由灰度均值來預測構造深度.這種方法只需要調用Matlab中的mean2（）函數就可實現.

本文采用室內模式試驗模擬細集料離析（F）、無離析（N）、低離析（L）、中等離析（M）及嚴重離析（H）五種不同集料離析程度的瀝青混合料.分別成型10塊無離析（N）AC20瀝青混合料車轍板和10塊不同程度集料離析（F，N，L，M，H各兩塊）的瀝青混合料車轍板，對車轍板進行圖像采集后，并用MATLAB軟件計算每張圖像的灰度均值，同時采用人工鋪砂法測量其構造深度TD，其結果見表1和表2.

運用統計學方法將無離析車轍板的TD和灰度均值以及不同離析程度車轍板的TD和灰度均值分別進行線性回歸，其結果如圖2所示.

從圖2（a）中可以看出，當瀝青混合料未產生集料離析時，車轍板圖像灰度均值與構造深度TD的線性相關系數達到0.726，說明混合料無集料離析條件下的圖像灰度均值與TD值具有較好的相關性.但從圖2（b）中發現，當瀝青混合料發生不同程度集料離析時，車轍板圖像灰度均值與構造深度TD的線性相關系數下降到0.515，雖然這兩個指標大致符合線性反比關系，但其相關系數小于0.6，Remero提出以相關系數[17-18]大于0.6作為儀器效能的界限值，因此在混合料產生不同程度集料離析時，由灰度均值來預測構造深度TD還是存在較大的誤差，因此本文認為用該指標判別瀝青路面的離析情況還需進一步研究.

1.3分形維數法

分形維數[19-20]法也被用來評價瀝青路面的離析.由于瀝青路面紋理構造復雜、極不規則，且不同離析程度瀝青混合料表面的復雜程度存在顯著的差異，可以采用分形維數對其表面紋理構造進行定量描述.本文采用基于FractalFox分維數計算軟件采用盒子分維法[21]提取路面紋理構造灰度圖像中的分形維數，定量描述不同離析程度瀝青混合料表面的粗糙程度，嘗試評價路面離析.

選用上文中的10塊不同程度集料離析的瀝青混合料車轍板圖像，對所有圖像進行灰度化處理后，分別導入FractalFox軟件計算各自的分形維數，計算結果如表3所示.

將表3的數據運用統計學的方法將D與TD進行線性回歸，結果如圖3所示.并繪制出D隨混合料離析逐漸增大而變化的柱狀圖，如圖4所示.

從圖3可以看到，瀝青混合料構造深度TD值與分形維數D值的線性相關系數為0.618 8，表明兩個指標呈現一定的相關規律，即隨著路面構造深度的增大，其分形維數呈減小趨勢，說明瀝青混合料粗集料離析越嚴重，其表面構造深度越大，而表面紋理構造的復雜程度降低了.

理論上F，N，L，M和H五種集料離析程度混合料的分維數大小排序應為：F>N>L>M>H，但根據圖4所示，雖然隨著混合料集料離析程度的增加，分形維數整體呈遞減趨勢，但是F，N，L三種離析程度混合料的分形維數的大小區別不明顯，L的分形維數反而最大，同時M和H的分形維數值也不易區分.究其原因，F，N，L表面紋理構造的復雜程度雖有差異，但該差異的跳躍性并不明顯，因而由此計算出的分形維數區分度不大，M，H同理.

雖然瀝青路面構造紋理的復雜程度與分形維數存在一定的相關性，但相關規律并不明顯，因此本文認為分形維數暫不適用于瀝青混合料集料離析的評價.

2宏觀構造寬度K值法

路面紋理構造可分為微觀構造和宏觀構造.微觀構造又稱微觀紋理，是路面集料表面的微小構造；宏觀構造又稱宏觀紋理，是路面集料之間構成的間隙.從三維立體角度分析，細集料離析混合料的宏觀構造豎向深度明顯小于粗集料離析混合料，無離析混合料介于兩者之間；從二維平面角度分析，細集料離析混合料的宏觀構造下凹區域的橫向寬度明顯小于粗集料離析混合料，同時無離析混合料介于兩者之間.因此，本文嘗試提取路面紋理中下凹區域即宏觀橫向寬度作為研究對象，通過建立宏觀構造的橫向寬度與豎向深度之間的關聯，從路面的二維圖像中提取三維信息，從而提出一種新的評價離析的方法.

以下是宏觀構造寬度K值的計算模型.

1）選取分割閾值：將灰度圖像轉換為二值圖像的關鍵就是分割閾值的選取，閾值的選取直接決定了二值圖像中白色區域面積的大小，本文根據數字圖像的最大類內、類間距離閾值判定準則，同時考慮路表宏觀構造的漸變特性，由灰度極大值和區域內的灰度均值得到將灰度圖像轉換成二值圖像的分割閾值.計算公式如下：

T=k×（Fmax—F均）

式中：T為分割閾值；Fmax為灰度極大值；F均為灰度均值；k為修正系數，取值范圍為5%～15%，本文取10%.

2）提取白色區域面積：在Matlab軟件中調用rgb2gray（）函數將灰度圖像轉換成二值圖像后，再調用bwarea（）函數快速讀取二值圖像中白色像素點個數X，即為白色區域總面積.

3）計算白色區域平均寬度K′值：調用bwmorph（I，'skel'，inf）函數將二值圖像骨架化，得到一幅骨架圖像（見圖5），該骨架線條是由各白色小區域最中間的像素點連成的，可以看成是白色面積區域的中心線，調用bwarea（）函數快速讀取骨架線條像素點個數Y，即為中心線的“長度”，根據面積除以長度得到寬度的思想，二值圖像中各白色區域平均寬度為K′=X/Y.

4）計算宏觀構造寬度K值：統計圖像內校準尺長度L范圍內對應的數字圖像像素數量M，由公式ε= L/ M得到數字圖像的像素當量，再由公式K=εK′得到路面宏觀構造寬度K值.

分別對10塊無離析（N）AC20瀝青混合料車轍板和10塊不同程度集料離析（F，N，L，M，H各兩塊）的瀝青混合料車轍板進行圖像采集，并計算每塊車轍板的K值，計算結果如表4和表5所示.

將表4和表5的數據運用統計學的方法將K與TD進行線性回歸，結果如圖6和圖7所示.并繪制出K值隨混合料離析逐漸增大而變化的柱狀圖，如圖8所示.

從圖6和圖7可知，當瀝青混合料無離析時，其構造深度TD值與構造寬度K值的線性相關系數R2為0.850，當瀝青混合料產生不同程度集料離析時，其TD值與K值的線性相關系數R2為0.868，說明K值與TD值具有較好的線性相關性，因此通過計算K值來預測路面構造深度TD值是可行的，同時研究表明，當路面產生不同程度集料離析時，K值與TD值的相關系數不降反升，表明K值與TD值的線性相關性不會隨混合料的集料離析而降低.

圖8是不同離析程度混合料的K值變化規律，從圖中可以發現，K值隨混合料粗集料離析程度的增加呈遞增趨勢，且存在明顯的規律性，大致規律為F

鑒于宏觀構造寬度K值的計算方法簡便，圖像采集過程操作簡單，且計算值與傳統的手工鋪砂法具有良好的線性相關性，能較好地表征瀝青路面集料的離析程度，因此推薦基于路面宏觀構造寬度K值的圖像處理模型作為評價瀝青路面非均勻性的主要無損檢測工具之一.

3基于K值評價瀝青混合料集料離析的方法

這里提出兩種利用K值評價瀝青路面集料離析的方法：

1）第1種方法是找出路面圖像K值與構造深度TD值的線性相關性，通過計算K值來預測TD值，最終還是用構造深度TD值的大小來評價路面集料的離析程度.

2）第2種方法是基于本文對瀝青混合料不同離析程度的室內模擬試驗，對不同離析程度的車轍板進行圖像采集并計算其K值，整合出不同離析程度的K值劃分標準.本文通過室內離析混合料的模擬試驗，整合的K值劃分標準如表6所示.

4結論

本文分析比較了幾種數字圖像技術評價瀝青路面的離析可能性的方法，研究發現：

1）由灰度矩陣擬合構造曲面模型在評價瀝青混合料離析中，得到的三維數據點擬合的空間曲面方程的相關系數不到0.2，較難保證精度要求.

2）灰度均值法，在混合料無離析條件下的圖像灰度均值與TD值具有較好的相關性，但當瀝青混合料發生不同程度離析時，車轍板圖像灰度均值與構造深度TD的線性相關系數下降到0.515.

3）瀝青路面構造紋理的復雜程度與分形維數D存在一定的相關性，但相關規律并不明顯，因此本文認為分形維數暫不適用于瀝青混合料離析的評價.

4）K值隨混合料離析程度的增加呈遞增趨勢，存在明顯的規律性，能較好地表征不同離析程度路面的紋理構造特征，因此可用構造深度K值來評價瀝青路面的離析情況.

參考文獻

[1]韓立志，梁增潔，彭余華，等.AC25級配離析的分形評價標準及預測方法[J].長安大學學報：自然科學版，2014，34（6）：27-33.

HAN Lizhi， LIANG Zengjie， PENG Yuhua， et al. Frctal evaluation criteria and prediction method of AC25 gradition segregation[J]. Journal of Changan University：Natural Science，2014， 34（6）： 27-33.（In Chinese）

[2]MEEGODA J N ， HETTIARACHCHI C H， BANDARA N， et al. Correlation of surface texture， segregation，and measuremento f air voids[J]. Proc XIICSMFE Sanfranscisco， 2002，3：19-26.

[3]李想，梁乃興，趙毅.基于數字圖像技術瀝青路面集料離析的研究[J]. 中外公路，2014，34（4）： 90-93.

LI Xiang， LIANG Naixing， ZHAO Yi. Research on digital image technology of asphalt pavement based on aggregate segregation[J].Journal of China & Foreign Highway， 2014， 34（4）： 90-93.（In Chinese）

[4]包秀寧，張肖寧.路面顆粒材料離析性的評價方法[J].華南理工大學學報：自然科學版，2010，38（3）：31-36.

BAO Xiuning， ZHANG Xiaoning. Segregability evaluation of granular pavement materials [J]. Journal of South China University of Technology： Natural Science， 2010， 38（3）： 31-36.（In Chinese）

[5]馬榮貴，王建鋒，李平.瀝青路面構造深度精確檢測方法研究[J].科學技術與工程，2014，14（8）：265-268.

MA Ronggui，WANG Jianfeng， LI Ping.Research on high precision measurement of pavement texture depth[J].Science Technology and Engineering，2014， 14（8）： 265-268.（In Chinese）

[6]JTG E60-2008公路路基路面現場測試規程[S]. 北京：人民交通出版社，2008：90-92.

JTG E60-2008Field test methods of subgrade and pavement for highway engineering[S].Beijing： China Communications Press， 2008：90-92. （In Chinese）

[7]王端宜，李維杰，張肖寧.用數字圖像技術評價和測量瀝青路表面構造深度[J].華南理工大學學報：自然科學版，2004，32（2）：42-45.

WANG Duanyi， LI Weijie， ZHANG Xiaoning. Evaluation and measurement of asphalt pavement surface texture depth with digital image technique[J]. Journal of South China University of Technology： Natural Science， 2004， 32（2）： 42-45.（In Chinese）

[8]MASAD E， MUHUNTHAN B， SHASHIDHAR N， et al.Quantifying laboratory compaction effects on the internal structure of asphalt concrete[J].Transportation Research Record：Journal of the Transportation Research Board， 1999，1691：179-185.

[9]WANG L B， WANG X，MOHAMMAD L，et al.Application of mixture theory in the evaluation of mechanical properties of asphalt concrete[J]. Journal of Materials in Civil Engineering， 2004，16（2） ： 167-174.

[10]燕海峰.基于非均勻性的瀝青路面施工質量控制與評價研究[D].西安：長安大學公路學院，2011：57-60.

YAN Haifeng. Research on construction quality control and evaluation of asphalt pavement based on nonuniformity[D]. Xian： Highway College of Chang'an University， 2011：57-60.（In Chinese）

[11]王端宜，吳文亮，張肖寧，等.基于數字圖像處理和有限元建模方法的瀝青混合料劈裂試驗數值模擬[J]. 吉林大學學報：工學版，2011，41（4）：968-973.

WANG Duanyi， WU Wenliang， ZHANG Xiaoning， et al. Numerical simulation of splitting test of asphalt mixture based on DIPFEM[J]. Journal of Jilin University： Engineering and Technology， 2011，41（4） ： 968-973.（In Chinese）

[12]張肖寧.瀝青混合料材料結構數字圖像技術的基本方法[J].華南理工大學學報：自然科學版，2012，40（10）：166-173.

ZHANG Xiaoning. Basic methods of digital image technology for material structure of asphalt mixtures [J]. Journal of South China University of Technology： Natural Science， 2012，40（10） ： 166-173.（In Chinese）

[13]彭勇，孫立軍，王元清，等.數字圖像處理在瀝青混合料均勻性評價中的應用[J].吉林大學學報：工學版，2007，37（2）：334-337.

PENG Yong，SUN Lijun，WANG Yuanqing， et al. Application of digital image processing in evaluating homogeneity of asphalt mixture [J]. Journal of Jilin University： Engineering and Technology， 2007，37（2）：334-337.（In Chinese）

[14]宋毅.應用分形理論描述瀝青混合料級配離析及其相關問題研究[D].廣州：華南理工大學土木與交通學院 ，2010：6-16.

SONG Yi. Study on grading segregation description of asphalt mixtures and related problems based on fractal theroy[D]. Guang Zhou： College of Civil Engineering and Transportation， South China University of Technology， 2010：6-16.（In Chinese）

[15]周吳軍.基于分形理論瀝青路表面離析判別方法研究[D]. 長沙：湖南大學土木工程學院，2009：28-44.

ZHOU Wujun. Research on construction quality control and evaluation of asphalt pavement based on nonuniformity[D]. Chang

sha： College of Civil Engineering， Hunan University， 2009：28-44. （In Chinese）

[16]劉思峰，方志耕，謝乃明.基于核和灰度的區間灰數運算法則[J].系統工程與電子技術，2010，32（2）：313-316.

LIU Sifeng， FANG Zhigeng， XIE Naiming. Algorithm rules of interval grey numbers based on the“Kernel” and the degree of greyness of grey numbers [J]. Systems Engineering and Electronics，2010，32（2） ： 313-316.（In Chinese）

[17]PENTLAND P. Fractalbased description of natural scenes[J]. IEEE Transactions on Pattern Analysis and Machine Intelligence， 1984，6（6）：661-674.

[18]MANDELBROT B B. The fractal geometry of nature[M]. San Francisco； W.H.Freeman and Company，1982： 361-366.

[19]王耀南，王紹源，毛建旭. 基于分形維數的圖像紋理分析[J]. 湖南大學學報： 自然科學版， 2006，33（5）：67-72.

WANG Yaonan，WANG Shaoyuan，MAO Jianxu. Image texture analysis based on fractal dimension [J]. Journal of Hunan University： Natural Sciences，2006，33（5） ： 67-72.（In Chinese）

[20]趙明華，戴杰，張玲，等. 基于分形理論的粉煤灰滲透率研究[J]. 湖南大學學報：自然科學版，2015，42（1）：75-80.

ZHAO Minghua， DAI Jie， ZHANG Ling， et al. Fractal theorybased study of the permeability of fly ash [J]. Journal of Hunan University： Natural Sciences，2015，42（1）：75-80.（In Chinese）

[21]趙海英，楊光俊，徐正光.圖像分形維數計算方法的比較[J].計算機系統應用，2010，20（3）：238-246.

ZHAO Haiying， YANG Guangjun， XU Zhengguang. Comparison of calculation methodsbased image fractal dimension [J]. Computer Systems & Applications， 2010，20（3） ： 238-246.（In Chinese）