三亩地 三亩地SAN MU DI · CODE DIARY
ARTICLE DETAIL

日记详情

真实记录编程学习的某一天,欢迎挑你感兴趣的翻一翻。

ENVI 5.6 保姆级教程:手把手教你处理 Landsat 8 遥感影像(从下载到预处理)

ENVI 5.6 保姆级教程:手把手教你处理 Landsat 8 遥感影像(从下载到预处理)

ENVI 5.6 零基础实战:Landsat 8 遥感影像处理全流程详解

1. 遥感影像处理入门基础

对于刚接触遥感影像处理的初学者来说,Landsat 8数据是一个理想的起点。作为美国地质调查局(USGS)和NASA联合运营的卫星,Landsat 8提供了免费且高质量的地球观测数据,覆盖了从可见光到热红外的多个波段。ENVI 5.6作为业界领先的遥感影像处理软件,其直观的界面和强大的功能使其成为学习和专业工作的首选工具。

在开始实际操作前,我们需要了解几个核心概念:

  • 空间分辨率:Landsat 8的多光谱波段空间分辨率为30米,全色波段为15米
  • 时间分辨率:Landsat 8每16天对同一地区进行一次观测
  • 光谱特性:包含11个波段,每个波段对地物特征的响应不同

提示:处理遥感影像时,建议使用英文路径和文件名,避免软件兼容性问题

2. 数据获取与准备

2.1 Landsat 8 数据下载

获取Landsat 8数据有多种途径,以下是最常用的三种方式:

  1. USGS EarthExplorer:官方数据源,提供最完整的数据集
  2. 地理空间数据云:中文界面,下载速度较快
  3. Google Earth Engine:适合批量处理和大范围数据获取

以地理空间数据云为例,详细下载步骤如下:

  1. 访问官网并注册账号
  2. 在"高级检索"中选择"Landsat 8 OLI_TIRS 卫星数字产品"
  3. 设置区域范围(可通过行政区或手动绘制)
  4. 选择时间范围(建议云量低于10%)
  5. 点击检索并选择合适的数据进行下载

2.2 数据解压与检查

下载完成后,通常会得到一个压缩包,解压后包含以下关键文件:

文件类型说明重要性
MTL.txt元数据文件★★★★★
B1-B11各波段数据★★★★★
QA质量评估波段★★★★
# 示例解压命令(Windows) tar -xzf LC08_L1TP_123032_20201020_20201105_01_T1.tar.gz

3. ENVI 5.6 基础操作

3.1 软件界面介绍

首次打开ENVI 5.6,主界面分为以下几个主要区域:

  • 工具栏:包含常用工具和功能快捷方式
  • 图层管理器:显示已加载的数据集
  • 视图窗口:影像显示区域
  • 工具箱:所有处理功能的入口

3.2 数据加载与显示

加载Landsat 8数据的正确方法:

  1. 点击"File" → "Open As" → "Optical Sensors" → "Landsat" → "GeoTIFF with Metadata"
  2. 选择MTL.txt文件
  3. 等待ENVI自动识别并加载所有波段

注意:直接加载MTL文件而非单个波段文件,可以确保ENVI正确识别波段关系和元数据

4. 预处理流程详解

4.1 辐射定标

辐射定标是将原始DN值转换为物理量的关键步骤。在ENVI中执行辐射定标:

  1. 打开"Radiometric Correction" → "Radiometric Calibration"
  2. 选择多光谱数据
  3. 设置参数:
    • 定标类型:Radiance
    • 输出格式:BIL(FLAASH兼容)
    • 数据类型:Float
  4. 点击"Apply FLAASH Settings"应用推荐参数
  5. 指定输出路径并运行
# 伪代码展示辐射定标计算过程 def radiometric_calibration(DN, gain, bias): """ DN: 原始数字值 gain: 增益系数(来自MTL文件) bias: 偏移量(来自MTL文件) """ return DN * gain + bias

4.2 大气校正

FLAASH是ENVI中最常用的大气校正工具,操作步骤如下:

  1. 计算区域平均高程(使用ENVI自带DEM数据)
  2. 打开"Radiometric Correction" → "FLAASH Atmospheric Correction"
  3. 关键参数设置:
    • 传感器类型:Landsat-8 OLI
    • 地面高程:从DEM计算获得
    • 大气模型:根据成像时间和纬度选择
    • 气溶胶模型:根据地区特性选择
  4. 运行校正并保存结果

大气模型选择参考表:

季节纬度范围推荐模型
冬季<23°Tropical
23°-45°Mid-Latitude Winter
>45°Sub-Arctic Winter
夏季<23°Tropical
23°-45°Mid-Latitude Summer
>45°Sub-Arctic Summer

4.3 图像镶嵌与裁剪

当研究区域跨越多个影像时,需要进行镶嵌处理:

  1. 打开"Mosaicking" → "Seamless Mosaic"
  2. 添加所有待镶嵌影像
  3. 设置颜色校正(直方图匹配)
  4. 生成接边线并设置羽化距离
  5. 选择输出参数并运行

裁剪操作则更为简单:

  1. 加载研究区域的矢量边界文件(如.shp)
  2. 使用"Subset Data from ROIs"工具
  3. 选择输入影像和边界文件
  4. 设置输出参数并执行

5. 常见问题与解决方案

在实际操作中,新手常会遇到以下问题:

  • 问题1:ENVI无法打开MTL文件

    • 检查文件路径是否包含中文
    • 确认下载的数据完整无损
  • 问题2:FLAASH运行失败

    • 确保输入数据是辐射亮度值(Radiance)
    • 检查高程单位是否正确(km)
    • 确认输出路径有写入权限
  • 问题3:镶嵌后色彩不一致

    • 尝试不同的直方图匹配方法
    • 调整羽化距离参数
    • 考虑手动设置接边线

6. 进阶技巧与优化建议

经过多次项目实践,我发现以下技巧能显著提高处理效率:

  1. 批量处理:使用ENVI的批处理功能或编写IDL脚本自动化重复任务
  2. 内存管理:对于大区域数据,适当增加ENVI内存分配(默认设置→首选项→内存)
  3. 质量控制:在处理每个步骤后,使用"Quick Stats"工具检查数据范围是否合理
  4. 元数据记录:使用ENVI的"Edit Header"功能添加处理历史记录

对于需要更高精度的项目,建议考虑:

  • 使用地面控制点进行几何精校正
  • 结合更高分辨率数据(如Sentinel-2)进行融合
  • 尝试不同的辐射传输模型(如6S)进行大气校正
← 返回列表