參考文獻 ： 基於 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 *****》 藍波段反演aod度敏感度非常高
原理》藍波段大氣反射較強，而地表反射較弱。將清晰天的地表反射率（即地表反射率庫）代入地表反射率公式，反演氣溶膠。
備註》當地表反射率<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作為雲去除閾值。
演算法入門練習記錄
沫璃邀請她的朋友參加週末的派對。沫璃買了3種顏色的氣球，現在她要有這些氣球來裝飾餐桌，每個餐桌只用恰好3個氣球裝飾，要求3個氣球的顏色不能完全一樣，可以是2種或者3種顏色。沫璃想知道這些氣球最多能裝飾多少張餐桌。分析 取出值最大的顏色，只要其他兩種顏色加起來達到總顏色的1 3，那麼答案是氣球總數 3...
演算法入門（一）
1.for迴圈的使用 其實簡單的for 初始化 判斷條件 調整 迴圈，一般都會用。tip break用於跳出for迴圈，continue用於跳出此次迴圈，進入下一次 for 死迴圈 如 int i,j 3,n 5 for i 0 ij 迴圈結束時，i 4 q1.求1000 9999中的aabb型別的...
演算法入門例項（一）
在n n的矩陣中，填入矩陣維度，然後形成乙個蛇形矩陣，例如，當填入維度為4時，形成的蛇形矩陣 10 11 12 1 9 16 13 2 8 15 14 3 7 6 5 4 可以想象對這個蛇形矩陣的繪製時，是在繪製一條路。筆的起點是在第一行的最後一列。筆的軌跡是下，下，下，左，左，左，上，上，上，右，...