圖像法快速測量色斑直徑和雨滴直徑

盧俊輝，楊志紅，鄭先科

（江漢大學 物理與信息工程學院，湖北 武漢 430056）

圖像法快速測量色斑直徑和雨滴直徑

盧俊輝，楊志紅，鄭先科

（江漢大學 物理與信息工程學院，湖北 武漢 430056）

采用圖像法處理濾紙色斑測量雨滴直徑。應用距離變換生成分水嶺，以分水嶺對濾紙色斑進行分割，獲取濾紙上色斑面積，計算色斑等效直徑 D，率定自然雨滴直徑 d=1.1～4.1 mm，與色斑直徑關系d=0.356 1×D0.7233。圖像法不受色斑毛刺影響，色斑等效直徑測量準確，工作效率高。

圖像處理；雨滴直徑；色斑直徑

0 引言

雨滴譜的觀測方法比較多，包括雷達觀測法［1］、雨滴譜儀觀測法［2］、濾紙色斑法［3］等，濾紙色斑法原理明確、設備簡單。

濾紙色斑法是基于雨滴在濾紙上形成的色斑面積與雨滴的體積成正比的假設。大量試驗分析表明，色斑直徑與雨滴直徑具有冪函數關系：

式中：D為色斑直徑，K和n分別為系數和指數，K、D和n與濾紙的種類有關。

測量濾紙色斑直徑，以往采用質量法或圖形描繪法。質量法是將雨滴色斑描繪在透明的硫酸紙上，在電子天平上稱出帶有色斑紙的質量，與單位紙的質量對比，計算出色斑直徑。圖形描繪法是采用autocad、CorelDraw等繪圖軟件描繪出色斑圖案，并用該軟件直接計算出單個色斑直徑［4-5］。但是，雨滴在濾紙上形成的色斑帶有形狀不規則的毛刺，色斑面積測量結果存在較大的誤差，測量工作相當復雜。

筆者采用Matlab圖像法［6］，根據色斑位置對濾紙色斑圖像進行分割，獲取每個色斑圖像面積，色斑本身的毛刺也包含在內，通過色斑面積計算色斑直徑。圖像處理法直接檢測色斑像素，色斑直徑測量精度好、工作效率高。

1 色斑形成

1.1 色斑材料

濾紙選用直徑15 cm定性中速濾紙，涂料選用曙紅水溶性染料和滑石粉的混合粉末，其質量比為1∶8，使用前將混合粉末薄薄地均勻涂在濾紙上，當雨滴落在濾紙上時，留下永久性的近似圓形色斑。

1.2 滴定針頭

滴定針頭選用不銹鋼點膠針頭，同時采用了水滴細分方法獲得更小的水滴直徑［7］，點膠針頭編號、針頭內徑及產生的水滴直徑見表1，其中26a～26l是對26號點膠針頭的細分編號。

表1 點膠針頭編號及水滴直徑dTab.1 Tube number and waterdrop diameterd

1.3 水滴色斑

為了獲得近似自然雨滴所形成的色斑，點膠針頭安裝在高度12 m的平臺上，用濾紙承接點膠針頭形成一系列直徑的水滴，獲得對應的色斑，所形成的色斑與實際雨滴所形成的色斑近似，共得到25張帶有色斑的濾紙，所對應的水滴直徑為1.1～4.1 mm。其中圖1（a）為直徑3.5 mm水滴從12 m高度落下形成的色斑，圖1（b）為自然雨滴形成的色斑。

圖1 水滴色斑和雨滴色斑Fig.1 Waterdrop stain and raindrop stain

2 圖像方法獲取色斑直徑

圖像法測量色斑直徑的過程，以圖1（a）中直徑3.5 mm水滴色斑直徑測量為例詳細說明。

2.1 色斑預處理

利用掃描儀以1∶1的比例掃描濾紙色斑，將濾紙圖像統一調整為500×500像素，刪除殘缺和重疊的色斑，適當調整圖像的對比度和飽和度，預處理后的濾紙色斑見圖2（a）。

2.2 色斑二值化

因濾紙上粉末涂抹不均勻，導致色斑顏色不均勻。為了準確獲取色斑圖像，采用二值化增強色斑顏色，色斑圖像二值化見圖2（b），色斑不均勻得到改善。

圖2 水滴色斑及其二值圖Fig.2 Waterdrop stain and its binary image

2.3 色斑形態學處理

有些色斑二值圖像存在缺陷，采用圖像學進行膨脹運算重構色斑，并對色斑進行腐蝕開、閉運算，色斑開運算后見圖3（a），色斑閉運算后見圖3（b），消除色斑存在的缺陷，生成了完整的實心色斑。

圖3 二值圖開運算、閉運算Fig.3 Open operation and close operation of binary image

2.4 色斑距離變換和分水嶺形成

根據色斑位置對圖像進行分割，需要找出色斑中心。色斑中心通過距離變換實現，距離變換結果如圖4（a），顏色越深，代表該像素越靠近色斑中心，根據色斑中心生成了8個圓盤，見圖4（b），其距離局部最小區域見圖4（c），再根據區域劃分在圖像內部形成分水嶺見圖4（d）。

圖4 二值圖距離變換和分水嶺Fig.4 Distance transformation and watershed of binary image

2.5 色斑分割

以分水嶺為邊界對色斑二值圖進行分割，將圖3（b）水滴色斑圖分割為8個區域，同時右上角2個粘連的水滴色斑被分割，8個色斑的像素值采用1～8進行標記，見圖5（a）。按照同樣方法將圖1（b）雨滴色斑圖分為44個區域，有部分粘連的色斑也被分割，即有44個雨滴，44個色斑的像素值采用1～44進行標記，見圖5（b）。

2.6 色斑直徑

通過像素標記，計算出分割的8個水滴色斑面積，色斑平均直徑29.3 mm。雨滴色斑最大直徑30.8 mm，雨滴色斑最小直徑28.1 mm。

圖5 水滴色斑和雨滴色斑分割Fig.5 Image segmentation of waterdrop and raindrop

3 雨滴色斑直徑識別

3.1 雨滴直徑與色斑等效直徑關系

按照上述步驟，對直徑d=1.1～4.1 mm水滴形成的色斑進行了圖像處理，分別獲得色斑對應的平均直徑D=4.6～29.3 mm，見表2。對表2水滴直徑和色斑直徑進行指數擬合，擬合關系見式（2），與文獻［3，8-10］的擬合關系對比見圖6，擬合關系效果很好。

表2 水滴直徑d和色斑直徑DTab.2 Waterdrop diameterdand stain diameterD /mm

圖6 色斑直徑和雨滴直徑擬合Fig.6 Image fitting of stain diameter and raindrop diameter

3.2 雨滴色斑直徑識別

將圖5（b）所測量的色斑直徑代入式（2），計算出圖5（b）中的最大雨滴直徑為d=4.2 mm，最小雨滴直徑d=0.4 mm，符合自然雨滴直徑范圍。

4 結語

采用點膠針頭，以12 m高度落下水滴形成的色斑為基準，采用圖像處理法測量濾紙色斑直徑，率定了水滴直徑d=1.1～4.1 mm與其色斑直徑的關系，并對自然雨滴色斑進行了測量，圖像法不受色斑毛刺影響，色斑直徑測量準確，工作效率高。

（

）

［1］趙振維，吳振森，沈廣德.利用雨滴尺寸分布數據確定雷達測雨參數［J］.電子學報，2000，28（3）：109-111.

［2］LOFFLER-MANG M，JOSS J.An optical disdrometer for measuring size and velocity of hydrometeors［J］.Atmos and Oceanic Technology，2000，17（2）：130-139.

［3］舒若杰，高建恩，吳普特，等.基于計算機繪圖軟件的雨滴譜測定新方法［J］.中國水土保持科學，2006，4（3）：65-69.

［4］李紅，任志遠，湯躍，等.噴頭噴灑雨滴粒徑測試的改進研究［J］.農業機械學報，2005，36（10）：50-53.

［5］廖煒，衛苗苗，黃煜煜.采用濾紙色斑法對雨滴直徑的研究［J］.武漢理工大學學報：交通科學與工程版，2008，32（6）：1165-1168.

［6］岡薩雷斯.數字圖像處理的MATLAB實現［M］.阮秋琦，譯.北京：清華大學出版社，2013.

［7］盧俊輝.用于液滴細分的同軸雙孔針頭、裝置及方法［P］.CN：201410428273.6，2014-08-27.

［8］徐向舟，張紅武，朱明東.雨滴粒徑的測量方法及其改進研究［J］.中國水土保持，2004（2）：22-24.

［9］竇葆璋，周佩華.雨滴的觀測和計算方法［J］.水土保持通報，1982，2（1）：44-47.

［10］陳振宇，葉全民，趙美香，等.噴霧霧滴濾紙率定法研究初探［J］.山西農業大學學報，1996，16（4）：422-425.

（責任編輯：陳 曠）

Measuring of Stain Diameter and Raindrop Diameter with Image Processing Method

LU Junhui，YANG Zhihong，ZHENG Xianke
（School of Physics and Information Engineering，Jianghan University，Wuhan 430056，Hubei，China）

Mainly introduces a kind of measuring methods for raindrop size.On the basis of the theory of distance transformation and image segmentation，the raindrop can be mapped to stain pixel，and its diameter is calculated easily.According to the experiments，the relation between the raindrop diameter（d）and the stain diameter（D）is d=0.356 1×D0.7233.This image processing method is highly accurate and efficient because it is not affected by splash burr.

image processing；raindrop diameter；stain diameter

P426.62；TP391.41

A

1673-0143（2015）02-0116-05

10.16389/j.cnki.cn42-1737/n.2015.02.004

2015-01-05

武漢市青年晨光計劃資助項目（2013071004010470）

盧俊輝（1974—），男，副教授，博士，研究方向：環境感知和環境識別。