Landsat简介及数据预处理

1、Landsat 8 简介及数据预处理OLI 陆地成像仪包括 9 个波段，空间分辨率为 30 米，其中包括一个 15 米的全色波段，成像宽幅为 185x185km。OLI 包括了 ETM+传感器所有的波段，为了避免大气吸收特征，OLI 对波段进行了重新调整，比较大的调整是 OLI Band5(0.8450.885 m)，排除了 0.825m 处水汽吸收特征；OLI 全色波段 Band8波段范围较窄，这种方式可以在全色图像上更好区分植被和无植被特征；此外，还有两个新增的波段：蓝色波段 (band 1; 0.4330.453 m) 主要应用海岸带观测，短波红外波段(band 9; 1.3601.390 m) 包括水汽强吸收特征可用于云检测；近红外 band5 和短波红外band9 与 MODIS 对应的波段接近。OLI 陆地成像仪 TM序号 波段 波段（um） 空间分辨率(m)波段（um） 空间分辨率(m)1 蓝色波段 0.4330.453 302 Blue 0.4500.515 30 0.4500.52 303 Green 0.5250.600 30 0.520.60 304 Red 0.

2、6300.680 30 0.630.69 305 Near IR 0.8450.885 30 0.760.90 306 SWIR 1.5601.660 30 1.551.75 307 2.1002.300 30 10.4012.50（热红外） 1208 全色 0.5000.680 15 2.082.35 309 短波红外 1.3601.390 3010 中心波长 10.9 微米10011 中心波长 12.0 微米100Landsat8 数据打开和辐射定标处理美国的 USGS（http:/glovis.usgs.gov/）网站提供最新的 Landsat8 数据下载，产品类型标示 L1GT，与之前的数据格式类似，每个波段以.tif 文件提供，元数据存放在_MTL.txt 文件中。Landsat8 增加了几个波段，详细信息浏览：http:/ ENVI5.0SP3 中非常容易打开 Landsat8 数据，如下：(1) 选择 File-Open ，选择_MTL.txt 文件打开。(2) ENVI 自动显示 RGB 显示真彩色图像，打开 Data Manager 对话框，可以看到ENVI 自动读取

3、元数据信息，包括中心波长信息、波段名称等。并将数据根据类型自动划分为三类。(3) 从文件信息中可以看到，热红外数据被重采样为 30 米分辨率，与可见光-近红外波段一致，全色为 15 米分辨率。图 1：Data Manager 对话框打开之后就可以很方便的进行其他处理，比如辐射定标、大气校正、融合等处理。下面使用 ENVI 下的通用定标工具进行 Landsat8 的辐射定标。(1) 选择 ToolBox/Radiometric Correction/Radiometric Calibration，选择可见光-近红外数据。(2) 在 Radiometric Calibration 面板中，可以选择定标类型：辐射亮度值和大气表观反射率。(3) 其他选项是方便用于 FLAASH 大气校正。(4) 选择文件名和路径输出(5) 如图 3 所示，得到大气表观反射率数据。图 2：Radiometric Calibration 面板图 3：大气表观反射率结果ENVI 下的 Landsat8 大气校正（初试）Landsat8 OLI 陆地成像仪比之前的 TM/ETM+多了两个波段， 0.4330.453

4、m 和1.3601.390 m ，怎么多的波段对于地表反演更加有利。ENVI5.1 直接支持Landsat8 的大气校正（2013 下半年发布），利用 ENVI5.1 提供的 Landsat8 波谱响应函数在 ENVI5.0SP3 下也能完成大气校正。大气校正之前，启动 ENVI Classic ，设置 preferences-Miscellaneous:Cache Size：2048( 最大内存 75%)Image Tile Size ： 100 （推荐 1-4M ）注：电脑内存为 8g ， 64 位操作系统保存后重启 ENVI5 。波谱响应函数文件下载： http:/ OLI 和 TIRS 两个传感器第一步：辐射定标选择 File-Open ，选择 _MTL.txt 文件打开。(2)ENVI 自动显示 RGB 显示真彩色图像，打开 Data Manager 对话框，可以看到ENVI 自动读取元数据信息，包括中心波长信息、波段名称等。并将数据根据类型自动划分为三类。(1)选择 ToolBox/Radiometric Correction/Radiometric Calibration

5、 ，选择可见光 -近红外数据。(2)在 Radiometric Calibration 面板中定标类型（ Calibration Type ）：辐射亮度值（ Radiance ）输出储存顺序 （ Output Interleave ）： BIL输出数据类型： Float单击 FLAASH Settings 按钮，自动获取辐射亮度单位转换系数 Scale Factor ：0.1其他选项是方便用于 FLAASH 大气校正。(3) 选择文件名和路径输出图 1 ： Radiometric Calibration 面板第二步： FLAASH 大气校正选择 Toolbox/Radiometric Correction/Atmospheric Correction Module/FLAASH Atmospheric Correction ，打开 FLAASH 大气校正工具。（1）文件输入与输出信息项目单击 Input Radiance Image 按钮，选择上一步准备好的辐射亮度值数据LC81230322013132LGN02_rad.dat 。在 Radiance Scale Factors 对话

6、框中选择 Use single scale factor for all bands （ Single scale factor ： 1 ），在辐射定标中对单位进行了转换。单击 Output Reflectance File 按钮选择输出文件名和路径。（2）传感器与图像目标信息l Lat ： 40 19 39.46 ， Lon ： 116 42 2.98 （ FLAASH 自动获取）l Sensor Type ： UNKONWN-MSIl Ground Elevation （ km ）： 0.043 （从相应区域的 DEM 获得平均值）l Flight Date ： 2013-05-12 Flight Time ： 02:55:26注：在右边图层管理器中， 单击右键选择 View Metadata ，在 Metadata viewer中浏览 time 可以看到飞行时间图 2 ：图像成像时间查看（3） 大气模型（ Atmospheric Model ）： Sub-Arctic Summer （ 5 月份 纬度： 40-50 ）（4） 气溶胶模型（ Aerosol Model ）： Urb

7、an（5） 气溶胶反演（ Aerosol Retrieval ）： 2-Band （ K-T ）（6） 初始能见度（ Initial Visibility ）： 40 。图 3 ： FLAASH 基本参数设置多光谱设置（ Multispectral Settings ）l Defaults 下拉框： Over-Land Retrieval Standard （ 660 ： 2100 ）。l Filter Function File ：选择 ldcm_oli.sli 波谱响应文件图 4 ：多光谱设置（8）高级设置（ Advanced Settings ）：tile 设置为 100M ，其余按照默认设置。（9）单击 Apply 按钮，执行 FLAASH 。图 5 ：估算能见度、水汽柱结果第三步：浏览结果打开大气校正结果，浏览植被波谱曲线如下，大致可以看出大气校正后消除了大气散射的影响。图 6 ：大气校正后的植被波谱曲线同时发现 1.3601.390 m 波段数据大气校正之后结果全部为 0 ，单独打开这个波段的原始文件 LC81230322013132LGN02_B9.TIF 或者在 ENV

8、I 中标识为Cirrus （1.3730 ）波段，发现这个波段的图像噪声非常大，主要用于识别卷云，辐射定标后的值在集中在 0.011697- 0.023395 。图 7 ：卷云 Cirrus （ 1.3730 ）波段图像讨论 ldcm_oli.sli 波谱响应文件包括了可见光 - 红外，全色 9 个波段的响应，而我们大气校正使用了前 8 个波段，另外卷云 Cirrus 波段噪声比较大，是否将Cirrus 、Pan 两个波段的响应函数去除，同时将两个波段从图像中移除，只针对 7 个波段进行大气校正，精度是否会更高？图 8 ： oli 波谱响应函数2、 TM 数据辐射定标ENVI basic tools preprocessing calibration utilities Landsat calibration，弹出如下对话框，图 3：图 3 辐射定标参数设置对话框3、 储存顺序调整Flassh 大气校正对于波段存储的要求为：BIL，BIP 格式，上述计算得到的存储方式为 BSQ，在此进行波段存储顺序的转化，具体操作如下：ENVI basic tools convert data (BS

9、Q ,BIL ,BIP)图 4 存放顺序转换4、 Flaash 校正参数设置大气校正的前期准备工作完毕，现在进行校正参数的设置：ENVI basic tools preprocessing calibration utilities FLAASH，弹出对话图 5 FLAASH 参数界面设置图 6 多光谱设置对对话框根据上述图中的参数设置，然后点击 ok，运行 flaash 大气校正。Landsat8 移除卷云 Cirrus 波段的大气校正测试在“ ENVI 下的 Landsat8 大气校正（初试 ）”文章最后提出了一个讨论，这里根据这个设想做出了另外一个结果。结果分析显示，两种方法得到的结果基本一致。第一步：重新制作波谱响应函数（1）启动 ENVI classic ，选择 Window-Start New Plot Window 。（2）在 ENVI Plot Window 窗口中，选择 File-Input Data Spectral Library，打开 ldcm_oli.sli 波谱响应文件。（3）如下图所示选择 7 个波段的波谱响应函数，选择 File-Save plot as- Spectral Library ，按照默认参数保存为 .sli 文件。 图 1 ：选择 7 个波段的波谱响应文件第二步：大气校正（1）使用 layer stacking 工具将辐射定标后的文件保存为 7 个波段的文件，也就是去除卷云 Cirrus 波段。（2） 打开 FLAASH 工具进行大气校正。第三步：浏览结果分别对两种结果进行统计，如下图所示，每个波段的均值和方差相差非常小（个位数以内），折合 01 反射率在 10 -3 范

《Landsat简介及数据预处理》由会员平***分享，可在线阅读，更多相关《Landsat简介及数据预处理》请在金锄头文库上搜索。