LANDSAT8影像数据岩性制图1、数据准备.PDFVIP

北京天图科技有限公司 LANDSAT8 影像数据岩性制图 1、数据准备 1.1 数据下载 从USGS 官网或者地理空间数据云按照研究区区域范围下载质量好的影像数据。 本文研究原始数据： LC81200302016264LGN00.tar.gz （产品级别L1T ） 1.2 数据导入 Landsat 数据通常都是以压缩包（.tar.gz）形式进行存储，且压缩包内各个数据波段以单 波段影像文件的方式来进行存储，为了影像处理与分析，有时候需要将多个多光谱的单波段 文件组合为一个多波段影像文件。此操作也可通过Import Data 工具来实现。 （1）点击Manage Data | Import Data ，打开如下导入数据的界面： （2 ）在Format 文件类型处选择Landsat-7 or Landsat-8 from USGS，在Input File(*.gz) 处加载tar.gz 压缩包文件，在Output File(*.img)处设置输出IMG 文件路径及文件名，点击 OK 按钮，弹出下面的对话框： 1 北京天图科技有限公司 （3 ）取消勾选第二项至最后一项数据导入，仅合成多光谱波段（1-7 ），点击OK 按钮 完成数据导入。 波段合成后影像：lc81200302016264lgn00-msi.img 1.3 数据查看 （1）在视图窗口点击右键，点击Open Raster Layer 加载波段合成后的多光谱影像， 进行影像浏览； 2 北京天图科技有限公司 （2 ）查看影像元数据信息 1）使用文本编辑器（UltraEdit 或写字板）打开元数据文件： LC81200302016264LGN00_MTL.txt 2 ）查看基本信息：影像获取时间、太阳高度角、辐射定标参数 （增益、偏移） 1.4.cal 校准文件制作 基于影像元数据中查询到的辐射定标系数制作 ATCOR 模块进行大气校正所需的辐射 定标文件： （1）在C:\Program Files \Hexagon \ERDAS IMAGINE 2015 \etc\atcor\cal\landsat8_oli 路 径下打开模板定标文件landsat8_oli_template.cal ，使用文本编辑器打开； （2 ）由于元数据文件中辐射定标系数单位为W m2 sr1 m1 ，与.cal 文件中的默认 单位有0.1 的换算关系，需将辐射定标系数文件中的数字按照下面的公式计算后填入。 （3 ）另存文件为且放在同一个安装路径下。 3 北京天图科技有限公司 1.5 数据裁切 由于ATCOR 模块去云去雾仅支持8000*8000 以下的数据量，且有时数据量过大会出现 算法失败，因此本文对原始影像（7831*7951 ）进行裁切，仅选取右下角一小部分包含云雾 的区域进行示例。 （1）在视图窗口点击右键，点击Open Raster Layer 加载波段合成后的多光谱影像； （2 ）点击Raster | Drawing | Polygon 绘制.aoi 范围文件； （3 ）点击Raster | SubsetChip | Create Subset Image 图标，在弹出的Subset 对话框 中点击AOI…按钮，打开Choose AOI 对话框，选择Viewer 选项加载视图窗口已绘制的 AOI 文件，点击OK 按钮关闭Choose AOI 对话框，点击OK 按钮实现裁切。 裁