1 概述 
　　每年春季，我國(guó)北方大部分地區(qū)都會(huì)受到不同程度的沙塵天氣的影響，一方面，沙塵引起的能見(jiàn)度降低與大顆粒沙塵的沉降嚴(yán)重影響了運(yùn)輸業(yè)、農(nóng)業(yè)等行業(yè)的正常發(fā)展，對(duì)于國(guó)民經(jīng)濟(jì)造成極大的直接和間接的損失；另一方面，懸浮在空氣中粒徑小于100μm的總懸浮顆粒物會(huì)對(duì)人們的健康帶來(lái)不利的影響。 
　　國(guó)內(nèi)對(duì)沙塵天氣的檢測(cè)當(dāng)前以遙感為主、地面站點(diǎn)檢測(cè)為輔。比較而言，使用衛(wèi)星傳感器獲取的遙感影像具有檢測(cè)范圍廣且時(shí)效性強(qiáng)的特點(diǎn)，以搭載MODIS的Aqua和Terra衛(wèi)星為例，它們可以以1～2天為周期獲取數(shù)據(jù)，這樣便可以及時(shí)發(fā)現(xiàn)沙塵源地并將其規(guī)模和運(yùn)動(dòng)方向標(biāo)記出來(lái)，達(dá)到動(dòng)態(tài)、省時(shí)省力檢測(cè)大范圍沙塵天氣的目的。隨著遙感技術(shù)的發(fā)展，衛(wèi)星傳感器更加廣泛地應(yīng)用于環(huán)境管理與污染檢測(cè)等方向。 
　　2 數(shù)據(jù)源 
　　2.1 MODIS數(shù)據(jù) 
　　MODIS是Terra與Aqua上搭載的十分重要的地球觀測(cè)傳感器，兩顆衛(wèi)星上的MODIS傳感器可以每1～2天對(duì)地球觀測(cè)一次，同時(shí)獲取36個(gè)波段（0.4～14.3μm）的數(shù)據(jù)，包括16個(gè)熱紅外波段以及20個(gè)可見(jiàn)光與近紅外波段，這些數(shù)據(jù)常用于反演水汽含量、地表溫度、氣溶膠參數(shù)等，可以幫助研究人員理解陸地、海洋以及低層大氣的動(dòng)態(tài)變化進(jìn)程，并有助于構(gòu)建可以預(yù)測(cè)變化的地球系統(tǒng)模型。 
　　2.2 臭氧檢測(cè)儀（OMI）的數(shù)據(jù) 
　　2.2.1 Aura/OMI傳感器。Aura衛(wèi)星亦隸屬于EOS系統(tǒng)，該衛(wèi)星被設(shè)計(jì)用于獲取臭氧、氣溶膠和大氣中關(guān)鍵氣體的定量信息，因此，Aura搭載了四個(gè)相關(guān)的傳感器，分別是高分辨動(dòng)力發(fā)聲器、微波分叉發(fā)聲器、臭氧層觀測(cè)儀和對(duì)流層放射光譜儀。其中，OMI可以對(duì)264～504nm范圍的波段進(jìn)行觀測(cè)，其光譜分辨率在0.42～0.63nm，利用OMI獲取的數(shù)據(jù)可以對(duì)大氣中的NO2、SO2、BrO、HCHO以及氣溶膠進(jìn)行量測(cè)，同時(shí)使用OMI還可以區(qū)分不同的氣溶膠類(lèi)型，如沙塵和硫化物。 
　　2.2.2 紫外波段氣溶膠指數(shù)（UVAI）。大氣氣溶膠可以通過(guò)吸收作用對(duì)紫外輻射產(chǎn)生影響，研究表明，沙塵的增加使達(dá)到地面的紫外輻射量減少，而這個(gè)過(guò)程又可以對(duì)對(duì)流層中臭氧光化學(xué)反應(yīng)產(chǎn)生影響，因此可以使用臭氧觀測(cè)儀收集的數(shù)據(jù)衍生出反映氣溶膠變化的紫外波段氣溶膠指數(shù)，可以認(rèn)為UVAI提供了氣溶膠對(duì)紫外波段吸收作用的量化指標(biāo)。Torres等人的研究表明紫外觀測(cè)對(duì)氣溶膠具有高度的敏感性，并能有效降低地表反射率對(duì)觀測(cè)的影響。在實(shí)際應(yīng)用中，常常用于對(duì)煤煙或沙塵的檢測(cè)。 
　　在獲取UVAI時(shí)，將通過(guò)氣溶膠觀測(cè)到的光譜輻射率記為I360，而將無(wú)氣溶膠時(shí)觀測(cè)到的光譜輻射率記為I331，分別計(jì)算二者實(shí)測(cè)值與理論值的比值，并計(jì)算差值。使用等式可以表示為： 
　　3 算法原理與算法描述 
　　3.1 沙塵檢測(cè)算法原理以及亮度溫度的計(jì)算 
　　3.1.1 亮度溫度差值/比值算法。本文進(jìn)行亮度溫度差值和比值運(yùn)算時(shí)選擇MODIS 29、31和32波段，這三個(gè)波段分別處于8～9μm和10～12μm光譜區(qū)間，這樣選擇的原因在于氣溶膠顆粒在這一光譜區(qū)間內(nèi)其光譜特征會(huì)隨著波長(zhǎng)的改變而產(chǎn)生較大的波動(dòng)，通過(guò)對(duì)不同的影像進(jìn)行對(duì)比便可以有效識(shí)別不同的氣溶膠類(lèi)型。不僅如此，大氣對(duì)該光譜區(qū)間的干擾往往是小的，這樣得到的影像便可以反映氣溶膠的實(shí)際情況。因此，基于單一波段影像數(shù)據(jù)，科學(xué)家成功實(shí)現(xiàn)了對(duì)大氣中火山灰的觀測(cè)。而且隨著高光譜數(shù)據(jù)的應(yīng)用，研究人員利用窄波段光譜數(shù)據(jù)實(shí)現(xiàn)了對(duì)沙塵等氣溶膠更高精度的觀測(cè)。然而單一波段的影像不僅會(huì)受到地面輻射的干擾，同時(shí)氣體與氣溶膠還對(duì)輻射具有不可忽略的吸收作用?？梢岳脗鞲衅鹘邮盏降募t外波段的輻射亮度值公式對(duì)上述內(nèi)容進(jìn)行描述： 
　　式中：Bλ為普朗克方程；Ta表示大氣中氣溶膠層的溫度；Ts表示地表溫度；δλ為氣溶膠光學(xué)厚度。 
　　從等式中可以看出，即使當(dāng)氣溶膠層足夠厚，使地面輻射的干擾可以忽略不計(jì)，即 ，所接收到的輻射亮度依然受到氣體與氣溶膠本身的影響，而正如前文所述，在8～9μm和10～12μm內(nèi)，氣體對(duì)輻射的干擾會(huì)達(dá)到小，在這種情況下，就只需要考慮氣溶膠本身的吸收作用。 
　　3.1.2 Dust*-parameter（D*）算法。Dust*-parameter（D*）算法是基于上文論述的BTD11-12和BTD8.6-11兩種差值算法得到的，其表達(dá)式為： 
　　式中：C和E分別為BTD11-12、BTD8.6-11的溫度位移，本文取值分別為0.09與15.0。 
　　通過(guò)前文所述，可以得知BTD11-12與BTD8.6-11均可以對(duì)以石英為主要成分的沙塵進(jìn)行有效觀測(cè)，除此以外，BTD8.6-11還對(duì)硅質(zhì)礦物非常敏感，因?yàn)楣栀|(zhì)礦物擁有較寬的殘余輻射帶，相較于11μm，可以吸收更多在8.6μm處的輻射，造成BTD8.6-11的值為負(fù)值；而云層與此相反，因?yàn)樵茖訒?huì)更多地吸收11μm處的輻射，這樣BTD8.6-11值即為正值?；诖?，如果沙塵中含有硅質(zhì)礦物，D*算法就可以有效進(jìn)行識(shí)別，同時(shí)取自然指數(shù)可以對(duì)沙塵信號(hào)進(jìn)行放大。 
　　3.1.3 MODIS亮度溫度計(jì)算方法。本文使用普朗克函數(shù)對(duì)亮度溫度進(jìn)行求解，在進(jìn)行計(jì)算之前，需要將MODIS影像的像元值轉(zhuǎn)換為對(duì)應(yīng)的熱輻射強(qiáng)度值，衛(wèi)星接收到的熱輻射強(qiáng)度轉(zhuǎn)換方法如下：sr-1·μm-1；NDi為第i波段影像的像元值；DRi與DRSi則是對(duì)應(yīng)波段的輻射常量。 
　　接下來(lái)，利用普朗克函數(shù)對(duì)各波段亮度溫度計(jì)算，亮度溫度Ti計(jì)算公式為： 
　　式中：λi是對(duì)應(yīng)波段的有效波長(zhǎng)，在本文中使用的是29、31和32波段的有效波長(zhǎng)分別為8.6μm、11μm和12μm。C1和C2均為光譜常量，它們的值分別是C1=119104356W·m-2·sr-1·μm4、C2=1.4387685×104μm·K。在計(jì)算時(shí)為簡(jiǎn)便直觀可以將C2/λi與C1/λi5看作一個(gè)整體。 
　　實(shí)際上，亮度溫度是描述一般地物的“等效”溫度參數(shù)即在一定的波段范圍內(nèi)，這一參數(shù)指的是當(dāng)一般地物與黑體具有相等的輻射亮度時(shí)，黑體所具有的溫度。它表示地物的輻射量與具有多大溫度的黑體的輻射相當(dāng)，由于一般地物的比輻射率ε均小于1，故地物的亮度溫度往往小于它對(duì)應(yīng)的實(shí)際溫度。從影像中可以看到計(jì)算得到的亮度溫度的取值在220K和300K之間，比同季節(jié)正常的地表溫度低一些，根據(jù)上述結(jié)論，這個(gè)計(jì)算結(jié)果是合理的。 
　　3.2 云層的處理 
　　3.2.1 決定系數(shù)（R2）去云。通過(guò)計(jì)算BT7.3與BT11間的決定系數(shù)來(lái)確定高緯度云層的地理位置從而將其去除，當(dāng)決定系數(shù)趨近于1，說(shuō)明觀測(cè)區(qū)域是云層的可能性極大；當(dāng)這個(gè)數(shù)值遠(yuǎn)離1，則該區(qū)域應(yīng)為無(wú)云區(qū)域。這個(gè)方法的關(guān)鍵在于尋找一個(gè)合適的閾值來(lái)界定云層與非云層區(qū)域，如果閾值選擇偏小，可能會(huì)將有用的信息劃為云層而遭到掩蓋；如果偏大，則有部分高緯度云層會(huì)有所遺漏。在本實(shí)驗(yàn)中，首先選取同一地區(qū)有云與無(wú)云時(shí)刻的兩幅影像，并按上述方法計(jì)算其28、31波段的亮度溫度與決定系數(shù)，根據(jù)決定系數(shù)影像上像元值的直方圖分布決定該閾值的取值。 
　　3.2.2 BT11去云。31波段被設(shè)計(jì)用于獲取云層或地表溫度，而云層的溫度往往低于后者。因此從11μm波段的亮度溫度影像上便可以對(duì)云層覆蓋區(qū)域進(jìn)行判別。本文首先統(tǒng)計(jì)出BT11影像上的大值Tmax，為了排除噪聲對(duì)統(tǒng)計(jì)結(jié)果的影響，選擇大的五個(gè)像元值并計(jì)算其平均數(shù)作為T(mén)max，同時(shí)根據(jù)相關(guān)文獻(xiàn)，可以獲知云層在11μm處的大亮度溫度約為273K，而對(duì)應(yīng)的地表溫度為288.2K，其比值為0.95，因此本文將溫度大值的0.95倍作為云層的閾值溫度，即T云=0.95×Tmax。那么，BT11影像上像元值大于T云且小于Tmax的區(qū)域可作為無(wú)云區(qū)域；同樣的，當(dāng)像元值小于T云，則相應(yīng)區(qū)域可判定為有云區(qū)域。 
　　4 結(jié)語(yǔ) 
　　本文主要依據(jù)紅外波段對(duì)沙塵粒徑、光學(xué)折射系數(shù)等屬性的敏感性展開(kāi)沙塵檢測(cè)算法的研究，在構(gòu)建算法時(shí)，主要考慮了云層對(duì)沙塵數(shù)據(jù)的掩蓋以及地面輻射、大氣吸收散射作用等因素對(duì)沙塵信息的干擾。 
　　云層不僅會(huì)對(duì)沙塵數(shù)據(jù)掩蓋，同時(shí)細(xì)小的卷云有可能會(huì)與提取的沙塵信息混淆造成精度下降，因此本文通過(guò)雙重云層掩膜算法對(duì)云層去除，介紹了亮度溫度差值算法（BTD）、亮度溫度比值算法（RBT）和Dust*算法的沙塵檢測(cè)算法，考慮到云層的影響，又介紹了決定系數(shù)（R2）去云和BT11去云方法，終獲得沙塵檢測(cè)結(jié)果。本文提到的檢測(cè)方法可以實(shí)現(xiàn)一定精度的沙塵檢測(cè)，但是隨著沙塵越來(lái)越嚴(yán)重，對(duì)沙塵的檢測(cè)及產(chǎn)生原因的研究還需要更加深入。