全天空成像儀云量計(jì)算方法的改進(jìn)

周文君，牛生杰，許瀟鋒

(1．南京信息工程大學(xué)中國氣象局大氣物理與大氣環(huán)境重點(diǎn)開放實(shí)驗(yàn)室，江蘇南京210044;2．江蘇省鹽城市氣象局，江蘇鹽城224000)

0 引言

云是表征天氣、氣候特點(diǎn)的重要因素之一，也是大氣動(dòng)力、熱力、水分輸送過程綜合作用的外在表現(xiàn)。云對(duì)地氣系統(tǒng)輻射平衡影響很大，其輻射強(qiáng)迫對(duì)地球氣候變化有著非常重要的作用(陳洪濱，1997;Long and Thomas，2000;Calbo and Pages，2002;Long，2004)。因此，對(duì)云的探測極為重要。目前，已經(jīng)有一系列的探測手段對(duì)云的各個(gè)要素進(jìn)行觀測分析，例如:飛機(jī)穿云觀測云中微物理特征，衛(wèi)星遙感為云的分布及微物理特征量的研究提供了重要依據(jù)，以及地面對(duì)云的目測或使用地基測云儀測量云量、云高等(陳渭民等，1997;宋慶利，2003;Martins et al．，2003)。

隨著科技的不斷進(jìn)步，地基測云儀的技術(shù)也得到了快速發(fā)展。美國大氣輻射研究機(jī)構(gòu)為了定量研究云量對(duì)輻射的影響，與國家海洋大氣科研所以及氣候能源工作實(shí)驗(yàn)室合作研發(fā)了半球天空成像儀(hemispheric sky imager，HSI)。該儀器可以在白天拍攝全天空160°范圍內(nèi)的圖像。Long and Deluisi(1998)詳細(xì)地介紹了 HSI儀器的觀測原理和計(jì)算方法，并對(duì)HSI的云量觀測資料進(jìn)行了分析，結(jié)果表明:HSI儀器的觀測結(jié)果比衛(wèi)星和雷達(dá)的觀測結(jié)果更精確。之后，在美國YANKEE公司的合作之下，HSI被發(fā)展為TSI-440(total sky imager 440，全天成像儀)。該儀器在圖像處理技術(shù)及構(gòu)建方面均取得了進(jìn)一步發(fā)展，是目前使用較為廣泛的地基測云儀。

TSI-440儀器可以自動(dòng)進(jìn)行全天空云量的持續(xù)性觀測，時(shí)空分辨率較高，得到的云量計(jì)算結(jié)果較人工目測、衛(wèi)星探測的結(jié)果更精確(Hodges，1998)。Long et al．(2001)使用獨(dú)立的分析方法將TSI儀器的觀測結(jié)果與WSI儀器(whole sky imager)觀測結(jié)果進(jìn)行對(duì)比分析，發(fā)現(xiàn)兩者相關(guān)性較好。

此外，Sabburg and Long(2004)和 Sabburg and Parisi(2008)使用TSI-440儀器觀測云量，并根據(jù)云量的變化，分析云對(duì)紫外輻射增強(qiáng)的作用。Calbo and Sabburg(2008)采用TSI以及全天空相機(jī)同時(shí)拍攝天空，對(duì)兩張圖片中的同一塊云體采用統(tǒng)計(jì)分析、傅里葉變換等獲取圖像信息，對(duì)晴空、波狀云、卷云、層狀云和積狀云這五種情況進(jìn)行識(shí)別，準(zhǔn)確率達(dá)到了76%。Kassianov et al．(2005)等使用 TSI和 HSI兩臺(tái)儀器，根據(jù)雙基站測量云高的原理進(jìn)行試驗(yàn)分析，得到云高計(jì)算公式，并且將分析結(jié)果與微脈沖激光雷達(dá)實(shí)測云高進(jìn)行對(duì)比，結(jié)果表明:TSI儀器對(duì)于平均云底高度低于3.5 km的單層云體的云高測量是可行的。按照地面觀測自動(dòng)化的發(fā)展要求，國內(nèi)也逐步開展地基測云儀的利用和研制，霍娟和呂達(dá)仁(2006)利用地基測云儀分析云量對(duì)輻射變化的影響，并以此反演氣溶膠的變化特征(霍娟和呂仁達(dá)，2005;霍娟，2007)。傅德勝和王新芝(1995)、孫學(xué)金等(2009a，2009b)、孫曉剛等(2008)使用紅外測云儀提取云圖特征，進(jìn)行云狀分類。高太長等(2010)、秦健等(2011)總結(jié)了目前地基測云儀的研究現(xiàn)狀，詳細(xì)分析了全天空云圖的獲取、鏡頭保護(hù)、云點(diǎn)識(shí)別、云量計(jì)算等多項(xiàng)測云技術(shù)問題。

雖然，TSI測云儀已經(jīng)得到了廣泛應(yīng)用，計(jì)算結(jié)果也較為精確，但是在不同的天空情況下，仍然存在一定的云量計(jì)算誤差。Pfister and McKenzie(2003)通過分析TSI云量和地表輻射量變化，指出:如果太陽被薄云遮掩但依然清晰可見時(shí)，很可能出現(xiàn)輻射增強(qiáng)的現(xiàn)象;此外，在觀測過程中，該儀器對(duì)太陽周邊云量的探測以及太陽是否被云體遮掩的判斷準(zhǔn)確度還需要進(jìn)一步提高。

綜上所述，TSI的云量計(jì)算存在一定誤差，還需進(jìn)一步研究。本文主要利用太湖地區(qū)2008年5—10月的TSI-440觀測資料展開云量計(jì)算誤差的分析。首先分析不同能見度條件下天空?qǐng)D像的成像特征，然后分析陰天圖像、復(fù)雜天空(多云)圖像的紅藍(lán)比值分布情況及云量計(jì)算誤差，并對(duì)減小誤差的方法進(jìn)行了初步研究。

1 TSI基本工作原理及資料來源

TSI儀器主要的部件有:照相機(jī)、帶有加熱裝置的半球鏡面、鏡面上方的遮光帶以及下方的電子設(shè)備系統(tǒng)(圖1)。該儀器可以在白天自動(dòng)進(jìn)行全天空云量的持續(xù)監(jiān)測，其工作原理是:通過儀器上方的照相機(jī)垂直向下拍攝帶有加熱裝置的半球鏡面，得到當(dāng)時(shí)天空所呈現(xiàn)的圖像，并將拍攝得到的圖像自動(dòng)存儲(chǔ)到計(jì)算機(jī)上進(jìn)行云量計(jì)算和處理。

大氣、云粒子對(duì)可見光具有不同的散射原理:大氣分子對(duì)可見光的散射與λ4成反比，如果天空為晴空，對(duì)藍(lán)光波段的散射遠(yuǎn)遠(yuǎn)大于對(duì)紅光波段的散射，因此晴空呈現(xiàn)藍(lán)色;而云粒子對(duì)可見光的散射，在不同的波段散射的程度是相當(dāng)?shù)?，因此云體呈現(xiàn)白色。根據(jù)不同的散射原理，可以從圖像信息值(RGB值)中區(qū)分云與藍(lán)天。自然界中的可見光都是按照一定的比例混合而得到，任意一種顏色都可以用紅(R)、綠(G)、藍(lán)(B)三色混合而成，圖像中的每個(gè)像素點(diǎn)均含有RGB三個(gè)亮度值，而RGB值在一定程度上反映了紅(中心波段為650 nm)、綠(中心波段為570 nm)、藍(lán)(中心波段為450 nm)這三個(gè)波段的輻射強(qiáng)度。因此，首先提取圖像中的信息值(RGB值)，計(jì)算每個(gè)像素點(diǎn)的紅藍(lán)比值(定義為:R/B)，然后設(shè)定合理的閾值來判別像元點(diǎn)是否為云點(diǎn)，這種方法被稱為閾值識(shí)別法(Morris，2005)。

圖1 全天空成像儀TSI-440(a)及其組成部分(b)Fig．1 (a)Total sky imager 440 and(b)its components

本文所用資料取自南京信息工程大學(xué)太湖觀測站(120.31°E，31.58°N)。該站自 2008 年 4 月—2009年4月在太湖地區(qū)進(jìn)行自動(dòng)觀測，具備多種氣象觀測儀器，如:全天空成像儀(TSI-440)、微脈沖激光雷達(dá)、微波輻射計(jì)、多通道旋轉(zhuǎn)光度計(jì)等，為云量及輻射的研究提供了全面的數(shù)據(jù)資料。本文為了研究TSI云量計(jì)算誤差，采用了2008年5—10月共計(jì)6個(gè)月的TSI資料以及無錫站的常規(guī)觀測資料。TSI資料包括:原始天空?qǐng)D像、處理后圖像以及相應(yīng)的云量計(jì)算結(jié)果。圖像分辨率為352×288，是24位全彩色圖像，圖像采集率為1 min/張，所拍攝的天空范圍是可視張角160°，其工作時(shí)間平均為07:00—17:00(北京時(shí)間，下同)，儀器的具體工作時(shí)間與太陽高度角有關(guān)，當(dāng)太陽高度角大于3°時(shí)，儀器自動(dòng)開始獲取圖片。

2 TSI資料分析

通常人工觀測時(shí)將天空劃分為十等分來估算云量大小，文中的云量均按照0～1計(jì)算。根據(jù)分析圖像特征發(fā)現(xiàn):當(dāng)天空完全被掩蔽，但云體仍顯明亮?xí)r，這種云較高，如:厚薄不均的蔽光高層云或透光高層云。反之，若云呈現(xiàn)灰白或灰黑色，天空陰沉，這種云則較低，如:蔽光層積云、雨層云等。文中所定義的薄云是高云(云高大于5 km)且云層較薄，例如毛卷云、鉤卷云、卷層云等，在圖像中此類云體與非云大氣的霧霾粒子難以區(qū)分，是云量計(jì)算中的難點(diǎn)。而不透光云體是指云層較厚，可以完全遮擋太陽直射光的中低層云，成像時(shí)可能出現(xiàn)陰影區(qū)域。

2.1 不同能見度條件下的圖像分析

一般情況下，當(dāng)氣溶膠、霧霾不嚴(yán)重時(shí)，人眼所見的晴空應(yīng)為藍(lán)色，同樣相機(jī)成像也是如此。但是，如果能見度較低，天空將呈現(xiàn)灰白色。圖2是4種不同能見度條件下的晴空?qǐng)D像及儀器自帶的處理結(jié)果，可以看出:隨著能見度的減小，原始天空?qǐng)D像的晴空區(qū)域逐漸由深藍(lán)色變?yōu)樗{(lán)色進(jìn)而變?yōu)榛野咨?，而儀器自帶的處理結(jié)果中云量區(qū)域隨著能見度的減小而增大。主要是由于氣溶膠前向散射作用使得太陽周邊像素點(diǎn)飽和，能見度越小，受太陽影響的區(qū)域越大。

能見度是影響大氣成像的重要因素，圖2可以說明大氣成像與能見度密切相關(guān)。能見度不同直接影響成像效果及后期的分析處理工作，因此首先分析不同能見度條件下圖像的成像特征。圖3中給出了3種不同能見度條件下晴空?qǐng)D像紅藍(lán)比值的分布情況，可以看出:隨著能見度的減小，紅藍(lán)比值有明顯增大的趨勢，同時(shí)，比值為1(R=B=255的飽和像素點(diǎn))的像元數(shù)隨著能見度的減小而增加?；谀芤姸葘?duì)天空成像的影響，本文中所選擇的天空?qǐng)D像均為能見度大于5 km的圖像。

同時(shí)，從晴空?qǐng)D像中也可以看出:在晴空天氣情況下，由于氣溶膠前向散射作用，會(huì)產(chǎn)生一定的云量計(jì)算誤差。針對(duì)晴空?qǐng)D像所產(chǎn)生的云量誤差，Long(2010)提出了云量改進(jìn)方案。該方案是根據(jù)天空其他區(qū)域的云量分布情況及變化特征來估算太陽圈及近地平區(qū)域的云量。經(jīng)與目測結(jié)果對(duì)比，發(fā)現(xiàn)該方案可以有效地減小晴空天氣情況下全天空成像儀的云量計(jì)算誤差。

圖225 km(a、b)、15 km(c、d)、8 km(e、f)、3 km(g、h)能見度條件下的晴空?qǐng)D像(a、e、c、g)及儀器自帶(b、f、d、h)的處理結(jié)果(藍(lán)色、白色和灰色區(qū)域分別代表晴空、不透光云和薄云;黑色條帶表示遮光帶，遮光帶上圓點(diǎn)表示太陽所在位置，其中黃點(diǎn)代表太陽未被云體掩蔽，白點(diǎn)表示太陽被云體掩蔽)Fig．2 (a，e，c，g)The clear sky images and(b，f，d，h)the corresponding TSI processed images at different visibilities(In the processed images，the blue regions represent cloudless and the write and gray regions stand for opaque cloud and thin cloud respectively．The black strip is shading belt．The dot on the shading belt denotes the position of the sun．The yellow dot represents sunny days without cloud and the white dot denotes cloudy days) a，b．25 km;c，d．15 km;e，f．8 km;g，h．3 km

圖3 不同能見度條件下晴空?qǐng)D像的紅藍(lán)比值分布情況Fig．3 The red blue ratio distributions of the clear sky images at different visibilities

2.2 陰天圖像紅藍(lán)比值分析及云量計(jì)算

不同天空狀況，云的成像特征并不相同，例如:陰天圖像比晴空?qǐng)D像明顯偏暗。這主要是因?yàn)殛幪炖锾栔鄙涔獗辉仆耆趽?，圖像的亮度整體偏小。TSI儀器在進(jìn)行圖像處理時(shí)，對(duì)所有的圖像均采用統(tǒng)一閾值，一般情況下，紅藍(lán)比值小于0.7的為晴空點(diǎn)，比值大于0.8的為云點(diǎn)，介于0.7～0.8之間的為薄云點(diǎn)，此閾值對(duì)晴空?qǐng)D像或亮度較高的圖像比較適用，而對(duì)于亮度偏暗的圖像會(huì)造成一定的云量計(jì)算誤差。圖4中是2008年10月4日13:10的陰天圖像及儀器自帶的處理結(jié)果，地面觀測資料顯示當(dāng)天全天云量均為10成，即云量為1。而儀器的結(jié)果為不透光云量0.332，薄云量0.441，晴空區(qū)域0.227，這與地面觀測資料及目測識(shí)別圖像的結(jié)果有較大差異。

文中選擇了500幅陰天圖像做RGB值直方圖統(tǒng)計(jì)分析，RGB值取值范圍在0～255之間，將其分為51檔(以5為間隔單位)，統(tǒng)計(jì)每幅圖像天空范圍像素點(diǎn)的RGB值分布情況。結(jié)果發(fā)現(xiàn):圖像中的紅色像元值(R)多集中于125～240之間，綠色像元值(G)多集中于150～245之間，藍(lán)色像元值(B)多集中于185～255之間。而晴空?qǐng)D像的R值多集中于80～200之間，G值多集中于150～220之間，B值多集中于230～255之間，太陽周邊的像素點(diǎn)均為R=B=G=255的飽和像素點(diǎn)。由此可見:陰天圖像的R值比晴空?qǐng)D像的R值偏大，而藍(lán)色像元值B值相對(duì)較小，且陰天圖像中的飽和像素點(diǎn)較少，說明太陽完全被云體遮掩。

圖5是圖4陰天圖像對(duì)應(yīng)的紅藍(lán)比值分布情況，可以看出:此時(shí)像元點(diǎn)紅藍(lán)比值多集中于0.65～0.9之間，且在紅藍(lán)比值為1處并無像元點(diǎn)，即無飽和像素點(diǎn)。說明此時(shí)太陽已被完全遮掩。容易看出:若對(duì)該圖采用儀器的統(tǒng)一閾值(0.7;0.8)，將會(huì)造成較大的云量計(jì)算誤差。

圖4 2008年10月4日13:10的陰天圖像(a)、儀器自帶程序的處理后圖像(b)及改進(jìn)閾值后的圖像(c)Fig．4 (a)The overcast sky image，(b)the corresponding TSI processed image and(c)the image after setting the new threshold at 13:10 BST on October 4，2008

圖5 2008年10月4日13:10的陰天圖像的紅藍(lán)比值的分布Fig．5 The red blue ratio distribution of overcast sky images at 13:10 BST on October 4，2008

選擇200幅陰天圖像進(jìn)行紅藍(lán)比值概率分布的分析，采用直方圖統(tǒng)計(jì)分析法，將天空像素點(diǎn)的紅藍(lán)比值范圍0.3～1.0分為35檔(以0.02為間隔單位)，首先對(duì)每幅陰天圖像計(jì)算每個(gè)單位檔上的像素點(diǎn)個(gè)數(shù)及概率分布值，然后將這200組數(shù)據(jù)每個(gè)單位檔做平均值計(jì)算，得到像素點(diǎn)的平均概率分布。同時(shí)，分別計(jì)算這200組數(shù)據(jù)與平均概率分布之間的標(biāo)準(zhǔn)偏差，結(jié)果發(fā)現(xiàn):所有數(shù)據(jù)的標(biāo)準(zhǔn)偏差均在0.03～0.06之間，標(biāo)準(zhǔn)偏差較小。說明像素點(diǎn)的平均概率分布可以反映陰天圖像紅藍(lán)比值的分布情況。此外，為了進(jìn)一步檢驗(yàn)平均概率分布的代表性，另選取擇了200幅陰天圖像，分析紅藍(lán)比值的分布情況，并計(jì)算每幅圖像的紅藍(lán)比值與平均概率分布的相關(guān)系數(shù)，共計(jì)200個(gè)相關(guān)系數(shù)。其中，22.7%的相關(guān)系數(shù)大于0.8，45.5%的相關(guān)系數(shù)介于0.7～0.8之間，還有27.2%的相關(guān)系數(shù)介于0.6～0.7之間，即相關(guān)系數(shù)大于0.6的數(shù)據(jù)占了整個(gè)樣本數(shù)量的95%。其余一部分?jǐn)?shù)據(jù)與平均概率分布的相關(guān)性不是非常好，通過分析這部分?jǐn)?shù)據(jù)的紅藍(lán)比值分布發(fā)現(xiàn):比值多數(shù)集中于0.8～1.0之間，觀察圖像資料發(fā)現(xiàn)這部分資料均為天空輻亮度較大的圖像資料。

圖6 陰天圖像紅藍(lán)比值的平均概率分布Fig．6 The mean probability distribution of red blue ratio on cloudy days

圖6給出了陰天圖像紅藍(lán)比值的平均概率分布，可以看出，陰天圖像比值集中于0.62～1.0之間，峰值位于0.7～0.8之間。因此，可以判定紅藍(lán)比值大于0.62的像素點(diǎn)均為云點(diǎn)，比值介于0.3～0.62的像素點(diǎn)均為晴空點(diǎn)，不透光云點(diǎn)閾值設(shè)定為0.66(即比值介于0.66～1.0像素點(diǎn)為云點(diǎn))，而0.62～0.66之間的像素點(diǎn)可能是薄云點(diǎn)，也可能是圖像中的陰影區(qū)域，暫時(shí)無法定論，需要根據(jù)不同的天空?qǐng)D像進(jìn)行分析。根據(jù)上述設(shè)定的閾值，重新計(jì)算后得到了圖4a中陰天圖像對(duì)應(yīng)的計(jì)算結(jié)果及改進(jìn)閾值之后的圖像(圖4c)，其計(jì)算結(jié)果為:不透光云量1，總云量1，與實(shí)際觀測結(jié)果一致。

由于陰天時(shí)輻亮度較小，使得圖像中陰影區(qū)域被誤識(shí)別為薄云或晴空，表1詳細(xì)地給出了多個(gè)陰天圖像的原始結(jié)果及改進(jìn)閾值之后的結(jié)果，可以看出:原始結(jié)果高估了薄云量，誤將不透光云量以薄云量計(jì)算，且總云量也存在一定的計(jì)算誤差;改進(jìn)閾值后得到的結(jié)果，薄云量的計(jì)算誤差減小，不透光云量及總云量更接近實(shí)際觀測資料?？梢姡槍?duì)陰天圖像重新設(shè)定閾值可以有效地減小因輻亮度不同而產(chǎn)生的薄云計(jì)算誤差，提高了云量計(jì)算的精度。

表1 2008年10月陰天圖像的原始云量結(jié)果及改進(jìn)閾值后的云量結(jié)果Table 1 The raw result and new result after setting the new threshold of cloud cover in the overcast sky image in October 2008

2.3 復(fù)雜天空?qǐng)D像分析

本文所指的復(fù)雜天空?qǐng)D像是云量大于0.2的多云圖像，在圖像分析過程中發(fā)現(xiàn)，有較多的圖像成像時(shí)出現(xiàn)陰影區(qū)域，這些區(qū)域在圖像處理過程中被作為薄云計(jì)算，造成了云量計(jì)算誤差，此類情況與陰天時(shí)出現(xiàn)誤差的情況類似。圖7a是一幅藍(lán)天白云分界清晰的天空?qǐng)D像，從儀器自帶的處理結(jié)果中(圖7b)可以看出:天空中云體較厚的部分(圖中灰色區(qū)域)在圖像處理過程中被作為薄云計(jì)算;而太陽未被云體遮掩，因此周邊像素點(diǎn)飽和的區(qū)域，被識(shí)別為不透光云。圖7c顯示了該圖像的紅藍(lán)比值分布情況，可以看出:圖中晴空區(qū)域比值最小(藍(lán)色區(qū)域);云層較高、較明亮的云體比值偏大，如太陽周邊云體(紅色區(qū)域);而云體較厚、較暗的部分紅藍(lán)比值相對(duì)較小，如圖中的黃色區(qū)域(對(duì)應(yīng)圖7b的灰色區(qū)域)。由于此類圖像的云量計(jì)算誤差與陰天情況類似，因此重新設(shè)定的閾值與陰天圖像相同，比值為0.3～0.62的像素點(diǎn)為晴空點(diǎn)，0.66～1.0的像素點(diǎn)為不透光云點(diǎn)，而介于0.62～0.66之間的像素點(diǎn)為薄云點(diǎn)。

圖7 原始天空?qǐng)D像(a)、儀器自帶的處理圖像(b)及原始天空?qǐng)D像與儀器自帶的處理圖像比值(c)Fig．7 (a)The raw image，(b)the corresponding TSI processed image and(c)the red blue ratio distribution image

圖8是2008年10月16日14:07的天空?qǐng)D像、儀器自帶的處理圖像及改進(jìn)閾值后得到的圖像。從原始圖像中，可以清楚地分辨出紅色框內(nèi)區(qū)域?yàn)椴煌腹庠企w，但由于云層較低，云體較厚，呈現(xiàn)灰色，在原始處理過程中被作為薄云計(jì)算，原始計(jì)算結(jié)果為:薄云量0.267，不透光云量0.239，總云量為0.506。通過目測識(shí)別可知紅色框內(nèi)區(qū)域?yàn)椴煌腹庠企w，被誤識(shí)別為薄云體。利用新閾值(晴空0.62，云0.66)對(duì)圖像處理后，得的計(jì)算結(jié)果為:薄云量0.042，不透光云量為0.602。地面觀測資料顯示:總云量為6成，即0.6，低云量為0.6，更改閾值之后的計(jì)算結(jié)果與實(shí)際觀測資料更為接近。由此可見，采用新閾值減小了儀器薄云量的計(jì)算誤差，同時(shí)也提高了總云量的計(jì)算精度。

圖82008年10月16日14:07的原始天空?qǐng)D像(a)、處理后圖像(b)和改進(jìn)閾值后的圖像(c)Fig．8 (a)The raw image，(b)the corresponding processed image and(c)the simulated image after setting the new threshold at 14:07 BST on October 16，2008

結(jié)合太湖觀測站的TSI資料及臺(tái)站觀測資料，選出每日08:00、11:00及14:00的觀測資料共計(jì)450組數(shù)據(jù)，將不同天氣下的云量進(jìn)行誤差統(tǒng)計(jì)分析。由于本文所采用的新閾值是對(duì)多云及陰天情況，因此對(duì)云量為0.3～1.0的8種情況進(jìn)行誤差計(jì)算。圖9給出了不同情況下儀器自帶處理結(jié)果、改進(jìn)閾值后的結(jié)果分別與目測云量的絕對(duì)誤差情況，可以看出，儀器自帶的處理結(jié)果在云量大于等于0.7時(shí)，絕對(duì)誤差大于0.1，其中陰天時(shí)誤差達(dá)到最大。改進(jìn)閾值之后得到的計(jì)算結(jié)果在多云及陰天時(shí)，絕對(duì)誤差均小于0.08，說明重新設(shè)定閾值可以有效地減小多云及陰天情況下的云量計(jì)算誤差。

3 結(jié)論與討論

利用全天空成像儀代替人工觀測云量，在多數(shù)情況下效果良好。但是，由于天空成像特征與天氣變化情況、能見度等密切相關(guān)，如隨著能見度的減小，晴空?qǐng)D像的紅藍(lán)比值逐漸增大，所以導(dǎo)致全天空成像儀自帶程序在復(fù)雜天空狀況下云量計(jì)算出現(xiàn)較大誤差。

圖9 不同天氣情況下的云量絕對(duì)誤差Fig．9 The absolute differences under different weather conditions

在分析太湖觀測站TSI資料的過程中發(fā)現(xiàn):該儀器在晴空、陰天及云體較厚的情況下易產(chǎn)生云量計(jì)算誤差。針對(duì)陰天情況，文中對(duì)多幅天空?qǐng)D像進(jìn)行紅藍(lán)比值統(tǒng)計(jì)分析，結(jié)果發(fā)現(xiàn):比值主要集中于0.62～1.0之間。根據(jù)圖像的成像特征，將晴空點(diǎn)閾值設(shè)定為0.62(＜0.62判為晴空)，云點(diǎn)閾值設(shè)定為0.66(＞0.66判為云)，采用新閾值得到的云量結(jié)果與常規(guī)氣象觀測結(jié)果更為接近。同時(shí)，對(duì)天空?qǐng)D像中存在陰影區(qū)域的其他圖像也使用新閾值，得到的云量結(jié)果比原結(jié)果更準(zhǔn)確，減小了因輻亮度不同而造成的云量計(jì)算誤差。

在圖像分析過程中，還發(fā)現(xiàn)圖像的紅藍(lán)比值分布情況與太陽高度角、太陽方位角也密切相關(guān)。因此，要解決以上這些問題，應(yīng)根據(jù)不同的天氣情況、像素點(diǎn)的位置、太陽中心點(diǎn)位置等信息設(shè)定不同的閾值，建立一套完整的云點(diǎn)閾值識(shí)別法。另外，也可以嘗試運(yùn)用其他計(jì)算途徑來獲取云量信息，例如二維傅里葉變換、統(tǒng)計(jì)分析法或圖像處理技術(shù)中的邊緣提取法等進(jìn)行云量的計(jì)算，進(jìn)一步改善計(jì)算結(jié)果。

陳洪濱．1997．關(guān)于云和云天大氣對(duì)太陽輻射的吸收異常［J］．大氣科學(xué)，21(6):750-757．

陳渭民，邊多，郁凡．1997．由衛(wèi)星資料估算晴空太陽直接輻射和散射輻射［J］．南京氣象學(xué)院學(xué)報(bào)，20(3):326-333．

傅德勝，王新芝．1995．云圖紋理特征的抽取與云的自動(dòng)分類［J］．南京氣象學(xué)院學(xué)報(bào)，18(4):530-535．

高太長，劉磊，趙世軍，等．2010．全天空測云技術(shù)現(xiàn)狀及發(fā)展［J］．應(yīng)用氣象學(xué)報(bào)，21(1):100-109．

霍娟．2007．地基全天空成像系統(tǒng)云與氣溶反演及其應(yīng)用研究［D］．北京:中國科學(xué)研究院．

霍娟，呂仁達(dá)．2005．利用全天空數(shù)字圖像對(duì)北京上空云況分布特征的試驗(yàn)分析［J］．氣象科學(xué)，25(3):238-243．

霍娟，呂達(dá)仁．2006．晴空與有云大氣輻射分布的數(shù)值模擬及其對(duì)全天空?qǐng)D像云識(shí)別的應(yīng)用［J］．氣象學(xué)報(bào)，64(1):31-38．

秦健，劉磊，高太長，等．2011．基于統(tǒng)計(jì)方法的地基器測云天識(shí)別［J］．氣象科學(xué)，31(5):604-612．

宋慶利．2003．利用衛(wèi)星資料研究云對(duì)地面凈輻射的影響［D］．南京:南京氣象學(xué)院．

孫曉剛，孫學(xué)金，牛珍聰，等．2008．全天空云圖獲取的一種方式及算法實(shí)現(xiàn)［J］．氣象科學(xué)，28(3):338-341．

孫學(xué)金，劉磊，高太長，等．2009a．基于模糊紋理光譜的全天空紅外圖像云分類［J］．應(yīng)用氣象學(xué)報(bào)，20(2):157-163．

孫學(xué)金，劉磊，高太長，等．2009b．基于LBP算法的全天空紅外云圖的分類研究［J］．大氣科學(xué)學(xué)報(bào)，32(4):490-497．

Calbo J，Pages D．2002．Empirical studies of cloud effects on UV radiation:A Review［J］．J Geophys Res，43(2):1-28．

Calbo J，Sabburg J M．2008．Feature extraction from whole sky groundbased images for cloud-type recognition ［J］．Amer Meteor Soc，25:3-14．

Hodges G．1998．A Comparison of fractional cloud cover determination by a sky image with visual estimation by trained national weather service observers［EB/OL］．［2012-10-10］．http://www．noaa．gov/．

Kassianov E，Long C N，Christy J．2005．Cloud-base-height estimation from paired ground-based hemispherical observations［J］．J Appl Meteor，44(8):1221-1233．

Long C N．2004．The shortwave clear-sky detection and fitting algorithm:Algorithm operational details and explanations［R］//Atmospheric radiation measurement program．Tech Rep ARM TR-004．

Long C N．2010．Correcting for circumsolar and near-horizon errors in sky cover retrievals from sky images［J］．J Atmos Soc，4(1):45-52．

Long C N，Deluisi J J．1998．Development of an automated hemispheric sky imager for cloud fraction retrievals［C］//10th symposium on meteorological observations and instrumentation．Arizona:Amer Meteor Soc:171-174．

Long C N，Thomas P．2000．Identification of clear skies from broadband dynamometer measurements and calculation of down-welling shortwave cloud effects［J］．J Geophys Res，105(12):15609-15626．

Long C N，Slater D W，Tooman T．2001．Total sky imager(TSI)model 880 status and testing results［R］//Atmospheric radiation measurement program．Tech Rep ARM TR-004．

Martins F R，Souza M P，Pereira E B．2003．Comparative study of satellite and ground techniques for cloud cover determination［J］．Space Res，32(11):2275-2280．

Morris V R．2005．Total sky imager hand book［EB/OL］．［2012-10-10］．http://www．a(chǎn)rm．gov/．

Pfister G，McKenzie R L．2003．Cloud coverage based on all-sky imaging and its impact on surface solar irradiance［J］．Amer Meteor Soc，42:1421-1434．

Sabburg J M，Long C N．2004．Improved sky imaging for studies of enhanced UV irradiance［J］．Atmos Chem Phys，4:2543-2552．

Sabburg J M，Parisi A V．2008．Cloud observation for the statistical evaluation of the UV index at Toowoomba，Australia［J］．Int J Biometeorol，52:159-166．