于一凡，潘 軍，邢立新，蔣立軍，孟 濤，韓曉靜，周彩彩

(吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026)

短波紅外波段高溫目標(biāo)識別的可行性分析

于一凡，潘 軍，邢立新，蔣立軍，孟 濤，韓曉靜，周彩彩

(吉林大學(xué)地球探測科學(xué)與技術(shù)學(xué)院，長春 130026)

目前主要采用熱紅外遙感數(shù)據(jù)識別高溫目標(biāo)，研究發(fā)現(xiàn)短波紅外數(shù)據(jù)也具有較好的高溫目標(biāo)識別能力。在高溫目標(biāo)混合像元中，高溫地物的溫度和面積是待求解的未知數(shù)，也是決定混合像元波譜特性的重要參數(shù)。首先基于能量守恒原理建立了地表高溫目標(biāo)混合像元的輻射能量方程，然后進(jìn)行方程參數(shù)的敏感性分析。研究結(jié)果表明，高溫目標(biāo)面積百分比和常溫地物反射率對高溫地物的溫度和面積反演最為敏感。以山西保德和陜西府谷研究區(qū)ETM+第7波段的數(shù)據(jù)為例進(jìn)行了高溫目標(biāo)識別研究，所識別的火點(diǎn)像元輻射通量密度為背景的1.36～4.76倍。經(jīng)野外驗(yàn)證，用馬氏距離法識別高溫目標(biāo)的精度達(dá)到88%，表明短波紅外波段的遙感數(shù)據(jù)可用于高溫目標(biāo)的識別。

短波紅外;高溫目標(biāo);常溫地物;可行性分析

0 引言

在高溫目標(biāo)遙感識別研究中多采用中、熱紅外數(shù)據(jù)［1－2］，但受其空間分辨率的限制，難以實(shí)現(xiàn)對局地高溫目標(biāo)地物異常的識別。與中、熱紅外數(shù)據(jù)相比，短波紅外數(shù)據(jù)也具有高溫目標(biāo)識別優(yōu)勢［3］，而且在高溫目標(biāo)識別中具有更高的精度［4］。近年來，一些學(xué)者對基于短波紅外和中、熱紅外數(shù)據(jù)的火點(diǎn)識別方法進(jìn)行了比較研究［5－6］。因此，對短波紅外數(shù)據(jù)高溫目標(biāo)識別的可能性進(jìn)行研究就更有理論意義和實(shí)用價(jià)值。在遙感圖像中，高溫與常溫地物的混合像元輻射能量主要包括4部分，即混合像元內(nèi)高溫目標(biāo)的發(fā)射輻射能量、高溫目標(biāo)的反射輻射能量、常溫地物的反射輻射能量和常溫地物的發(fā)射輻射能量?？梢詮纳鲜鲚椛涮匦灾謱Χ滩t外波段高溫目標(biāo)識別進(jìn)行可行性分析。本文通過對高溫目標(biāo)混合像元輻射能量方程及相關(guān)參數(shù)進(jìn)行分析得到各參量的敏感性，在此基礎(chǔ)上說明短波紅外數(shù)據(jù)可以用于高溫目標(biāo)的識別。

1 地物輻射特性及原理

1.1 高溫目標(biāo)混合像元輻射能量方程

利用短波紅外數(shù)據(jù)對高溫目標(biāo)混合像元的能量組成進(jìn)行分析，確定混合像元的輻射能量包含常溫地物和高溫目標(biāo)的反射與發(fā)射能量(圖1)，即

式中:M為混合像元的輻射通量密度;M1為混合像元內(nèi)高溫目標(biāo)的發(fā)射輻射通量密度;M2為混合像元內(nèi)常溫地物的反射輻射通量密度;M3為混合像元內(nèi)高溫目標(biāo)的反射輻射通量密度;M4為混合像元內(nèi)常溫地物的發(fā)射輻射通量密度;S為混合像元內(nèi)高溫目標(biāo)的面積百分比。

圖1 高溫目標(biāo)混合像元短波紅外輻射傳輸能量組成Fig.1 Short wave－infrared radiation transmission energy of high temperature target mixed pixel

1.2 高溫目標(biāo)混合像元輻射原理

據(jù)Stefan－Boltzmann定律和Wien’s位移定律可得，在混合像元中，地表高溫目標(biāo)的發(fā)射能量為高溫目標(biāo)的發(fā)射率與普朗克函數(shù)的乘積，即

式中:ε為發(fā)射率;λ為波長;B為溫度為T的黑體在波長為λ時(shí)的黑體普朗克函數(shù);h為普朗克常數(shù)，h=6.63×10－34Js;c為真空中的光速，c=3 ×108m/s;k為波爾茲曼常數(shù)，k=1.38 ×10－23J/K。

地表常溫地物的反射能量為

式中:Tθ為大氣透過率;E為穿過大氣后的太陽輻照度;E0為大氣上界太陽輻照度;θ為太陽天頂角;ρ為常溫地物的反射率;為日地天文單位距離。

地表高溫目標(biāo)的反射能量為

地表常溫地物的發(fā)射能量為

式中T'為常溫，T'=300 K。

地表混合像元輻射能量為

式中ρ0為混合像元的視反射率。

上述公式中的參數(shù)ε，Tθ，ρ，S和T均為未知數(shù)，因此需要對其進(jìn)行敏感性分析。

2 研究區(qū)概況

研究區(qū)位于山西省保德縣與陜西省府谷縣的交界處，E110°22'～ 111°19'，N38°39'～ 39°35'之間。保德和府谷分別位于黃河的兩岸，區(qū)內(nèi)的石炭系、二疊系和三疊系地層中均含有煤層，煤炭資源較為豐富。

2002年前后，該區(qū)廣泛分布著土法煉焦?fàn)t、金屬冶煉廠等多種高溫目標(biāo)，因此選擇2002年7月14日獲取的Landsat7 ETM+數(shù)據(jù)作為本文方法研究的基礎(chǔ)數(shù)據(jù)。對ETM+數(shù)據(jù)進(jìn)行了輻射定標(biāo)和大氣校正等一系列預(yù)處理，并通過裁剪得到研究區(qū)圖像(圖2)。

圖2 研究區(qū)ETM+7(R)4(G)2(B)假彩色合成圖像Fig.2 False color image of study area composed of ETM+7(R)4(G)2(B)

研究區(qū)內(nèi)有居民地、道路、林地、耕地(山地或低地)、河流、河漫灘和高溫目標(biāo)共7類地物。對這7類地物進(jìn)行樣本選取并檢驗(yàn)其訓(xùn)練樣本的J－M距離和散度，檢驗(yàn)結(jié)果表明，訓(xùn)練樣本的可分離性可滿足應(yīng)用要求。

太陽天頂角θ在ETM+圖像的頭文件中可以得到，θ=27.237839 6°; 日地天文單位距離可由下式得出，即

式中JD為儒略日(Julian day)，以儒略日計(jì)日是為便于計(jì)算年代相隔久遠(yuǎn)或不同歷法的2個(gè)事件所間隔的日數(shù)。經(jīng)查歷法表，2002年7月14日的JD=195，ds=1.016566 781，故=1.033408 02。

3 參數(shù)的敏感性分析

高溫目標(biāo)混合像元輻射能量方程中存在ε，Tθ，ρ和S等4個(gè)未確定參數(shù)，故需要對其與溫度T之間的敏感性做出分析。為便于進(jìn)一步對短波紅外波段中的高溫混合像元進(jìn)行敏感性分析，本文對上述4個(gè)參數(shù)各取一個(gè)適當(dāng)?shù)亩ㄖ?。各參量與溫度T間的關(guān)系為

式中:ρ0為混合像元視反射率，ρ0=0.454;E0(用于計(jì)算E)為從中國遙感衛(wèi)星地面站獲取的ETM+7波段太陽光譜輻射量［7］，E0=82.1 W·m－2·μm－1;λ為 ETM+7 波段標(biāo)稱波長，λ =2.208 μm。

近年來的研究發(fā)現(xiàn)，同一物體在不同溫度下的發(fā)射率變化不大［8］，故本文采用常溫的煤、炭和焦炭的發(fā)射率代替高溫目標(biāo)的發(fā)射率。本文使用近紅外礦物分析儀在野外不同位置采集的煤、炭和焦炭的光譜數(shù)據(jù)作為高溫目標(biāo)的反射率。光譜儀測量的波譜范圍為1300～2500 nm，波長掃描間隔為2 nm和4 nm。發(fā)射率等于1與反射率的差。研究中選取了39個(gè)樣品，測得的光譜反射率在ETM+7波段的均值為0.094437，因此高溫目標(biāo)的發(fā)射率為0.905563。

本文使用MODTRAN［9］來獲取大氣的平均透過率。根據(jù)研究區(qū)的具體情況，大氣廓線選擇中緯度夏季，波數(shù)為4529 cm－1(即ETM+7波段的波長)，得到Tθ=0.851231。對于常溫地物的反射率ρ，由于高溫地物多數(shù)在沖溝附近，所以為了更接近真實(shí)值，混合像元中常溫地物的反射率取沖溝反射率的平均值，得到其反射率為0.181373。研究區(qū)高溫目標(biāo)多為土法煉焦?fàn)t，煉焦?fàn)t的大小為10 m×3 m，每個(gè)煉焦廠中大約有5～7個(gè)煉焦?fàn)t，估算時(shí)取煉焦?fàn)t個(gè)數(shù)的均值6個(gè)。敏感性分析取ε，Tθ，ρ和S中任意3個(gè)參量為定值，剩下的一個(gè)變量在一定區(qū)間內(nèi)以適當(dāng)?shù)牟介L變動(dòng)，從而獲取該變量與溫度T之間的關(guān)系，得出如下結(jié)論:

1)根據(jù)所測光譜大致范圍，高溫目標(biāo)發(fā)射率ε以0.01的步長在0.8～1的范圍內(nèi)變動(dòng)，在此范圍內(nèi)變化會引起5 K的溫度浮動(dòng)(圖3)。

圖3 高溫目標(biāo)的發(fā)射率ε與T的關(guān)系Fig.3 Relationship between ε and T

2)通過調(diào)整MODTRAN的輸入?yún)?shù)，得到大氣透過率Tθ以0.005的步長在0.8～1的范圍內(nèi)的變動(dòng)，在此范圍內(nèi)變化會引起12 K的溫度浮動(dòng)(圖4)。

圖4 大氣透過率Tθ與T的關(guān)系Fig.4 Relationship between Tθand T

3)據(jù)Pan等［10］對ETM+7波段中典型地物的數(shù)據(jù)統(tǒng)計(jì)，耕地和居民地的反射率為0.08～0.18，黃土的反射率為0.11～0.34，裸地的反射率為0.19～0.26。根據(jù)研究區(qū)中的地物類別，常溫地物取反射率ρ以0.01的步長在0.1～0.3范圍內(nèi)變動(dòng)，在此范圍內(nèi)的變化會引起31 K的溫度浮動(dòng)(圖5)。

圖5 常溫地物的反射率ρ與T的關(guān)系Fig.5 Relationship between normal atmospheric temperature feature reflectiuity ρ and T

4)根據(jù)煉焦廠中煉焦?fàn)t的個(gè)數(shù)，高溫目標(biāo)面積比S以0.01的步長在0.1～0.3范圍內(nèi)變動(dòng)，在此范圍內(nèi)的變化會引起38 K的溫度浮動(dòng)(圖6)。

圖6 高溫目標(biāo)的面積百分比S與T的關(guān)系Fig.6 Relationship between high temperature torget S and T

綜上所述，4個(gè)參量對T的影響程度從高到低依次為:高溫目標(biāo)面積比S，常溫地物反射率ρ，大氣透過率Tθ，高溫目標(biāo)發(fā)射率ε。

4 可行性分析

根據(jù)式(2)，波長取ETM+7波段的標(biāo)稱波長，溫度取以20 K為步長從300～1080 K的變化范圍，計(jì)算得到地表高溫目標(biāo)發(fā)射的輻射通量密度(表1，λ =2.208 μm)。

表1 不同溫度地表的短波紅外發(fā)射輻射通量密度Tab.1 Shortwave infrared emission radiation flux density of different surface temperature

結(jié)合式(3)與式(7)，得出不同地物在ETM+7波段的反射輻射通量密度(表2，λ=2.208 μm)。

表2 不同地物短波紅外反射輻射通量密度Tab.2Shortwave infrared reflection radiation flux density of different surface features(W·m－2·μm－1)

對高溫目標(biāo)和常溫地物組成的混合像元而言，其輻射能量為反射與發(fā)射的能量綜合。將各參量分別帶入式(4)(5)中，得到 M3=5.14 W·m－2·μm－1相對于600 K的高溫目標(biāo)發(fā)射輻射較小，為了方便進(jìn)一步分析，可忽略M3項(xiàng); 而M4=0.0048 W·m－2·μm－1相對于300 K的常溫地物平均輻射通量密度13.89 W·m－2·μm－1，亦可被忽略。故得到短波紅外波段高溫目標(biāo)的輻射通量密度為

令M1=bM2，b為倍數(shù);并假設(shè)可識別的臨界條件為M≥2M2，則得到

整理后得到

式(11)取等號，由此可確定短波紅外波段可識別高溫目標(biāo)的最低溫度和最小面積(圖7)，落在A區(qū)范圍內(nèi)的高溫目標(biāo)可以被識別出來，而落在B區(qū)范圍內(nèi)的高溫目標(biāo)不能夠被識別出來。

圖7 高溫目標(biāo)識別的臨界面積和溫度Fig.7 Critical S and T for recognition of high temperature target

5 實(shí)際驗(yàn)證

在實(shí)際應(yīng)用中，林火溫度一般在600 K以上［11］，土法煉焦的煉焦過程溫度為950～1200℃(相當(dāng)于950+273.15 ～1200+273.15 K)。因此取S=1和土法煉焦最高溫度T=1470 K作為極值，對圖7曲線取幾個(gè)典型值，得到高溫地物可被識別的臨界值(表3)?？梢钥吹娇勺R別的最低溫度為525 K，可被識別的最小面積為像元的3/10000。

表3 高溫地物可被識別的臨界值Tab.3 Critical value for recognition of high temperature features

對研究區(qū)ETM+數(shù)據(jù)的1545 840個(gè)像元用馬氏距離法［12－13］計(jì)算，在 F 分布［14］下，以信度 α =0.05及自由度 m=6和 n－m=1545 834獲取了289個(gè)高溫像元，經(jīng)野外驗(yàn)證其識別精度達(dá)到88%。對獲取的289個(gè)高溫像元及背景像元進(jìn)行統(tǒng)計(jì)分析，以步長0.2、倍數(shù)從 1.36 ～4.76統(tǒng)計(jì)得到相應(yīng)的頻數(shù)(圖8)。

圖8 高溫目標(biāo)像元與常溫背景輻射通量密度間的倍數(shù)與頻數(shù)關(guān)系Fig.8 Relationship between multiple and frequency of radiation flux densities of high temperature target and background

其中高溫像元輻射通量密度大于1.5倍常溫像元的達(dá)到99%，大于2倍常溫像元的達(dá)到93%。當(dāng)高溫目標(biāo)充滿整個(gè)像元時(shí)(即S=1)，則有M=M1，所需各參量都取均值計(jì)算，得到的溫度為526 K，與上述結(jié)果相吻合。

6 結(jié)果分析

高溫目標(biāo)所占比例S越高，b值越小，則可識別的臨界最小溫度T就越低;溫度T越高，b值越大，則高溫目標(biāo)所占比例S的最小面積就越小。也就是說，當(dāng)高溫目標(biāo)充滿整個(gè)像元，可識別的最低溫度為525 K;當(dāng)高溫目標(biāo)的面積比分別為1/10和1/100時(shí)，可識別的最低溫度分別為615 K和775 K。當(dāng)高溫目標(biāo)的面積比小到占1個(gè)像元的3/10000時(shí)，盡管整個(gè)像元的能量明顯降低，但只要溫度達(dá)到1470 K就可被識別出來。此時(shí)，盡管高溫目標(biāo)輻射峰值的波長不在探測波段范圍內(nèi)，但由于混合像元的輻射能量高于常溫地物的輻射通量密度，故高溫目標(biāo)仍可被有效地識別。當(dāng)溫度較高的高溫目標(biāo)達(dá)到飽和時(shí)，傳感器接收的高溫目標(biāo)值達(dá)到最大，超出所能測定的極限，此時(shí)高溫目標(biāo)一定能被識別出來，只是在溫度計(jì)算時(shí)不能算出其真實(shí)溫度。因此，用短波紅外波段遙感數(shù)據(jù)識別高溫目標(biāo)是可行的。

7 結(jié)論

本文基于能量守恒原理建立了地表高溫目標(biāo)混合像元的輻射能量方程，在進(jìn)行方程參數(shù)的敏感性分析基礎(chǔ)上得出以下結(jié)論:

1)高溫目標(biāo)的溫度取決于ε，Tθ，ρ和S等4個(gè)參量，其中高溫目標(biāo)面積比S和常溫地物反射率ρ對高溫地物的溫度反演最為敏感，是溫度反演精度的決定性因素。

2)輻射能量方程中，高溫目標(biāo)的發(fā)射輻射能量和常溫地物的反射輻射能量是方程最重要的組成部分，由它們與總能量間的關(guān)系可得出S，b和T等3個(gè)參量之間的關(guān)系，進(jìn)一步可得到能被識別的面積和溫度范圍。

3)實(shí)際應(yīng)用中，在土法煉焦溫度范圍內(nèi)，得到高溫目標(biāo)可被識別的最小面積為像元的3/10000，可被識別的最低溫度為525 K。通過對研究區(qū)數(shù)據(jù)實(shí)際驗(yàn)證，得到當(dāng)高溫目標(biāo)充滿整個(gè)像元時(shí)的溫度為526 K，與上述結(jié)果相吻合。

4)經(jīng)野外驗(yàn)證，馬氏距離法應(yīng)用于短波紅外，其識別精度達(dá)到88%，表明短波紅外波段可用于高溫目標(biāo)的識別。

［1］肖 霞.基于類間方差的MODIS森林火災(zāi)監(jiān)測方法研究［D］.合肥:中國科學(xué)技術(shù)大學(xué)，2010.Xiao X.The study of an improved method for forest fire detection based on variance between － class［D］.Hefei:University of Science and Technology of China，2010.

［2］趙 彬，趙文吉，潘 軍，等.NOAA－AVHRR數(shù)據(jù)在吉林省東部林火信息提取中的應(yīng)用［J］.國土資源遙感，2010(1):77－80.Zhao B，Zhao W J，Pan J，et al.The application of NOAA －AVHRR to forest fire detection in eastern Jilin Province［J］.Remote Sensing for Land and Resources，2010(1):77 －80.

［3］劉安濤.地表熱異常遙感溫度反演［D］.長春:吉林大學(xué)，2007.Liu A T.Land surface temperature retrieval of the hot abnormity on remote sensing［D］.Changchun:Jilin University，2007.

［4］朱亞靜，邢立新，潘 軍，等.短波紅外遙感高溫地物目標(biāo)識別方法研究［J］.遙感信息，2011(6):33 －36，41.Zhu Y J，Xing L X，Pan J，et al.Method of identifying high － temperature target using shortwave infrared remote sensing data［J］.Remote Sensing Information，2011(6):33 －36，41.

［5］Barducci A，Guzzi D，Marcoionni P，et al.Comparison of fire temperature retrieved from SWIR and TIR hyperspectral data［J］.Infrared Physics and Technology，2004，46(1/2):1 －9.

［6］Matheson P S，Malheson D S.Comparison of fire temperature and fractional area modeled from SWIR，MIR，and TIR multispectral and SWIR hyperspectral airborne data［J］.Remote Sensing of Environment，2011，115(3):876 －886.

［7］馮鐘葵.關(guān)于地物輻射值的計(jì)算問題［J］.中國遙感衛(wèi)星地面站:用戶簡訊，2002(49):1.Feng Z K.Estimation problem about the surface features radiation value［J］.Remote Sensing Satellite Earth Station of China:User News Brief，2002(49):1.

［8］Chen F L，Zhuang J L，Xu X R，et al.The concept of effective emissivity of nonisothermal mixed pixel and its test［J］.Chinese Science Bulletin，2000，45(9):788 －795.

［9］Berk A，Anderson G P，Acharya P K，et al.Modtran4 user’s manual［M］.Air Force Research Laboratory，1991:11 －53.

［10］Pan J，Xing L X，Wen J C，et al.Inversion method study on short wave infrared remote sensing data high temperature surface feature temperature［C］//Processing of the 2nd International Congress on Image and Signa Processing，Tainijin:IEEE，2009:1 －4.

［11］劉 誠，李亞軍，趙長海，等.氣象衛(wèi)星亞像元火點(diǎn)面積和亮溫估算方法［J］.應(yīng)用氣象學(xué)報(bào)，2004，15(3):273 －280.Liu C，Li Y J，Zhao C H，et al.The method of evaluating sub － pixel size and temperature of fire spot in AVHRR data［J］.Journal of Applied Meteorolgical Science，2004，15(3):273 －280.

［12］孔祥生，苗 放，劉鴻福，等.遙感技術(shù)在監(jiān)測和評價(jià)土法煉焦污染源中的應(yīng)用［J］.成都理工大學(xué)學(xué)報(bào):自然科學(xué)版，2005，32(1):92－96.Kong X S，Miao F，Liu H F，et al.Application research of remote sensing technique in monitoring pollution source by digenous coke－ production［J］.Journal of Chengdu University of Technology:Science and Technology Edition，2005，32(1):92 －96.

［13］紀(jì)宏金.地球化學(xué)背景與異常劃分的多元方法［J］.長春地質(zhì)學(xué)院學(xué)報(bào)，1988，18(3):311 －320.Ji H J.Multivariate method for distinguishing geochemical background from geochemical anomalies［J］.Journal of Changchun University of Earth Science，1988，18(3):311 －320.

［14］紀(jì)宏金.多元正態(tài)總體假設(shè)檢驗(yàn)在礦化帶識別中的應(yīng)用［J］.長春地質(zhì)學(xué)院學(xué)報(bào)，1991，21(3):321 －326.Ji H J.An application of test of hypothesis of mulitivariate normal population to mineralization zone recogniztion［J］.Journal of Changchun University of Earth Science，1991，21(3):321 －326.

Feasibility analysis of shortwave infrared band for recognition of high temperature target

YU Yifan，PAN Jun，XING Lixin，JIANG Lijun，MENG Tao，HAN Xiaojing，ZHOU Caicai
(College of Geoexploration Science and Technology，Jilin University，Changchun 130026，China)

At present，high temperature target recognition mainly uses thermal infrared remote sensing data.In this study，the authors have found that shortwave infrared band has better recognition capability than thermal infrared band.In the mixed pixel of high temperature targets，temperature and area of high temperature objects are unknown.They are the key parameters that can determine the spectral character of mixed pixels.Based on the constant energy principle，the authors formulated the radiation energy equation for the mixed pixel of high temperature targets on the Earth’s surface.The results of the sensibility analysis for the equation parameters show that the area percentage of high temperature targets and the reflection of the normal temperature targets are most sensitive to the invertion of the temperature and the area of high temperature targets.ETM+7 data obtained in Baode of Shanxi and Fugu of Shaanxi were used for study of high temperature target recognition.The radiation flux density of the recognized fire is about 1.36 to 4.76 times that of the background value.Field verification shows that Mahalanobis method has the precision of 88%，suggesting that shortwave infrared band can be used to recognize high temperature targets.

short wave infrared;high temperature target;normal temperature object;feasibility analysis

TP 75

A

1001－070X(2014)01－0025－06

2013－03－18;

2013－04－21

高等學(xué)校博士學(xué)科點(diǎn)專項(xiàng)科研基金新教師類資助課題項(xiàng)目(編號:20110061120067)資助。

10.6046/gtzyyg.2014.01.05

于一凡，潘軍，邢立新，等.短波紅外波段高溫目標(biāo)識別的可行性分析［J］.國土資源遙感，2014，26(1):25－30.(Yu Y F，Pan J，Xing L X，et al.Feasibility analysis of shortwave infrared band for recognition of high temperature target［J］.Remote Sensing for Land and Resources，2014，26(1):25 －30.)

于一凡(1988－)，女，碩士研究生，研究方向?yàn)檫b感數(shù)據(jù)處理及異點(diǎn)識別。Email:yuyifanjy@163.com。

潘 軍(1971－)，男，副教授，主要從事遙感與地理信息系統(tǒng)教學(xué)和科研工作。Email:panj@jlu.edu.cn。

(責(zé)任編輯:劉心季)