《环境小卫星多光谱数据预处理方法》由会员分享，可在线阅读，更多相关《环境小卫星多光谱数据预处理方法（14页珍藏版）》请在金锄头文库上搜索。
1、01Lwdul:SPOT聞NOQid-fllrd询口 rWivwftiVfert-5GetiEyr-|iEREAVHRHDaWFEEQSEfiOi.EJW15ATM(KiZARTDSAT-LDK40ATRDMCDM3 NOAkjFORMOSAPZOISAT-iRifiisEyeHJ-ULBrR&da*-MillpyDighd EEuiinnU9G&b-UU L&ARn-EP STAJft-ietieri： Format?环境小卫星多光谱数据数据预处理【裁剪、配准、FLAASH精确大气校正】第一步：安装环境小卫星数据处理补丁使用 HJ-1 数据读取补丁蚪 envi_hjlalb_tQoL.sa拷
2、贝该 sav 文件到 ENVI 安装目录的 save_add 目录下* 心 Eselis ENVIS. tlassi-c save_add启动 ENVI-File-Open External File 有 HJ-1A1B 选项即为安装成功9沁刃门即火1 Fe | 3icTcc5 Classificrtion rarcfwi 匚He Spedyal Mappen Infi&oe 冃怙口pri VhcIet FihOpen Female Fls-Ope” Ew1e*nJ Ah pn Prniain Fl.EdkENd HeaderGener-ale Tm Elsga Virvirda Rb- *；
3、Impad frem IDL Vdabk- Export Ie D0L Varina CiDmpik ULModuhIDL CUT Ulittinfican 口aadcr|i LiteCh却*eOipMtc号mian ba ScriptLvM-M Suiip So*ipiHeslcrp Mspiaj- 6口EfMVl Quiw MmHE|rrENV1 Lq Man&ge*C1os# 4J AIikPmFvnnEnExit第二步：数据读取和定标主菜单-File-Open External File-HJ-1A1B，打开环境小卫星数据读取补丁。在HJ-1A/1B Tools V3.3面板中，选择C
4、CD，点击Input File输入基于遥感图像的圆柏生长过程研究HJ-1数据84817T文件夹中的.xml文件，点击Output Path设置数据的输出路径，勾选“Calibration “Layer Stacking两个选项，单击Apply按钮。第三步：工程区裁剪由于整景数据范围非常大，本专题工程区只是其中一小部分，在进行大气校正之前，先 将相关研究区域裁剪出来。(1) 打开上一步生成的文件：HJ1A-CCD2-21-72-20120915-L20000848171_Calbrated_LayerStacking(2 ) 主菜单-File-Save File As-ENVI Standard
5、，弹出 New File Builder 面板，( 3) 在 New File Builder 面板中， 单击 Import File， 弹出的 Create New File Input File 面板，( 4) 在 Create New File Input File 面板中，选中 Select Input File 列表中的裁剪数据， 单击 Spatial Subset 按钮，(5 )在 Select Spatial Subset 面板中，单击 Image，弹出 Subset by Image 对话框，（6） 在 Subset by Image 对话框中， 按住鼠标左键拖动图像中的红色矩
6、形框确定裁剪区域，裁剪出相关的区域，单击 OK，（7） 在 Select Spatial Subset 面板中，可以看到裁剪区域信息，单击 OK，（8） 在 Create New File Input File 对话框中，单击 OK，（9 ）在New File Builder，单击Choose设置输出文件名“裁剪”及路径，单击0K。完成相关区域的裁剪。IE：| CaxcfllNev-Ate Builderg Cr-eate Mew Ffe mpyt FilmEHVI Filii Eor FiR.BLrcmt STj-trfliJMis Filts ?DutpuL Ktiul t to Jbp
7、aceInporL JaliirdaT1 ima:. . . | | D.l.b訂Stfiplui D Id 15T2：L tfS 15T2LTull Sue : 06.573.&3q bytei Subtil SdK ffil. 573.631 LylesS-dbsat ty Imggc： D i tplsy 俎 |StfiplL 5G3Cancel Select Spatial Eubt珀h: HJ-CCie-21 -72-201 MG-LSOMla4!Ih苗：LSTFl I 13701 tFloabLne Foint】Li naiIt. 13?0L HL 3T0LSubs*!： Ua咤Q
8、 Subset by Image11 町 | Fi | aiMHEYF | ScicCLlSprctril Eui=t q/iiTdSuuoy Typ*. UnlnioxitBjf-L Order Hd3 L (Znll)FrcijeEtioDi Ol, Irxne 4T Horth find. 3D Utu-Dk-Din VE-B4Vmlmflbi 0. ITS to 0 83VppK!： Ltfl. CriL4t： L,DaacrLp-La au* )U1 A. CCQ2 D kt x, Gr tup =LEVH2, SEKhtDitft =ZD1MW-|.5 03：3B：42.DT, =
9、4tflllLt4ll = HJi. spuwcoTd = CCD2, FmElamtMu = 朝.53Q. durLhtmu.lhILnlLrai = 313 TLO. DB-JIII|第四步：大气校正环境小卫星提供了波谱响应函数，以文本形式提供，第一列表示波长（nm ）,后面 四列分别表示4个波段对应波长的波谱响应值。需要制作波谱曲线来描述波谱响应函数， 用于大气矫正。1)制作波谱曲线(1 )主菜单 WindowTStart New Plot Window ,打开 ENVITPIot Window 面板，在波 谱绘制窗口中，选择FileTInput DataTASCII，导入” HJ1A-
10、CCD2-光谱响应函数.txt”文本 文件如图所示自动将第一列作为X轴 后面4列作为Y轴 波长单位选择nanometers， 单击 OK 。 In put ASCH Fi le立 T *Zi.y-J.-. Flic： C : V:er ack i nD e sktA *l：aLon 皿5. Eovs ： 4114COO.OOE+DO0. OOEI-OOE401O.OOE+DO0. OOEtOO4C2O.OOE+DO0. OOEI-OOE4D3O.OOE+DO0. ODEi-OO404O.OOE+DO0. OOE4-OCIE如一p Zl卜I Ajci s C olunri LSelect Y X
11、jiis CoLwim.5VaveleriEth Uriit生 lanomeT Scile Factor 1- UULUDUCKCaiLLftl(2 )如图，在绘制窗口生成了 4条曲线，选择EditTData Parameters，编辑每条线的名称为b1，b2，b3，b4，便于区分。T一口启0,4Q.2File Edit Options Plot_Function Hdp40)5006007 DO SOO 90010c| EMVI Plot Window(3 )选择 FileTSave Plot AsTSpectral Library，在 Output Plot to Spectral Lib
12、rary 面板 中，单击 Select All Items，单击 OK。(4 )在Output Plot to Spectral Library面板中，有输出曲线相关参数设置，按默认， 选择输出路径和文件名，单击OK ,将波谱曲线保存为波谱库文件：环境1A星CCD2光 谱响应.sli。2)FLAASH 大气矫正(1)数据准备FLAASH 对图像文件有以下几个要求： 数据是经过定标后的辐射亮度(辐射率)数据，单位是：(pW )/( cm2*nm*sr ) 数据带有中心波长(wavelenth )值，如果是高光谱还必须有波段宽度(FWHM ), 这两个参数都可以通过编辑头文件信息输入( Edit
13、Header) 。 数据类型：支持四种数据类型：浮点型(floating )、长整型(long integer )、整 型(integer)和无符号整型(unsigned int)。数据存储类型：ENVI标准栅格格式文件， 且是 BIP 或者 BIL。 波谱范围： 4002500nm。本次用的环境小卫星经过以上处理，已经定标为W*mA(-2)*srA(-1)*umA(-1)单位、 浮点型的辐射率数据，有中心波长信息，下面将 BSQ 格式转成 BIL 格式。选择主菜单Basic ToolsTConvert Data(BSQ,BIL,BIP)，选择已经经过定标和配准的数据 裁剪img，在 Convert File Parameters 中，Output Interleave 选择 BIL，选择 Convert In Place : yes，单击 OK。Q? Convert File Pa meters1召）Inpmt Interlave :
