记录第一次ENVI实操（辐射校正+大气校正）
一、研究区域裁剪
首先我们要有两期不同年份同一地区的遥感影像（可以通过这里选取数据集进行下载）
根据个人喜好我首先会选择先进行影像的裁剪，尽可能地将图像放在研究区域内
-
“File”中打开影像
-
“File”——“save as”——“save as …（ENVI）”
-
在弹出来的窗口中选择头文件影像
-
点击右侧的箭头进行选择研究区域的裁剪
这张图里在右边黑框上边“第三个图标”可以选择希望作为底图的影像直接进行裁剪，我这里选择了用剪切过的2018年的影像来裁剪新的
选取完成后选择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”一栏选择“LandSat 8”（下载的数据是哪个平台就选几）
9. 根据之前的操作记录填“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
END.结束啦
有看到的客官的话欢迎批评指正，第一次写没有经验感谢大家！
1666
被折叠的 条评论
为什么被折叠?
到【灌水乐园】发言
