虽然我不专业但这是我一直想写的东西从这里起个头,记录一下ENVI遥感图像处理等基础(遥感专业的东西不要问我)。上周末开始一直在加班就焦人的很。让我安心借调吧。
公众号原文链接,欢迎关注
ENVI下Landsat8辐射定标与大气校正不专业的Landsat8辐射定标与大气校正https://mp.weixin.qq.com/s/rQd3Wm7MoHairamNrPPdRA
一、Landsat8简介
首先了解一下Landsat-8的各个波段情况,Landsat-8卫星上携带两个传感器,分别是OLI陆地成像仪和TIRS热红外传感器。卫星一共有11个波段。
波段名称 | 波长范围(μm) | 分辨率(m) |
Band 1 Coastal -海岸波段 | 0.43 – 0.45 | 30 |
Band 2 Blue - 蓝波段 | 0.45 – 0.51 | 30 |
Band 3 Green - 绿波段 | 0.53 – 0.59 | 30 |
Band 4 Red - 红波段 | 0.64 – 0.67 | 30 |
Band 5 NIR - 近红外波段 | 0.85 – 0.88 | 30 |
Band 6 SWIR 1 - 短波红外1 | 1.57 – 1.65 | 30 |
Band 7 SWIR 2 - 短波红外2 | 2.11 – 2.29 | 30 |
Band 8 Pan - 全色波段 | 0.50 – 0.68 | 15 |
Band 9 Cirrus - 卷云波段 | 1.36 – 1.38 | 30 |
Band 10 TIRS 1 - 热红外1 | 10.6 – 11.19 | 100 |
Band 11 TIRS 2 - 热红外2 | 11.5 – 12.51 | 100 |
而不同的波段组合也有不同的作用,Landsat8常见的波段组合有以下这些。
自然色(Natural Color) | 4 3 2 |
假彩色(False Color (urban)) | 7 6 4 |
红外彩色(Color Infrared (vegetation)) | 5 4 3 |
农业(Agriculture) | 6 5 2 |
大气渗透(Atmospheric Penetration) | 7 6 5 |
健康植被(Healthy Vegetation) | 5 6 2 |
陆地/水体(Land/Water) | 5 6 4 |
大气去除(Natural With Atmospheric Removal) | 7 5 3 |
短波红外(Shortwave Infrared) | 7 5 4 |
植被分析(Vegetation Analysis) | 6 5 4 |
(资料介绍来自USGS官网和ENVI技术殿堂博客)
二、数据获取与简介:
数据来源为地理空间数据云Landsat8 OLI_TIRS影像数据。剑阁县矢量边界。考虑到之后想统一拿一套数据进行操作,经过一段时间的上班摸鱼,最后选择了我在四川经历的故事比较多的一个地方—剑门关!位于四川省广元市剑阁县。也是因为只用下载一景便覆盖了这个县,不用费事再去拼接影像。
接着看一下下载下来的数据,Landsat8命名是有规则的,像这样:
1、文件命名:例如我下载下来的压缩包名称是LC81290382019223LGN00
LC8:landsat-8卫星
129:WRS path
038:WRS row
2019:影像获取年份
223:影像获取日期,从1月1为1算起,即223即8月11日
LGN:接站代码
00:产品级别
2、数据命名:打开里面的数据会有这样格式的
LC08_L1TP_129038_20190811_20190820_01_T1_ANG(投影信息)或MTL(影像元数据)或BQA(影像质量信息)或B1-B11(波段),ENVI中打开数据就是打开MTL文件。
LC08:Landsat-8卫星
L1TP(Level 1 Precision Terrain (Corrected)):表示L1级数据,TP表示该数据已经进行地形校正和几何校正。
129038:获取的数据地区编号在WRS-2参考系统中的条带号为129,行编号为38。
第一个时间20190811即数据获取时间,第二个时间20190820即数据处理时间。
既然说起WRS和条带号,就需要了解一下是什么。目前WRS有两个系统,分别为WRS-1(1983年之前)和WRS-2(1983年之后)。Landsat1-3使用WRS-1,Landsat4、5、7、8使用WRS-2。了解WRS的目的主要是为了能够更方便的定位影像。例如:我们想下载某一地区的Landsat影像数据,可以根据该区域所在的Path/Row来定位影像区域。对着图找一下,129038是不是所下载的区域。另外USGS官网提供行带号矢量文件的下载,可以加载至Arcgis进行查看。
我国范围内WRS-2 Path/Row与行政区划、经纬度对照图
三、辐射定标
辐射定标是将DN值转换为辐射亮度值、反射率值和温度等物理量,而DN值本身就是一些无量纲的数字,没有什么物理意义,就是数字而已,遥感图像通常是用DN值记录信息的,当进行遥感定量化分析时候就需要进行转换。
打开【Radiometric Correction】中的ENVI中辐射定标工具【Radiometric Calibration】,选择MultiSpectral多光谱数据。定标类型选择辐射亮度值(Radiance),输出格式为BIL,因为后续FLAASH大气校正的输入数据类型为BIL。输出数据类型为Float,系数为0.1。或者直接点击Apply FLAASH Settings按钮设置以上参数。点击ok,等待一会,右图为辐射定标后的影像。
四、影像裁剪
这一步就是把整景影像按照研究区域裁剪出来,如果不想多做一步裁剪的话,在上一步中有个Mask掩膜,这里不做,仅做一下裁剪。
首先,将剑阁县边界的shp格式转为evf格式的,打开ENVI Classic(经典版本那个)。点击【File】—【Open Vector File】打开有我们shp边界数据的文件夹,选择shp格式的那个文件,因为影像是UTM投影,矢量数据需要保持一致。
点击ok,进行下一步,选中这个文件,【Load Selected】,就可以看到左边这个边界是不是就出来了。然后关掉这些界面即可。
打开本来用的ENVI5.3,【File】—【open】找到保存边界的文件夹,点击加载即可,看,这样就进来了。至于为什么要转换为evf格式。
使用【Subsat Data from ROIs】工具,选择辐射定标后的影像,点ok。
下一步,Mask pixels outside of ROI选择YES,不然裁剪的时候是按照边界的外接矩形裁剪的。点击ok。裁剪后的影像如右图。
五、大气校正
大气校正的目的是消除大气和光照等因素对地物反射的影响,获得地物反射率、辐射率、地表温度等真实物理模型参数,用来消除大气中水蒸气、氧气、二氧化碳、甲烷和臭氧等对地物反射的影响,消除大气分子和气溶胶散射的影响。
Radiometric Correction—Atmospheric Correction Module—FLAASH Atmospheric Correction。这个界面有几个需要设置的地方:
1、文件输入与输出项目
【Input Radiance Image】 文件输入就是上面辐射定标裁剪后的影像。
【Output Reflectance File 】输出信息项目文件路径最好保持一致。
选择输入文件后会弹出个对话框,选第二个。
(2)传感器与图像信息
Lat和Lon都是自动获取的不用管。
传感器类型(Sensor Type):多光谱数据(MultiSpectral)—Landsat8 OLI。
图像区域平均海拔(Ground Elevation): 0.68km,ENVI里自带90m的高程数据,如下图,找到软件所在的路径,可以通过统计工具统计高程,可以看到平均高程是684.7米,或者可以直接百度百科查一下也行。
成像时间(Flight Date),右键图层,View Metadata可查,2019-8-11-3:33:06,也可以从MTL文件中找到。
(3)大气模型(Atmospheric Model):Tropical,依据图像纬度和8月成像时间。
Latitude | Jan. | March | May | July | Sept. | Nov. |
80 | SAW | SAW | SAW | MLW | MLW | SAW |
70 | SAW | SAW | MLW | MLW | MLW | SAW |
60 | MLW | MLW | MLW | SAS | SAS | MLW |
50 | MLW | MLW | SAS | SAS | SAS | SAS |
40 | SAS | SAS | SAS | MLS | MLS | SAS |
30 | MLS | MLS | MLS | T | T | MLS |
20 | T | T | T | T | T | T |
10 | T | T | T | T | T | T |
0 | T | T | T | T | T | T |
-10 | T | T | T | T | T | T |
-20 | T | T | T | MLS | MLS | T |
-30 | MLS | MLS | MLS | MLS | MLS | MLS |
-40 | SAS | SAS | SAS | SAS | SAS | SAS |
-50 | SAS | SAS | SAS | MLW | MLW | SAS |
-60 | MLW | MLW | MLW | MLW | MLW | MLW |
-70 | MLW | MLW | MLW | MLW | MLW | MLW |
-80 | MLW | MLW | MLW | SAW | MLW | MLW |
(4)气溶胶模型(Aerosol Model):Rural郊区
(5)气溶胶反演(Aerosol Retrieval): 2-Band(K-T)图里选错了。这个反演方法的话我查了下是必须要有短波红外才能反演,哨兵和Landsat都有。
(6)多光谱设置(Multispectral Settings): 右边那个设置,Defaults下拉框选第一个。点击ok。
其它都保持默认即可,点击Apply,稍微一等,时间稍微会长一点点,最后得到大气校正后的影像(假彩色显示)。
做的时候出现了报错,总是提示内存资源不足,百度的方法都试了没啥用,直接重新装了遍软件就好了....CSDN有一篇博客总结了大气校正常见的报错,如碰上可去搜一下看看。
END
如果在实验过程中有什么地方遇到问题可以后台留言,竭力解答,大概率不会。需要练习数据也请后台留言。如果可以的话,希望能够点个关注转发和分享,点个在看并且点个赞~~
欢迎关注公众号【一点规划】