基于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功能可以自动化这个流程后续教程会详细介绍。