6階數(shù)據(jù)特征的雹云預測方法
【專利摘要】6階數(shù)據(jù)特征的雹云預測方法屬于統(tǒng)計學領(lǐng)域，尤其屬于圖像分析領(lǐng)域?，F(xiàn)有對雹云的預測往往是依賴于經(jīng)驗進行判別，這樣的方法存在很多不足。為了克服現(xiàn)有技術(shù)的不足，本發(fā)明公開：一種預測冰雹的方法，其特征是，所述的方法包括以下步驟：第一步，利用云圖獲得dBZ，第二步，把RGB圖像轉(zhuǎn)化為灰度圖像，第三步，灰度圖像的反射率分成15個量級，第四步，計算每個量級的比重，第五步，計算出每個量級的均值x1、方差x2、峰度x3、角二階矩x4、熵x5、相關(guān)性x6，第六步，利用公式：y＝148.50246·x1+44.33304·x2?213.01248·x3?123.00231·x4+22.56971·x5+3457·x6y0＝41.16641故有若：y＞y0,則X∈G1，否則X∈G2，G1為降雹云集合，G2為無雹云集合。本發(fā)明可以預測雹云。
【專利說明】
6階數(shù)據(jù)特征的雹云預測方法
技術(shù)領(lǐng)域
[0001] 本發(fā)明屬于統(tǒng)計學領(lǐng)域，尤其屬于圖像分析領(lǐng)域。
【背景技術(shù)】
[0002] 冰雹作為一種強對流天氣，其特點是空間尺度小、生命史短、突發(fā)性強、發(fā)展演變 迅速，其預報難度是眾所周知的。經(jīng)過"十五"、"十一五"的發(fā)展和建設(shè)，新疆氣象防雹技術(shù) 水平有了很大提升，但與國民經(jīng)濟和社會發(fā)展日益增長的需求相比，仍然存在一定不足，主 要表現(xiàn)為:冰雹天氣預測的準確率和精細化水平有待提高，特別是對局部冰雹天氣的臨近 (0~3小時）和短時（3~12小時)預報能力亟待增強。"十三五"期間是新疆社會穩(wěn)定和經(jīng)濟 發(fā)展的關(guān)鍵時期，迫切需要對新疆冰雹天氣提供更為精細化的監(jiān)測和預報。而對雹云識別 業(yè)務(wù)基礎(chǔ)主要來源于天氣雷達，實時觀測對流云團的雷達圖像，圖像涵蓋了云團生命發(fā)展 變化的信息，通過對該云團的雷達圖像進行實時處理，提前預測到該云團是否降雹對于防 雹減災有重要意義。
[0003] 現(xiàn)有對雹云的預測往往是依賴于經(jīng)驗進行判別，這樣的方法存在很多不足。

【發(fā)明內(nèi)容】

[0004] 為了克服現(xiàn)有技術(shù)的不足，本發(fā)明公開：
[0005] 1，6階數(shù)據(jù)特征的雹云預測方法，其特征是，所述的方法包括以下步驟：
[0006] 第一步，采集雷達RGB云圖，利用云圖獲得dBZ，
[0007] 第二步，把第一步所得的RGB圖像轉(zhuǎn)化為灰度圖像，
[0008] 第三步，把第二步所得的灰度圖像的反射率分成15個量級，
[0009] 第四步，計算出每個量級在圖像中的比重，
[0010]第五步，計算出每個量級的均值xl、方差χ2、峰度x3、角二階矩x4、熵x5、相關(guān)性x6，
[0011] 第六步，把第五步所得的均值X1、方差x2、峰度x3、角二階矩x4、熵x5、相關(guān)性x6代 入公式：
[0012] y = 148.50246 · χι+44·33304 · X2-213.01248 · X3-123.00231 · X4+22.56971 · X5+ 3457 · X6
[0013]故有若:y>y。，則XeGi否則XGG2，GA降雹云集合，62為無雹云集合。
[0014] 2,其特征是：所述 y〇 = 41.16641。
[0015] 3,其特征是:還包括：
[0016] 利用公式：
[0018] 得到 y〇。
[0019] 降雹云集合Gi中的平均值，降雹云集合G2中的平均值，m為降雹云集合Gi中 的元素個數(shù)，Π 2為無雹云集合G2中的元素個數(shù)。
[0020]反射率按照5為單位分開，(^2值81對應(yīng)的灰度值計算為S1-2.5到81+2.5的灰度平 均值。
[0021 ] 4，其特征是:所述的反射率、量級、RGB值、灰度值對應(yīng)的關(guān)系為：
[0022] 反射率為:-5,對應(yīng)的量級為：1，對應(yīng)的RGB值為:-201，201，201，對應(yīng)的灰度值為： 201；
[0023] 反射率為:0，對應(yīng)的量級為：2，對應(yīng)的RGB值為:-118，118，118，對應(yīng)的灰度值為： 118；
[0024] 反射率為：5，對應(yīng)的量級為：3，對應(yīng)的RGB值為:-255，170，170，對應(yīng)的灰度值為： 196；
[0025] 反射率為：10，對應(yīng)的量級為:4，對應(yīng)的RGB值為:-238，140，140，對應(yīng)的灰度值為： 170；
[0026] 反射率為：15，對應(yīng)的量級為:5，對應(yīng)的RGB值為:-201，112，112，對應(yīng)的灰度值為： 139；
[0027] 反射率為：20，對應(yīng)的量級為：6，對應(yīng)的RGB值為：0，251，144，對應(yīng)的灰度值為： 164；
[0028] 反射率為:25，對應(yīng)的量級為:7，對應(yīng)的RGB值為:0，187，0，對應(yīng)的灰度值為：110;
[0029] 反射率為:30，對應(yīng)的量級為:8，對應(yīng)的RGB值為:-255，255，112，對應(yīng)的灰度值為： 239；
[0030] 反射率為：35，對應(yīng)的量級為：9，對應(yīng)的RGB值為：208，208，96，對應(yīng)的灰度值為： 195；
[0031] 反射率為:40，對應(yīng)的量級為：10，對應(yīng)的RGB值為：255，96，96，對應(yīng)的灰度值為： 144；
[0032] 反射率為:45，對應(yīng)的量級為:11，對應(yīng)的RGB值為:218，0，0，對應(yīng)的灰度值為:65;
[0033] 反射率為:50，對應(yīng)的量級為12，對應(yīng)的RGB值為：174，0，0，對應(yīng)的灰度值為:52;
[0034] 反射率為:55，對應(yīng)的量級為：13，對應(yīng)的RGB值為:0，0，255，對應(yīng)的灰度值為:29;
[0035] 反射率為：60，對應(yīng)的量級為：14，對應(yīng)的RGB值為：-160，255，255，對應(yīng)的灰度值 為:227;
[0036] 反射率為：65，對應(yīng)的量級為：15，對應(yīng)的RGB值為：231，0，255，對應(yīng)的灰度值為： 98 〇
[0037] 5,其特征是：由RGB與灰度值的換算關(guān)系式
[0038] Gray = R X 0.2999+G X 0.587+B X 0.114
[0039] 獲得對應(yīng)的灰度值。該方法可以提高雹云預報的準確率，降低誤報率。有效的提高 了識別雹云圖像準確率，為冰雹預報提供參考。
【具體實施方式】：
[0040] 1強度統(tǒng)計
[0041] 雷達反射率圖像是基于RGB彩色空間的，以基本反射率因子圖的圖例為基礎(chǔ)，利用 RGB彩色空間的成色原理，讀取回波反射圖。即將反射率為-5dBZ-65dBZ，分為15種相應(yīng)的 量級。其中不同量級所對應(yīng)的RGB值如表所示，由RGB與灰度值的換算關(guān)系式（1)，可知對應(yīng) 的灰度值。
[0042] Gray = R X 0.2999+G X 0.587+B X 0.114( 1)
[0043] 表1用來說明反射率與RGB值、灰度值對應(yīng)關(guān)系。
[0044] 表1:
[0047] 雷達反射率圖像的不同顏色代表的是不同的反射強度，即可以通過統(tǒng)計不同顏色 在圖像中所占的比重，反應(yīng)出圖像中回波強度的特征，不同的強度所占的比重可以由概率 表示，計算如式⑵所示。如一幅大小為MXN的圖像，則含有MXN個像素，反射率強度分為15 個量級，這不同強度在圖像中的比重為：
[0048] Pi = si/MXN i = l,2,---,15 (2)
[0049] 其中，Sl為第i個量級的顏色在圖像上的像素個數(shù)。
[0050] 2、六階雹云識別的模型設(shè)計
[0051]根據(jù)氣象學的理論可知，強對流的極端天氣，只可能出現(xiàn)在反射強度大于40dBZ的 區(qū)域。故只選取量級為10-15的，按照以下方法計算得到雹云和非雹云的四個一階統(tǒng)計測 度值。
[0052]均值：

[0070] m為降雹云集合Gi中的元素個數(shù)，n2為無雹云集合G2中的元素個數(shù)。
[0071]判別函數(shù)為：
[0072] y = 148.50246 · χι+44·33304 · X2-213.01248 · X3-123.00231 · X4+22.56971 · X5+ 3457 · X6 (3)
[0073]
[0074] 判別臨界點為:y〇 = 41.16641，且;^(I> $721，故有若:y>yQ則XGGi，否貝iJXeGsjA 降雹云集合，62為無雹云集合。
[0075] 3、模型檢驗
[0076]利用本發(fā)明所述的方法，發(fā)明人進行了如下驗證，以證明本發(fā)明對預測冰雹有較 高的準確率。
[0077]數(shù)據(jù)選取2009年間石河子、沙灣地區(qū)及阿克蘇地區(qū)的不同于訓練樣本的數(shù)據(jù)，作 為待測樣本表2。由于待測樣本的分類已知，則利用此判別函數(shù)進行判別，并將判別結(jié)果與 真實結(jié)果進行比較，即可知道判別函數(shù)對于反射率圖像判別的準確率。結(jié)果如表3所示。
[0078]表2各地區(qū)的反射率圖像強度的一階統(tǒng)計特征-待測樣本
[0083]由表2和表3可知，此判別模型的判別準確率很好，且該模型對于雹云的識別是比 較準確的，準確率為100 %。
【主權(quán)項】
1.6階數(shù)據(jù)特征的富云預測方法，其特征是，所述的方法包括w下步驟： 第一步，采集雷達RGB云圖，利用云圖獲得地Z， 第二步，把第一步所得的RGB圖像轉(zhuǎn)化為灰度圖像， 第Ξ步，把第二步所得的灰度圖像的反射率分成15個量級， 第四步，計算出每個量級在圖像中的比重， 第五步，計算出每個量級的均值X1、方差x2、峰度x3、角二階矩x4、賭巧、相關(guān)性x6， 第六步，把第五步所得的均值XI、方差x2、峰度x3、角二階矩x4、賭巧、相關(guān)性x6代入公 式： 7 = 148.50246 · X1+44.33304 · X2-213.01248 · X3-123.00231 · X4+22.56971 · X已+ :3457 · X6 故有若:y >y〇，貝化E Gi，否貝化e G2，G功降富云集合，G勸無富云集合。2. 如權(quán)利要求1所述的6階數(shù)據(jù)特征的富云預測方法，其特征是:所述y〇 = 41.16641。3. 如權(quán)利要求1所述的6階數(shù)據(jù)特征的富云預測方法，其特征是:還包括： 利用公式：/1降富云集合G沖的平均值，J0)降富云集合G2中的平均值，m為降富云集合Gi中的元 素個數(shù)，Π 2為無富云集合G2中的元素個數(shù)， 得到y(tǒng)O。4. 如權(quán)利要求1所述的6階數(shù)據(jù)特征的富云預測方法，其特征是:所述的反射率、量級、 RGB值、灰度值對應(yīng)的關(guān)系為： 反射率為:-5，對應(yīng)的量級為：1，對應(yīng)的RGB值為:-201，201，201，對應(yīng)的灰度值為:201; 反射率為:0，對應(yīng)的量級為:2，對應(yīng)的RGB值為:-118，118，118，對應(yīng)的灰度值為：118; 反射率為:5，對應(yīng)的量級為:3，對應(yīng)的RGB值為:-255，170，170，對應(yīng)的灰度值為：196; 反射率為：10，對應(yīng)的量級為:4，對應(yīng)的RGB值為:-238，140，140，對應(yīng)的灰度值為：170; 反射率為：15，對應(yīng)的量級為:5，對應(yīng)的RGB值為:-201，112，112，對應(yīng)的灰度值為：139; 反射率為:20，對應(yīng)的量級為:6，對應(yīng)的RGB值為:0，251，144，對應(yīng)的灰度值為：164; 反射率為:25，對應(yīng)的量級為:7，對應(yīng)的RGB值為:0，187，0，對應(yīng)的灰度值為：110; 反射率為:30，對應(yīng)的量級為:8，對應(yīng)的RGB值為:-255，255，112，對應(yīng)的灰度值為:239; 反射率為:35，對應(yīng)的量級為:9，對應(yīng)的RGB值為:208，208，96，對應(yīng)的灰度值為：195; 反射率為:40，對應(yīng)的量級為：10，對應(yīng)的RGB值為:255，96，96，對應(yīng)的灰度值為：144; 反射率為:45，對應(yīng)的量級為：11，對應(yīng)的RGB值為:218，0，0，對應(yīng)的灰度值為:65; 反射率為:50，對應(yīng)的量級為12，對應(yīng)的RGB值為：174，0，0，對應(yīng)的灰度值為:52; 反射率為:55，對應(yīng)的量級為：13，對應(yīng)的RGB值為:0，0，255，對應(yīng)的灰度值為:29; 反射率為：60，對應(yīng)的量級為：14，對應(yīng)的RGB值為：-160，255，255，對應(yīng)的灰度值為： 227； 反射率為:65，對應(yīng)的量級為：15，對應(yīng)的RGB值為:231，0，255，對應(yīng)的灰度值為:98。5.如權(quán)利要求1所述的6階數(shù)據(jù)特征的富云預測方法，其特征是：由RGB與灰度值的換算 關(guān)系式 Gray = R X 0.2999+G X 0.587+B X 0.114，獲得對應(yīng)的灰度值。
【文檔編號】G06K9/00GK106096546SQ201610405806
【公開日】2016年11月9日
【申請日】2016年6月8日
【發(fā)明人】徐文霞, 李國東, 吳晨瑜, 錢斯祺
【申請人】李國東