【ENVI课业练习记录】-EW帮帮网

记录第一次ENVI实操（辐射校正+大气校正）

一、研究区域裁剪

首先我们要有两期不同年份同一地区的遥感影像（可以通过这里选取数据集进行下载）

根据个人喜好我首先会选择先进行影像的裁剪，尽可能地将图像放在研究区域内

“File”中打开影像
“File”——“save as”——“save as …（ENVI）”
在弹出来的窗口中选择头文件影像
点击右侧的箭头进行选择研究区域的裁剪

选取完成后选择OK 完成裁剪文件的命名（不要有中文）就完成裁剪啦

二、辐射定标

这一步我是在裁剪过图像之后做的

打开已经裁剪过的图像
在右侧工具栏找到“Radiometric correction”——“Radiometric Calibration”
选择要进行辐射定标的影像，出现如下图所示的工具框
选择“OK”后，在“Output Interleave”选项下拉选择：BIL，“Scale Factor”输入0.10，选择文件输出路径，点击OK，此时完成2018年裁剪影像的辐射校正；用同样的步骤对2020年的裁剪影像进行辐射校正
要做的内容都定标完成后辐射定标就结束了

三、大气校正

先进行高程计算：“File”—“Open World Data”—“Elevation”，打开高程地图
在右侧工具栏双击“Statistics”—“Compute Band Statistics”，在出现的工具框中选择刚才打开的高程地图，点击“Stats Subset”
在出现的工具框中点击File
选择2018年辐射校正完的图像（自己做的话就是自己已经辐射校正完的就行），点击“OK”
出现的“Compute Statistics Parameters”工具框中点击“OK”
记录下地面高程值34.043573，但这时的数值单位是米，换算成千米就是0.034（后期会用到这个千米的高程值）
在原始解压的数据文件中，打开2018年图像的MTL文本文件（即头文件，记事本打开），查找遥感图像的时间并记录（后来我发现选上LandSat 8 好像自动出现时间）
在右侧工具栏找到“Radiometric Correction”–“Atmospheric Correction Module”–“FLAASH Atmospheric Correction”
在工具栏“Input Radiance Image”选择2018年辐射校正完的文件，下图选择“Use single seale feator for all bands”。选择输出路径，在“Sensor Type”一栏选择2018年
根据之前的操作记录填“Flight Date”栏的时间和“Ground Elevation”栏的高程值
点击Help，根据图像的纬度值和时间选出对应“Atmospheric Model”的选项
在Multispectral Settings"中点击“Defaults”在下拉选项中选择第一项“Over-Land Retrieval standard(660:2100nm)”，点击“OK”
点击“Apply”运行处理
处理成功后会有这个界面

四、几何校正

要有两期辐射校正完的影像

打开经典版ENVI
我这里以年份为示例：
打开2020年大气校正后的影像，波段选择RGB打开：即近红外，红光，绿光
打开2018年大气校正后的影像，波段选择RGB打开：即近红外，红光，绿光
在工具条选择Map—Registration—Select GCPs Image to Image.
以2020年的影像作为参照影像，2018年的影像作为被校正影像。
合理进行布点选择：
要求同名点选取，不能有三个点在一条直线上，在像元角点进行选择
Goto可以找到该点在图像上的位置，On/Off可以打开和关闭该点，Delete可以删除该点；找到至少20个点，总误差RMS Error小于0.1且点在图像上均匀分布
保存几何校正图像Options—Warp File（as Image to Map）…–2018dq—OK
在出现的“Registration Parameters”工具框“Degree”选项改为：2，“Resampling”选择三次卷积内插法“Cubic Convolution”，选择保存路径，点击ok