记录第一次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.结束啦
有看到的客官的话欢迎批评指正,第一次写没有经验感谢大家!