基于ENVI的遥感影像处理实战——以Landsat8数据为例(上)
1. Landsat8数据与ENVI环境准备
第一次接触遥感影像处理时,我被那些专业术语吓得不轻。直到用ENVI处理了Landsat8数据才发现,只要掌握正确方法,遥感图像预处理就像用美图秀秀修照片一样直观。Landsat8这颗卫星每天能拍740张地球照片,每张覆盖185×185公里范围,特别适合区域环境监测。
要获取数据,我推荐从地理空间数据云平台下载。打开浏览器输入网址,在高级检索中选择"Landsat 8 OLI_TIRS"数据集。这里有个实用技巧:时间范围建议选5-6月或9-10月,这两个时段云量少且植被特征明显。去年我做武汉东湖水质分析时,选了6月影像,处理效果比冬季数据清晰30%以上。
下载的压缩包解压后会看到多个文件,关键要找到带"MTL.txt"后缀的元数据文件。这个文本文件就像相片的EXIF信息,记录了传感器参数、拍摄时间等关键数据。有次我误删了这个文件,导致后续处理全部报错,不得不重新下载。
安装ENVI时要注意,建议选择5.3以上版本,这个系列对Landsat8支持最完善。首次启动记得勾选"Extensions"里的"Atmospheric Correction"模块,后面做FLAASH大气校正时会用到。我电脑配置是i7处理器+16G内存,处理单景影像大约需要15分钟,配置低的话可能需要更久。
2. 辐射定标实战详解
打开ENVI时,新手常犯的错误是直接双击图像文件。正确做法是:点击File → Open As → Optical Sensors → Landsat → GeoTIFF with Metadata,然后选择MTL文件。这样能自动识别各波段组合,我第一次操作时没注意这点,结果波段顺序全乱了。
进行2%线性拉伸显示是个好习惯。在图层管理器右键选择"Stretch Type" → "Linear 2%",能让图像细节更突出。去年处理内蒙古草原影像时,没做拉伸前根本看不清植被分布,调整后连放牧路径都清晰可见。
进入辐射定标环节,Radiometric Calibration工具里有几个关键参数:
- 定标类型选Radiance(辐射亮度值)
- 输出格式必须选BIL
- 数据类型选Float
- 点击Apply FLAASH Settings自动填充参数
有次我手误选了反射率输出,导致后续大气校正完全失败。建议输出文件命名时加上"_rad"后缀,比如"Hubei_202006_rad.dat",这样后期管理更方便。定标完成后,可以用Quick Stats工具查看数值范围,正常情况应该在0-100之间。
3. 大气校正关键技巧
FLAASH校正是最易出错的环节,我整理了三个避坑要点:
- 高程数据必须准确:先用GMTED2010.jp2计算区域平均高程,误差超过50米就会影响结果
- 大气模型选择要看拍摄时间和纬度:9月份31°N左右选Tropical模型
- 气溶胶模型根据地表类型选:城市区域用Urban,森林用Rural
参数设置界面有个隐藏技巧:在Multispectral Settings里,把Default改为"Over-Land Retrieval Standard"。去年做鄱阳湖湿地分析时,这个设置让水体分类精度提升了12%。输出文件建议命名为"flaash_result.dat",方便后续步骤调用。
校正完成后,可以对比波段4(近红外)的直方图变化。成功的大气校正会使直方图向左偏移,同时保留更多细节。有次我发现校正后图像反而变模糊,检查发现是气溶胶模型选错了。
4. 图像镶嵌与裁剪实战
当研究区跨越多景影像时,镶嵌就是必备技能。我总结的镶嵌四步法:
- 在Seamless Mosaic工具中添加所有校正后的影像
- 设置Data Ignore Value为0(忽略黑边)
- 颜色校正选Histogram Matching
- 接边线羽化距离设10-15像素
去年拼接三峡库区影像时,我发现接边线自动生成有时会穿过建筑物。这时可以切到Manual模式手动调整,虽然耗时但效果更好。输出时记得选Cubic Convolution重采样方法,比最近邻法平滑得多。
裁剪环节最重要的是ROI文件准备。建议在ArcGIS里先处理好行政边界shp文件,属性表里确保没有非法字符。在Subset Data from ROIs工具中,Mask pixels outside of ROI一定要选Yes,否则会保留黑边。有个项目我忘了设置这项,结果后续分类时多出20%的噪声像元。
5. 常见问题排查手册
遇到报错别慌,我整理了这些年的排错经验:
- 中文路径报错:把整个工程移到英文目录下
- 内存不足:在ENVI首选项里把Cache Size调到最大
- 波段错乱:检查是否通过MTL文件打开数据
- 结果异常:逐步验证每个环节的输出直方图
特别提醒:处理大范围影像时,建议分块处理。去年做全省土地利用分类时,我先把影像按县市裁剪,最后再合并结果,比直接处理全省影像快3倍。ENVI的Batch Processing功能可以自动化这个流程,后续教程会详细介绍。
