本發(fā)明涉涉及對地觀測技術(shù)領(lǐng)域,特別是涉及一種以大氣吸收波段11.73μm為參考波段,對熱紅外高光譜遙感數(shù)據(jù)進行在軌光譜定標方法。
背景技術(shù):傳感器成像環(huán)境變化所導致的光電器件變形、機械振動所導致的光電元件錯位以及光電元件本身的日久老化使得傳感器每個波段的中心波長和半高全寬(Fullwidthathalfmaximum,F(xiàn)WHM)隨成像環(huán)境變化會發(fā)生較大的系統(tǒng)性漂移。因此,實驗室光譜定標所獲取的中心波長和半高全寬并不能完全代表航空遙感數(shù)據(jù)獲取時傳感器的中心波長和半高全寬。尤其是傳感器各個波段的中心波長極容易受到光電器件變形和機械振動的影響而產(chǎn)生偏移,進而影響發(fā)射率和溫度的反演精度。尤其是在大氣吸收波段附近的發(fā)射率反演精度,它很容易受到中心波長偏移的影響。目前,對于高光譜遙感數(shù)據(jù)的光譜定標研究主要集中于基于場景的可見光、近紅外和短波紅外譜段的高光譜遙感數(shù)據(jù)光譜定標?;趫鼍暗墓庾V定標原理主要是假設傳感器的光譜響應函數(shù)(Spectralresponsefunction,SRF)為高斯函數(shù),比較模擬數(shù)據(jù)與原始數(shù)據(jù)在氧氣吸收波段(760nm)、水汽吸收波段(940nm,1140nm)和二氧化碳吸收波段(2060nm)及其附近光譜的差異,并調(diào)整模擬數(shù)據(jù)的中心波長和FWHM,直至二者差異最小,最終得到準確的中心波長和FWHM。按照最終的目標不同,這些光譜定標方法可以分為三類,一是在對比過程中,假設SRF和FWHM不變(使用實驗室定標結(jié)果),只求正確的中心波長;二是假設SRF為高斯函數(shù),同時改變中心波長和FWHM,最終求得正確的中心波長和FWHM;三是假設中心波長和FWHM改變很小,選擇不同類型的SRF得到模擬光譜,并與原始數(shù)據(jù)比較,最終求得最優(yōu)的SRF。除了直接比較模擬數(shù)據(jù)與原始數(shù)據(jù)之外,也有學者提出通過調(diào)整中心波長和FWHM,計算大氣輻射傳輸方程中的參數(shù),并反演反射率,直至所反演的反射率曲線中尖峰消失,得到正確的中心波長和FWHM。。有學者嘗試分別用8-13μm內(nèi)21個水汽吸收波段對SEBASS(Spatiallyenhancedbroadbandarrayspectrographsystem)熱紅外(7.3-13.6μm,128波段)高光譜成像儀獲取的航空熱紅外高光譜遙感數(shù)據(jù)進行了光譜定標,光譜定標所用大氣模擬程序為FASCODE。光譜定標結(jié)果顯示,采用基于場景的光譜定標方法對中心波長的變化更敏感,對FWHM的變化相對不敏感,SEBASS數(shù)據(jù)的中心波長偏移在±0.5個波段寬度內(nèi)??傮w來說,由于熱紅外高光譜傳感器較少,實際獲取的數(shù)據(jù)也較少,導致了對于熱紅外高光譜遙感數(shù)據(jù)的光譜定標的關(guān)注相對較少,還沒有實用化的在軌光譜定標技術(shù)。隨著對地觀測對熱紅外高光譜數(shù)據(jù)的迫切需求,也必將需要實用化的熱紅外高光譜遙感數(shù)據(jù)在軌光譜定標技術(shù)。
技術(shù)實現(xiàn)要素:針對現(xiàn)有的技術(shù)空白和缺點,本發(fā)明所要解決的技術(shù)問題是提供一種準確性好,速度快,易操作的熱紅外高光譜數(shù)據(jù)在軌定標方法。為了解決上述技術(shù)問題,本發(fā)明提供的一種熱紅外高光譜遙感數(shù)據(jù)在軌光譜定標方法,其特點是,光譜定標方法的步驟如下:1)將熱紅外高光譜遙感數(shù)據(jù)獲取時的大氣參數(shù)、成像參數(shù)輸入MODTRAN5.3大氣輻射傳輸模擬軟件得到模擬的1nm超高分辨率入瞳輻射亮度光譜。輸入的大氣參數(shù)主要有大氣模型、水汽含量、氣溶膠模型、氣象視距等。輸入的成像參數(shù)主要有傳感器高度、觀測天頂角、觀測方位角、地表高程、地表反射率(ρs)、地表溫度(Ts)等。一般對于熱紅外遙感來說,輸出的1nm分辨率的輻射亮度光譜所覆蓋譜段范圍為7.0-14.0μm;2)以實驗室光譜定標獲得的中心波長λi和FWHMΔλi為初始值,不斷調(diào)整中心波長(λi+δλi)和FWHM(Fullwidthathalfmaximum,F(xiàn)WHM)(Δλi+δΔλi)對超高分辨率入瞳輻射亮度光譜進行卷積運算,得到模擬的傳感器入瞳輻射亮度光譜。模擬的入瞳輻射亮度計算方法為:其中,Li(δλi,δΔλi)為波段i的入瞳輻射亮度;F(λ)為超高分辨率輻射亮度光譜;Si(λ,λi,Δλi,δλi,δΔλi)為波段i的光譜響應函數(shù),Si(λ,λi,Δλi,δλi,δΔλi)的計算方法為:其中,λ為波長,λi為一般設置為實驗室內(nèi)光譜定標得到的中心波長;△λi一般設置為實驗室內(nèi)光譜定標得到的FWHM。3)選擇模擬的大氣吸收波段(11.73μm)及其附近入...