趙曉莉,苑 躍,黃曉龍,劉相全

(1.四川省氣象災害防御技術中心，四川 成都 610072；2.四川省氣象探測數(shù)據(jù)中心，四川 成都 610072；3.重慶眾仁科技有限公司，重慶 408000)

0 引言

自20世紀60年代初，全國逐步開展風向風速自動觀測，主要使用EL型電接風向風速計和達因式風向風速觀測儀器，并以自記紙的形式保存了近40 a的風向風速歷史資料[1]。風自記紙是寶貴的客觀表征風向風速日變化的唯一實測的基礎資料，進行數(shù)字化處理后，這些資料對于研究我國風向風速特征研究具有重要意義[2]；將為開發(fā)利用風能資源、預防風災以及城市建設、環(huán)境保護等提供權(quán)威的基礎氣象科學數(shù)據(jù)[3-9]。

全國已完成EL電接風向風速自記紙的圖像掃描工作，數(shù)字化自記紙圖像中數(shù)據(jù)的自動識別將成為亟待解決的問題。本文根據(jù)《地面氣象觀測規(guī)范》(2003版)[10]中對EL型電接風自記紙整理規(guī)范要求，研究EL型電接風自記紙風向數(shù)字化識別方法，并對該方法開發(fā)的風向提取軟件的風向識別效果進行檢驗。

1 EL型電接風向風速計

1.1 設備介紹

EL型電接風向風速計是一種由風速風向傳感器、指示器和記錄器3部分組成的風向風速有線遙測儀器。記錄器由8個風向電磁鐵、1個風速電磁鐵、自記鐘、自記筆、筆擋、充放電線路等部分組成[11-14]。每兩個磁鐵為一組控制一支記錄筆，當控制風向記錄時刻的繼電器吸合時，與風向標所在方位的方位塊相對應的電磁鐵吸動，控制相應的記錄筆記下風向標所在方位。

1.2 EL型電接風向風速自記紙

EL型電接風向風速自記紙上的數(shù)據(jù)是記錄器記錄所得，分為兩個部分：上部為風向，由4根平行線分割成8份，按照從上到下的順序分別表示S、N、NE、SW、E、W、SE、NW 8個風向，剩余8個風向(NNE、ENE、ESE、SSE、SSW、WSW、WNW、NNW)則通過記錄器在同一時間兩個記錄筆同時畫下的風向線組合得到；下部為風速，按空氣行程200 m電接1次，風速自記筆相應跳動1次來記錄，正點前10 min內(nèi)的風速根據(jù)跡線通過自記紙上平分格線的格數(shù)(1格相當于1.0 m/s)計算[15-22]。

1.3 自記紙風向的整理

根據(jù)《地面氣象觀測規(guī)范》(2003版)對EL型電接風向風速計自記整理的要求，自記紙從各正點前10 min內(nèi)的5次風向記錄中挑取出現(xiàn)次數(shù)最多的。如最多風向有兩個出現(xiàn)次數(shù)相同，應舍去最左面的1次劃線，而在其余4次劃線中來挑?。蝗羧杂袃蓚€風向相同，再舍去左面的1次劃線，按右面的3次劃線來挑取。如5次劃線均為不同方向，則以最右面的1次劃線的方向作為該時記錄。

2 EL型電接風向風速計自記紙圖像預處理

風自記紙經(jīng)掃描后形成的自記紙數(shù)字圖像會有輕微的傾斜或變形，需要通過圖像預處理糾正這些錯誤，以保證數(shù)據(jù)的精度。通過局部投影取極值，找長直線，即對應邊界線的準確位置；再根據(jù)找出長直線的傾斜角，獲取圖像的傾斜角，并對圖像進行傾斜校正。

根據(jù)業(yè)務要求，觀測員每天需對風自記紙進行整理，并將人工從自記紙上讀取的數(shù)據(jù)記錄在風自記紙上，這樣勢必造成人工字跡影響風向、風速跡線的情況。因此，需要刪除干擾跡線的人工字跡，以保證后續(xù)跡線自動識別的準確性，一般通過對紅(R)、綠(G)、藍(B) 3個顏色通道閾值的設置區(qū)分人工記錄的鉛筆字跡與風向、風速自記跡線。

3 風向跡線與底圖的分離

風向風速自記紙掃描要求24位彩色，標準掃描分辨率為200 DPI，為了提高圖片質(zhì)量，掃描分辨率一般大于200 DPI。計算機圖像處理中，圖像的色彩是通過對紅(R)、綠(G)、藍(B) 3個顏色通道的變化以及它們相互之間的疊加來得到各式各樣的顏色。自記紙圖像中，底色和坐標的紅色分量較大，綠色和藍色分量較小，而風向、風速跡線的紅色分量相對較小，綠色和藍色分量較大；綜合效果看，曲線與底色、坐標色的灰度差異不如單純的紅色差異顯著。因此，區(qū)分底色與風向、風速跡線采用紅色分量作為閾值條件。一般情況下，以小于等于200判定為風向、風速跡線，大于200判定為底色和坐標背景信息，以此來區(qū)分自記紙底圖坐標與跡線，最終提取風向跡線，如圖1所示。

4 風向數(shù)據(jù)提取

根據(jù)風向線自記紙的特征將風向線數(shù)據(jù)提取分成4個子圖單獨處理，每個子圖處理包含2個步驟：①根據(jù)風向線與底圖分離后風向線點集自動獲取風向線所在位置；②人機交互。

4.1 自動跟蹤風向線

根據(jù)記錄器每2.5 min記錄一次瞬時風向的規(guī)則，設有風向為1，無風向為0。將風向軸線上下1/3軸線距矩形區(qū)域(如圖2所示)是否有豎直投影作為風向的判別依據(jù)。若同一時次風向軸線上部、下部均有投影，則根據(jù)投影值的大小進行判別，投影值大認為有風向記為1，對應時次投影值相對小的則記為0；若同一時次風向軸線上部、下部均無投影，則對應時次均記0。以圖2風向線子圖中14時30分風向跡線為例，以風向軸為分界線，分別計算該時刻風向軸上部和下部風向跡線的像素點之和，而后求差值，結(jié)果小于0，說明下部投影值大于上部，故下部風向應記1，上部風向記0。

將4個風向線子圖同一時刻風向判別結(jié)果(即8×1的數(shù)組)，作為該時刻瞬時風向記錄。以某站17時正點前10 min風向自記為例(如圖3所示)，綠色框內(nèi)的8×1的數(shù)組為17時正點的瞬時風向數(shù)組，利用風向組合規(guī)則(如表1所示)，判斷該時刻瞬時風向為ENE。按照上述方法判斷17時正點前10 min內(nèi)的5次瞬時風向均為ENE，故17時10 min風向為ENE。

圖1 風向跡線提取前后對比圖Fig.1 Contrast map before and after wind direction trace extraction(a, Part of wind self-recording paper, b, Recognition of the wind direction trace)

圖2 風向線子圖處理區(qū)域示意圖Fig.2 Subgraph of wind direction line processing area diagram

圖3 瞬時風向判別示意圖Fig.3 Diagram of instantaneous wind direction discrimination

表1 瞬時風向識別對照表Tab.1 look-up tab of Instantaneous wind direction recognition

4.2 人機交互

人機交互是通過人工手動操作的方式修正自動跟蹤錯誤的風向線的過程，包括對風向線進行增加、刪除及修改。通過鼠標點擊風向記錄線的方式進行個別風向線增加和刪除。人工修改大面積風向線通過繪制等距風向線實現(xiàn)：鼠標框選需要修改的風向線大致范圍，在鼠標點擊附近搜尋第一根和最后一根風向線位置，根據(jù)風向線間最大長度以及搜尋的第一根風向線和最后一根風向線之間的長度綜合計算風向線平均間距，等間距繪制新的風向線用于替換修改區(qū)域內(nèi)原有風向線。

5 風向識別效果檢驗

根據(jù)上述EL型電接風自記紙風向數(shù)字化識別方法研發(fā)了EL型電接風自記紙數(shù)字化軟件，并隨機對11個臺站100張EL型風速風向自記紙(如表2所示)共計2 400個風向正點值進行了效果測試。

表2 測試臺站及測試EL型風自記紙數(shù)量Tab.2 Test station and EL wind self-recording paper quantity

為了驗證風向識別方法的可用性，將軟件識別的正點10 min風向與全國地面氣象資料A6文件中的數(shù)據(jù)進行對比，其中，A6文件中的正點10 min風向為人工整理EL型風自記紙讀取的記錄。風向完全一致認為軟件提取風向數(shù)據(jù)正確，不一致則需要人工核查風自記紙，判斷是否是由于人工讀數(shù)錯誤造成不一致。如表3所示，風向一致率較高，軟件提取風向數(shù)據(jù)正確率達97.2%；人工讀取錯誤導致的不一致占2.7%；軟件讀數(shù)錯誤僅占0.2%。軟件讀取風向錯誤主要與圖像中跡線的清晰程度有較大的關系，圖像跡線越清晰軟件讀取正確率越高，且正確率優(yōu)于人工讀取。以綿陽站(56 196)為例，該站自記紙圖像較為清晰，風向線整齊有序，測試該站1個月的風自記紙圖像，僅有15個風向與A6文件不一致，經(jīng)核查這15個風向均為人工讀數(shù)錯誤。

表3 測試結(jié)果與A6文件對比Tab.3 Comparisons between test results and A6 files

6 小結(jié)

本文主要介紹軟件風向識別算法，算法以《地面氣象觀測規(guī)范》(2003版)中EL型電接風自記紙整理規(guī)范為依據(jù)，參考人工讀取EL型電接風向風速自記紙方式建立。利用該識別算法開發(fā)的EL型電接風自記紙數(shù)字化軟件，從測試結(jié)果可見，基于風向識別算法開發(fā)的EL型電接風自記紙數(shù)字化軟件風向識別正確率達97.2%，表明本文研究的風向識別方法可以應用到EL型電接風自記紙數(shù)字化處理中。