定量/高光谱遥感之——大气校正（二）
2.6 FLAASH大氣校正
FLAASH是基于MODTRAN4+輻射傳輸模型，MODTRAN模型是由進行大氣校正算法研究的領先者SpectralSpectral Sciences, Inc和美國空軍實驗室（Air Force Research Laboratory）共同研發(fā)。Exelis?VIS公司負責集成和GUI設計。
1.??????FLAASH特點
l?支持傳感器種類多，包括多光譜的ASTER，AVHRR，GeoEye-1，IKONOS，IRS，Landsat，MODIS，SeaWiFS，SPOT，QuickBird，RapidEye等，高光譜HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)、AISA等。可以通過自定義波譜響應函數(shù)支持更多的傳感器。工程化應用價值比較明顯。
l?FLAASH采用了MODTRAN4+輻射傳輸模型，該算法精度高。任何有關影像的標準MODTRAN大氣模型和氣溶膠類型都可以直接使用。
l?通過影像像素光譜上的特征來估計大氣的屬性，不依賴遙感成像時同步測量的大氣參數(shù)數(shù)據(jù)。
l?可以有效地去除水蒸氣/氣溶膠散射效應，同時基于像素級的校正，矯正目標像元和鄰近像元交叉輻射的“鄰近效應”。
l?對由于人為抑止而導致波譜噪聲進行光譜平滑處理。作為結(jié)果，除了真實地表反射率外，還可以得到整幅圖像內(nèi)的能見度、卷云與薄云的分類影像、水氣含量數(shù)據(jù)。
2.??????使用ENVI大氣校正模塊
ENVI大氣校正模塊的使用主要又以下7個方面組成：1、輸入文件準備，2、基本參數(shù)設置，3、多光譜數(shù)據(jù)參數(shù)設置，4、高光譜數(shù)據(jù)參數(shù)設置，5、高級設置，6、輸出文件，7、處理結(jié)果。下面介紹這7個方面內(nèi)容。
(一)?輸入文件準備
1)????????支持傳感器類型
高光譜包括：HyMAP、AVIRIS、CASI、HYDICE、HYPERION(EO-1)、AISA等；多光譜包括：ASTER，AVHRR，GeoEye-1，IKONOS，IRS，Landsat，MODIS，SeaWiFS，SPOT，QuickBird，RapidEye等，
航空：860nm-1135nm波長范圍
2)????????數(shù)據(jù)是經(jīng)過定標后的輻射亮度（輻射率）數(shù)據(jù)，單位是：（μW）/（cm2*nm*sr）。
3)????????數(shù)據(jù)帶有中心波長（wavelenth）值，如果是高光譜還必須有波段寬度（FWHM）,這兩個參數(shù)都可以通過編輯頭文件信息輸入（Edit Header）。
4)????????數(shù)據(jù)類型
支持四種數(shù)據(jù)類型：浮點型（floating）、長整型(long integer )、整型（integer）和無符號整型(unsigned int)。數(shù)據(jù)存儲類型：ENVI標準柵格格式文件，且是BIP或者BIL。
5)????????波譜范圍：flaash能夠做的數(shù)據(jù)光譜范圍是0.4－2500μm。
(二)?基本參數(shù)設置
1)????????輸入文件及輸出路徑設置，如圖6所示。
2)????????傳感器基本信息設置
選擇傳感器類型，成像中心點經(jīng)緯度，成像時間，高度信息（成像區(qū)域和傳感器飛行高度），傳感器這些都可以在數(shù)據(jù)自帶信息文件里獲得，需要注意的是南半球和西半球要用負值。
?
圖6基本參數(shù)設置
3)????????大氣模型（Atmospheric Model）
提供6種大氣模型，如表2所示，可根據(jù)緯度和成像季節(jié)對照表3查找對應的大氣模型。
?
| Model Atmosphere | Water Vapor (std atm-cm) | Water Vapor (g/cm2) | Surface Air Temperature |
| Sub-Arctic Winter (SAW) | 518 | 0.42 | -16 °C or 3 °F |
| Mid-Latitude Winter (MLW) | 1060 | 0.85 | -1 °C or 30 °F |
| U.S. Standard (US) | 1762 | 1.42 | 15 °C or 59 ° |
| Sub-Arctic Summer (SAS) | 2589 | 2.08 | 14 °C or 57 ° |
| Mid-Latitude Summer (MLS) | 3636 | 2.92 | 21 °C or 70 ° |
| Tropical (T) | 5119 | 4.11 | 27 °C or 80 ° |
??????????????????????????????????表2六種標準的大氣模型
| Latitude (°N) | Jan. | March | May | July | Sept. | Nov. |
| 80 | SAW | SAW | SAW | MLW | MLW | SAW |
| 70 | SAW | SAW | MLW | MLW | MLW | SAW |
| 60 | MLW | MLW | MLW | SAS | SAS | MLW |
| 50 | MLW | MLW | SAS | SAS | SAS | SAS |
| 40 | SAS | SAS | SAS | MLS | MLS | SAS |
| 30 | MLS | MLS | MLS | T | T | MLS |
| 20 | T | T | T | T | T | T |
| 10 | T | T | T | T | T | T |
| 0 | T | T | T | T | T | T |
| -10 | T | T | T | T | T | T |
| -20 | T | T | T | MLS | MLS | T |
| -30 | MLS | MLS | MLS | MLS | MLS | MLS |
| -40 | SAS | SAS | SAS | SAS | SAS | SAS |
| -50 | SAS | SAS | SAS | MLW | MLW | SAS |
| -60 | MLW | MLW | MLW | MLW | MLW | MLW |
| -70 | MLW | MLW | MLW | MLW | MLW | MLW |
| -80 | MLW | MLW | MLW | SAW | MLW | MLW |
表3數(shù)據(jù)經(jīng)緯度與獲取時間對應的大氣模型
4)????????水氣反演Water Retrieval
水氣反演設置，采用兩種方式對水氣進行反演：
a)????????利用水氣反演模型恢復影像中每個像元的水氣量
使用水氣反演模型，數(shù)據(jù)必須具有15nm以上波譜分辨率，且至少覆蓋以下波譜范圍之一：
l?1050-1210 nm (對應1135 nm)
l?870-1020 nm (對應940 nm)
l?770-870 nm (對應820 nm)
對于大多數(shù)多光譜傳感器，水氣反演默認顯示的是NO，因為大多數(shù)傳感器沒有適當?shù)牟ǘ蝸硌a償水氣的影響。
b)????????單一的水氣因數(shù)用于整體影像，默認是1，
對于多光譜數(shù)據(jù)使用水氣反演模型，可以在多光譜設置中手動設置水氣波段
5)????????氣溶膠模型（Aerosol Model）
a)????????提供四種標準MODTRAN氣溶膠模型
Rural（鄉(xiāng)村）、Urban（城市）、Maritime（海洋）、Tropospheric（對流層，能見度在40km以上）
b)????????兩種氣溶膠反演方法
2-Band（K-T）方法(類似模糊減少法)，如果沒有找到適應的黑值（一般是陰影區(qū)或者水體），系統(tǒng)將采用能見度值來計算；所以即使選擇了該選項也要給能見度。
選擇None，采用能見度值參與氣溶膠反演，能見度值大約參考值參見表4
| 天氣條件 | 能見度 |
| 晴朗 | 40 to 100 km |
| 中等霧、陰霾 | 20 to 30 km |
| 厚霧、陰霾 | 15 km 或者更少 |
表4天氣條件與能見度對照表
6)????????光譜打磨(高光譜) Spectral Polishing
對高光譜數(shù)據(jù)的光譜進行打磨處理，對波譜曲線進行微調(diào)，使波譜曲線更加近似于真實地物的波譜曲線。
7)????????重新定標波長（Recalibratingthe Input Wavelengths）
???各個波段的中心波長對大氣校正結(jié)果影像比較大。
(三)?多光譜數(shù)據(jù)參數(shù)設置
???當基本設置里設置了水氣反演以及氣溶膠反演，相應的要在此設置相關參數(shù)（單擊Multispectral Settings）。有兩種設置方式：文件方式和圖形方式，一般選擇圖形方式。
1)????????水氣反演模型參數(shù)
| Water Retrieval | 1135 nm | absorption | 1117 - 1143 nm |
| reference upper wing | 1184 - 1210 nm | ||
| reference lower wing | 1050 - 1067 nm | ||
| 940 nm | absorption | 935 - 955 nm | |
| reference upper wing | 870 - 890 nm | ||
| reference lower wing | 995 - 1020 nm | ||
| 820 nm | absorption | 810 - 830 nm | |
| reference upper wing | 850 - 870 nm | ||
| reference lower wing | 770 - 790 nm | ||
| Aerosol Retrieval | KT upper | 2100 - 2250 nm | |
| KT lower | 640 - 680 nm | ||
| Cloud Masking | cirrus clouds | 1367 - 1383 nm | |
表5推薦使用波長與波段對應關系
2)????????氣溶膠模型參數(shù)設置（用氣溶膠模型要求數(shù)據(jù)波段覆蓋660nm和2100nm波譜），具體波段選擇可以從表5所列中選擇。
?
圖7多光譜數(shù)據(jù)參數(shù)設置
(四)?高光譜數(shù)據(jù)參數(shù)設置
3)????????自動選擇通道定義(推薦)
4)????????設置通道定義
?
圖8高光譜設置
(五)?高級設置
?
圖9高級參數(shù)設置
在高級設置面板里可以設置以下幾個參數(shù)：
1)????????光譜定義文件：內(nèi)置AVIRIS、HYMAP、HYDICE、HYPERION、CASI、AISA。
2)????????氣溶膠厚度系數(shù)：用于技術鄰域效應范圍。一般值為1~2km。
3)????????CO2混合比率：默認為390ppm。
4)????????使用領域糾正（Use Adjacency Correction）。
5)????????使用以前的MODTRAN模型計算結(jié)果。
6)????????設置MODTRAN模型的光譜分辨率（推薦值5 cm-1)。
7)????????設置MODTRAN多散射模型。
提供三種模型供選擇Isaacs，DISORT和Scaled DISORT。默認是Scaled DISORT和streams為8，這種模型對于小于1000nm具有較高的精度；
8)????????天頂角\方位角（針對非星下點傳感器）。
9)????????輸出反射率縮放系數(shù)（Output Reflectance Scale Factor）：為了降低結(jié)果儲存空間，默認反射率乘于10000。
(六)?輸出文件
經(jīng)過大氣校正模塊處理，可以得到以下結(jié)果：
l?表面反射率影像
l?水氣含量數(shù)據(jù)
l?云圖
l?日志文件
l?FLAASH大氣校正模板文件
(七)?處理結(jié)果
此模塊的波譜還原精度高，使用大氣校正模塊進行大氣校正，能高保真地恢復地物波譜信息。
?
圖10校正前后效果圖
3.?????????常見問題
??? 1、如圖11所示錯誤，為因為沒有設置輸出反射率文件名，解決方法是單擊Output Reflectance File按鈕，選擇反射率數(shù)據(jù)輸出目錄及文件名。
?
?
2、在做modis數(shù)據(jù)時候，不能使用領域糾正（Use Adjacency Correction）。
3、當數(shù)據(jù)有地理坐標時候，像元大小要也參數(shù)中設置保持一致。如你tm數(shù)據(jù)是25米的，那么選擇Landsat時候，像元大小默認是30米，這個時候如果改成25米會報錯。
4、得到的結(jié)果有負值，這個不好解釋，可以理解為誤差或者其他因素。或者對太陽光強吸收區(qū)域，如水體。
更多錯誤詳細參考：http://blog.sina.com.cn/s/blog_764b1e9d0100pvrk.html
總結(jié)
以上是生活随笔為你收集整理的定量/高光谱遥感之——大气校正（二）的全部內(nèi)容，希望文章能夠幫你解決所遇到的問題。
- 上一篇： 苹果三代耳机_苹果准备三款AirPods
- 下一篇： Lilu：Rails Mockup驱动开