定量/高光谱遥感之——传感器定标
我們常用影像的像元值大多是經過量化的、無量綱的DN值，而進行遙感定量化分析時，常用到輻射亮度值、反射率值、溫度值等物理量。傳感器定標就是要獲得這些物理量的過程。本專題敘述的主要是衛星傳感器的定標，航空傳感器有很多相似地方。
本專題包括以下內容：
l?傳感器定標概念
l?傳感器定標類型
l? ENVI下的傳感器定標
1傳感器定標概念
傳感器定標很多地方又名為輻射定標，嚴格意義上講，輻射定標是傳感器定標的一部分內容。以下是國內的定義，如趙英時等《遙感應用分析原理與方法》上描述：定標是將遙感器所得的測量值變換為絕對亮度或變換為與地表反射率、表面溫度等物理量有關的相對值的處理過程。或者說，遙感器定標就是建立遙感器每個探測器輸出值與該探測器對應的實際地物輻射亮度之間的定量關系；建立遙感傳感器的數字量化輸出值DN與其所對應視場中輻射亮度值之間的定量關系（陳述彭）。輻射亮度的典型的單位為：W/cm2.μm.sr（瓦特/平方厘米.微米.球面度）。
我們總結以上的定義，通俗的說法：傳感器定標就是將圖像的數字量化值（DN）轉化為輻射亮度值或者反射率或者表面溫度等物理量的處理過程。其中反射率又分為大氣外層表觀反射率和地表實際反射率，后者又屬于大氣校正的范疇，有的時候也會將大氣校正納入傳感器定標的一種途徑。
2傳感器定標類型（定標參數測量）
傳感器定標可分為絕對定標和相對定標。絕對定標是獲取圖像上目標物的絕對輻射值等物理量；相對定標是將圖像目標物輻射量歸一化某個值范圍內，比如以其他數據作為基準。
傳感器定標可分為三個階段或者說三個方面內容：①發射前的實驗室定標；②基于星載定標器的星上定標；③發射后的定標（場地定標）。
一、???????????實驗室定標
在遙感器發射之前對其進行的波長位置、輻射精度、光譜特性等進行精確測量，也就是實驗室定標。它一般包含兩部分內容：
l?光譜定標
確定遙感傳感器每個波段的中心波長和帶寬，以及光譜響應函數。
l?輻射定標
在模擬太空環境的實驗室中，建立傳感器輸出的量化值（DN）與傳感器入瞳處的輻射亮度之間的模型，一般用線性模型表示。
Y=aX+b
a和b就是我們常說的偏移和增益。
也有些傳感器使用其他表達式，如環境小衛星中的元數據說明文件中的這段話明確說明了定標公式、參數值、定標單位信息：
(gain1,Fielddata,L=DN/g+L0,W*m^(-2)*sr^(-1)*um^(-1))B1:g0.5782, L0 3.4608, B2:g 0.5087, L0 5.8769, B3:g 0.6825, L0 8.0069, B4:g 0.6468,L0 8.8583
二、???????????星上定標
有些衛星載有輻射定標源、定標光學系統，在成像時實時、連續的進行定標。
三、???????????場地定標
場地定標指的是遙感器處于正常運行條件下，選擇輻射定標場地。一般選擇沙漠地區，它的光譜響應穩定，如利比亞沙漠用于定標AVHRR、北非沙漠定標SPOT影像、敦煌西戈壁沙漠定標CBERS影像、美國的白沙導彈靶場常用于高分辨率圖像的定標。通過選擇典型的均勻穩定目標，用精密儀器進行地面同步測量感器過頂時的大氣環境參量和地物反射率，利用遙感方程，建立圖像與實際地物間的數學關系，得到定標參數以完成精確的傳感器定標。
基本原理：在遙感器飛越輻射定標場地上空時，在定標場地選擇偌干個像元區，測量傳感器對應的地物各波段光譜反射率和大氣環境參量（大氣氣溶膠光學厚度，大氣中水，臭氧含量等）等參量，并利用大氣輻射傳輸模型等手段求解傳感器入瞳處各光譜帶的輻射亮度，最后確定它與傳感器輸出的數字量化值之間的數量關系，求解定標系數，并估算定標不確定性。
主要方法有：
反射率法：在衛星過頂時同步測量地面目標反射率因子和大氣光學參量（如大氣光學厚度、大氣柱水汽含量等）然后利用大氣輻射傳輸模型計算出遙感器入瞳處輻射亮度值。具有較高的精度。
輻亮度法：采用經過嚴格光譜與輻射標定的輻射計，通過航空平臺實現與衛星遙感器觀測幾何相似的同步測量，把機載輻射計測量的輻射度作為已知量，去標定飛行中遙感器的輻射量，從而實現衛星的標定。最后輻射校正系數的誤差以輻射計的定標誤差為主。
輻照度法：又稱改進的反射率法，利用地面測量的向下漫射與總輻射度值來確定衛星遙感器高度的表觀反射率，進而確定出遙感器入瞳處輻射亮度，。這種方法是使用解析近似方法來計算反射率，從而可大大縮減計算時間和計算復雜性。
因此，我們常用的定標參數，有使用實驗室定標的結果（如高分辨率傳感器QuickBird、WorldView-1等）；也有使用實驗室定標與星上定標相結合的參數（如NOAA、MSS等）；由于設備老化，Landsat TM5的定標參數有用實驗室定標的（2003年前），也有用經過場地定標的參數（2003年后）；也有些學者為了獲得更精確定標數據，使用場地定標的方法完成一些傳感器數據的定標過程。
3定標實現
上述傳感器定標的方法很多，看起來也非常復雜，其實我們對遙感影像進行定標時候，只需要知道這類影像的定標公式和參數即可，這些公式和參數一般在數據帶的元數據文件中，或者固定的公式和參數。現在很多軟件提供的特定傳感器的定標工具，如果沒有也可以使用波段運算的方式實現定標。
4 ENVI下的傳感器定標
ENVI支持很多數據的定標，包括ASTER、MODIS、AVHRR, MSS,TM ,QuickBird, WorldView-1,TIMS等。也可以根據定標參數利用BandMath工具很方便的完成定標。QuickBird, WorldView-1的定標比較簡單。或者使用Apply gaine and offset工具可以實現Y=aX+b公式的定標。
下面介紹在ENVI下的ASTER、MODIS、AVHRR, Landsat的定標過程。
3.1 ASTER與MODIS定標
對于ASTER L1A/L1B和MODIS 02級數據，在打開數據時會自動完成對數據的定標。
如圖1所示打開ASTER L1B的結果，在波段列表中，自動讀取各個波段的中心波長信息，并按照波段范圍信息（VNIR、SWIR、TIR）分組波段。其中VNIR、SWIR自動定標為輻射亮度，單位是：W/m2/sr/μm；TIR數據定標為大氣表觀溫度值，單位：開爾文。打開其中一個數據，瀏覽像元值，可以看到已經定標為浮點型的輻射亮度值。
?
圖1： ASTER L1B數據
如圖2為打開MODIS 02級1km數據，其中250米和500米的波段經過重采樣為1km加入這個數據集中。ENVI根據各個波段的中心波長信息定標為三個類型數據：反射率數據（Reflectance）、輻射亮度值數據（Radiance）和發射率數據（Emissive）。其中反射率和發射率為0~1無單位值，輻射亮度值單位是：W/m2/μm/sr。
?
圖2 MODIS 02級數據
如果打開原始的ASTER和MODIS的DN值數據，可以在ENVI主菜單中選擇File->preferences，切換到Miscellaneous面板，將Auto-Correct ASTER/MODIS項設置為NO。
3.2 Landsat數據定標
對于Landsat4/5數據可以手動選擇以下兩種定標公式：
式中:
? QCAL為原始量化的DN值
? LMINλ為QCAL = 0時的輻射亮度值
? LMAXλ為QCAL = QCALMAX時的輻射亮度值
??注：LMINλ和LMAXλ的值取自Chander, Markham, and Helder (2009)的研究成果。
? QCALMIN是最小量化定標像素值（與LMINλ類似）。取值如下：
??1：LPGS產品
??1：04 April2004之后的NLAPS產品
??0：04 April2004之前的NLAPS產品
??注：如果沒有元數據信息，QCALMIN取默認值1（TM和ETM+)）或者0 (MSS)。
? QCALMAX為最小量化定標像素值（與LMAXλ類似）。根據元數據信息取值為127, 254, 255。當缺少元數據時，QCALMAX取默認值：255 (TM和ETM+)或者127 (MSS)。作為結果的輻射亮度值的單位為(W/(cm2*sr*μm))。
定標參數使用Chander, Markham, and Helder (2009)的研究成果，其中LPGS和NLAPS分別是兩種數據處理系統得到的產品：the Level 1 Product Generation System (LPGS) 和the
National Land Archive Production System(NLAPS)，從2008年12月份開始，L7 ETM+ 和L5都是以LPGS系統處理，L4 TM和MSS以NLAPS系統處理。具體參數如下所示：
?
圖3：NLAPS處理和LPGS處理產品的MSS定標參數
?
圖4：NLAPS處理產品的TM4、5定標參數
說明：L5 TM在(1984–1991)間數據使用(LMAXλ=169, 333)定標參數。
?
圖5：LPGS處理產品的ETM+數據定標參數
同時ENVI還可以定標為表觀大氣反射率（ρp）：
?????
式中：
Lλ為輻射亮度值
d為天文單位的日地距離
ESUNλ為太陽表觀輻射率均值，對于Landsat7 ETM+，ENVI使用“the Landsat 7 Science Data Users Handbook”上記錄的參數；對于LandsatTM 4/5,ENVI使用Chander and Markham (2003)研究成果
θs是以度為單位的太陽高度角。
如果定標的Landsat數據帶有元數據，ENVI會自動識別定標參數。對于帶元數據文件的Landsat 7 GeoTIFF（*_MTL.txt），需要通過這種方式打開文件File→Open External File→Landsat →GeoTIFFwith Metadata。選擇Basic Tools→Preprocessing→Calibration Utilities→Landsat Calibration，可以打開定標工具。
?
?
圖6： Landsat 5 GeoTIFF（*_MTL.txt）定標界面
當遇到沒有元數據文件的數據時，需要手動輸入定標參數，定標參數可以通過NASA網站查詢或者是一些研究成果，由于不同處理系統或者數據經過重采樣后原始的DN會有所不同，因此定標參數也有很多套標準。如下網址為NASA公布的定標參數。
http://landsathandbook.gsfc.nasa.gov/handbook/handbook_htmls/chapter9/chapter9.html
下面表1和表2是Landsat5的其他版本定標參數供參考。
表1 Landsat TM4/5飛行前輻射定標系數（實驗室定標）
?
| Band | Gain | offset |
| 1 | 0.602 | -1.52 |
| 2 | 1.170 | -2.84 |
| 3 | 0.806 | -1.17 |
| 4 | 0.815 | -1.51 |
| 5 | 0.108 | -0.37 |
| 6 | 0.055 | 1.24 |
| 7 | 0.057 | -0.15 |
表2 Landsat TM 5輻射定標系數（Gyannesh Chander等（2003年））
?
| ? | 從1984年3月1號至2003年5月4號 | 2003年5月5號之后 | ||||||
| Band | LMIN | LMAX | Gain | offset | LMIN | LMAX | Gain | offset |
| 1 | -1.52 | 152.10 | 0.602431 | -1.52 | -1.52 | 193.0 | 0.762824 | -1.52 |
| 2 | -2.84 | 296.81 | 1.175100 | -2.84 | -2.84 | 365.0 | 1.442510 | -2.84 |
| 3 | -1.17 | 204.30 | 0.805765 | -1.17 | -1.17 | 264.0 | 1.039880 | -1.17 |
| 4 | -1.51 | 206.20 | 0.814549 | -1.51 | -1.51 | 221.0 | 0.872588 | -1.51 |
| 5 | -0.37 | 27.19 | 0.108078 | -0.37 | -0.37 | 30.2 | 0.119882 | -0.37 |
| 6 | 1.2378 | 15.303 | 0.055158 | 1.2378 | 1.2378 | 15.303 | 0.055158 | 1.2378 |
| 7 | -0.15 | 14.38 | 0.056980 | -0.15 | -0.15 | 16.5 | 0.065294 | -1.15 |
?
3.3 AVHRR數據定標
ENVI提供的AVHRR定標工具可以對來自NOAA-6、-7、-8、-9、-10、-11、-12、-14、-15、-16、-17、-18和-19衛星的AVHRR數據進行定標。波段1和2被定標為反射率，波段3、4和5被定標為亮溫值（單位：開爾文）。其中NOAA-6、-7、-8、-9、-10、-11選擇主模塊->Basic Tools（或Spectral）->Preprocessing-> Data-SpecificUtilities->AVHRR->Calibrate NOAA 6/7/8/9/10/11。界面比較簡單。
在SST技術方面，ENVI使用多通道海面溫度（Multi-Channel Sea Surface Temperature ——MCSST）計算算法（包括三種：Split-window、dual-window和triple-window），分別用于白天和夜間的數據。每種方法的計算公式如下所示。
1.???????NOAA-12,-14,-15：
l? Day MCSST Split
Ts = a0 + a1*band4 + a2* (band4 - band5) +a3* (band4 - band5) * (sec(φ) -1)
l? Night MCSST Split
Ts = a0 + a1*band4 + a2* (band4 - band5) +a3* (band4 - band5) * (sec(φ)-1)
l? Night MCSST Dual
Ts = a0 + a1*band4 + a2* (band3 - band4) +a3* (sec(φ) -1)
l? Night MCSST Triple
Ts = a0 + a1*band4 + a2* (band3 - band5) +a3* (sec(φ)-1)
2.???????NOAA-16,-17,-18：
l? Day MCSST Split
Ts = a0 + a1*band4 + a2*band5 + a3* (band4- band5) * (sec(φ) -1)
l? Night MCSST Split
Ts = a0 + a1*band4 + a2*band5 + a3* (band4- band5) * (sec(φ) -1)
l? Night MCSST Dual
Ts = a0 + a1*band4 + a2* (band3 - band4) +a3* (sec(φ) -1)
l? Night MCSST Triple
Ts = a0 + a1*band3 + a2*band4 + a3*band5 +a4* (band3 - band5) * (sec(φ) -1)
3.???????NOAA-19
l? Day MCSST Split
Ts = a0 + a1*band4 + a2(band4 - band5) + a3(band4- band5)(sec(φ) -1)
l? Night MCSSTSplit
Ts = a0 + a1*band4 + a2(band4 - band5) + a3(band4- band5)(sec(φ) -1)
l? Night MCSSTDual
Ts = a0 + a1*band4 + a2(band3 - band4) + a3(sec(φ)-1)
l? Night MCSSTTriple
Ts = a0 + a1*band4 + a2(band3 - band5) + a3(band3- band5)(sec(φ) -1)
4.???????NOAA MetOp-A FRAC：
l? Day MCSST Split
Ts = a0 + a1*band4 + a2* (band4 - band5) +a3* (band4 - band5) * (sec(φ) -1)
l? Night MCSST Split
Ts = a0 + a1*band4 + a2* (band4 - band5) +a3* (band4 - band5) * (sec(φ) -1)
l? Night MCSST Dual
Ts = a0 + a1*band4 + a2* (band3 - band4) +a3* (sec(φ) -1)
l? Night MCSST Triple
Ts = a0 + a1*band4 + a2* (band3 - band5) +a3*(band3 - band5) * (sec(φ) -1)
其中Sec(φ)是衛星高度角的正切值，Band3、4、5是定標后的亮溫值（單位:開爾文），“a”參數值是根據漂流浮標和熱帶太平洋固定浮標數據，利用回歸模型獲得。不同的衛星對應的參數不一樣，詳細參見表3~6。
?
表3 Day MCSST Split
| Satellite | a0 | a1 | a2 | a3 |
| NOAA-12 | -263.006 | 0.963563 | 2.579211 | 0.242598 |
| NOAA-14 | -278.43 | 1.017342 | 2.139588 | 0.779706 |
| NOAA-15 | -261.029735 | 0.959456 | 2.663579879 | 0.570613 |
| NOAA-16 | -273.77 | 3.301267 | -2.30195 | 0.628966 |
| NOAA-17 | -271.206 | 0.992818 | -2.49916 | 0.915103 |
| NOAA-18 | -280.43 | 1.02453 | 2.10044 | 0.0784059 |
| NOAA-19 | -278.74596 | 1.01922 | 1.72270 | 0.80263 |
| MetOp-A FRAC | -273.816 | 1.00255 | 2.39451 | 0.903773 |
表 4 Night MCSST Split
| Satellite | a0 | a1 | a2 | a3 |
| NOAA-12 | -263.94 | 0.967077 | 2.384376 | 0.480788 |
| NOAA-14 | -282.24 | 1.029088 | 2.275385 | 0.752567 |
| NOAA-15 | -271.3969724 | 0.993892 | 2.7523466369 | 0.662999 |
| NOAA-16 | -273.15 | -2.53655 | 3.5316 | 0.753291 |
| NOAA-17 | -276.59 | 1.01015 | -2.5815 | 1.000541 |
| NOAA-18 | -276.075 | 1.00841 | 2.23459 | 0.736946 |
| NOAA-19 | -277.71304 | 1.01432 | 1.91798 | 0.72064 |
| MetOp-A FRAC | -277.447 | 1.01377 | 2.52362 | 1.03056 |
?
?
表5 Night MCSST Dual
| Satellite | a0 | a1 | a2 | a3 |
| NOAA-12 | -279.846 | 1.031355 | 1.288548 | 2.265075 |
| NOAA-14 | -273.914 | 1.008751 | 1.409936 | 1.975581 |
| NOAA-15 | -283.5117285 | 1.041037 | 1.5875819344 | 1.67743 |
| NOAA-16 | -273.082 | 1.50825 | 1.00413 | 1.52452 |
| NOAA-17 | -276.603 | 1.01805 | 1.49789 | 1.96181 |
| NOAA-18 | -279.755 | 1.02958 | 0.0502887 | 1.78302 |
| NOAA-19 | -276.61174 | 1.01873 | 1.47374 | 1.88560 |
| MetOp-A FRAC | -273.235 | 1.00711 | 1.49927 | 1.88373 |
表6 Night MCSST Triple?
| Satellite | a0 | a1 | a2 | a3 | a4 |
| NOAA-12 | -271.971 | 1.000281 | 0.911173 | 1.710028 | ? |
| NOAA-14 | -275.364 | 1.010037 | 0.920822 | 1.760411 | ? |
| NOAA-15 | -276.7558563 | 1.015354 | 1.0635723508 | 1.294955 | ? |
| NOAA-16 | -271.763 | 0.733532 | 1.01684 | 0.344182 | -0.753123 |
| NOAA-17 | -274.622 | 1.00903 | 0.913248 | 0.440015 | ? |
| NOAA-18 | -274.398 | 1.00820 | 0.841674 | 0.377061 | ? |
| NOAA-19 | -275.24563 | 1.01084 | 0.81643 | 0.43235 | ? |
| MetOp-A FRAC | -273.044 | 1.00424 | 0.894349 | 0.508159 | ? |
5總結
對于一些新上天的傳感器，實驗室定標參數基本能滿足應用需求。當對精度要求較高時候，往往需要通過場地定標的方式提高定標的結果。對于一些服役時間較長的傳感器，也需要通過場地定標的方式來校對實驗室定標結果，最典型的是Landsat系列衛星。
值得注意的是，為保證精度，定標時的數據源往往是初級別的，因為定標是對DN進行運算，如果經過重采樣的數據其DN值會發生一定變化。
?
總結
以上是生活随笔為你收集整理的定量/高光谱遥感之——传感器定标的全部內容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 基于SpringBoot框架的网上购物商
- 下一篇： 《九年双11：互联网技术超级工程》免费下