地表反射率計(jì)算

1、黃石市地表反射率計(jì)算一、數(shù)據(jù)預(yù)處理1、打開：用ENVI5.1將黃石市2000年遙感影像數(shù)據(jù)的3,、4、5波段打開（1）用鼠標(biāo)左鍵雙擊ENVI5.1圖標(biāo)，打開ENVI5.1程序；（2）打開黃石市2000年遙感影像數(shù)據(jù)的3,、4、5波段。 FileOpen Image File選擇黃石市2000年遙感影像數(shù)據(jù)的3、4、5波段打開。2、 合成：對(duì)黃石市2000年遙感影像數(shù)據(jù)的3、4、5波段進(jìn)行合成。 Basic ToolsLayer StackingImport File選擇黃石市2000年遙感影像數(shù)據(jù)的3、4、5波段OkChoose命名（2000_band543_hecheng）打開Ok3、裁剪：

2、用黃石市邊界矢量數(shù)據(jù)裁剪合成后的2000黃石市遙感影像。（1）打開合成后的黃石市2000年遙感影像數(shù)據(jù)的3、4、5波段； FileOpen Image File選圖（2000_band543_hecheng）打開Ok（2） 打開黃石市邊界矢量數(shù)據(jù)；VectorOpen Vector File選圖（黃石市邊界范圍.evf)打開 備注：建立掩膜時(shí)一定要將2000_band543_hecheng和黃石市矢量邊界的影像打開。（3）以黃石市邊界矢量數(shù)據(jù)建立掩膜；Basic ToolsMaskingBulid MaskDisplay #1Options Import EVFS選圖（111）OkChoose

3、命名（2000_band543_hecheng_yanmo)打開Apply（4）應(yīng)用掩膜；Basic ToolsMaskingApply Mask2000_band543_hechengSelect Mask Bang2000_band543_hecheng_yanmoOkOkChoose命名（2000_band543_hecheng_clip)打開Ok （5）備注：若裁剪后Scroll窗口內(nèi)黑色背景面積太大可以進(jìn)行調(diào)整。Basic ToolsResize Data2000_band543_hecheng_clipSpatial SubsetROI/EVF選圖（EVF:111）OkOkOkCh

4、oose命名（2000_band543_hecheng_clip _xiugai)打開Ok2、 地表反射率計(jì)算1、 利用6S模型進(jìn)行大氣輻射校正，從而得到計(jì)算地表反射率的公式和三個(gè)參數(shù)（1） 運(yùn)行6S軟件 （2）打開黃石市遙感圖像頭文件和6S操作手冊(cè)查看要輸入的數(shù)據(jù) （3）輸入數(shù)據(jù)Geometrical conditions (幾何條件)igeom 0-7:7(在6S操作手冊(cè)p1第34頁(yè)對(duì)應(yīng)上TM影像，即Landset對(duì)應(yīng)的數(shù)字為7)時(shí)間、經(jīng)度、緯度（在黃石市頭文件）：時(shí)間均換算成以小時(shí)為單位Atmospheric model（大氣模型）idatm 0-8: 2（中緯度夏季大氣模式）（在6S操

5、作手冊(cè)p1第35頁(yè)）Aerosol model(type)（氣溶膠類型） （在6S操作手冊(cè)p1第36頁(yè)）iaer0-12:3 Aerosol model (concentration)（氣溶膠濃度）the visibility:16 (輸入光學(xué)厚度)The altitude of target （ 目標(biāo)的高度）xps: -0.02 (目標(biāo)高度0.02km)The sensor altitude （ 傳感器的高度）xpp: -1000 (傳感器在衛(wèi)星高度)The spectral conditions （ 波段狀況）iwave25-30:27 （25即TM1,26即TM2，27即TM3，28即T

6、M4，29即TM5,30即TM6。但是tm6是熱紅外波段，不適用于我們的此時(shí)地表反射率的計(jì)算。）Ground reflectance (type)（地面反射率）inhomo: 0 (地面均勻)idirec: 0 (地面無方向影響)igroun: 0 (反射率不隨波長(zhǎng)變化)ro: 0 (地面波段反射率0)Atmospheric correction mode （激活大氣訂正的方式）rapp: -0.2 （表觀反射率)（4）書寫方程完成后將形成的文件重命名為。波段4的操作同上，只是在輸入The spectral conditions （ 波段狀況）時(shí)改為28即可，得到。打開文件，根據(jù)這兩個(gè)公式y(tǒng)=

7、xa*(measured radiance)-xb; acr=y/(1.+xc*y) ,輸入的是影像所測(cè)的輻射亮度Radiance，得到的是反射率影像。將這兩個(gè)公式合并為一個(gè)公式：y=xa*(measured radiance)-xb;acr=y/(1+xc*y)合并后：acr=xa*(measured radiance)-xb/1+xc*xa*(measured radiance)-xb2、 地表反射率計(jì)算打開ENVI軟件Basic ToolBand Math輸入(0.00336*b1-0.04143)/1+0.07179*(0.00336*b1-0.04143)Add to List（若添加成功說明公式正確）Ok選2000_band543_hecheng_clip_xiugai的第三波段Cho