參考文獻 ：         基於 HJ 衛星的北京市氣溶膠光學厚度遙感反演      冉凡坤

深藍演算法反演氣溶膠光學厚度的依據>>>

在亮像元地表，藍光波段的地表反射率較小，只有可見光波段的 1/2～1/4，可以通過對藍光波段地表反射率的模擬實現對氣溶膠光學厚度的反演。同時，深藍波段的分子散射和吸收都遠遠的低於其它波段，氣溶膠垂直廓線的影響對深藍波段氣溶膠光學厚度反演的影響甚微。

反演演算法>>>
1. 計算衛星在  412、490 和 670nm 的表觀反射率；

2. 對衛星資料進行高程修正，資料基於數字高程；

3. 基於藍波段的反射率變化和吸收氣溶膠指數對雲進行提取；

4. 建立412nm、490nm 和 670nm 的地表反射率庫，地表反射率取的是晴朗無雲條件下的地表反射率；

5. 基於地表反射率和輻射傳輸模型 6S，通過建立查詢表得到氣溶膠光學厚度；

敏感度測試>>>

當地表反射率較小時（< 0.2），衛星觀測到的表觀反射率與氣溶膠光學厚度有著很好的線性相關性，且敏感度很大。

但隨著地表反射率的增大（>0.2），衛星觀測到的表觀反射率與氣溶膠光學厚度的相關性幾乎為 0。

因此當利用深藍演算法進行氣溶膠光學厚度反演時，必須保證地表反射率不高於 0.1。

藍波段地表反射率的分佈比較集中，80%介於 0～0.1 之間，60%小於 0.05，只有 10%高於 0.2 

參考文獻 ：         基於深藍演算法的 ＨＪ－１ＣＣＤ資料快速大氣校正模型        王中挺 等

原理>>>>藍波段大氣反射較強，而地表反射較弱。將清晰天的地表反射率（即地表反射率庫）代入地表反射率公式，反演氣溶膠。

備註>>>>當地表反射率<0.1，反射率隨AOD有較好的線性關係； 

當地表反射率>0.2，基本無變化。                 =====》去雲閾值的設立依據

MODIS地表反射率產品>>>>包括250  m單日合成、500 m與1 km單日合成、250 m 8日合成、500 m 8日合成、CMG網格全球日合成(解析度為0.05°)等5種形式

                                                  為消除地表反射率的短期變化，選擇8日合成產品進行地表反射率庫的構建。該產品採用正弦投影，將全球分為10°×10°的小塊存放。

                                                   將MODIS的8日合成地表反射率產品進行波段提取、重投影、重取樣、影象剪下和鑲嵌等處理得到地表反射率庫。

                                       &&&對地表反射率資料進行綜合統計，分析藍波段（MODIS第3波段）地表反射率在不同值域區間的分佈情況。

                                       &&&得到藍波段地表反射率累積百分比。====》藍波段地表反射率整體較小。

                                                                                          ====》在地表反射率小於0.1時，表觀反射率隨著AOD有著較好的線性關係；

                                                                                                       在地表反射率大於0.2時，表觀反射率隨著AOD基本沒有變化。

                                                                                                       因此，當地表反射率小於0.1時，本演算法能夠較好的執行。

查詢表>>>>1152個觀測幾何，12個太陽天頂角(0—66°，間隔6°)、6個觀測天頂角(0—35°，間隔為7°)、16個太陽與衛星之間的相對方位角(0—180°，間隔12°)；

                     大氣氣溶膠模式引數假設為大陸型氣溶膠，並設立6個大氣氣溶膠光學厚度值(0.0001、0.25、0.5、1、1.5和1.95)；

                     波段自定義為HJ-1  CCD相機第一波段的響應函式；暫時不考慮地表反射，反射率設為0。

？？？？如果不使用HJ-1CCD資料，如何確定自定義波段的響應函式？？？？

HJ-1CCD觀測幾何的轉換>>>.

輻射定標及雲去除>>>>首先，從CCD資料的輔助xml檔案中讀取輻射定標係數，將DN值轉換為表觀輻亮度；

                                         然後，根據上一步生成的太陽天頂角，，將獲得的表觀輻亮度轉換為表觀反射率。

            雲的去除則是通過閾值實現。研究表明(Acker-man  等，2010)，在紅光波段衛星觀測到的雲的反射率一般大於0.2====》選取0.2作為雲去除閾值。