记录学习
一、数据介绍
二、数据下载
1.地理空间数据云(地理空间数据云 (gscloud.cn))
2.USGS官网(https://earthexplorer.usgs.gov/EarthExplorer (usgs.gov))
三、预处理
Landsat 8-OLI 数据是在地理空间数据云里下载的,这些数据是已经几何校正过的
融合,辐射校正、大气校正、裁剪
3.1影像融合
将8波段30 m分辨率的多光谱数据和15 m的全色数据进行融合,从而得到15 m分辨率的多光谱数据
1、启动ENVI 5.3,选择File——Open,选择MTL.txt文件打开;选择“B8.TIF”文件打开。
2、在工具箱(Toolbox)中,点击Image Sharpening,双击Gram-Schmidt Pan Sharpening打开对话框。在对话框中选择低分辨率的多光谱数据文件,点击OK,再选择高分辨全色波段数据文件,点击OK。
3、在Pan Sharpening Parameters参数面板如图下图,选择传感器类型为:Landsat8_oli,重采样方法选择Cubic Convolution,设置输出路径和文件名。
3.2辐射定标
辐射定标是将DN值转换为辐射亮度值、地表反射率等相关物理量的相对值的过程
第一步:在ENVI 5.3的工具箱中查找工具:Radiometric Correction——Radiometric Calibration双击,选择要校正的多谱数据(而本次试验选择以上融合的多光谱),“LC81200332020291LGN00_rh.dat”进行辐射定标。
第二步:选择参数对话框中点击:【Apply Flaash Setting】后会自动读取元数据中的正射参数,也可以自己输入参数,主要参数如下图:选择多光谱数据。定标类型选择辐射亮度值(Radiance),输出格式为BIL,因为后续FLAASH大气校正的输入数据类型为BIL。输出数据类型为Float,系数为0.1。或者直接点击应用FLAASH设置按钮设置以上参数。确定输出路径后点击确定。
3.3大气校正
打开Radiometric Correction>Atmospheric Correction Module>FLAASH Atmospheric Correction工具,导入经过辐射定标的数据,因为在辐射定标中已经进行了单位换算,所以在弹出的Radiance Scale Factors中选择第二项
Lat和Lon都是自动获取的不用管。
传感器类型(Sensor Type):MultiSpectral(多光谱数据)》Landsat8 OLI。
气溶胶模型:渭-库绿洲(除6-9月用MLS(中纬度夏季),其他用SAS
查看成像时间(Flight Date)
方法一:在Layer Manager(图层管理器)窗口右击原始影像,按照提示选择View Metadata(查看元数据),点击Time查看成像时间
方法二:从MTL文件中找到。具体名称:DATE_ACQUIRED = 2013-05-12;SCENE_CENTER_TIME = 02:55:26.6336980Z;
查看平均高程值 :
方法一:由于FLAASH大气校正需要影像区域的平均高程,所以可以使用ENVI自带全球高程数据进行计算。选择File》Open World Data》Elevation,打开ENVI自带的全球900mDEM数据(GMTED2010.jp2数据)。再将需要计算高程数据的影像打开(需要带坐标信息),可以放大看到影像数据叠加在DEM数据上。打开工具箱(Toolbox)中的Statistics——Compute Statistics,在输入文件对话框中选择GMTED2010.jp2数据,再单击Stats Subset,单击File,选择需要统计高程信息对应的图像,然后OK。之后的参数均设置为默认。最后得到统计的平均高程信息,可以看到影像平均高程为1690.69m。
方法二:使用FLAASH Setting Guide插件计算高程(ENVI扩展工具:FLAASH 参数设置向导 - ENVI-IDL技术殿堂 - 博客园 (cnblogs.com))
在ToolBox中,启动\Extensions\FLAASH Setting Guide,双击打开
Input Radiance Raster:输入影像(FLAASH大气校正要求输入数据为Radiance),
Maritime Raster:输入影像是否位于海上(用于判断气溶胶反演方法),默认为No
设置多光谱参数,进入Multispectral Settings,将Default设置为Over-Land Retrieval Standard(600:2100)。其他参数按照默认设置,最后点击Apply进行大气校正。
图 FLAASH 大气校正结果中获取的水体波谱曲线(左-校正前,右-校正后)
显示大气校正结果图像,查看像元值,可以看到像元值扩大 10000 后,值在几百到几千不等。如果要得到 0-1 范围内的反射率数据,可以使用 BandMath 除以 10000.0。
图 FLAASH 大气校正结果中获取的植被波谱曲线(左-校正前,右-校正后)
大气校正前与校正后图像对比如图12,上图为校正前,下图为校正后,可以看到校正后的图像亮度有所增强。同时对比校正前后同一像素点植被的波谱曲线,右图为校正后波谱曲线,左方为校正前,可以看到校正后植被的波谱曲线趋势符合真实植被波谱。
3.4 裁剪
将研究区矢量文件加载进来(注:至于行政边界线的矢量文件如何获取, 行政边界矢量数据下载教程),加载进来后,在右面工具箱搜索“Subset Data from ROIs,双击打开,选择要裁剪的影像然后进入如下界面,这里是选择ROI区域,也就是裁剪区。
“Mask pixel outside of ROI?”
选择yes,选yes之后意思就是按照不规则边界区域进行裁剪,
选择no就是根据ROI 的那个矩形框进行裁剪,裁剪出来是个矩形。
“Mask Background”背景值设置为零,意思是黑色。