landsat8数据预处理

最新推荐文章于 2024-04-02 17:01:51 发布

今天也要加油吖！

最新推荐文章于 2024-04-02 17:01:51 发布

阅读量4.4k

点赞数 24

文章标签：学习方法课程设计笔记经验分享

本文链接：https://blog.csdn.net/YYQX212/article/details/134926800

版权

记录学习

一、数据介绍

二、数据下载

1.地理空间数据云（地理空间数据云 (gscloud.cn)）

2.USGS官网（https://earthexplorer.usgs.gov/EarthExplorer (usgs.gov)）

三、预处理

Landsat 8-OLI 数据是在地理空间数据云里下载的，这些数据是已经几何校正过的

融合，辐射校正、大气校正、裁剪

3.1影像融合

将8波段30 m分辨率的多光谱数据和15 m的全色数据进行融合，从而得到15 m分辨率的多光谱数据

1、启动ENVI 5.3，选择File——Open，选择MTL.txt文件打开；选择“B8.TIF”文件打开。

2、在工具箱（Toolbox）中，点击Image Sharpening，双击Gram-Schmidt Pan Sharpening打开对话框。在对话框中选择低分辨率的多光谱数据文件，点击OK，再选择高分辨全色波段数据文件，点击OK。

3、在Pan Sharpening Parameters参数面板如图下图，选择传感器类型为：Landsat8_oli，重采样方法选择Cubic Convolution，设置输出路径和文件名。

3.2辐射定标

辐射定标是将DN值转换为辐射亮度值、地表反射率等相关物理量的相对值的过程

第一步：在ENVI 5.3的工具箱中查找工具：Radiometric Correction——Radiometric Calibration双击，选择要校正的多谱数据（而本次试验选择以上融合的多光谱），“LC81200332020291LGN00_rh.dat”进行辐射定标。

第二步：选择参数对话框中点击：【Apply Flaash Setting】后会自动读取元数据中的正射参数，也可以自己输入参数，主要参数如下图：选择多光谱数据。定标类型选择辐射亮度值（Radiance），输出格式为BIL，因为后续FLAASH大气校正的输入数据类型为BIL。输出数据类型为Float，系数为0.1。或者直接点击应用FLAASH设置按钮设置以上参数。确定输出路径后点击确定。

3.3大气校正

打开Radiometric Correction>Atmospheric Correction Module>FLAASH Atmospheric Correction工具，导入经过辐射定标的数据，因为在辐射定标中已经进行了单位换算，所以在弹出的Radiance Scale Factors中选择第二项

Lat和Lon都是自动获取的不用管。

传感器类型（Sensor Type）：MultiSpectral（多光谱数据）》Landsat8 OLI。

气溶胶模型：渭-库绿洲（除6-9月用MLS（中纬度夏季），其他用SAS

查看成像时间（Flight Date）

方法一：在Layer Manager(图层管理器）窗口右击原始影像，按照提示选择View Metadata(查看元数据），点击Time查看成像时间

方法二：从MTL文件中找到。具体名称：DATE_ACQUIRED = 2013-05-12；SCENE_CENTER_TIME = 02:55:26.6336980Z；

查看平均高程值：

方法一：由于FLAASH大气校正需要影像区域的平均高程，所以可以使用ENVI自带全球高程数据进行计算。选择File》Open World Data》Elevation，打开ENVI自带的全球900mDEM数据（GMTED2010.jp2数据）。再将需要计算高程数据的影像打开（需要带坐标信息），可以放大看到影像数据叠加在DEM数据上。打开工具箱（Toolbox）中的Statistics——Compute Statistics，在输入文件对话框中选择GMTED2010.jp2数据，再单击Stats Subset，单击File，选择需要统计高程信息对应的图像，然后OK。之后的参数均设置为默认。最后得到统计的平均高程信息，可以看到影像平均高程为1690.69m。