地表反射率計(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波段。File Open Image File 選擇黃石市2000年遙感影像數(shù)據(jù)的3、4、5 波段t打開。Q Available Bdnd, -File Options79LI a 7211 (Jytft) BQD im =2、合成：對(duì)黃石市2000年遙感影像數(shù)據(jù)的3、4、5波段進(jìn)行合成Basic Tools Layer Stacking Import File 選擇

2、黃石市 2000年遙感影像數(shù)據(jù)的 3、4、5 波段01 Choose命名(2000_band543_hecheng 打開OkZont IJI * K 5 Sti Zcnt.X Pl JI el S 1 Ke -U. JU .U JU1.11.1It Lar ET Pi k1 Sir*i30 OOCO.IOCIOLayer亡Idng ParametersITiiits .Hwe 1 11FiLt. mt習(xí)Vfriin二TiftC有也巧市邊界芯&0L, 432 bylsEacuj dib.1SJl疋 t. A1.JL !_D vl w-el A2A L暗Dmtput Proj ecti on JT

3、ew,.,Art i tr aryGmbgraphic Lqt/LorLStttiWAD ZT，St 亀t曽 f Itfrt (HAD 03) Ar tent - 工 m* j. Ar(entint - Znut Z Ar耳*n七in豈- 2ojta 3Naarest Hal ghbor S3)Arnii- 2on*1ArctxitintZomEAr*ntin7orn3D譏則wITniti. F.Zikt 別=車 HSZcikt.Outpiit Resml t toFil電MemoryEnt er O utput Filsnun Chocisft1 Pik1 Six Ineluiva rtil

4、 th* fil aJBLEmcIusiv rUlfompLists fil 口#rlp 30 00000000MtY Pixtl Silt 0 00000000MlariEieaJftpliAg Naarest N&i ghbor v Scroll (0.03ZS6) 一(5) 備注：若裁剪后Scroll窗口內(nèi)黑色背景面積太大可以進(jìn)行調(diào)整Basic Tools 宀Resize Data2000_band543_hecheng_clip 宀SpatialSubsetROI/EVi選圖(EVF:111) OkOkOkChoose命名(2000_ba nd543_heche ng_clip _xi

5、ugai) 打開Ok鬥 Scroll (0.08132) 一2、地表反射率計(jì)算1、利用6S模型進(jìn)行大氣輻射校正，從而得到計(jì)算地表反射率的公式和三個(gè) 參數(shù)(1)運(yùn)行6S軟件，津們応，好燦乍I AST垂躊 I壬西師密*飪2 - 65斟6S.EXE(2) 打開黃石市遙感圖像頭文件和 6S操作手冊(cè)查看要輸入的數(shù)據(jù)LniZZOaOOUOBJCOOMTLttt TMO. )(3) 輸入數(shù)據(jù)】F:PrOgr arm F&_eacuteab16.!EJCuj 即ICEOM Is In t bv ranv 陽” gwin - 7IM LAHDSflT)=Month - 18cu CunlMersdl t ine

6、hh.ddd - 2.43 m1qi = 1142 xlat ( I at i tilde of= 29 T1frTHOSPHERC MODEL Idatn- 2ftEROSCL HODEiL (TVPE.laer ,v us= ILXFS i_H to Exprc sE the alt itudr 0f target. xps in lm - 一XPP ik to eMpvflsc thft cencnk1 alttiuda. xpp -10W:虹ER THE EPECTBflL COHDITJOHS：iwdiAe Cl to 34* see intraduct Lon for 27CMOU

7、NHD R1FLKCTAHC1 iMBH flH WKHMHKMIHH Hi HBJHUUMiMLHHInhonu C0-1 - Gidiret: Iditec 二 9* no ireclon-A 1 effete ； iflkrec = L P h4sionaL effect idirec = fl lOroun = a ro flJMPP is a pn*Heter 11de to 4Ct lv4te tnosphvc Curre-c t Ion Rode +即h 7L JGeometrical conditions （幾何條件）igeom 0-7:7（ 在6S操作手冊(cè)pl第34頁對(duì)應(yīng)上T

8、M影像，即Landset對(duì)應(yīng)的數(shù)字 為7）筋r-lsr CE.皆mc-a忖 F邑二筆,GEOMETRICAL CONDIIIONS：I GtOH is in the t*ange 0-7, igeoti） = 7時(shí)間、經(jīng)度、緯度（在黃石市頭文件）：時(shí)間均換算成以小時(shí)為單位ELEV1TION_SOURCE 工 PLE2000 0UTPUT_F0RHA7 = GEOTIFF EFHEMERIS_.TFE 二DEFINITIVE*SPACECRAFT. ID = LAMBSATS SENSOR.ID 二 *TI*SENSOR_MOLE = SAM*VRS-PATH = 122陽 S_RW 匸 039

9、|DATEACXUIRED - 200006-18 SCBKE_CEWTER_TIIE = C2：26：3287304402 CORNER UL 丄萄PRODUCT - 31.26164CORNER UL LON PRODUCT 二 114,C7111bblionJ*CORKERUR丄占T_PRODUCT = 31.曲428 COR HER _UR _LOW_PRODnCT = 116.56290 CORNER丄L丄AT_PRODUCT = 2- 3i221 CX)RNER_Ll_LOri_FKODUCT 二 114.12851 CORNER LR LAT PRODUCT 二 29. 3424

10、2MH -IMCORNER_LR_LOK_PKODUCT 二 116-5714STM ；month nonth- 6Jday (lay) - 18tu 2-43xlon (lony itiLde of 114.2model （大氣模型）Atmosphericidatm 0-8:l-at = 29_12 （中緯度夏季大氣模式）（在6S操作手冊(cè)pl第35頁）NZ 9fiTMOPH ERIC MODEL:Aerosol model（type）（氣溶膠類型）（在6S操作手冊(cè)pl第36頁）iaer0-12:3 3L JrHuriZ o ncrAEFOSCL MODEL .ider （8-12= 3 Ae

11、rosol model （c oncen trati on）（氣溶膠濃度）the visibility:16（輸入光學(xué) 厚度）ftEROSOL MODEL COHCENTRATION?u -The altitudeof target （目標(biāo)的高度）xps: -0.02（目標(biāo)高度=0.02km）XPS is Cq xprcss Vhc= altitude of target. xps （. r km e 0 ,The sensor altitude （傳感器的高度）xpp: -1000（傳感器在衛(wèi)星高度）XPP s：七 expiress it he sensiOk* alt; liude B x

12、pp -1H00 The spectralconditions（波段狀況）iwave25-30:27（25即TM1,26即叮M2 27即TM3 28即TM4 29即TM5,30即TM6但是tm6是熱紅外波段，不適用于我們的此時(shí)地表反射率的計(jì)算。）ENTEK THE GFECTRfiL CONDITIONS：IviaUeto M4. see xnicclucit lun For de t-a ils） = 2V- 30 Grou nd reflecta nee （type）（地面反射率） in homo: 0（地面均勻）idirec: 0（地面無方向影響）igroun: 0（反射率不隨波長變 化

13、）ro: 0（地面波段反射率=0）GRQUNRD REFLECTfiNCE inhorw n 以idii-ec : idlrecn d dit-nct Lonl aFFec t : idir-ec = 1R （ms dHrct ioral eF f ctdii-AC = B= H）*0 = Atmospheric correction mode （激活大氣訂正的方式）rapp: -0.2（表觀反射率）BAPP is a paranet er a. 1 lows to act iute atnosphec correction node. rapp -H.2 I I （4）書寫方程完成后將形成的

14、文件重命名為波段4的操作同上，只是在輸入The spectral conditions （波段狀況）時(shí) 改為 28即可，得到 re#lertflncDUt 。打開慳文件，根據(jù)這兩個(gè)公式y(tǒng)=xa*（measured radiance）- xb; acr=y/（1.+xc*y）,輸入的是影像所測的輻射亮度Radianee，得到的是反射率影像。3tm3pheric correction resultInput apparent Teflectancemeasured radiance v/iL2/sr/*lcatiosphericil ly correct cd reflect aneBS. 509

15、.24icoefficients xa xb xc*y=a* (Aeaeured. radianceJ-tfe; acr=y/(1+ +ec*y).QC33E：,041434* * 扌*+ * * 4*4+!：+ + + 卓卄* +4+W+t+ +4+4+ * 卓 +r+t+ +4+將這兩個(gè)公式合 并為一個(gè)公式:y=xa*(measured radiance)-xb;acr=y/(1+xc*y)合并后：acr=xa*(measured radia nce)-xb/1+xc*xa*(measuredradia nce)-xb2、地表反射率計(jì)算打開 ENV 軟件Basic Tool Ba nd M

16、ath輸入(0.00336*b1-0.04143)/1+0.07179*(0.00336*b1-0.04143)AddtoList （若添加成功說明公式正確）Ok選2000_band543_hecheng_clip_xiugai 的第三波段Choose命名為 tm3_M表 反射率打開OkBand MathPrenous Band lath Expressions:(0.00336*bl-0 O4143)/1W OT17 (0* 00336*11-0 04143SaveClearDelateEnter s rxpression!00336*bl-0. 04143)/HO 07173 00336

17、1-0 04143)Add to ListOKCancel HelpVariables to Bands PairirgsExp: (0.00336*bl-0.04143)/l+0.07179*(0.00336*bl-0.Vari ables used in expression：Bl 辰si (H&sk(Band 1 :LT5122039Z000170BJCOO_B：3le Bands Li st)0 _b and543 _h 邑 c h 電cl i P_x i u gaiResize (Mask (Layer (Band 1:LT51220392000170BJC00J5. TI：Resi

18、ze (Mask (Lajr (Band 1 ：LT51220392000170BJC00_B4. TI：Reside(Layer (Bmd 1 :LT5122039200017OB兀00_B3 TIMap InfoMap Vtri a*ble to Input FileSpatial SubsetFtjII SceneOutput Result to (:FilelikmryEnter Output Filename Choose CompressOKQueue CancelJtelp Cl earOn Band Math (0+00336*b1-0.04143)-. - 口4Fife Overlay Enhance Tools Window| File3、相同步驟操作最后得到tm4地表反射率 波段4的數(shù)據(jù)aiijs：