ASTER数据处理.docxVIP

ASTER数据 ?ASTER 数据有14个波段组成，1、2、3波段为可见光/近红外波段，3B 波段为后视成像波段，是卫星飞过去了几十秒后对先前垂直成像区域的重新成像，3N 波段 3B 波段是一样的波段范围成像，只不过 3N是垂直成像3B为后视成像，3N 波段与3B波段组成立体像对用于 ASTER 立体测图生成 DEM。用于制作调绘片我们采用 1、2、3 波段。1、2、3波段的空间分辨率为 15m，4-9 波段为短波红外波段空间分辨率为30m，10-14 波段为热红外波段空间分辨率为90m 。? 在 ENVI 里处理相对比较方便，可以直接读取，然后选择自己想要的波段另存文件。 利用ENVI的ASTER数据处理 1.用ENVI打开ASTER数据，观察数据取值，如果是整形则进入下一步；否则，若数据出现浮点型，那么在ENVI主菜单File-Preferences-Miscellaneous当中把Auto-Correct ASTER/MODIS选项改为No，然后重新打开，ASTER数据就变成了整形。ENVI自动对这两种数据进行辐射转换，但是参数获得来源不明，还是自己手动设置放心一些。如果在浮点型数据基础上再做一次辐射转换，后果可想而知（本人阵亡于此） 2.利用CrossTalk 3.0做串扰纠正； 3.将红外波段30米分辨率的数据向可见光-近红外配准，重采样 4.为每个波段输入波长、增益、偏差； 5.应用增益偏差参数进行辐射转换； 6.大气纠正（FLAASH） 7.裁剪 FLAASH模块的大气校正1.1　FLAASH模块简介FLAASH是由世界一流的光学成像研究所－波谱科学研究所（Spectral Sciences）在美国空气动力实验室支持下开发的大气校正模块。波谱科学研究所在1989年大气辐射传输模型开发初期就广泛从事MODTRAN的研究工作，已成为大气辐射传输模型开发过程中不可缺少的一员。FLAASH适用于高光谱遥感数据（如HyMap，AVIRIS，HYIDCE，HYPERION，Probe-1，CASI和AISA）和多光谱遥感数据（如陆地资源卫星，SPOT，IRS和ASTER）的大气校正。当遥感数据中包含合适的波段时，用FLAASH还可以反演水气、气溶胶等参数。ENVI中大气校正模型FLAASH，是高光谱辐射能量影像反射率反演的首选大气校正模型。FLAASH能够精确补偿大气影响，其适用的波长范围包括可见光至近红外及短波红外，最大波长范围为3μm。其他的大气校正模型是计算方法基于查找表（Look-up Table）、利用插值方法计算，而FLAASH是直接移植了modtran4中的辐射传输计算方法。用户可以选取代表研究区的大气模型和气溶胶类型，并且对每景影像，Modtran都有独特的解决方案。1.2　ASTER数据预处理ASTER　L1B数据是记录是DN（Digital Number）值，而基于FLAASH大气校正过程中，需要的是辐射能量值。因此，需要对ASTER　L1B数据辐射定标，即把无量纲的DN值转换成有量纲的分辐辐射亮度值的过程（式1），Radiance=gain*DN+offset　　　　　　　　　　　　　（式1）其中，gain是增益,offset是偏差。经辐射定标后，得到天顶辐射能量值，其量纲为W/(m2.sr.um)。ASTER数据多以HDF格式储存，利用ENVI软件中Baisc Tools-Preprocessing-Data-Specific Utilities-View HDF Global Attribute功能，读取相应ASTER　HDF文件中的增益、偏差、成像时间和中心点坐标信息，增益和偏差信息见表1。根据ASTER数据特点［6］，需要对短波红外波段数据重采样，以使其与可见光－近红外波段影像的像元数相同。然后，把可见光－近红外波段与短波红外段按波段顺序合并成新的文件。利用ENVI软件，编辑新文件的头文件。主要输入波长值见文献6、波长单位、增益、偏差和中心波长半极值宽度（FWHM，Full Width Half Maximum）。中心波长半极值宽度可以近似用每个波段的带宽代替，如果波长单位为微米，应该转成纳米单位。然后，利用ENVI软件中Baisc Tools-Preprocessing-General-Purpose Utilities-Apply Gain and Offset功能下，生成辐射能量数据。由于ENVI软件保存的数据格式通常为波段顺序格式（BSQ），而FLAASH大气校正模块使用的是波段逐行交叉顺序（BIL）或波段逐像元交叉顺序（BIP）格式的文件，因此，ASTER数据预处理的最后一步是把BSQ格式的文件转换成BIL或BIP格式的文件.FLAASH模块参数设置1.3.1　尺